ZHANG Hongxue, LIU Weiqun. Evolution Mechanism of Shale Gas Reservoirs Permeability under Thermal-Fluid-Solid Coupling[J]. Chinese Journal of High Pressure Physics, 2023, 37(3): 035303. doi: 10.11858/gywlxb.20230615
Citation: HUANG Min, ZHU Benhao, XIAO Gesheng, QIAO Li. Simulation on Deformation Damage and Strain Rate Effect of Nb3Sn Composite Superconductors under Cycling Load at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2024, 38(2): 024201. doi: 10.11858/gywlxb.20230755

Simulation on Deformation Damage and Strain Rate Effect of Nb3Sn Composite Superconductors under Cycling Load at Extreme Low Temperature

doi: 10.11858/gywlxb.20230755
  • Received Date: 17 Oct 2023
  • Rev Recd Date: 25 Nov 2023
  • Accepted Date: 25 Dec 2023
  • Available Online: 21 Jan 2024
  • Issue Publish Date: 09 Apr 2024
  • The study on damage and fracture of superconducting Nb3Sn under cyclic loading is an indispensable part of understanding the origin of the irreversible strain limit in Nb3Sn. This paper uses molecular dynamics simulation to investigate the fracture and deformation damage behavior of polycrystalline and single crystal Nb3Sn/Nb composite materials under cyclic loading at extremely low temperatures. The effects of strain rate on crack initiation and growth were carefully analyzed in both polycrystalline and single crystal Nb3Sn/Nb composite materials. The results indicate that slip occurs in single crystal Nb3Sn/Nb composite materials after cyclic loading. When the local stress at the slip band intersection exceeds the material strength, microcracks initiate at the slip band intersection, leading to fracture failure of the composite material. In contrast, the failure of polycrystalline Nb3Sn/Nb composite materials is due to the inability of stress at grain boundaries to relax under cyclic loading, which leads to the initiation of microcracks at the grain boundaries and intergranular fracture of the composite material. The analysis of the different damage, fracture, and failure mechanisms of polycrystalline and single crystal Nb3Sn/Nb composite materials at different strain rates shows that the fracture is brittle at low strain rates. As the strain rate rises, the number of slip bands in the single crystal Nb3Sn layer increases, enhancing the toughness of the single crystal Nb3Sn/Nb composite material. Conversely, the influence of grain boundaries on material strength decreases in polycrystalline materials as the strain rate increases. Moreover, polycrystalline Nb3Sn/Nb composite materials exhibit significant residual strength after local fracture of Nb3Sn at high strain rates. The research results will contribute to a better understanding of the damage evolution process of Nb3Sn/Nb composite materials under cyclic loading and offer theoretical guidance for optimizing material performance.

     

  • 我国页岩气可采资源量为2.508×1013 m3[1],四川盆地具有丰富的页岩气资源,已成为我国页岩气开发的主战场。四川盆地的页岩气储层埋深较大,其中龙马溪组的地层埋深普遍在4 000~4 500 m之间,埋深较浅处约为3 000~4 000 m[2]。由于地下的地温梯度和静水压力梯度分别为0.03 K/m和0.01 MPa/m[3],所以龙马溪组的地层温度范围为378~423 K(地表温度按288 K计算)。Gensterblum等[3]测定了甲烷在303、318、333和348 K(对应的储层深度分别为500、1 000、1 500和2 000 m)时的吸附等温线,发现甲烷的吸附能力与温度和孔隙压力有关。在不同温度(373~873 K)和不同压力(0.1~15.0 MPa)条件下,Geng等[4]对油页岩进行了热解反应实验发现,随着温度的升高,孔隙率逐渐增大,裂隙的最大张开度也随之增大;最大孔隙率为室温下的13.52倍,裂隙总数为813条,裂隙最大孔径为0.383 mm。对于页岩储层,裂隙是气体流动的主要通道,气体解吸引起的体应变对页岩气的开采有着显著的影响。因此,建立热-流-固耦合作用下页岩的渗透率模型对准确预测深部页岩气的渗流规律尤为重要。

    目前,国内外学者主要通过理论分析和室内实验相结合的方法研究温度对页岩渗透率的影响。利用岩石三轴渗流仪,赵瑜等[5]研究了孔压和温度对页岩渗流的影响机理。利用自制的煤岩稳态实验装置,Wang等[6]测试了我国无烟煤在不同气体压力和温度下的渗透率。吴迪等[7]利用三轴渗流实验装置,研究了孔隙压力、体积应力和温度对页岩中CH4气体渗流规律的影响。张道川等[8]采用多场耦合作用下不同相态CO2致裂驱替CH4实验装置,研究了有效应力、孔隙压力、温度以及吸附膨胀等因素对裂隙页岩渗透特性的影响。基于典型S&D模型,李波波等[9]建立了考虑温度和恒定外应力条件下的渗透率模型。Wang等[10]采用高温三轴渗透率实验设备研究了现场应力条件下油页岩的各向异性渗透率。Schwartz等[11]认为在恒定孔隙压力和温度变化下,页岩渗透率的演化规律可以用热解吸应变和页岩的热膨胀应变来描述,由此提出了热解吸渗透率增强模型。Teng等[12]提出了变温度下煤气相互作用的热敏感性渗透率模型,该模型能够解释热膨胀、热压裂、基质吸附和裂隙表面的热挥发对煤渗透率的影响。在围压10~25 MPa、温度298~343 K、注入压力6~12 MPa作用下,Ju等[13]测试了煤样的CO2渗透率,建立了表征裂隙煤CO2渗透率的力学模型。Sinha等[14]利用稳态渗透率测试设备研究了应力、压力、温度、吸附作用和气体类型对完整页岩岩样渗透率的影响。

    综上所述,国内外学者主要通过实验研究温度、围压、孔压和吸附等因素对页岩渗透率的影响,能够解释上述因素(特别是温度)对页岩渗透率影响的解析模型凤毛麟角。为此本研究将深部页岩储层简化为多孔弹性介质模型,考虑温度和吸附作用对渗透率的影响,推导页岩气开采过程中页岩储层的有效应力-渗透率模型,阐述围压、孔隙压力、有效应力、热解吸和温度对页岩渗透率的影响机理。然后,结合多孔介质弹性理论(考虑温度和解/吸附对变形的影响)建立页岩储层在单轴应变条件下的渗透率解析模型,利用实验数据验证该模型的有效性。最后,通过该模型分析不同温度下页岩渗透率的演化规律。

    页岩的孔隙率可以表示为

    ϕ=VpVb
    (1)

    式中:Vp为页岩孔隙的体积,m3Vb为页岩的体积,m3。对式(1)微分,获得

    dϕϕ=dεp+dεb
    (2)

    式中:dεp为孔隙的体应变增量,dεb为页岩的体应变增量。

    考虑有效应力、气体的解吸附和温度对页岩变形的影响,dεb和dεp分别为

    dεb=dVbVb=1K(dσαdp)dεs3αTdT
    (3)
    dεp=dVpVp=1Kpdσ(1Kp1Ks)dpdεs
    (4)

    式中:dVb和dVp分别为页岩的体积增量和孔隙的体积增量,m3KKpKs分别为页岩的体积模量、孔隙体积模量和固体颗粒的体积模量,Pa;dσ为平均法向应力增量,Pa;dp为孔隙压力增量,Pa;dεs为吸附引起的体应变增量;α为Biot系数,α=1−K/KsαT为页岩的线胀系数,K−1;dT为温度增量,K。

    联立式(2)、式(3)和式(4),得到

    dϕϕ=(1K1Kp)(dσdp)3αTdT
    (5)

    对式(5)积分

    ϕϕ0=exp{(1K1Kp)[(σσ0)(pp0)]3αT(TT0)}
    (6)

    式中:ϕ0为页岩的初始孔隙率;σσ0分别为平均法向应力和初始平均法向应力,Pa;pp0分别为孔隙压力和初始孔隙压力,Pa;TT0分别为页岩气储层的温度和初始温度,K。

    页岩气储层的渗透率和孔隙率的关系[15]

    kk0=(ϕϕ0)3
    (7)

    式中:kk0分别为页岩气储层的渗透率和初始渗透率,m2

    将式(6)代入式(7),考虑页岩的体积模量远大于其孔隙体积模量,得到考虑温度、吸附作用的有效应力-渗透率模型

    kk0=exp[3Kp(σeσe0)9αT(TT0)]
    (8)

    式中:σeσe0分别为有效应力和初始有效应力,单位Pa,σe=σpσe0=σ0p0

    该模型第1项为有效应力变化引起的体应变增量,第2项为温度(储层深度)变化引起的体应变增量。页岩气开采过程中储层渗透率的演化本质是储层的应变演化影响孔隙体积,该渗透率模型综合考虑了有效应力、温度对储层渗透率的影响。

    当不考虑温度对页岩气储层渗透率的影响时,式(8)可简化为文献[16-17]推导的煤的有效应力-渗透率模型

    kk0=exp[3Cp(σeσe0)]
    (9)

    式中:Cp为孔隙压缩系数,单位Pa−1Cp=1/Kp

    将页岩气储层简化为线弹性多孔介质,随着温度的升高,孔隙率和裂隙张开度均随之增大[4]。因此,考虑温度作用,规定压应力和压应变为正,页岩气储层的本构关系

    σij=λεkkδij+2Gεij+αpδij+Kεsδij+3KαTTδij
    (10)

    式中:σij为应力分量,单位Pa;λ为拉梅常数,单位Pa,λ=Eμ/[(1+μ)(1−2μ)];εkk为体应变;δij为Kronecker符号;G为剪切模量,单位Pa,G=E/[2(1+μ)];E为岩体的弹性模量,单位Pa;μ为泊松比;εij为应变分量;α为Biot系数,α=1−K/KsK为体积模量,单位Pa,K=E/[3(1−2μ)];Ks为固体颗粒体积模量,单位Pa。

    当考虑温度作用时,气体吸附/解吸引起的体应变利用下式计算[11]

    εs=εLb0(T/TT0T0)1/2exp[Q/QRT(RT)]p1+b0(T/TT0T0)1/2exp[Q/QRT(RT)]p
    (11)

    式中:εL为Langmuir体应变;b∞0为在无限大压力下的亲和常数,单位Pa−1b0=1pL1exp[Q/Q(RT0)(RT0)]pL为Langmuir压力,单位Pa;Q为吸附热,单位J/mol;R为普适气体常数,单位J/(mol·K)。

    联立式(8)、式(10)和式(11),得到页岩气储层在单轴应变条件下的热解吸渗透率模型

    kk0=exp{3Kp[1+μ3(1μ)(pp0)2E9(1μ)(εLb0(T/TT0T0)1/2exp(Q/QRT(RT))p1+b0(T/TT0T0)1/2exp(Q/QRT(RT))pεLb0exp(Q/QRT0(RT0))p01+b0exp(Q/QRT0(RT0))p0)]}exp{3Kp[2E3(1μ)+3Kp]αT(TT0)}
    (12)

    式中,Kp=ϕ0K[18]ϕ0=2b0/2b0a0a0[19-20]

    文献[11]测试了Marcellus裂隙页岩和完整页岩岩样的渗透率比随温度的变化,分别见表1表2。实验所用气体为CH4,实验过程中保持围压(24 MPa)和孔隙压力(6 MPa)不变。

    表  1  裂隙页岩的渗透率比
    Table  1.  Permeability ratio of fractured shale
    Temperature/Kk/k0 of fractured shale
    2951.000
    3031.044
    3131.113
    3231.211
    下载: 导出CSV 
    | 显示表格
    表  2  完整页岩的渗透率比
    Table  2.  Permeability ratio of intact shale
    Temperature/Kk/k0 of intact shaleTemperature/Kk/k0 of intact shale
    307.11.000332.71.095
    308.00.978333.31.045
    313.81.051337.71.058
    321.11.037337.91.061
    321.21.008338.01.084
    325.91.085
    下载: 导出CSV 
    | 显示表格

    将式(12)热解吸渗透率模型的理论解和不同温度下Marcellus裂隙页岩的实验数据进行对比,当模型中的参数T0=295 K[11]εL=1.7×10−3[18]b∞0=1.83×10−4 MPa−1[11]Q=1.7×104 J/mol[11]p=p0=6 MPa、R=8.31 J/(mol·K)[11]μ=0.35、E=32.75 GPa[21]Kp=24.26 MPa、αT=1.8×10−6 K−1时,该模型的理论解与Marcellus裂隙页岩的实验数据对比如图1所示。

    图  1  不同温度下热解吸模型解析解与裂隙页岩的实验数据对比
    Figure  1.  Comparison of analytical results of thermal-sorptive model and experimental data offractured shale under different temperatures

    图1中,当温度从295 K升高到323 K时,Marcellus裂隙页岩的渗透率增大了21.1%。由于黏土和有机质的解吸以及岩桥的热膨胀,裂隙页岩的渗透率比值随着温度的升高而逐渐增大。由此可见,该热解吸渗透率模型的理论解与Marcellus裂隙页岩的实验数据基本一致。

    将式(12)的热解吸渗透率模型理论解与不同温度下Marcellus完整页岩(含天然裂隙)的实验数据进行对比,当模型中T0=295 K[11]εL=1.7×10−3[18]b∞0=1.83×10−4 MPa−1[11]Q=1.7×104 J/mol[11]p=p0=6 MPa、R=8.31 J/(mol·K)[11]μ=0.35,E=32.75 GPa[21]Kp=43.53 MPa、αT=1.8×10−6 K−1时,该模型的理论解与Marcellus完整页岩的实验数据对比见图2

    图  2  不同温度下热解吸模型解析解与完整页岩的实验数据对比
    Figure  2.  Comparison of analytical results of thermal-sorptive model and experimental data of intact shaleunder different temperatures

    图2可知,随着温度的升高,完整页岩的渗透率比值呈反弹趋势,且整体逐渐增大,这是由于热解吸引起的孔隙膨胀量大于固体颗粒膨胀引起的孔隙闭合量。该渗透率模型的解析解能够与不同温度下Marcellus完整页岩的实验数据较好地吻合。

    由式(12)计算得到,页岩储层在单轴应变条件下的渗透率随温度和孔隙压力的演化规律如图3所示,储层初始温度为300 K,初始孔压为6 MPa。

    图  3  渗透率随温度和孔压的演化规律
    Figure  3.  Evolution of permeability withtemperature and pore pressure

    图3所示,当储层温度介于300~420 K范围内且保持不变时,随着孔隙压力的增大,渗透率比值先下降然后反弹,这是因为温度不变时,影响渗透率的因素主要为有效应力和气体吸附/解吸,随着孔隙压力的增大,有效应力减小,孔隙体积变大,导致渗透率增大;另一方面,孔隙压力增大,气体吸附量增大,孔隙体积减小,导致渗透率减小。孔隙压力增大前期,有效应力引起的渗透率增大量小于气体吸附引起的渗透率减小量,导致渗透率下降;孔隙压力增大后期,有效应力引起的渗透率增大量大于气体吸附引起的渗透率减小量,从而导致渗透率反弹。温度越高,渗透率比值下降越平缓,反弹越急剧。

    当储层温度介于420~440 K范围内,且保持不变时,随着孔压的增大,渗透率逐渐增大。这是因为温度较高时,随着孔隙压力的增大,有效应力引起的渗透率增大量始终大于气体吸附引起的渗透率减小量,导致渗透率逐渐增大。

    图3还可以看出,当孔隙压力保持不变时,随着温度升高,渗透率比值先增大然后逐渐减小;孔隙压力越大,渗透率比值开始下降时的温度越大。这是因为当孔隙压力不变时,影响渗透率的主要因素为气体热解吸/热吸附和岩体热膨胀,随着温度升高,气体热解吸量增大,孔隙体积增大,导致渗透率增大;另一方面,随着温度升高,页岩固体颗粒膨胀,孔隙体积缩小,导致渗透率减小。温度升高前期,热解吸引起的渗透率增大量大于热膨胀引起的渗透率减小量,从而导致渗透率增大;而温度升高后期,热解吸引起的渗透率增大量小于热膨胀引起的渗透率减小量,导致渗透率减小。

    不同孔隙压力作用下,页岩气储层渗透率随温度变化的演化规律如图4所示。当孔隙压力保持不变(5、15 和25 MPa)时,随着温度升高,渗透率比值先增大,然后逐渐减小,孔隙压力越小,热解吸引起的渗透率增大的速率越大,导致渗透率比值开始减小时的温度越小;当渗透率比值开始减小以后,孔隙压力越小,热膨胀引起的渗透率减小的速率越大,从而导致渗透率比值减小的越急剧。

    图  4  不同孔压下渗透率随温度的演化规律
    Figure  4.  Evolution of permeability with temperatureunder different pore pressures

    当孔隙压力保持35和45 MPa不变时,随着温度升高,渗透率比值轻微下降,然后反弹,最后逐渐下降。孔隙压力较小时,反弹现象较为明显。从图4可以看出,当温度在300~440 K范围内变化,孔隙压力为5 MPa时,渗透率比值变化范围为0.7~1.2。而孔隙压力为5和40 MPa时,渗透率比值在0.9~1.1范围内波动。由此可见,孔隙压力越大,温度对储层渗透率的影响越小。

    不同温度作用下,页岩气储层渗透率随孔压的演化规律如图5所示。当储层温度保持不变(318、348、378和408 K,对应的储层深度分别为1000、2000、3000和4000 m)时,随着孔隙压力的减小,渗透率比值先减小然后逐渐增大,其演化规律呈反弹现象。储层温度(储层深度)越大,渗透率的反弹现象越不明显,当储层深度达到3000 m时,渗透率基本不反弹。

    图  5  不同温度下渗透率随孔压的演化规律
    Figure  5.  Evolution of permeability with porepressure under different temperatures

    当储层温度均保持不变时,影响储层渗透率的主要因素为有效应力和孔压降低引起的气体解吸。当孔隙压力下降时,一方面吸附于有机质颗粒表面和孔隙表面的气体解吸,孔隙体积和裂隙张开度增大,导致储层渗透率增大;另一方面,孔隙压力下降导致储层有效应力增大,孔隙和裂隙被压缩,导致储层渗透率减小。当气体解吸引起的基质收缩量小于有效应力增大引起的孔隙和裂隙的压缩量时,渗透率下降;当气体解吸引起的基质收缩量大于有效应力增大引起的孔隙和裂隙压缩量时,渗透率增大,这也正是恒温作用下储层渗透率随着孔隙压力的下降呈现反弹现象的原因。

    在式(12)所示的热解吸渗透率模型中,泊松比、孔隙体积模量、热膨胀系数和Langmuir体应变对储层渗透率的影响较为显著,下面将分别讨论恒压、恒温作用下上述参数对储层渗透率的影响。

    当孔隙压力保持10 MPa不变时,泊松比对储层渗透率的影响如图6所示。可以看出,泊松比不变时,随着温度升高,储层渗透率比值先增大然后减小,呈现“倒U形”演化规律。泊松比越大,渗透率比值的梯度越大。

    图  6  恒压下泊松比对渗透率的影响
    Figure  6.  Effect of Poisson’s ratio on permeabilityunder constant pore pressure

    究其原因,储层孔隙压力不变,热解吸/吸附和热膨胀是影响页岩储层渗透率的主要因素。当温度升高时,页岩岩体膨胀导致储层孔隙和裂隙体积压缩,使得储层渗透率下降;另外,由于温度升高,页岩的吸附能力逐渐降低[22],从而导致储层孔隙体积逐渐增大,裂隙张开度随之增大[4],储层渗透率也随之增大。综上所述,热解吸和热膨胀对页岩储层渗透率的影响是相反的。温度升高初期,热解吸引起的渗透率变化量大于热膨胀引起的渗透率变化量;温度升高后期,热解吸引起的渗透率变化量小于热膨胀引起的渗透率变化量,从而导致恒压下渗透率随温度呈现“倒U形”变化。泊松比越大,热解吸引起的渗透率变化量越大,热膨胀引起的渗透率变化量越小,因此渗透率比值梯度随之增大。

    当孔隙压力保持10 MPa不变时,孔隙体积模量对储层渗透率的影响见图7。如图7所示,当孔隙压力恒定时,储层渗透率比值随温度升高仍然呈现“倒U形”的演化规律,而且随着孔隙体积模量的增大,储层渗透率比值梯度减小。这是因为孔隙体积模量越大,热膨胀引起的渗透率变化量越大,而热解吸引起的渗透率变化量越小,从而导致渗透率比值梯度随孔隙体积模量的增大而减小。

    图  7  恒压下孔隙体积模量对渗透率的影响
    Figure  7.  Effect of pore volume modulus onpermeability under constant pore pressure

    当储层的孔隙压力不变时,页岩的线胀系数对储层渗透率的影响见图8。如图8所示,当线胀系数小于临界值(渗透率随温度的演化规律不再呈现“倒U形”时的线胀系数)时,渗透率比值随温度的演化规律呈现“倒U形”现象,可通过式(12)计算出该临界值;当线胀系数大于该临界值时,渗透率比值随温度的升高而逐渐下降,不再呈现“倒U形”。在温升初期,线胀系数越大,渗透率比值梯度越大;在温升后期,渗透率比值梯度随线胀系数的增大而减小。这是因为孔隙压力不变时,只有气体的热解吸和岩体的热膨胀影响储层渗透率,线胀系数越大,页岩岩体因温度升高而引起的体积膨胀量越大,挤压的孔隙体积和裂隙体积越大,渗透率下降的越快。当热膨胀引起的渗透率变化量总是大于热解吸引起的渗透率变化量时,渗透率比值随温度的升高呈现逐渐下降的趋势。

    图  8  恒压下线胀系数对渗透率的影响
    Figure  8.  Evolution of coefficient of linear expansionon permeability under constant pore pressure

    恒压作用下朗缪尔体应变对储层渗透率的影响见图9。如图9所示,当朗缪尔体应变大于临界值(渗透率随温度的演化规律不呈现“倒U形”时的朗缪尔体应变)时,储层渗透率随温度的变化曲线呈“倒U形”趋势,可以通过式(12)计算出此临界值;当朗缪尔体应变小于临界值时,渗透率随着温度的升高而逐渐降低。在温升初期,朗缪尔体应变越小,渗透率比值梯度越大;而在温升后期,朗缪尔体应变越小,渗透率比值梯度越小。这是因为朗缪尔体应变发生变化时,热解吸引起的体应变随之变化,热膨胀引起的体应变不发生变化,因此当升温过程中热解吸引起的渗透率变化量始终小于热膨胀引起的渗透率变化量时,渗透率比值随温度的升高逐渐下降。朗缪尔体应变越大,温升初期热解吸引起的渗透率变化量越大,此渗透率比值梯度越小,温升后期热解吸引起的渗透率变化量越小,渗透率梯度比值越大。

    图  9  恒压下Langmuir体应变对渗透率的影响
    Figure  9.  Evolution of Langmuir volume strain onpermeability under constant pore pressure

    当储层温度不变时,泊松比对储层渗透率的影响见图10。如图10所示,恒温下储层渗透率随孔隙压力的下降,先逐渐降低然后反弹逐渐增大。泊松比越大,有效应力变化引起的孔隙和裂隙的体积变化越大,从而导致渗透率变化量越大;泊松比越大,气体解吸引起的孔隙和裂隙的体积变化越小,进而导致渗透率变化量越小。因此,泊松比越大,渗透率比值梯度越大。

    图  10  恒温下泊松比对渗透率的影响
    Figure  10.  Effect of Poisson’s ratio on permeabilityunder constant temperature

    当储层温度不变时,孔隙体积模量对储层渗透率的影响见图11。如图11所示,当储层温度保持不变时,渗透率比值随孔隙压力的减小均呈现先变小再反弹增大的趋势。孔隙体积模量越大,孔隙抵抗变形的能力越强,相同的有效应力变化量引起的孔隙体积变化量越小,所以有效应力变化引起的渗透率变化量越小。另外,孔隙体积模量越大,气体解吸引起的渗透率变化量越小。因此,孔隙体积模量越大,渗透率比值梯度越小。

    图  11  恒温下孔隙体积模量对渗透率的影响
    Figure  11.  Effect of pore volume modulus onpermeability under constant temperature

    当储层温度不变时,朗缪尔体应变对渗透率的影响见图12。如图12所示,朗缪尔体应变越大,渗透率随孔隙压力变化呈现的先变小再反弹现象越明显,随着朗缪尔体应变的减小,反弹现象逐渐消失。这是因为有效应力引起的渗透率变化量与朗缪尔体应变无关,而气体解吸引起的渗透率变化量取决于朗缪尔体应变。朗缪尔体应变越小,气体解吸引起的渗透率变化量越小,从而导致在孔隙压力的下降过程中有效应力变化引起的渗透率变化量(减小)始终大于气体解吸引起的渗透率变化量(增大)。

    图  12  恒温下Langmuir体应变对渗透率的影响
    Figure  12.  Effect of Langmuir volume strain onpermeability under constant temperature

    (1) 考虑温度、气体压力、有效应力和热解吸对储层渗透率的影响,推导了页岩气储层的有效应力-渗透率模型,该模型能解释有效应力不变时渗透率随孔隙压力的演化机制。

    (2) 建立了储层在单轴应变条件下的热解吸渗透率模型,该模型能准确地阐述Marcellus裂隙页岩和完整页岩岩样的渗透率随温度的演化机制。

    (3) 热解吸/吸附、热膨胀和孔隙压力是影响页岩气储层渗透率的主要因素。当孔隙压力恒定时,渗透率随温度呈“倒U形”演化规律。孔隙压力越大,温度对渗透率的影响越小;当温度恒定时,渗透率随孔压变化呈“U形”演化规律,温度(对应储层深度)越高,渗透率随孔压下降的反弹现象越不明显。

    (4) 储层孔压恒定时,泊松比越大,渗透率比值梯度越大;孔隙体积模量越大,渗透率比值梯度越小;当线胀系数大于临界值时,渗透率比值随温度的演化规律不呈现“倒U形”;当朗缪尔体应变小于临界值时,渗透率比值随温度的演化规律也不呈现“倒U形”。

    (5) 储层温度恒定时,泊松比越大,渗透率比值梯度越大;孔隙体积模量越大,渗透率比值梯度越小。当朗缪尔体应变小于临界值时,渗透率比值随孔压的变化不呈现“U形”演化规律。

  • [1]
    ZHANG R, GAO P F, WANG X Z. Strain dependence of critical superconducting properties of Nb3Sn with different intrinsic strains based on a semi-phenomenological approach [J]. Cryogenics, 2017, 86: 30–37. doi: 10.1016/j.cryogenics.2017.07.007
    [2]
    DEVRED A, BACKBIER I, BESSETTE D, et al. Status of ITER conductor development and production [J]. IEEE Transactions on Applied Superconductivity, 2012, 22(3): 4804909. doi: 10.1109/TASC.2012.2182980
    [3]
    KIYOSHI T, MATSUMOTO S, KOSUGE M, et al. Superconducting inserts in high-field solenoids [J]. IEEE Transactions on Applied Superconductivity, 2002, 12(1): 470–475. doi: 10.1109/TASC.2002.1018445
    [4]
    DEVRED A, BACKBIER I, BESSETTE D, et al. Challenges and status of ITER conductor production [J]. Superconductor Science and Technology, 2014, 27(4): 044001. doi: 10.1088/0953-2048/27/4/044001
    [5]
    BOTTURA L, DE RIJK G, ROSSI L, et al. Advanced accelerator magnets for upgrading the LHC [J]. IEEE Transactions on Applied Superconductivity, 2012, 22(3): 4002008. doi: 10.1109/TASC.2012.2186109
    [6]
    EISENSTATT L R, WRIGHT R N. Plastic deformation in polycrystalline Nb3Sn [J]. Metallurgical Transactions A, 1980, 11(7): 1131–1138. doi: 10.1007/BF02668137
    [7]
    CLARK J B, WRIGHT R N. Laboratory extrusion of Nb3Sn [J]. Metallurgical Transactions A, 1983, 14(11): 2295–2299. doi: 10.1007/BF02663304
    [8]
    OCHIAI S, OSAMURA K. Influence of cyclic loading at room temperature on the critical current at 4.2 K of Nb3Sn superconducting composite wire [J]. Cryogenics, 1992, 32(6): 584–590. doi: 10.1016/0011-2275(92)90045-C
    [9]
    WESSEL W A J, NIJHUIS A, ILYIN Y, et al. A novel “test arrangement for strain influence on strands” (TARSIS): mechanical and electrical testing of ITER Nb3Sn strands [J]. AIP Conference Proceedings, 2004, 711(1): 466–473.
    [10]
    NIJHUIS A, VAN DEN EIJNDEN N C, ILYIN Y, et al. Impact of spatial periodic bending and load cycling on the critical current of a Nb3Sn strand [J]. Superconductor Science and Technology, 2005, 18(12): S273–S283. doi: 10.1088/0953-2048/18/12/009
    [11]
    VAN DEN EIJNDEN N C, NIJHUIS A, ILYIN Y, et al. Axial tensile stress-strain characterization of ITER model coil type Nb3Sn strands in TARSIS [J]. Superconductor Science and Technology, 2005, 18(11): 1523–1532. doi: 10.1088/0953-2048/18/11/020
    [12]
    NIJHUIS A, MIYOSHI Y, JEWELL M C, et al. Systematic study on filament fracture distribution in ITER Nb3Sn strands [J]. IEEE Transactions on Applied Superconductivity, 2009, 19(3): 2628–2632. doi: 10.1109/TASC.2009.2018082
    [13]
    MITCHELL N. Analysis of the effect of Nb3Sn strand bending on CICC superconductor performance [J]. Cryogenics, 2002, 42(5): 311–325. doi: 10.1016/S0011-2275(02)00041-3
    [14]
    SHETH M K, LEE P, MCRAE D M, et al. Procedures for evaluating filament cracking during fatigue testing of Nb3Sn strand [J]. AIP Conference Proceedings, 2012, 1435(1): 201–208.
    [15]
    SHEN F Z, ZHANG H C, HUANG C J, et al. Experimental study on strain sensitivity of internal-tin Nb3Sn superconducting strand based on non-destructive technology [J]. Physica C: Superconductivity and Its Applications, 2021, 584: 1353784. doi: 10.1016/j.physc.2020.1353784
    [16]
    JIANG L, SU X Y, SHEN L Y, et al. Damage behavior of Nb3Sn/Cu superconducting strand at room temperature under asymmetric strain cycling [J]. Fusion Engineering and Design, 2021, 172: 112869. doi: 10.1016/j.fusengdes.2021.112869
    [17]
    TABIN J, SKOCZEŃ B, BIELSKI J. Discontinuous plastic flow in superconducting multifilament composites [J]. International Journal of Solids and Structures, 2020, 202: 12–27. doi: 10.1016/j.ijsolstr.2020.05.033
    [18]
    QIULI SUN E. Multi-scale nonlinear stress analysis of Nb3Sn superconducting accelerator magnets [J]. Superconductor Science and Technology, 2022, 35(4): 045019. doi: 10.1088/1361-6668/ac5a11
    [19]
    MITCHELL N. Finite element simulations of elasto-plastic processes in Nb3Sn strands [J]. Cryogenics, 2005, 45(7): 501–515. doi: 10.1016/j.cryogenics.2005.06.003
    [20]
    WANG X, LI Y X, GAO Y W. Mechanical behaviors of multi-filament twist superconducting strand under tensile and cyclic loading [J]. Cryogenics, 2016, 73: 14–24. doi: 10.1016/j.cryogenics.2015.11.002
    [21]
    JIANG L, ZHANG X Y, ZHOU Y H. Nonlinear static and dynamic mechanical behaviors of Nb3Sn superconducting composite wire: experiment and analysis [J]. Acta Mechanica Sinica, 2023, 39(3): 122322. doi: 10.1007/s10409-022-22322-x
    [22]
    LEE J, POSEN S, MAO Z G, et al. Atomic-scale analyses of Nb3Sn on Nb prepared by vapor diffusion for superconducting radiofrequency cavity applications: a correlative study [J]. Superconductor Science and Technology, 2019, 32(2): 024001. doi: 10.1088/1361-6668/aaf268
    [23]
    ZHANG Y, ASHCRAFT R, MENDELEV M I, et al. Experimental and molecular dynamics simulation study of structure of liquid and amorphous Ni62Nb38 alloy [J]. The Journal of Chemical Physics, 2016, 145(20): 204505. doi: 10.1063/1.4968212
    [24]
    KO W S, KIM D H, KWON Y J, et al. Atomistic simulations of pure tin based on a new modified embedded-atom method interatomic potential [J]. Metals, 2018, 8(11): 900. doi: 10.3390/met8110900
    [25]
    CHUDINOV V G, GOGOLIN V P, GOSHCHITSKII B N, et al. Simulation of collision cascades in intermetallic Nb3Sn compounds [J]. Physica Status Solidi (A), 1981, 67(1): 61–67. doi: 10.1002/pssa.2210670103
    [26]
    LUTSKO J F. Stress and elastic constants in anisotropic solids: molecular dynamics techniques [J]. Journal of Applied Physics, 1988, 64(3): 1152–1154. doi: 10.1063/1.341877
    [27]
    PENG Y X, KNIGHT C, BLOOD P, et al. Extending parallel scalability of LAMMPS and multiscale reactive molecular simulations [C]//Proceedings of the 1st Conference of the Extreme Science and Engineering Discovery Environment: Bridging from the eXtreme to the Campus and Beyond. Chicago: ACM, 2012: 37.
    [28]
    PAPADIMITRIOU I, UTTON C, TSAKIROPOULOS P. Ab initio investigation of the intermetallics in the Nb-Sn binary system [J]. Acta Materialia, 2015, 86: 23–33. doi: 10.1016/j.actamat.2014.12.017
    [29]
    SUNDARESWARI M, RAMASUBRAMANIAN S, RAJAGOPALAN M. Elastic and thermodynamical properties of A15 Nb3X (X = Al, Ga, In, Sn and Sb) compounds-first principles DFT study [J]. Solid State Communications, 2010, 150(41/42): 2057–2060.
    [30]
    ZHANG R, GAO P F, WANG X Z, et al. First-principles study on elastic and superconducting properties of Nb3Sn and Nb3Al under hydrostatic pressure [J]. AIP Advances, 2015, 5(10): 107233. doi: 10.1063/1.4935099
    [31]
    DE MARZI G, MORICI L, MUZZI L, et al. Strain sensitivity and superconducting properties of Nb3Sn from first principles calculations [J]. Journal of Physics: Condensed Matter, 2013, 25(13): 135702. doi: 10.1088/0953-8984/25/13/135702
    [32]
    ESHELBY J D. Dynamical theory of crystal lattices: Max Born and Kun Huang: Oxford University Press, 1954, pp. xii+420. 50s. net. [J]. Journal of the Mechanics and Physics of Solids, 1955, 3(3): 231.
    [33]
    SCHEUERLEIN C, DI MICHIEL M, BUTA F, et al. Stress distribution and lattice distortions in Nb3Sn multifilament wires under uniaxial tensile loading at 4.2 K [J]. Superconductor Science and Technology, 2014, 27(4): 044021. doi: 10.1088/0953-2048/27/4/044021
    [34]
    时海芳, 任鑫. 材料力学性能 [M]. 2版. 北京: 北京大学出版社, 2015: 181–183.

    SHI H F, REN X. Mechanical properties of materials [M]. 2nd ed. Beijing: Peking University Press, 2015: 181–183.
    [35]
    LIANG L W, WANG Y J, CHEN Y, et al. Dislocation nucleation and evolution at the ferrite-cementite interface under cyclic loadings [J]. Acta Materialia, 2020, 186: 267–277. doi: 10.1016/j.actamat.2019.12.052
    [36]
    ESSMANN U, MUGHRABI H. Annihilation of dislocations during tensile and cyclic deformation and limits of dislocation densities [J]. Philosophical Magazine A, 1979, 40(6): 731–756. doi: 10.1080/01418617908234871
    [37]
    ZOTOV N, GRABOWSKI B. Molecular dynamics simulations of screw dislocation mobility in bcc Nb [J]. Modelling and Simulation in Materials Science and Engineering, 2021, 29(8): 085007. doi: 10.1088/1361-651X/ac2b02
    [38]
    SPARKS G, MAAß R. Shapes and velocity relaxation of dislocation avalanches in Au and Nb microcrystals [J]. Acta Materialia, 2018, 152: 86–95. doi: 10.1016/j.actamat.2018.04.007
    [39]
    PADILLA II H A, BOYCE B L. A review of fatigue behavior in nanocrystalline metals [J]. Experimental Mechanics, 2010, 50(1): 5–23. doi: 10.1007/s11340-009-9301-2
    [40]
    KUMAR K S, SURESH S, CHISHOLM M F, et al. Deformation of electrodeposited nanocrystalline nickel [J]. Acta Materialia, 2003, 51(2): 387–405. doi: 10.1016/S1359-6454(02)00421-4
    [41]
    OVID’KO I A, SHEINERMAN A G. Triple junction nanocracks in fatigued nanocrystalline materials [J]. Reviews on Advanced Materials Science, 2004, 7(1): 61–66.
    [42]
    OVID’KO I A, SHEINERMAN A G. Grain size effect on crack blunting in nanocrystalline materials [J]. Scripta Materialia, 2009, 60(8): 627–630. doi: 10.1016/j.scriptamat.2008.12.028
    [43]
    SUENAGA M, JANSEN W. Chemical compositions at and near the grain boundaries in bronze-processed superconducting Nb3Sn [J]. Applied Physics Letters, 1983, 43(8): 791–793. doi: 10.1063/1.94457
    [44]
    OCHIAI S, OSAMURA K, UEHARA T. Grain size and its relation to tensile strength of Nb3Sn compound in bronze-processed multi-filamentary superconducting materials [J]. Journal of Materials Science, 1987, 22(6): 2163–2168. doi: 10.1007/BF01132954
    [45]
    OH S H, JEONG Y J, NA S H, et al. Atomic behavior of Ti in A15 Nb3Sn and its effects on diffusional growth of Nb3Sn layer [J]. Journal of Alloys and Compounds, 2023, 957: 170438. doi: 10.1016/j.jallcom.2023.170438
    [46]
    OCHIAI S, OSAMURA K, MAEKAWA M. Comparison of mechanical and superconducting properties of titanium-added Nb3Sn composite wire with those of non-added ones [J]. Superconductor Science and Technology, 1991, 4(6): 262–269. doi: 10.1088/0953-2048/4/6/009
  • Relative Articles

    [1]YANG Xiangyang, WU Dun, ZHU Youlin, LI Junguo, ZHANG Ruizhi, ZHANG Jian, LUO Guoqiang. Molecular Dynamics Simulation Study on Spallation Failure of [100] Single Crystal Aluminum under Different Waveform Loadings[J]. Chinese Journal of High Pressure Physics, 2024, 38(3): 030106. doi: 10.11858/gywlxb.20240786
    [2]CHANG Lijun, HUANG Xingyuan, YUAN Shenglin, CAI Zhihua. Mechanical Properties and Failure Analysis of UHMWPE Fiber Composite Laminates under Compressive Loading[J]. Chinese Journal of High Pressure Physics, 2023, 37(1): 014102. doi: 10.11858/gywlxb.20220633
    [3]NIE Feiqing, MA Ruiqiang, LI Zhiqiang. Compressive Properties of Ice Containing Cotton at Low Strain Rates[J]. Chinese Journal of High Pressure Physics, 2023, 37(3): 034104. doi: 10.11858/gywlxb.20230608
    [4]WANG Haoyang, WEI Ying, QIAO Li. Simulation of Dynamic Crack Propagation in Superconducting Nb3Sn at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884
    [5]WANG Jianan, WU Bao, HE Anmin, WU Fengchao, WANG Pei, WU Heng’an. Research Progress on Dynamic Damage and Failure of Metal Materials under Shock Loading with Molecular Dynamics Simulation[J]. Chinese Journal of High Pressure Physics, 2021, 35(4): 040101. doi: 10.11858/gywlxb.20210747
    [6]WEN Yanbo, HUANG Ruiyuan, LI Ping, MA Jian, XIAO Kaitao. Damage Evolution Equation of Concrete Materials at High Temperatures and High Strain Rates[J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 024103. doi: 10.11858/gywlxb.20200617
    [7]SHI Zhentian, YANG Xujia, WANG Haoyang, QIAO Li. Superconducting Transition of Nb3Sn Single Crystal under High-Pressure[J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 021102. doi: 10.11858/gywlxb.20200615
    [8]WANG Yonghuan, XU Peng, FAN Zhiqiang, WANG Zhuangzhuang. Mechanical Characteristics and Quasi-Static Compression Deformation Mechanism of Open-Cell Aluminum Foam with Spherical Cells[J]. Chinese Journal of High Pressure Physics, 2019, 33(1): 014201. doi: 10.11858/gywlxb.20180532
    [9]SUN Xiaobo, GAO Yubo, XU Peng. Failure and Fracture Characteristics of Al2O3 Ceramics under Impact Loading[J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 054202. doi: 10.11858/gywlxb.20180695
    [10]WU Qijian, ZHI Xudong. Strain Rate Effect of GFRP-Reinforced Circular Steel Tube under Low-Velocity Impact[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 044203. doi: 10.11858/gywlxb.20180653
    [11]GAO Guang-Fa. Hardening Effect of the Strain Rate on the Dynamic Tensile Strength of the Plain Concrete[J]. Chinese Journal of High Pressure Physics, 2017, 31(5): 593-602. doi: 10.11858/gywlxb.2017.05.013
    [12]GAO Guang-Fa. Effect of Strain-Rate Hardening on Dynamic Compressive Strength of Plain Concrete[J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 261-270. doi: 10.11858/gywlxb.2017.03.007
    [13]PANG Bao-Jun, YANG Zhen-Qi, WANG Li-Wen, CHI Run-Qiang. Dynamic Compression Properties and Constitutive Model with Strain Rate Effect of Rubber Material[J]. Chinese Journal of High Pressure Physics, 2011, 25(5): 407-415 . doi: 10.11858/gywlxb.2011.05.005
    [14]WANG Wei-Bin, YIN Jian-Zhong. Molecular Dynamics Simulation of CO2-Expanded Liquids[J]. Chinese Journal of High Pressure Physics, 2010, 24(4): 271-278 . doi: 10.11858/gywlxb.2010.04.006
    [15]QIN Kun, YANG Li-Ming, HU Shi-Sheng. Study on Strain Rate Sensitivities of FCC Metals Based on Multi-Scale Analysis[J]. Chinese Journal of High Pressure Physics, 2009, 23(3): 196-202 . doi: 10.11858/gywlxb.2009.03.006
    [16]LIU Yao-Dong, YU Ji-Lin, ZHENG Zhi-Jun. Effect of Inertia on the Dynamic Behavior of Cellular Metal[J]. Chinese Journal of High Pressure Physics, 2008, 22(2): 118-124 . doi: 10.11858/gywlxb.2008.02.002
    [17]LIU Hong, WEI Dong-Qing, ZHAO Ji-Jun, GUO Yong-Xin, GONG Zi-Zheng. Molecular Dynamic Study of Liquid Nitromethane under High-Pressure[J]. Chinese Journal of High Pressure Physics, 2004, 18(4): 319-327 . doi: 10.11858/gywlxb.2004.04.006
    [18]ZHANG Dai-Yu, LIU Fu-Sheng, LI Xi-Jun, JING Fu-Qian. Shock Temperatures of Porous Iron from MD Simulations[J]. Chinese Journal of High Pressure Physics, 2003, 17(1): 16-21 . doi: 10.11858/gywlxb.2003.01.003
    [19]CAI Ling-Cang, CHEN Qi-Feng, JING Fu-Qian, CHEN Dong-Quan. Study of Equation-of-State of Dense Helum by Molecular Dynamics Simulation[J]. Chinese Journal of High Pressure Physics, 1999, 13(4): 255-261 . doi: 10.11858/gywlxb.1999.04.003
    [20]YU Wan-Rui. Molecular Dynamics Studies of Material Behavior at High Strain Rates[J]. Chinese Journal of High Pressure Physics, 1989, 3(2): 143-147 . doi: 10.11858/gywlxb.1989.02.006
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(27)  / Tables(1)

    Article Metrics

    Article views(198) PDF downloads(37) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return