• Scopus数据库收录期刊
  • 中国科技核心期刊
  • 全国中文核心期刊
  • 美国工程索引(EI)收录期刊

二阶Stokes波作用下多孔弹性海床动力响应解析解

王哲, 周博豪, 张智卿, 李西斌

王哲, 周博豪, 张智卿, 李西斌. 二阶Stokes波作用下多孔弹性海床动力响应解析解[J]. 岩土工程学报, 2024, 46(9): 1791-1799. DOI: 10.11779/CJGE20230504
引用本文: 王哲, 周博豪, 张智卿, 李西斌. 二阶Stokes波作用下多孔弹性海床动力响应解析解[J]. 岩土工程学报, 2024, 46(9): 1791-1799. DOI: 10.11779/CJGE20230504
WANG Zhe, ZHOU Bohao, ZHANG Zhiqing, LI Xibin. Analytical solutions for dynamic response in poroelastic seabed under second-order Stokes waves[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(9): 1791-1799. DOI: 10.11779/CJGE20230504
Citation: WANG Zhe, ZHOU Bohao, ZHANG Zhiqing, LI Xibin. Analytical solutions for dynamic response in poroelastic seabed under second-order Stokes waves[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(9): 1791-1799. DOI: 10.11779/CJGE20230504

二阶Stokes波作用下多孔弹性海床动力响应解析解  English Version

基金项目: 

浙江省自然科学基金重点项目 LHZ21E090001

国家自然科学基金项目 52178367

详细信息
    作者简介:

    王哲(1978—),男,博士,教授,主要从事岩土工程的教学和科研工作。E-mail: wangzsd@zjut.edu.cn

    通讯作者:

    李西斌, E-mail: ytulxb@zafu.edu.cn

  • 中图分类号: TU43

Analytical solutions for dynamic response in poroelastic seabed under second-order Stokes waves

  • 摘要: 基于二阶Stokes波理论,研究了非线性海洋波作用下多孔弹性海床的动力响应问题。利用Biot固结理论以及Verruijt储存方程,在笛卡尔坐标下建立了海床的控制方程。通过将海洋波函数拓展到复数域,采用严格的数学推导求得了海床动力响应的解析解,并将该解与既有解进行对比验证。最后分析了波浪和海床特征参数对海床有效应力、剪切应力和孔压分布的影响,结果表明波浪和海床的参数对海床动力响应具有显著的影响。其中波浪周期和水深对二阶Stokes波的二阶项有较大影响,渗透系数和剪切模量则分别影响海床动力响应的变化速率和幅值。
    Abstract: Based on the theory of the second-order Stokes waves, the dynamic response of a poroelastic seabed induced by nonlinear ocean waves is investigated. The governing equations for the seabed are established in Cartesian coordinates using the Biot's general consolidation theory and the Verruijt's storage equation. The analytical solutions for the dynamic response of the seabed are obtained by extending the ocean wave function and the field quantities in the seabed to the complex domain using a rigorous mathematical derivation. The correctness and reliability of the analytical solutions are verified by comparing with the existing solutions. Finally, the effects of the characteristic parameters of the waves and seabed on the vertical effective stress, horizontal effective stress, shear stress and pore pressure distribution of the seabed are analyzed. The results show that the characteristic parameters of the waves and seabed have a significant effect on the dynamic response of the seabed. The wave period and water depth have a significant impact on the second-order term of the nonlinear waves, while the permeability and shear modulus respectively affect the rate of change and amplitude.
  • CO2+O2原地浸出过程中,溶浸液在砂岩内渗流,受到砂岩孔隙率的影响。同时溶浸液和砂岩接触发生化学反应,将U4+转化成U6+稳定地存在溶浸液内,反应侵蚀使砂岩内部不断产生孔隙、裂隙,影响着物理渗流过程。低渗透砂岩孔隙率低,溶浸液和砂岩接触的面积小,化学反应铀浸出量较低,开采就比较困难。CO2+O2浸出开采低渗透砂岩的关键在于提高低渗透砂岩的孔隙率,增强砂岩的渗透性[1]。如何有效地提高低渗透砂岩的孔隙率,使不适合CO2+O2地浸开采的低渗透砂岩矿层达到适合地浸开采所需的渗透率,对于提高浸出效果具有重要的实际经济意义。

    为提高浸出效果,世界各国学者采用了爆破爆炸[2]、表面活性剂[3-4]、振动波的方式。其中,爆破爆炸是利用炸药的爆破冲击波和应力波压碎岩体并形成裂隙[5-7],从而增加低渗透砂岩层的渗透性。表面活性剂是通过降低溶浸液表面张力[3],增加砂岩润湿性[4],使溶浸液极易进入微小孔隙、裂隙,从而提高溶浸液在低渗透砂岩内的渗透率。振动波有超声波振动和低频机械振动两种形式,超声波振动主要原理为:超声波引起空化气泡振动,使固体成分快速地溶解于溶浸剂,溶浸液和固体颗粒充分接触反应,从而加快了浸出。目前研究得出超声波能增大锗[8]、铜[9]、镍[10]、金[11]、银[12]、稀土[13]、锌[14]、铀[15-16]的浸出率,提高浸出速度。但是超声波浸出研究针对散体颗粒,未见针对块体岩石的报道。低频机械振动对铀浸出方面的研究目前仅有Makaryuk[17]进行了报道,在野外现场应用低频机械振动改善了砂岩层的颗粒堵塞,但是颗粒堵塞可能是人为钻井引起,无法直接说明振动改变了砂岩自身的孔隙率和渗透率,也未针对CO2+O2浸出条件下进行研究。因此,有必要将低频机械振动引入低渗透砂岩CO2+O2浸出,针对浸出过程中低频机械振动对砂岩的孔隙率、渗透率影响展开研究。

    本文从物理和化学角度出发,建立浸出过程中低频振动含铀砂岩孔隙率方程及渗透率表达式。并利用试验对渗透率表达式进行检验,验证方程和表达式所反映的孔隙率及渗透率变化规律的正确性,从而为CO2+O2原地浸出低渗透砂岩提供理论依据。

    依据孔隙率定义[18-19]

    φ=VPVB, (1)

    式中,VP为砂岩孔隙体积(m3),VB为砂岩总体积(m3)。

    CO2+O2振动浸出过程中,砂岩受到物理和化学侵蚀作用。不断的振动使砂岩受到物理机械损伤,而产生振动孔隙体积。同时振动使溶浸液加速流动,从而加速溶浸液和砂岩内部微孔裂隙的接触,砂岩微孔裂隙不断受到溶浸液涌入,并发生化学反应,产生化学侵蚀孔隙体积。因此砂岩孔隙体积是由振动孔隙体积、化学侵蚀孔隙体积、初始孔隙体积三部分构成。由于振动使砂岩骨架变得疏松膨胀,溶浸液侵蚀使砂岩骨架体积变小,因此孔隙率可表示为[19]

    φ=11φ01+εv[1+ΔVSσVchVS0]=11φ01+εv[1+ΔVVS0VchVB(1φ0)] (2)

    式中 εv为砂岩体积应变;φ0为振动作用下砂岩初始孔隙率;ΔV为外力振动引起的骨架体积变化量(m3);Vch为化学侵蚀孔隙体积(m3);VS0为砂岩初始骨架体积(m3)。

    由式(2)可知,找出化学侵蚀孔隙体积和振动产生的膨胀应变,就可以构建振动孔隙率方程。

    (1)理想状态下次生孔隙体积

    CO2+O2浸出铀的过程中,主要发生的化学反应为[20-21]

    2UO2+O2+4HCO=3=2[UO2(CO3)2]2+2HO2 (3)

    反应过程中溶浸液呈弱酸性,依据文献[22],溶液中Ca2+、Mg2+增多引起的孔隙率变化主要是由于CaCO3、CaMg(CO3)2溶解而产生,化学反应为

    CaCO3+CO2+H2O=Ca2++2HCO3, (4)
    CaMg(CO3)2+2CO2+2H2O=Ca2++Mg2++4HCO3 (5)

    由式(3)~(5)可知:1 mol的UO2与碳酸氢根反应后生成1 mol [UO2(CO3)2]2-,1 mol的CaCO3与二氧化碳、水反应后生成1 mol Ca2+,1 mol的CaMg(CO3)2化学反应生成1 mol Ca2+和1 mol Mg2+

    单位物质的量的矿物在理想状态下反应后的次生孔隙体积可由下式计算[23]

    Vch=V0ji=1Vi=N0Mr0ρ0ji=1NiMriρi (6)

    式中 Vch为理想状态下反应中某矿物发生反应形成的次生孔隙体积(m3);V0为反应矿物的初始体积,m3Vi为生成矿物的体积(m3);N0为反应物的单位物质的量(g/mol);Mr0为反应物的分子量;ρ0为反应物的密度(kg/m3);Ni为单位反应物物质的量对应生成物的物质的量(g/mol);Mri为第i个生成沉淀物的分子量;ρi为第i个生成沉淀物的密度(kg/m3);j为反应中生成沉淀物的总数。

    CO2+O2浸出条件下,CO23含量非常低,溶液中的碳酸盐全部以HCO3形式存在[20],因此可认为式(3)~(5)反应生成物中均没有直接生成沉淀物质,生成的Ca2+、Mg2+均以离子态被溶解到溶液内。则式(6)可简化为

    Vch=N0Mr0ρ0 (7)

    将CaCO3、UO2、CaMg(CO3)2的分子量、密度、反应单位物质的量统计于表13内,由式(7)计算可得出各个反应的次生孔隙体积。

    表  1  CaCO3反应产生的次生孔隙体积
    Table  1.  Secondary pore volume produced by CaCO3 reaction
    类别单位/mol分子量密度/(g·cm-3)体积/cm3次生孔隙体积/cm3
    CaCO311002.7136.936.9
    下载: 导出CSV 
    | 显示表格
    表  2  含铀矿物反应产生的次生孔隙体积
    Table  2.  Secondary pore volume produced by uranium-bearing mineral reaction
    类别单位/mol分子量密度/(g·cm-3)体积/cm3次生孔隙体积/cm3
    UO2127010.9624.6424.64
    下载: 导出CSV 
    | 显示表格
    表  3  CaMg(CO3)2反应产生的次生孔隙体积
    Table  3.  Secondary pore volume produced by CaMg(CO3)2 reaction
    类别单位/mol分子量密度/(g·cm-3)体积/cm3次生孔隙体积/cm3
    CaMg(CO3)211842.8564.5664.56
    下载: 导出CSV 
    | 显示表格

    (2)实际形成次生孔隙体积

    实际反应后,溶液中的离子不可能以整数物质的量状态形式存在,实际参加反应的反应物物质的量数应根据反应后溶液中实测离子来计算。

    由式(3),(5)可看出:1 mol UO2与1 mol CaMg(CO3)2分别对应生成1 mol [UO2(CO3)2]2-和1 mol Mg2+,CaMg(CO3)2和CaCO3发生反应后均有Ca2+产生。UO2或CaMg(CO3)2溶解生成的次生孔隙体积可由式(6)变化为

    Vch1=Vrc1Mr1ρ1 (8)

    式中 Vch1为UO2或CaMg(CO3)2溶解生成的次生孔隙体积(m3);Vr为反应后溶液的体积(m3);c1为反应后实测U6+或Mg2+浓度(mol/L);Mr1为UO2或CaMg(CO3)2的分子量;ρ1为UO2或CaMg(CO3)2的密度(kg/m3)。

    CaCO3溶解生成的次生孔隙体积可由式(6)变化为

    Vch2=Vr(c2c1)Mr2ρ2 (9)

    式中 Vch2为CaCO3溶解生成的次生孔隙体积(m3);c1为反应后实测溶液中Mg2+浓度(mol/L);c2为反应后实测溶液中Ca2+浓度(mol/L);Mr2为CaCO3分子量;ρ2为CaCO3密度(kg/m3)。

    实际反应后生成的总次生孔隙体积可表示为

    Vch=Vch1+Vch2=Vrc1Mr1ρ1+Vr(c2c1)Mr2ρ2 (10)

    由式(10)得到了化学侵蚀孔隙体积,只需再找出振动孔隙应变,就能得到振动孔隙率方程。

    (1)振动应力分析

    振动过程中,应力波在砂岩中以类正弦波传播并不断的衰减。当砂岩不受振动时,砂岩所受到的应力为0,则半个周期内的均值应力可表示为

    F=Tz20F0sin(2πωt)dt=F0πω (11)

    式中 ω为振动源的振动频率(Hz);F0为振动传播到砂岩时的力幅值(N);Tz为振动周期(s);t为砂岩受振动作用的时间(s)。

    一个周期内均值应力和半个周期内的均值应力大小相等,则t时刻岩体受到的应力σf可表示为

    σf=1TzTz20F0sin(2πωt)dttTz=F0ωtπ (12)

    应力波传播过程中振动作用时间t内,单位体积的砂岩所受到的振动应力可表示为[24-25]

    σf=1Vss00F0ωtπeηsds=ρsF0ωtmsπη(1eηs0) (13)

    式中 s0为振动在砂岩内传播的距离(m);η为应力波的衰减率;Vs为应力波影响到的砂岩体积(m3);ms为应力波影响到的砂岩质量(kg);ρs为砂岩视密度(kg/m3)。

    (2)损伤分析

    应力波在砂岩内传播过程中,损伤度D可表示为[26-27]

    D=1v2Dv2 (14)

    式中 vD为应力波在损伤试件中的传播速度(m/s);v为应力波在均质试件中的传播速度(m/s)。

    由于地震波、机械波、固体中的声波、超声波均为应力波,则波速衰减速率可表示为[28-29]

    η=vvDv (15)

    将式(14)代入式(15)可得

    η=11D (16)

    损伤度D又可表示为[30]

    D=D0+mεnv (17)

    其中,

    n=EmdEEmd, (18)
    m=(1D0)(n+1)εnd,max, (19)
    Emd=σf,maxεd,max (20)

    式中 Emd为损伤砂岩的峰值弹性模量(Pa);σf,max为受损砂岩受到的最大应力值(N);εd,max为受损砂岩的最大应变值;D0为砂岩的初始损伤;E为砂岩的初始弹性模量(Pa)。

    实际振动过程中,砂岩所受的最大振动应力值为未发生衰减的应力值;最大应变值为未衰减的振动应力作用下,物理应变值与化学侵蚀产生的化学应变值之和。对单位体积砂岩而言,应变的最大值可表示为

    εd,max=εch,max+εf,max=VchVC+σf,maxEv=VchVC+1EvVss00F0ωtπds=VchVC+3(12ν)ρsF0ωts0Emsπ (21)

    式中 εch,max为单位体积的化学应变最大值;εf,max为单位体积的物理应变最大值;VC为参加反应的砂岩反应前总体积(m3);Ev为砂岩体积模量(Pa);ν为砂岩泊松比。

    此时,振动作用下式(20)可变化为

    Emd=σf,maxεd,max=1VchVCσf,max+1Ev=EvVCσf,maxVchEv+VCσf,max (22)

    损伤计算过程中,正应变和负应变均使砂岩受到损伤,两者不能相互抵消,因此,在式(17)中,t时刻砂岩体积应变εvt时刻砂岩化学应变值与物理应变值之和,即

    εv=εch+εf=VchtVCtch+σfEv=VchtVCtch+F0ωtEvπ (23)

    式中 εch为化学应变值;εf为物理应变值;tch为化学反应的总时间(s)。

    将式(23)代入式(17),可得损伤表达式:

    D=D0+mεnv=D0+m(VchtVCtch+F0ωtEvπ)n (24)

    式(24)将化学侵蚀产生的应变、机械振动应变同时考虑,是物理和化学共同作用下的损伤表达式。

    (3)振动应变分析

    将式(24)代入式(16)可得衰减率的表达式:

    η=11D=11D0m(VchtVCtch+F0ωtEvπ)n (25)

    衰减振动力产生的应变可表示为

    ΔVSσVS0=σfEv=ρsF0ωtmsπη(1eηs0)Ev=3ρs(12ν)F0ωt{1e[11D0m(VchtVCtch+F0ωtEvπ)n]s0}πmsE[11D0m(VchtVCtch+F0ωtEvπ)n] (26)

    将式(10),(26)代入式(2)可得出振动浸出过程中砂岩孔隙率方程:

    (27)

    该方程以孔隙率定义为基础,考虑浸出过程的孔隙体积由振动引起的孔隙体积、化学侵蚀产生的孔隙体积、初始孔隙体积三部分构成。从化学反应方程式出发,分析实际反应的过程,得出化学侵蚀产生的孔隙体积。将机械振动引入,建立机械振动和化学侵蚀共同作用下,受损伤砂岩的应变表达式,从而建立浸出过程中振动砂岩孔隙率方程。该方程将物理振动损伤和化学损伤相结合,更能体现出浸出过程中低频机械振动和化学侵蚀对砂岩孔隙率的影响。为更加清晰地表达孔隙率方程所含各个表达式之间关系,将方程表示为

    φ=11φ01+εv[1+ΔVSσVS0VchVB(1φ0)] , (28a)
    ΔVSσVS0=ρsF0ωt(1eηs0)msπηEv , (28b)
    Vch=Vrc1Mr1ρ1+Vr(c2c1)Mr2ρ2 , (28c)
    η=11D0m(VchtVCtch+F0ωtEvπ)n , (28d)
    n=EmdEEmd  , (28e)
    m=(1D0)(n+1)εnd,max  , (28f)
    Emd=EvVCσf,maxVchEv+VCσf,max  , (28g)
    εv=VchtVCtch+F0ωtEvπ    (28h)

    为进一步获得振动条件下渗透率变化规律,必须找出渗透率和孔隙率之间的关系,其关系可表示为[31-32]

    k=k0(φφ0)3, (29)

    式中,k为振动作用下砂岩渗透率(D),k0为砂岩初始渗透率(D)。

    将式(27)代入式(29)即可获得振动浸出过程中砂岩渗透率表达式。

    目前尚未将低频振动应用于CO2+O2原地浸出低渗透砂岩,实际浸出过程中仅能测得未受到振动的砂岩孔隙率和渗透率。为检验建立的孔隙率方程及渗透率表达式是否能够正确反映低渗透砂岩孔隙率和渗透率变化规律,必须采用试验的方法进行检验。

    由于试验采用的试样孔隙率较小,现有仪器测量误差大,并且需要找出试样渗透率和孔隙率两者的变化规律。因此,试验采用渗透率变化来反应砂岩孔隙率的变化规律,通过测量振动前后试样的渗透率变化反映孔隙率变化。

    实际浸出过程中,由于顶板已经把纵向压力分解,处于含水层内的疏松砂岩仅仅受到水压作用,并且试样取样过程中会受到不同程度的振动影响。因此,试验采用压制的试样模拟砂岩原始的赋存状态,在无围压条件下进行振动试验。

    试验装置主要有伺服压力机、振动浸出试验系统、渗透仪。伺服压力机是为了压制试样,振动浸出试验系统是主要的反应场所,渗透仪可检测振动前后试样孔隙率变化规律。振动浸出试验系统如图1所示。

    图  1  振动浸出试验系统结构图
    Figure  1.  Structure of vibration leaching experimental system

    为正确反映浸出过程中振动对低渗透砂岩渗透率的影响,相同振动时间时,选用同一钻孔取得的原矿砂岩破碎,加入方解石、蒸馏水均匀混合,压制成砂岩试样。这样可确保相同振动时间选用的试验试样相同。不同振动时间时,选用不同钻孔的原矿砂岩压制,试验试样存在一定差异性,因此,试验可区分不同振动频率对试样渗透率的影响,同时也可以尽可能地还原地下实际浸出状况。将同一频率作用下试样浸出过程,看作地下不同深度的砂岩分层浸出过程,将后一个时间的渗透率累计相加前面所有时间的渗透率,可以更好地反映浸出实际中的同一频率作用下,渗透率随振动时间的变化规律。钻取的原铀矿样如图2所示,图中保鲜膜包裹的单个矿样为同一钻孔取得,不同钻孔取得的原矿样长短不一。振动试验方案如表4所示。

    图  2  钻取的原铀矿样
    Figure  2.  Drilled raw uranium ore samples
    表  4  振动试验方案
    Table  4.  Experimental scheme
    压样编号采用原矿振动时间/h
    10 Hz20 Hz30 Hz
    A1A2A3W-78
    B1B2B3W-1010
    C1C2C3W-412
    D1D2D3U-214
    E1E2E3U-1616
    F1F2F3U-518
    G1G2G3U-3220
    H1H2H3W-922
    下载: 导出CSV 
    | 显示表格

    试验过程主要有制备试样、振动浸出、振动前后渗透率测量4个过程。其中制备试样是试验的准备阶段,振动浸出是试验的主要阶段,渗透率测量是数据采集阶段。具体的试验过程如下:

    (1)制备试样

    将从某矿610 m深度左右钻取的低品位铀矿石破碎,筛取30~100目的矿石颗粒312 g,与800目的方解石104 g、蒸馏水30 g充分拌匀混合,装入压制模具。开启伺服压力机,以2.4 kN/s的加速度加载至204 kN后,保压20 min制得直径50 mm,高100 mm的圆柱体砂岩试件。压制成型的试样如图3所示。

    图  3  压制的砂岩试样
    Figure  3.  Pressed sandstone sample

    (2)振动前后渗透率测定

    为得到振动后试样渗透率变化规律,必须测定试样振动前后的渗透率。振动后减去振动前渗透率,可得到试样渗透率增量值,同时渗透率增量值可反映孔隙率的变化规律。

    振动前与振动后的渗透率测试方法相同,其测试过程为:将试样侧面包裹保鲜膜后装入渗透仪岩心夹持器,并抽真空。恒压作用下将试样饱和水,恒流驱替试样中的饱和水,直至夹持器前端压力表恒定,记录试样前后两端的压力值、泵流速,运用达西定律计算渗透率。

    (3)振动浸出试验

    将试样放置于通体开孔的内胆中,沿内胆侧面排间距14.13 mm,孔间隔6 mm,孔径4 mm的开孔进行烫孔。然后将内胆放置于反应釜内,抽真空5 min,向反应釜内注入1700 ml蒸馏水,并加热至30℃后充入1 MPa的氧气。调节减压阀向反应釜内继续充入二氧化碳,直至釜内气压达到2 MPa。打开扫频信号发生器、功率放大器,使激振器在相同力幅、不同频率下推动反应釜振动。达到试验时间后,关闭功率放大器及加热阀、打开通大气阀,将测试系统内的气体排空。打开反应釜取出试样,采集反应釜内的液体。

    10,20,30 Hz三个不同频率作用下,测得试样初始渗透率的变化范围是0.05~1.98 md,属于低渗透砂岩。振动后减去振动前试样的渗透率差如图4所示。

    图  4  低频振动作用下测定的砂岩渗透率增量变化图
    Figure  4.  Variation of measured sandstone permeability increment under low-frequency vibration

    图4中数据显示:在试验时间内,10 Hz渗透率增量变化范围为0.098~0.982 md,20 Hz渗透率增量变化范围为0.138~1.194 md,30 Hz渗透率增量变化范围为0.267~3.265 md。相对初始渗透率0.05~1.98 md,30 Hz渗透率是0 Hz的2.65~6.34倍,20 Hz渗透率是0 Hz的1.6~3.76倍,10 Hz渗透率是0 Hz的1.49~2.96倍;并且当振动时间相同时,振动频率越大,渗透率增量就越大;当振动频率相同时,渗透率增量随着振动时间呈波浪状变化趋势,但其数值均大于0。

    溶浸液在砂岩内渗流,不断侵蚀砂岩,使砂岩孔隙发生变化;同时振动作用使砂岩内部不断产生新的孔隙,侵蚀程度进一步增大。

    为分析砂岩侵蚀程度,选取振动22 h后,砂岩试样的外观进行对比,如图5所示。图5(a)~5(c)为不同频率振动后的试样顶底面外观。对比图中试样外观变化,可看出振动频率越大,试样被侵蚀程度越大,从而孔隙率、渗透率越大。

    图  5  试样顶面、底面外观
    Figure  5.  Appearance of top surfaces and bottom surfaces of samples

    试样初始孔隙率采用饱和及水中称重浮力法测定。由于试样体积较小,且孔隙率的测定需要较高的精度,为了减小误差,在计算过程中初始孔隙率采用实测试样孔隙率的平均值。试验参数见表5。化学反应后溶液中Ca2+、U6+、Mg2+的浓度值见表68

    表  5  试验参数
    Table  5.  Experimental parameters
    砂岩初始弹性模量E/Pa砂岩泊松比ν 砂岩反应前总体积VC/m3化学反应的总时间tch驱动力的幅值F0/N砂岩初始孔隙率φ0应力波影响到的砂岩质量m/kg砂岩的视密度ρs/(kg·m-3)振动在砂岩内传播距离S0/m砂岩初始损伤D0振动频率ω/Hz试验温度T/℃初始渗透率k0/md
    4.35×10100.251.9635×10-4数值与振动时间相同试验实测值为10.01324*0.4292.1849×103由于砂岩试样直径仅为0.05 m,计算过程中取其值为0.05 m压制试样认为其损伤值为010,20,30300.2994*
    注:*为取平均值。
    下载: 导出CSV 
    | 显示表格
    表  6  Ca2+浓度
    Table  6.  Ca2+ concentration
    振动时间/h振动频率/ Hz
    102030
    8116.23117.28132.26
    10108.22144.61212.42
    1284.17131.83272.54
    14108.22114.92124.25
    16328.66384.95533.06
    18198.40213.58224.45
    20320.88325.37345.39
    22247.53266.58298.56
    下载: 导出CSV 
    | 显示表格
    表  7  U6+浓度
    Table  7.  U6+ concentration
    振动时间/h振动频率/Hz
    102030
    80.0490.2390.645
    100.8090.9300.991
    120.1410.2400.400
    140.3870.5120.692
    161.3501.6092.360
    180.1290.1710.214
    200.1290.1840.214
    220.3100.4710.538
    下载: 导出CSV 
    | 显示表格
    表  8  Mg2+浓度
    Table  8.  Mg2+ concentration
    振动时间/h振动频率/Hz
    102030
    82.433.524.86
    104.865.567.29
    122.432.693.65
    144.865.296.08
    166.857.107.29
    185.365.686.08
    203.283.563.63
    222.983.093.56
    下载: 导出CSV 
    | 显示表格

    表58中的数据代入式(27),(29),计算得出的砂岩孔隙率及渗透率增量的变化规律如图6,7所示。图6中曲线显示:10 Hz孔隙率变化范围为0.0158~0.0204,20 Hz孔隙率变化范围为0.0168~0.0226,30 Hz孔隙率变化范围为0.0178~0.0265。

    图  6  方程计算砂岩孔隙率φ变化曲线
    Figure  6.  Variation curves of sandstone porosity φ calculated by equation
    图  7  方程计算砂岩渗透率增量∆k变化曲线
    Figure  7.  Variation curves of sandstone permeability increment ∆k calculated by equation

    图7中曲线显示振动后试样渗透率增量均大于零,10 Hz渗透率增量变化范围为0.2108~0.7564,20 Hz渗透率增量变化范围为0.3075~1.1956,30 Hz渗透率增量变化范围为0.4302~2.113。

    对比计算渗透率增量值和实测渗透率增量值可得出:计算值和实测值范围基本一致。对比图4,6,7可知:试验渗透率增量、计算渗透率增量、计算孔隙率变化规律一致,均为振动时间相同时,振动频率增加,值增加;振动频率相同时,振动时间增加,值呈现波浪状变化。其原因主要如下:

    由于试样制作过程中,对于相同振动时间,采用了同一钻孔取得的原矿砂岩颗粒,与方解石、蒸馏水均匀混合,制作出相同试样。不同振动时间选用不同钻孔的原矿砂岩压制,因此,不同振动时间采用的试样存在差异性,试样内物质成分有所差别,导致同一试验条件下,振动时间增加,渗透率、孔隙率增大程度不同,从而呈现波浪状变化。

    为更好地反映实际的浸出过程,将同一频率作用下试样浸出看作地下不同深度的砂岩分层浸出过程,将试验实测的渗透率增量进行累计相加,如图8所示。并将计算的渗透率增量累计相加,如图9所示。

    图  8  实测砂岩渗透率累计增量变化图
    Figure  8.  Variation of measured sandstone permeability cumulative increment
    图  9  计算砂岩渗透率累计增量变化图
    Figure  9.  Variation of calculated sandstone permeability cumulative increment

    对比图8,9可知:实测值和计算值变化范围基本一致,实测值10 Hz变化范围为0.098~2.3751,20 Hz变化范围为0.138~5.6915,30 Hz变化范围为0.267~10.1104;计算值10 Hz变化范围为0.2212~3.6129,20 Hz变化范围为0.3075~5.7689,30 Hz变化范围为0.4302~9.2162。实测值和计算值变化规律完全相同,均呈现线性变化,振动时间相同,振动频率增加,渗透率变大;振动频率相同时,振动时间增加,渗透率线性增加。因此,由以上分析可得出所建立的方程能正确反映铀浸出过程中低频振动对低渗透砂岩试样孔隙率、渗透率的变化规律,并且能正确反映出渗透率的变化范围。

    为增加低渗透含铀砂岩渗透率,创新性的将低频机械振动引入CO2+O2浸出,得到以下3点结论。

    (1)与前人得到的孔隙率方程相比,建立的孔隙率方程将物理振动损伤和化学损伤相结合,更能体现出CO2+O2浸出铀的过程中,低频机械振动和化学侵蚀对砂岩孔隙率的影响。

    (2)利用自制的振动浸出试验系统进行试验,发现低频振动明显提高了砂岩渗透率,渗透率增量最大值达3.265 md。

    (3)孔隙率方程及渗透率表达式能正确地反映浸出过程中,低频机械振动影响低渗透砂岩孔隙率、渗透率的变化规律及渗透率的变化范围。

  • 图  1   波浪和多孔弹性海床之间的动力相互作用

    Figure  1.   Dynamic interaction between ocean waves and poroelastic seabed

    图  2   线性波作用下本文解析解与Hsu等的解之间的对比

    Figure  2.   Comparison of seabed response between present solution and solution by Hsu et al. under linear waves

    图  3   二阶斯托克斯波作用下本文解析解与Zhang等的解之间的对比

    Figure  3.   Comparison of seabed response between present solution and solution by Zhang et al. under second-order Stokes waves

    图  4   二阶Stokes波作用下波浪周期T对海床动力响应的影响

    Figure  4.   Influences of wave period on dynamic response of seabed under second-order Stokes waves

    图  5   二阶Stokes波作用下水深d对海床动力响应的影响

    Figure  5.   Influences of water depth on dynamic response of seabed under second-order Stokes waves

    图  6   二阶Stokes波作用下波高H对海床动力响应的影响

    Figure  6.   Influences of seabed thickness on dynamic response of seabed under second-order Stokes waves

    图  7   二阶Stokes波作用下渗透系数K对海床动力响应的影响

    Figure  7.   Influences of permeability on dynamic response of seabed under second-order Stokes waves

    图  8   二阶Stokes波作用下剪切模量G对海床动力响应的影响

    Figure  8.   Influences of shear modulus on dynamic response of seabed under second-order Stokes waves

    表  1   多孔弹性海床和波浪的基本参数

    Table  1   Basic poroelastic properties and wave conditions

    波浪和海床土参数 数值
    波的周期T 10 s
    波高H 5 m
    水深d 0.125L*
    海床厚度h 20 m
    孔隙流体的实际体积弹性模量Kw 2×109 Pa
    孔隙流体密度ρf 1000 kg/m3
    孔隙率n 0.35
    泊松比μ 0.33
    渗透系数K 10-4 m/s
    剪切模量G 1×107 Pa
    海床土饱和度Sr 0.975
    *表中L代表波长。
    下载: 导出CSV
  • [1]

    LI X B, ZHANG Z Q, PAN E N. Wave-induced dynamic response in a transversely isotropic and multilayered poroelastic seabed[J]. Soil Dynamics and Earthquake Engineering, 2020, 139: 106365. doi: 10.1016/j.soildyn.2020.106365

    [2]

    JENG D S. Review of liquefaction around marine structures by B mutlu sumer[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2015, 141(5): 07515001. doi: 10.1061/(ASCE)WW.1943-5460.0000296

    [3]

    PUTNAM J A. Loss of wave energy due to percolation in a permeable sea bottom[J]. Transactions, American Geophysical Union, 1949, 30(3): 349. doi: 10.1029/TR030i003p00349

    [4]

    BIOT M A. Theory of propagation of elastic waves in a fluid-saturated porous solid: Ⅰ low-frequency range[J]. The Journal of the Acoustical Society of America, 1956, 28(2): 168-178. doi: 10.1121/1.1908239

    [5]

    SEKIGUCHI H, KITA K, OKAMOTO O. Response of poro-elastoplastic beds to standing waves[J]. Soils and Foundations, 1995, 35(3): 31-42. doi: 10.3208/sandf.35.31

    [6]

    BIOT M A. General theory of three-dimensional consolidation[J]. Journal of Applied Physics, 1941, 12(2): 155-164. doi: 10.1063/1.1712886

    [7]

    BIOT M A. Mechanics of deformation and acoustic propagation in porous media[J]. Journal of Applied Physics, 1962, 33(4): 1482-1498. doi: 10.1063/1.1728759

    [8]

    YAMAMOTO T, KONING H L, SELLMEIJER H, et al. On the response of a poro-elastic bed to water waves[J]. Journal of Fluid Mechanics, 1978, 87(1): 193-206. doi: 10.1017/S0022112078003006

    [9]

    MADSEN O S. Wave-induced pore pressures and effective stresses in a porous bed[J]. Géotechnique, 1978, 28(4): 377-393. doi: 10.1680/geot.1978.28.4.377

    [10]

    VERRUIJT A. Elastic storage of aquifers[M]// Flow through Porous Media. DE WIEST R J M. New York: Academic Press, 1969.

    [11]

    OKUSA S. Wave-induced stresses in unsaturated submarine sediments[J]. Géotechnique, 1985, 35(4): 517-532. doi: 10.1680/geot.1985.35.4.517

    [12]

    MEI C C, FODA M A. Wave-induced responses in a fluid-filled poro-elastic solid with a free surface: a boundary layer theory[J]. Geophysical Journal International, 1981, 66(3): 597-631. doi: 10.1111/j.1365-246X.1981.tb04892.x

    [13]

    HSU J R C, JENG D S. Wave-induced soil response in an unsaturated anisotropic seabed of finite thickness[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1994, 18(11): 785-807. doi: 10.1002/nag.1610181104

    [14]

    ZIENKIEWICZ O C, CHANG C T, BETTESS P. Drained, undrained, consolidating and dynamic behaviour assumptions in soils[J]. Géotechnique, 1980, 30(4): 385-395. doi: 10.1680/geot.1980.30.4.385

    [15]

    ULKER M B C, RAHMAN M S, JENG D S. Wave-induced response of seabed: various formulations and their applicability[J]. Applied Ocean Research, 2009, 31(1): 12-24. doi: 10.1016/j.apor.2009.03.003

    [16]

    LAUTON G, PATTIARATCHI C B, LENTINI C A D. Observations of breaking internal tides on the australian north west shelf edge[J]. Frontiers in Marine Science, 2021, 8: 629372. doi: 10.3389/fmars.2021.629372

    [17] 王忠涛, 栾茂田, 刘占阁, 等. 浅水区波浪非线性效应对砂质海床动力响应的影响[J]. 海洋工程, 2005, 23(1): 41-46. doi: 10.3969/j.issn.1005-9865.2005.01.007

    WANG Zhongtao, LUAN Maotian, LIU Zhange, et al. Studies on the effect of non-linear wave loading in shallow water region on dynamic response of sandy seabed[J]. The Ocean Engineering, 2005, 23(1): 41-46. (in Chinese) doi: 10.3969/j.issn.1005-9865.2005.01.007

    [18] 郭秀军, 朱大伟, 孟庆生, 等. 波浪作用下黄河口多层粉质土海床动力响应特征差异性分析[J]. 岩土工程学报, 2012, 34(12): 2270-2276. http://cge.nhri.cn/article/id/14954

    GUO Xiujun, ZHU Dawei, MENG Qingsheng, et al. Differences in dynamic response characteristics of multi-layer silty seabed under waves in Yellow River Estuary[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2270-2276. (in Chinese) http://cge.nhri.cn/article/id/14954

    [19]

    ZHOU X L, XU B, WANG J H, et al. An analytical solution for wave-induced seabed response in a multi-layered poro-elastic seabed[J]. Ocean Engineering, 2011, 38(1): 119-129. doi: 10.1016/j.oceaneng.2010.10.003

    [20]

    ZHANG Z Q, ZHOU B H, LI X B, et al. Second-order Stokes wave-induced dynamic response and instantaneous liquefaction in a transversely isotropic and multilayered poroelastic seabed[J]. Frontiers in Marine Science, 2022, 9: 1082337. doi: 10.3389/fmars.2022.1082337

    [21] 周晓智, 陈育民, 刘汉龙. 驻波作用下有限厚度海床动应力路径特性研究[J]. 岩土工程学报, 2018, 40(5): 890-899. doi: 10.11779/CJGE201805014

    ZHOU Xiaozhi, CHEN Yumin, LIU Hanlong. Study on characteristics of dynamic stress path of finite-thickness seabed under standing waves[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 890-899. (in Chinese) doi: 10.11779/CJGE201805014

    [22]

    LEMÉHAUTÉ B. An Introduction to Hydrodynamics and Water Waves[M]. New York: Springer-Verlag, 1976.

  • 期刊类型引用(5)

    1. 郭旭辉,朱鸿鹄,吴冰,高宇新,胡乐乐,曹鼎峰. 基于人工神经网络的黄土含水率光纤被动感测技术研究. 岩土力学. 2025(02): 653-664 . 百度学术
    2. 高磊,袁泽,王勤,高明军. 基于OFDR技术的土体含水率模型研究. 工程地质学报. 2025(01): 96-105 . 百度学术
    3. 严辉,林沛元. 深圳市岩溶地层标准贯入击数神经网络模型. 地质科技通报. 2025(02): 305-321 . 百度学术
    4. 李杰,朱鸿鹄,吴冰,刘喜凤,王家琛,曹鼎峰,施斌. 下蜀土降雨入渗光纤感测及渗透系数估算研究. 工程地质学报. 2024(02): 601-611 . 百度学术
    5. 朱鸿鹄. 工程地质界面:从多元表征到演化机理. 地质科技通报. 2023(01): 1-19 . 百度学术

    其他类型引用(0)

图(8)  /  表(1)
计量
  • 文章访问数:  618
  • HTML全文浏览量:  70
  • PDF下载量:  128
  • 被引次数: 5
出版历程
  • 收稿日期:  2023-06-04
  • 网络出版日期:  2024-03-24
  • 刊出日期:  2024-08-31

目录

/

返回文章
返回