Competitive adsorption characteristics and diffusion process of CO2−CH4 on mineral surface: A case study of the 2nd Section of Shanxi Formationin Yan 'an Gas Field
-
摘要:研究目的
CO2−EGR技术能提高天然气采收率，同时将CO2永久封存于地下，助力实现碳中和目标。CO2−CH4在纳米孔隙中竞争吸附和扩散是增采和封存的关键机理。
研究方法本文以鄂尔多斯盆地延安气田山2储层为靶区，利用分子动力学（MD）和巨正则系统蒙特卡罗（GCMC）方法建立模型，研究储层温压条件下CO2−CH4混合气体在关键矿物（石英和伊利石）纳米基质孔隙中的竞争吸附规律，分析CO2自扩散系数与温压关系。
研究结果（1）在等温（353.15 K）变压（5.9~17.7 MPa）和等压（11.8 MPa）变温（313.15~373.15 K）条件下，石英和伊利石对CO2吸附能力大于CH4，CO2−CH4在石英孔隙中的竞争吸附选择性大于伊利石孔隙；（2）等温（353.15 K）变压（5.9~17.7 MPa）和等压（11.8 MPa）变温（313.15~373.15 K）条件下, CO2−CH4在石英和伊利石孔隙中的竞争吸附选择性分别随压力、温度的增加而降低；（3）低压（5.9 MPa）高温（373.15 K）条件下，CO2在CO2−CH4−石英和CO2−CH4−伊利石系统的流动和扩散效率更高。
结论石英和伊利石对CO2吸附量更高，置换CH4能力更高，CO2封存效果更好。
创新点:阐明了纳米孔隙中CO2−CH4的竞争吸附和扩散规律；指出了延安气田山西组2段储层条件下CO2自扩散特性。
Abstract:This paper is the result of environmental geological survey engineering.
ObjectiveThe CO2−Enhanced Gas Recovery (CO2−EGR) technology significantly augments natural gas extraction efficiency while concurrently facilitating the permanent subsurface sequestration of CO2. This dual benefit substantially aids in achieving carbon neutrality goals. The mechanisms pivotal to enhanced recovery and storage include the competitive adsorption and diffusion of CO2−CH4 within nanopores.
MethodsThis study focuses on the 2ndsection of Shanxi formation in the Yan'an Gas Field located in the Ordos Basin. Using molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) methods, a model was established toinvestigate the competitive adsorption behaviors of CO2−CH4 mixed gases in the nanoporous matrices of key minerals, specifically quartz and illite, under reservoir−specific temperature and pressure conditions. Additionally, the study analyzes the correlation between the self−diffusion coefficient of CO2 and the variabilities in temperature and pressure.
ResultsThe study yields several findings: (1) At an isothermal condition of 353.15 K and varying pressures from 5.9 to 17.7 MPa, both quartz and illite exhibit heightened adsorptive capacities for CO2 in comparison to CH4. Furthermore, the competitive adsorption selectivity for CO2−CH4 is found to be greater in quartz pores than in illite pores. (2) Under similar isothermal conditions and at a constant pressure of 11.8 MPa with temperatures ranging from 313.15 K to 373.15 K, the competitive adsorption selectivity for CO2−CH4 in both quartz and illite pores is observed to diminish with increasing pressure and temperature. (3) Under conditions of low pressure (5.9 MPa) and high temperature (373.15 K), there is an enhancement in the mobility and diffusion efficiency of CO2 within both CO2−CH4−quartz and CO2−CH4−illite systems.
ConclusionQuartz and illite have higher CO2 adsorption capacity, greaterCH4displacement capacity, and better CO2 storage effect.
Highlights:The competitive adsorption and diffusion behavior of CO2−CH4 in nanopores is clarified, and the self-diffusion characteristics of CO2 under the reservoir conditions of the 2nd Section of the Shanxi Formation in the Yan’an Gas Field is identified.
-
1. 引 言
流速、流向及悬浮沉积物浓度对海洋物理、生物、化学要素的循环具有重要影响，是研究浅海及河口地区泥沙输运、物质运移的基本要素。同步获取研究区全剖面的流速、流向及悬浮沉积物浓度对于理解悬浮沉积物输运过程、预测营养盐及污染物的迁移和归宿等问题具有重要的意义。传统的测流方法常采用单点式流速仪，对测点不同深度流速进行逐点测量以获得整个剖面的流速；悬浮沉积物浓度的测量常采用取样过滤称重法，该方法不但费时费力，而且在时间和空间的分辨率有限。
声学技术在海洋观测领域的发展应用，使得海洋水文要素的测量有了质的改变，声学多普勒剖面流速仪（Acoustic Doppler Current Profiler，ADCP），能够快速获得某个完整剖面的流速、流向等参数，同时仪器同步记录的声学背向散射强度可间接计算所测海域悬浮沉积物浓度的分布情况，为海洋泥沙动力输运的研究提供可靠的数据支撑。声学法具有不采样、水体不受扰动、精度高等特点，已被广泛应用于海洋水文要素测量中（Simpson, 2001; Yorke and Oberg, 2002; Ha and Maa, 2008; Salehi and Strom, 2011;张英豪和赖锡军, 2014; Li et al., 2019）。船载走航式ADCP是指将ADCP搭载在调查船上，结合精准的导航系统，在船只走航过程中对所测海域的海流进行连续全剖面测量，具有不扰动流场、测量历时短、测量范围大、充分利用航时等特点，已被广泛应用到江河湖海的流速、流量测量中。如1997 年三峡工程大江截流和2002 年三峡工程明渠截流（杨丰和杨俊青，2004），以及广东北江飞来峡工程截流等（刘勇胜和黄程鹏，2014），杭州湾跨海大桥、港珠澳大桥等工程也均使用ADCP 进行了流速测量（刘彦祥，2016）。
祥云岛位于乐亭县东南部，京唐港与北港之间的沿海地带，岸滩分布有优质天然细沙，是著名的沿海旅游沙岛。祥云岛外海7~15 m等深线海域建有海洋牧场，占海面积1333 hm2。近年来，祥云岛受人为活动及海洋水动力条件影响，沿岸不同区域均出现了不同程度的侵淤现象（陈文超等，2016；汪翡翠等，2021）；与此同时，海洋牧场等近岸工程的建设带来的环境效应逐渐受到重视，如崔雷等（2017）通过耦合二维浅水水动力数值模型及物质输运数值模型，对京唐港填海工程中周边潮流场及悬浮沉积物运移进行模拟研究，并分析了对海域水质环境的影响。因此，对祥云岛海域流速、流向及悬浮沉积物浓度的快速同步调查，对祥云岛及其近海工程的生态环境具有重要意义。本文通过船载走航式ADCP在祥云岛海域进行现场海流、悬浮沉积物全剖面测量，为该海域海洋动力、泥沙输运及祥云岛岸线演化等海洋过程的研究提供支撑。
2. 现场调查
2.1 ADCP测流原理
ADCP基于多普勒原理（刘德铸，2010；张春海和董晓冰，2013），通过测定声波入射到海水中微颗粒背向散射在频率上的多普勒频移，从而得到不同水层水体的运动速度。所谓声学多普勒效应是指当物体相对于声源有径向位移时，导致观测到的声源发生改变的现象。多普勒频移与散射粒子速度之间的关系为：
Fd=2FVC 式中：Fd为声学多普勒频移；F为发射声波的频率；V为颗粒物沿声束方向的相对速度；C为声波在海水中的传播速度，式中“2”是因为ADCP发射并接收回波，所以多普勒频移加倍。假定水体中颗粒物的运动速度和水流速度一致，ADCP通过跟踪颗粒物的运动，测定水流相对于ADCP的速度，此速度扣除船速后，即可得到水流的绝对速度（杨锦坤等，2009）。
此外，由于散射ADCP发射波的颗粒物主要是悬浮沉积物，ADCP同步接收水中散射体的回波强度信号的大小能够间接反映水中悬浮沉积物浓度（杨惠丽等，2017）。利用ADCP背向散射强度反演悬浮沉积物含量是一种便捷，高效的方法，其非干扰式的测量特性更易获得真值。自Deines (1999)提出了声学多普勒剖面仪测得的背向散射强度与悬浮沉积物含量的关系后，许多学者也逐步开展了利用ADCP的背向散射强度反演悬浮沉积物含量的研究（Glenn et al., 1999; Hoitink, et al., 2005; Merckelbach et al., 2006; Gunawan et al., 2010）。已有研究表明源声波散射强度的大小与散射水体中的悬浮物浓度的大小成正比（Thorne and Campbell, 1992），在受潮流、潮汐作用强烈的近岸海域，利用回归公式反演的悬浮沉积物浓度与实测数据的平均误差小于20%（姬厚德等，2018）。对于有限海域，泥沙类型和粒径相对均一情况下，声学散射强度与泥沙浓度的关系可简化为公式：lg(C)=K1E+K2，其中C为悬浮沉积物浓度，E为回波强度，K1、K2为拟合因子，通过现场泥沙浓度和回波强度的回归分析来标定（Farmer, 1998）。
2.2 调查区域及测线布置
祥云岛长约13 km，呈NE－SW走向，岛屿面积20.68 km2，属于离岸沙坝岛。祥云岛东北段为潮汐通道，现已人为改造，西南端为大清河口；东部靠海侧建有简易的人工沙堤，中西部靠海侧残留低缓沙丘。潮间带的沙滩面积约为0.64 km2，潮滩平均坡度3°~8°。沉积物以浅黄色、黄棕色的细沙为主，平均粒径Mz介于2.06~2.34ϕ，分选系数介于0.24~0.47（陈文超等，2016；文明征等，2020；梅西等，2020）。
2019年7月25—28日，在祥云岛海域开展走航式ADCP测量。如图1所示本次测量共设计12条测线，其中测线L1~L10垂直于岸滩分布，平均测线长度约8.0 km；测线L11、L12平行于祥云岛分布，平均测线长度为17.5 km。由于祥云岛潮间带平缓，近岸段水深较浅，测量船容易发生触底的危险。为确保调查过程中设备和人员的安全，选择高潮位时间段进行垂直于岸线的测线测量。根据中国地质调查局天津地质调查中心在津冀沿海的潮汐站监测系统（李勇等，2020），各测线测量时间与潮位对应关系如图2所示。
![]()
图 2 测量期间的潮位数据（不同颜色显示相应调查测线的潮位信息）Figure 2. Tide level data during the survey period (Different colors show the tide level of the corresponding survey line)2.3 参数设置
走航式ADCP测量采用美国TRDI公司生产的瑞江牌（WHR600）声学多普勒流速剖面仪，具有4个换能器，主要技术指标有：工作频率600 kHz；换能器发射角度20°测速范围为±0.01~±10 m/s（水平方向），速度分辨率0.25 cm/s；测深范围：0.75~100 m。本次调查过程中ADCP换能器入水深度为1.10 m，波束角为20°；水流层数设置为50，层厚0.50 m，第一层中心位置为1.96 m。
3. 结果与分析
祥云岛位于半封闭的渤海内部，潮汐类型属于不规则半日潮，测量期间平均潮差0.79 m。测线L1、L3、L5~L7及L10~L12测量期间为涨潮期；测线L2、L4、L8及L9测量期为落潮期（图2）。
3.1 测量结果质量控制
现场流速测量过程中受风浪等环境的干扰，测量数据会出现明显的噪点。因此，数据分析前需要将这部分数据剔除。为控制测量数据质量，本次调查将满足回波强度（Echo Intensity）为50~200 counts，相关系数（correlation）为50~200，良好百分比检验参数大于50%的数据作为有效数据，其余数据予以剔除（杨锦坤等，2009）。通常认为ADCP单次回波所测的流速误差比较大，一般使用一段时间的平均流速来表示某一点或时刻的流速，ADCP流速的平均时间可根据具体情况而定（吴云帆等，2014）。本次测量数据中参与平均的呯（ping）集合数为20（测量时船速约2.1 m/s，相当于平均时长为1 min，平均距离约50 m）。
3.2 祥云岛海域流速特征
3.2.1 海流流向特征
走航式ADCP获得沿测线方向的全水深流速，为更直观地显示剖面流速流向特征，本文选取表、中、底不同深度的海流进行分析。如图3所示，测线L5不同深度水层（水深2.21 m、6.21 m、10.21 m）的航迹和流速矢量图结果表明，在涨、落潮阶段海流流向在不同深度水层的流向基本保持一致。对海流流向沿深度方向求均值，如图4所示：测线L1、L3、L5和L7为涨潮期（图4a），涨潮流为SW方向；测线L2、L4和L8为落潮期（图4b），落潮流为NE方向。因此，祥云岛外海海域流向具有明显的往复特征，海流流向以NE−SW为主，基本与祥云岛岸线平行。
![]()
图 3 测线L5不同深度水层航迹（红线所示）和流速矢量图底跟踪，ADCP通过接接收和处理由河底或海底的回波信号获得船速的方法Figure 3. Track (red line) and velocity vector graph of different water depths of line L5Bottom tracking, ADCP is a method of obtaining ship speed by receiving and processing echo signals from the river bottom or seafloor图2潮位信息显示，测线L9测量时间为平潮转落潮期，其全水深平均流向并没有表现为明显的NE方向（图4b），测线L9不同深度水层航迹和流速矢量结果如图5所示，在浅表层（2.96 m水深），近岸段流速方向为NE向，沿离岸方向流速方向为SW向，水深4.96 m水层流速矢量表明离岸段流向为SW向，水深9.96 m水层显示流速流向在整个测线上均表现为NE方向，与落潮流的方向保持一致。这表明在平潮转为落潮期间，海流转向始于底部海水，且海流转向与潮位变化具有一定的滞后性。为更清晰地描述这一问题，取测线L9不同等深线处海流流向数据绘制流向剖面分布图，如图6所示：相对水深（流速测量点水深/全水深）大于0.5时，海流流向为NE向，小于0.5是海流流向为SW向；测线L9由近岸向离岸方向测量，潮位逐渐由平潮转为落潮，沿等深线增加方向（离岸方向），流速转向的相对深度逐渐减小，再次证明该海域平潮转为落潮期间海流转向始于底层海水。
![]()
图 5 测线L9不同深度水层航迹（红线所示）和流速矢量图Figure 5. Track (red line) and velocity vector graph of different water depths of line L9![]()
图 6 测线L9不同等深线处海流流向剖面图Figure 6. Profile of current direction at different Isobaths of line L93.2.2 海流流速特征
本文选取涨落潮期间典型剖面流速幅值进行对比分析，测线L1、L3、L6和L10为涨潮期、测线L4和L8为落潮期，各测线流速幅值剖面分布如图7所示。涨潮期间，测线 L1位于祥云岛东北端，离岸方向距离海底约11 m高度处有一层相对稳定的高流速层，流速可达0.7~0.8 m/s，向祥云岛滩面方向水深逐渐变浅，水深小于11 m时，该高流速层减弱；测线L10位于祥云岛西南端，其流速特征与测线L1类似，流速沿向岸方向逐渐降低；测线L3流速幅值相对均匀，垂向和向岸方向流速变化不明显；测线L6流速幅值沿向岸方向呈现明显减小的趋势。结果表明：涨潮流流速沿向岸方向逐渐减小；相对于祥云岛东北和西南段，中部流速幅值整体偏小，测线L3和L6涨潮流离岸端流速可达到0.5 m/s以上。测线L4和L8位于祥云岛中部，其流速幅值明显小于测线L3和L6，祥云岛岸段落潮流速明显低于涨潮流速，而测线L8在近岸段的流速要明显大于离岸段。
测线L11与L12均位于涨潮期，涨潮过程中测线L12在时间上更接近于L1和L10（图2）。 以测线L12为例（图7）平行于祥云岛方向，涨潮期间祥云岛东北端存在明显的高流速层，高流速层位于距离海底约11 m处，这与测线L1的测量结果是一致的。为更直观地展现高流速层的分布形态，本文以测线L1为例，取测线L1中不同等深线（14.46 m、11.46 m和7.46 m）处的流速剖面进行对比，如图8所示。离岸方向等深线测量点数值越大，14.46 m处流速剖面结构表现为抛物线型，表明该处存在明显的高流速层，向岸方向，水深减小（7.46 m），流速剖面垂向分布趋于均匀，11.46 m处为过渡段，流速剖面表现为斜直线型。
3.3 祥云岛近岸悬浮沉积物分布
本次测量中未进行同步取样测试，未能建立背向散射强度与悬浮沉积物的回归公式。如前所述，ADCP记录的声学背向散射强度数据能够在一定程度上定性反映祥云岛周边海域的悬浮沉积物分布情况。本次测量中船载走航式ADCP记录了航测过程中的背向散射强度，以涨潮期间测线L1、L5、L10及L12为例，如图9所示。
涨潮期间祥云岛不同位置背向散射强度剖面的分布结果（图9）显示，祥云岛近底普遍存在一层高悬浮沉积物层，且向岸方向悬浮沉积物浓度逐渐增加。图10为测线L1中不同等深线处（14.46 m、11.46 m和7.46 m）的背向散射强度分布对比图，更直观地反映了悬浮沉积物向岸方向增加的现象。测线L12平行于祥云岛方向，沿NE方向在近底逐渐出现一层悬浮沉积物层，高悬浮沉积物浓度层的出现与流速分布（图7，测线L12）具有很好的一致性；水体中悬浮沉积物的维持需要一定的流速，流速幅值为0.30~0.60 m/s时，水体中更容易形成一层稳定的悬浮沉积物层。沿祥云岛岸线由东北向西南方向，流速降低相应的悬浮沉积物浓度也降低，表明悬浮沉积物在向西南方向输运过程中发生沉积。
![]()
图 10 测线L1不同等深线处背向散射强度Figure 10. Backscatter intensity at different isobaths of line L13.4 悬浮沉积物输运对岸线蚀退的影响
祥云岛是滦河早期从大清河口、长河口入海时建造的沙坝（高善明，1981），自滦河迁徙改道北移入海之后，陆源泥沙来源断绝，加之京唐港的修建，港口防沙堤不断向海延伸，截断了沿岸泥沙的补给，之后祥云岛岸线急剧退化。研究表明1987—2000年，祥云岛岸线呈现出北部侵蚀、南部淤积的状态，岸线变化速率分别为−9.20~+10.40 m/a；2000年以后，北部约1500 m砂质岸线转变为港口围堤，被京唐港扩建、祥云湾码头占用，在沿岸流作用下祥云岛东北段泥沙向西南段运移，引起西南段淤积（邢容容等，2017）。邱若峰等（2010）指出祥云岛西南端的金沙岛岸段，自1987年以来向海拓展22 m，该岸段向海推进的一个重要原因在于人为的吹沙造陆，使得该岸段泥沙补给充足（陈文超等，2016）。 本文测量结果显示祥云岛近海海域涨潮流速（SW向）大于落潮流速（NE向）且祥云岛东北端悬浮沉积物大于西南段悬浮沉积物。这表明目前祥云岛近岸海域的悬浮沉积物整体沿NE−SW方向进行运移，与目前祥云岛岸线的侵蚀现状是一致的。
此外，祥云岛海洋牧场于2007年开始申报建设，牧场位于京唐港西南的祥云湾外海海域，水深7~15 m等深线范围内。海洋牧场建设在其规划海域内投放了包括水泥沉箱、水泥铸块、花岗岩石礁、钢筋板框、船礁等在内的各类人工渔礁共计约20万m3（王湛，2018）。人工渔礁的投放必然对近岸海域的流场产生影响，天津地质调查中心2018年的模拟结果显示，人工渔礁的投放仅对投放区周围两倍面积范围内的流场有影响，并导致明显的侵蚀淤积现象
1 。但该结果尚缺乏现场实测数据支撑，有关海洋牧场建设对海域流场的影响，还需要后续海洋牧场建成后的持续跟踪调查。4. 结 论
利用走航式ADCP在祥云岛海域开展流速、流向及悬浮沉积物浓度等水文要素测量，该方法较传统的定点观测和取样测试，具有不扰动流场、测量历时短、测量范围大、充分利用航时等特点，效率大幅提升。测量结果显示：
（1）祥云岛海域海流以平行于岸线为主，具有明显的往复性，该区域涨潮流明显大于落潮流，涨潮流沿向岸方向逐渐减弱，祥云岛西南部局部落潮流存在向岸方向增强的现象。
（2）祥云岛海域在平潮向落潮过渡期间，海流转向始于底部海水，且海流转向与潮位变化具有一定的滞后性。
（3）测量期间背向散射强度数据表明，祥云岛近岸海域存在明显的高浓度层，且沿向岸方向浓度增加，平行于祥云岛岸线方向，悬浮沉积物浓度由东北向西南方向逐渐降低。
（4）祥云岛近岸海域在涨落潮流的作用下，悬浮沉积物沿祥云岛岸线由东北向西南运移，是导致祥云岛岸线在东北段侵蚀、西南段淤积的主要原因之一。
-
![]()
图 3 CO2−CH4在石英纳米孔隙的竞争吸附模拟图（绿色为CO2分子，紫色为CH4分子，采用联合原子力场）
Figure 3. Simulation of competitive adsorption of CO2−CH4 in quartz nanopores (green is CO2 molecule, purple is CH4 molecule)
![]()
图 4 不同目标压力Pt下的CO2（a）和CH4（b）在石英孔隙中的密度分布（红色虚线为±1 nm）
Figure 4. Density distribution of CO2 (a) and CH4 (b) in quartz pores under different target pressures Pt (red dashed line is ±1nm)
![]()
图 5 CO2−CH4 在石英和伊利石表面的吸附选择性与压力关系图
Figure 5. Relationship between adsorption selectivity and pressure of CO2−CH4 on quartz and illitesurface
![]()
图 6 不同温度下CO2(a)和CH4 (b)在石英表面的密度分布（红色虚线为±1 nm）
Figure 6. Density distribution of CO2 (a) and CH4 (b) on quartz surface at different temperatures (red dashed line is ±1 nm)
![]()
图 7 CO2−CH4在石英和伊利石表面的吸附选择性与温度的关系图
Figure 7. Relation between adsorption selectivity and temperature of CO2−CH4 on quartz and illite surface
![]()
图 8 CO2−CH4在伊利石孔隙的竞争吸附模拟（绿色为CO2分子，紫色为CH4分子，目标压力为11.8 MPa, 温度为353.15 K)
Figure 8. Simulation of competitive adsorption of CO2−CH4 in illite pores (green is CO2 molecule, purple is CH4 molecule, target pressure 11.8 MPa, temperature 353.15 K)
![]()
图 9 不同目标压力Pt 下CO2(a)和CH4 (b)在伊利表面的密度分布（红色虚线为±1 nm）
Figure 9. Density distribution of CO2 (a) and CH4 (b) on Erie surface under different target pressures Pt (red dashed line is ±1 nm)
![]()
图 10 在不同温度下的CO2 (a)和CH4 (b)在伊利石表面的密度分布（红色虚线为±1 nm）
Figure 10. Density distribution of CO2 (a) and CO2 (b) on illite surface at different temperatures (red dashed line is ±1 nm)
![]()
图 11 CO2在CO2−CH4−石英系统和CO2−CH4−伊利石系统中的自扩散系数与压力(a)和温度(b)的关系图
Figure 11. The relationship of self−diffusion coefficient of CO2 with pressure (a) and temperature (b) in CO2−CH4−quartz system and CO2−CH4−illite system
表 1 不同压力下CO2−CH4等摩尔化学势及计算压力和摩尔比（温度T = 353.15 K）
Table 1 Equimolar chemical potential, calculated pressure, and molar ratio of CO2−CH4 at different pressures (temperature T = 353.15 K)
目标压力/
MPaCO2化学势/
(Kcal/mol)CH4化学势/
(Kcal/mol)计算混合气体压力/
MPaCO2−CH4
摩尔比5.90 −9.4468 −8.2991 5.93 1.01 9.44 −9.1792 −8.0029 9.27 1.01 11.80 −9.0491 −7.8357 11.82 1.00 14.16 −8.9538 −7.716 14.13 1.02 17.70 −8.8396 −7.5482 18.07 1.01 
下载: 导出CSV
表 2 不同温度下CO2−CH4等摩尔化学势及计算压力和摩尔比（压力P=11.8 MPa）
Table 2 Equimolar chemical potential, calculated pressure, and molar ratio of CO2−CH4 at different temperatures (pressure P = 11.8 MPa)
温度/
KCO2化学势/
(Kcal/mol)CH4化学势/
(Kcal/mol)计算混合
气体压力/MPaCO2−CH4
摩尔比313.15 −6.7985 −7.9549 11.54 1.01 333.15 −7.3047 −8.4817 11.93 1.00 353.15 −9.0491 −7.8357 11.82 1.00 363.15 −8.0965 −9.3324 11.85 0.99 373.15 −8.3638 −9.6213 11.84 0.99
下载: 导出CSV
表 3 不同目标压力下混合气体压力在无孔隙（仅CO2和CH4）和孔隙中的对比
Table 3 Comparison of mixed gas pressures in no pores (only CO2 and CH4) and pores under different target pressures
目标压力/
MPa混合气体压力−
无孔隙/MPa混合气体压力−
石英孔隙/MPa混合气体压力−
伊利石孔隙/MPa5.90 5.93 7.22 7.78 9.44 9.27 11.88 9.48 11.80 11.82 15.06 14.68 14.16 14.13 18.39 16.61 17.70 18.07 23.86 21.78
下载: 导出CSV
表 4 不同温度下混合气体压力在无孔隙（仅CO2和CH4）和孔隙中的对比（目标压力为11.8 MPa）
Table 4 Comparison of mixed gas pressure in no pores (only CO2 and CO2) and pores at different temperatures (target pressure is 11.8 MPa)
温度
（K）混合气体压力−
无孔隙 (MPa)混合气体压力−
石英孔隙 (MPa)混合气体压力−
伊利石孔隙 (MPa)313.15 11.54 16.04 15.72 333.15 11.93 15.68 15.24 353.15 11.82 15.11 14.78 363.15 11.85 15.69 15.17 373.15 11.84 15.00 13.62
下载: 导出CSV
-
[1] Cygan R T, Liang J J, Kalinichev A G. 2004. Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field[J]. The Journal of Physical Chemistry B, 108: 1255−1266.
[2] Dong Dazhong, Zou Caineng, Li Jianzhong, Wang Shejiao, Li Xinjing, WangYuman, Li Denghua, Huang Jinliang. 2011. Shale gas resource potential and exploration and development prospect[J]. Geological Bulletin of China, 31(2): 324−336 (in Chinese with English abstract).
[3] Emami F S, Puddu V, Berry R J, Varshney V, Patwardhan S V, Perry C C, Heinz H. 2014. Force field and a surface model database for silica to simulate interfacial properties in atomic resolution[J]. Chemistry of Materials, 26: 2647−2658.
[4] Hamza A, Hussein I A, Al–Marri M J, Mahmoud M, Aparicio S. 2021b. CO2 enhanced gas recovery and sequestration in depleted gas reservoirs: A review[J]. Journal of Petroleum Science and Engineering, 196: 107685.
[5] Hamza A, Hussein I A, Al–Marri M J, Mahmoud M, Shawabkeh R. 2021a. Impact of clays on CO2 adsorption and enhanced gas recovery in sandstone reservoirs[J]. International Journal of Greenhouse Gas Control, 106: 103286.
[6] Hao Y, Yuan L, Li P, Zhao W, Li D, Lu D. 2018. Molecular simulations of methane adsorption behavior in illite nanopores considering basal and edge surfaces[J]. Energy & Fuels, 32: 4783−96.
[7] Harris J G, Yung K H. 1995. Carbon dioxide's liquid–vapor coexistence curve and critical properties as predicted by a simple molecular model[J]. The Journal of Physical Chemistry, 99: 12021−12024.
[8] Holmboe M. 2019. atom: A MATLAB package for manipulation of molecular systems[J]. Clays and Clay minerals, 67: 419−426.
[9] Humphrey W, Dalke A, Schulten K. 1996. VMD: Visual molecular dynamics[J]. Journal of Molecular Graphics, 14: 33−38.
[10] Jing Shasha. 2015. Molecular Simulation of CO2/CH4 Mass Transfer Process in Sandstone Micropores[D]. Chengdu: Southwest Petroleum University, 55 (in Chinese with English abstract).
[11] Ju Huijiao, Sun Wei, Yang Xipu, Han Zongyuan. 2011. Reservoir characteristics and main controlling factors of Shan–2 member, Yan 'an Area, Ordos Basin[J]. Fault–block Oil and Gas Field, 18(2): 142−145,157 (in Chinese with English abstract).
[12] Le T, Striolo A, Cole D R. 2015. CO2–C4H10 mixtures simulated in silica slit pores: Relation between structure and dynamics[J]. The Journal of Physical Chemistry C, 119: 15274−15284.
[13] Lei Qun, Wang Hongyan, Zhao Qun, Liu Dexun. 2008. Current situation and suggestions on exploration and development of unconventional oil and gas resources at home and abroad[J]. Natural Gas Industry, 28(12): 7−10 (in Chinese with English abstract).
[14] Lv Fangtao, Ning Zhengfu, Mu Zhongqi, Jia Zejiang, Liu Bei. 2023. Molecular simulation of methane flow in rough quartz nanopores[J]. Journal of Northeast Petroleum University, 47(5): 82−91 (in Chinese with English abstract).
[15] Oldenburg C M, Pruess K, Benson S M. 2001. Process modeling of CO2 Injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery[J]. Energy & Fuels, 15: 293−298.
[16] Plimpton S. 1995. Fast parallel algorithms for short–range molecular dynamics[J]. Journal of Computational Physics, 117: 1−19.
[17] Potoff J J, Siepmann J I. 2001. Vapor–liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen[J]. AIChE Journal, 47: 1676−1682.
[18] Shi L Z, Wang Z Z, Xing Z T, Meng S, Guo, S, Wu S M, Luo L Y. 2023. Geological characteristics of unconventional tight oil reservoir: A case study of Upper Cretaceous Qingshankou Formation, northern Songliao Basin, NE China[J]. China Geology, 6: 1−12.
[19] Song Zhengping, Zhang Bin, Kang Tianhe. 2018. Molecular simulation of competitive adsorption of CO2/CH4 in kaolinite based on adsorption site theory[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 37(4): 724–730 (in Chinese with English abstract).
[20] Sun H Y, Zhao H, Qi N, Qi X Q, Zhang K, Sun W C, Li Y. 2016. Mechanistic insight into the displacement of CH4 by CO2 in calcite slit–nanopores: the effect of competitive adsorption[J]. RSC Advances, 6: 104456−104462. doi: 10.1039/C6RA23456A
[21] Sun Haoyang. 2014. Study on Micro–mechanism of Carbon Dioxide Enhanced Shale Recovery[D]. Jinan: Shandong University, 1–103(in Chinese with English abstract).
[22] Sun Ying. 2021. Mechanism of CO2 Enhanced Oil recovery in shale Gas reservoirs based on competitive adsorption[D]. Dongying: China University of Petroleum (East China), 1–116 (in Chinese with English abstract).
[23] Wang Xiangzeng, Qiao Xiangyang, Mi Naizhe, Wang Ruogu. 2018. Beneficial development supporting technology of low permeability tight sandstone gas reservoir in Yan 'an Gas field[J]. Natural Gas industry, 38(11): 43−51 (in Chinese with English abstract).
[24] Wang Xiangzeng. 2016. Progress of unconventional gas exploration and development in Yanchang Petroleum Group[J]. Acta Petrolei Sinica, 37(1): 137−144 (in Chinese with English abstract).
[25] Watts R. 1996. Objectives of the U. S. DOE's research[J]. The Leading Edge, 15: 906−906.
[26] Yang Hu, Wang Jianmin. 2015. Study on tight sandstone reservoir and micro–pore throat characteristics of Shanxi Formation in Yanchang Gas field[J]. Journal of Xi 'an University of Science and Technology, 35(6): 755−762 (in Chinese with English abstract).
[27] Zhang Minghang, Guo Ping, Jiang Wei, Chen Hong. 2016. Study on adsorption characteristics of CO2/CH4 in illite[J]. World Science and Technology Research and Development, 38(5): 950−954 (in Chinese with English abstract).
[28] Zhang Xinmin, Zheng Hongqing, Li Gangyao. 2008. Research on fracturing technology of low permeability tight sandstone gas reservoir in Junggar Basin[J]. Xinjiang Oil and Gas, 4(S1): 69−72 (in Chinese with English abstract).
[29] 董大忠, 邹才能, 李建忠, 王社教, 李新景, 王玉满, 李登华, 黄金亮 2011. 页岩气资源潜力与勘探开发前景[J]. 地质通报, 31(2): 324–336. [30] 景莎莎. 2015. 砂岩微孔隙中CO2/CH4传质过程的分子模拟研究[D]. 成都: 西南石油大学, 1–55. [31] 琚惠姣, 孙卫, 杨希濮, 韩宗元 2011. 鄂尔多斯盆地延安地区山2段储层特征及其主控因素[J]. 断块油气田, 18(2): 142–145, 157. [32] 雷群, 王红岩, 赵群, 刘德勋. 2008. 国内外非常规油气资源勘探开发现状及建议[J]. 天然气工业, 28(12): 7−10. doi: 10.3787/j.issn.1000-0976.2008.12.003 [33] 吕方涛, 宁正福, 穆中奇, 贾泽江, 刘蓓. 2023. 粗糙石英纳米孔隙甲烷流动分子模拟[J]. 东北石油大学学报, 47(5): 82−91. doi: 10.3969/j.issn.2095-4107.2023.05.007 [34] 宋正平, 张彬, 康天合 2018. 基于吸附位理论的CO2/CH4在高岭石中竞争吸附的分子模拟[J]. 矿物岩石地球化学通报, 37(4): 724–730. [35] 孙浩洋. 2014. 二氧化碳提高页岩采收率的微观机制研究[D]. 济南: 山东大学, 1–103. [36] 孙仁远, 张云飞, 范坤坤, 史永宏, 杨世凯. 2015. 岩中黏土矿物吸附特性分子模拟[J]. 化工学报, 66(6): 2118−2122. [37] 孙莹. 2021. 基于竞争吸附的页岩气藏CO2提高采收率机理研究[D]. 东营: 中国石油大学(华东), 1–116. [38] 王香增. 2016. 延长石油集团非常规天然气勘探开发进展[J]. 石油学报, 37(1): 137−144. [39] 王香增, 乔向阳, 米乃哲, 王若谷. 2018. 延安气田低渗透致密砂岩气藏效益开发配套技术[J]. 天然气工业, 38(11): 43−51. doi: 10.3787/j.issn.1000-0976.2018.11.005 [40] 杨虎, 王建民 2015. 延长气田山西组致密砂岩储层及微观孔喉特征研究[J]. 西安科技大学学报, 35(6): 755–762. [41] 张明航, 郭平, 蒋炜, 陈红. 2016. CO2/CH4在伊利石中的吸附特性研究[J]. 世界科技研究与发展, 38(5): 950−954. [42] 张新民, 郑洪庆, 李纲要 2008. 准噶尔盆地低渗致密砂岩气藏压裂工艺技术研究[J]. 新疆石油天然气, 4(S1): 69–72.
计量
- 文章访问数: 775
- HTML全文浏览量: 173
- PDF下载量: 70
