深部咸含水层CO2注入的流体迁移模拟研究 ——以松辽盆地三肇凹陷为例

CO2地质封存是减少温室气体向大气排放的有效措施之一,而深部咸含水层CO2地质封存是目前可行的最有潜力的封存技术。先前研究表明,松辽盆地是一个潜在的封存场地。基于对松辽盆地地质资料的初步分析,选取三肇凹陷的姚家组1段和青山口组2、3段地层作为CO2的注入层,建立一个典型二维模型,研究CO2注入后的迁移规律。结果表明,CO2注入后会向上和侧向迁移,后期可能出现的对流作用能促进CO2的溶解。残留气体饱和度、注入层水平和垂直渗透率的比值对模拟结果影响最大。此外,储层中的薄页岩夹层有利于CO2的溶解,因此,在保证注入性和封存量的情况下,储层中低渗透性夹层是允许的。     

江汉盆地江陵凹陷二氧化碳地质封存数值模拟

通过建立江汉盆地江陵凹陷新沟嘴组储层和盖层的二维模型,对江陵凹陷进行了CO2地质封存数值模拟研究,分析了注入CO2的运移分布和溶解扩散情况,并对储盖层垂直方向的渗透率进行了敏感性分析.结果表明,所建的二维模型能够较好地描述CO2在储层中的运移过程和分布状况,而且在实际地质封存工作中必须要考虑一些地质上的因素如储盖层的渗透率、毛细管压力、岩石矿物组分丰度等的变化对封存效果的影响.     

注入速率对CO_2地质储存封存潜力的影响分析

国内某著名煤炭企业计划实施每年10万吨的CO2地质储存（CCS）项目,拟选了5组地层做为目标储层。但所选封存层平均渗透率在0.15~0.6mD,平均孔隙度在2~6%,属于低渗低孔地层,如不进行人工压裂提高注入层渗透率,要实现预定存储目标尚有一些困难。笔者在研究中发现,除对目标层进行一定的水裂酸化处理提高地层渗透特性可以显著提高注入性和存储能力外,CO2注入速率的变化对地层的封存能力和注入性也有明显影响。运用TOUGH2-ECO2N软件分别模拟了无水裂及水裂情况下8种不同注入速率下这些目标存储层的压力变化及CO2封存状态比例及理论最大封存能力。模拟结果表明使用水裂酸化方法对储层进行处理后,不仅可以使注入总量达到项目要求,还可使系统理论最大储存能力提高55%;并且在灌注过程中采用变速灌注方式,可以有效控制系统压力积聚,对将来实际灌注压力控制具有重要意义。     

岩盐沉淀对咸水层二氧化碳地质封存注入过程的影响:以江汉盆地为例

依托江汉盆地基础地质和水文地质数据,通过建立新沟嘴组一维、二维径向数值模型,对该储层进行CO2封存注入数值模拟研究,分析二氧化碳注入储层后,咸水层中发生的地球物理化学作用.以CO2的溶解作用和岩盐的沉淀作用为主,分析它们对储层压力和注入性的影响.模拟结果表明,二氧化碳注入储层后,近井周围会形成3个分区,即纯液相区、两相混相区、盐沉淀区;盐沉淀会导致压力积累严重,影响CO2可注入性;对盐度的敏感性分析表明,咸水含水层盐度越大,压力积累会越严重.因此,江汉盆地CO2地质封存应考虑盐度的影响,采取相关缓解压力积累的措施.     

数值模拟在CO2地质封存示范项目中的应用

CO2深部咸水层地质封存被认为是减缓温室效应的一种有效的工程技术手段。针对神华鄂尔多斯105 t/a CO2捕集与封存(CCS)示范项目,用数值模拟方法对CO2在地层中的运移过程进行了详细地刻画,分析了CO2的流动迁移、地层压力积聚过程及地层封存潜力。数值模型不但可以为工程的顺利进行提供技术支撑,而且可以节省人力财力。首先,根据实际监测数据对模拟参数进行校准,得到了合适的压力拟合曲线,确定了主要的水文地质参数。然后,对为期3 a的CO2续注工程进行预测,详细分析了CO2的晕扩散、溶解情况、地层压力变化情况、储层封存潜力等。得到如下结论：CO2在3 a内的最大迁移距离约为350 m;水裂可以有效提高CO2的注入性;隔离层能有效防止CO2逃逸。研究表明,尽管鄂尔多斯盆地属于低渗咸水层仍然能够有效安全地封存CO2。     

CO_2地质封存四维多分量地震监测技术进展

CCS技术是目前公认的快速减缓温室效应的最有效方法,CO2地质封存是CCS技术最核心的问题之一,监测CO2地质封存的安全性贯穿于CO2注入过程中与封存以后。四维地震监测技术是监测CO2是否泄漏、证实CO2封存安全性最有效的技术手段。常规四维地震技术通过对比CO2注入前后及注入不同阶段2次或者多次三维地震纵波振幅差异与旅行时差异,确定CO2在地下分布。而纵波振幅或旅行时差异是CO2饱和度与孔隙压力的综合反映,单纯的纵波信息难以区分饱和度与压力信息。目前,四维多分量地震监测技术的潜力并未挖掘,由于横波速度对于压力敏感,利用四维转换波信息监测CO2地质封存,可以识别注入CO2的压力分布范围。对于各向异性介质的储层,对比一次地震观测PS1,PS2旅行时、振幅差异与2次地震采集之间PS1,PS2旅行时、振幅差异,还可以有效确定注入CO2前与注入期间储层裂隙、裂缝的变化,以及储层与盖层的应力状态。四维多分量地震资料结合岩石物理资料和全波列测井资料,可以更准确地确定可能的CO2泄漏风险区域,更加可靠地评估CO2地质封存的安全性。     

页岩储层物性参数对CO2不同封存机制及封存量的影响

CO₂增强页岩气开采技术(CO₂-ESGR)一方面可提高CH4产量,另一方面又可实现CO₂地质封存。为了分析页岩储层物性参数对CO₂封存机制的影响,文章以鄂尔多斯盆地延长组页岩为研究对象,采用CMG-GEM软件建立双孔双渗均质模型,分析了CO₂-ESGR中页岩储层垂直渗透率与水平渗透率之比(Kv/Kh)、含水饱和度和孔隙度对不同CO₂封存机制封存量的影响;并设计了27组正交试验采用极差分析法比较了三种因素的影响程度。研究表明,Kv/Kh在0.1～1范围增大会增加不同CO₂封存机制的封存量,封存总量最大可增加69.96%,其中吸附封存量最大可增加97.96%,受到影响最大;含水饱和度在0～0.9范围增大引起CO₂封存总量先增加后减小,封存总量最大可减少67.12%,其中溶解封存量最大可减少83.35%,范围波动最大;页岩储层孔隙度在0.1～0.99范围增大会导致CO₂封存总量减少,封存总量最大...     

鄂尔多斯深部咸水层CO_2地质封存效果评价

中国首个陆上咸水层CO_2地质封存全流程示范项目于2010年正式实施。为更加清晰、准确地了解注入场地储层的注入性能和注入封存过程中可能遇到的潜在问题,基于场地储层结构和注入监测数据,采用储层多相流模拟软件TOUGH2-MP/ECO_2N对鄂尔多斯105 t/a CO_2注入1 620 m以深的特低渗砂岩咸水含水层进行数值模拟,对储层的压力积聚和CO_2羽体扩散的动态演化以及储层封存量进行评估。结果表明,所建立的模型比较准确地反映了实际注入过程和注入效果。3a注入引起的最大压力抬升小于15 MPa,CO_2在含水层中总体呈均匀扩散,CO_2注入地下3a和53a后,羽体在刘家沟储层中的横向迁移距离分别为550 m和700 m左右。在目前的统注方案下,CO_2主要封存层位在储层上部的刘家沟组(埋深为1 690~1 699 m),其吸气量占整个储层封存量的80%以上,储层吸气能力具有由浅到深变差的特征。53a模拟期内,进入泥岩盖层的CO_2总量不及注入总量的0.05%。     

超临界CO_2增强热卤水开采模型研究

面对全球气候变化和能源紧缺的巨大压力,CO2地质利用技术成为研究的焦点。为了实现CO2减排及资源化利用,提出了利用超临界CO2增强卤水提取和地热能开采的CO2地质利用集成系统的概念。通过建立CO2注入热卤水储层的质量平衡模型,分析了不同流量的卤水生产、CO2注入以及储层边界流对储层变化的影响,初步评估了该系统的CO2封存量、卤水提取量以及储层流体组成变化的时间尺度。研究表明,注入CO2提供了热卤水层的压力维持,促进卤水和地热资源的可持续开采,对于江陵凹陷研究区9×108 m3的储层有效体积,注入9.95×106 t CO2可提取17.12×106 t卤水,时间尺度超过30a。对于50~1 000kg/s的卤水生产速度,可以产生0.9~18.8 MW电力。同时,该技术增加了CO2的封存容量和效率,有利于CO2大规模安全封存,经济和环境效益显著。     

层状非均质性对CO_2在盖层中迁移泄漏规律的影响

盖层密封性能对CO_2地质封存工程的容量、经济性、安全性具有重要影响。盖层岩石具有明显的层状非均质特点,而这种层状非均质特点直接影响CO_2在盖层中的迁移与泄漏过程,从而影响盖层的密封性能。现有研究考虑了各向同性的非均质性对盖层密封性能的影响,但并未针对层状非均质性盖层展开相应的研究。建立了考虑孔隙度、渗透率和毛管压力参数变化关系的均质、非层状非均质和层状非均质3种盖层地质模型,并采用TOUGH2对统一边界条件的不同地质模型开展了数值模拟对比研究,分析了CO_2进入盖层、盖层内迁移和穿透盖层而形成缓慢泄漏的全过程。模拟结果表明:3种模型中,自由态CO_2饱和度分布与孔隙度、渗透率的分布基本一致;与均值模型相比,非层状非均质性会加速CO_2在盖层中的迁移,而层状非均质性将会阻碍CO_2在盖层中的迁移,增加穿透时间;与均质模型和非层状非均质模型相比,层状非均质模型将会大大降低CO_2的穿透盖层以后的泄漏速率和泄漏总量。针对建模和数值模拟过程中遇到的问题,提出了盖层地质建模过程中的一些局限性。其结果为更加客观和定量地评价CO_2在盖层内迁移规律及泄漏量奠定了一定的基础。     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

土压缩模量预测回归分析方法

在土压缩模量Es的分析计算中,常用的分析方法是根据土的e-p曲线试验资料建立公式求解,本文引入了数学上的回归分析方法来探讨Es与其它各项物理力学指标参数之间的相关关系,旨在找出对Es影响最显著的指标参数,并建立最优回归方程式,以此来预测Es。研究结果表明:该方法通用性强,应用范围广,亦可推广应用于其它相关领域。     

Causes of salinization of the Gulf of Taganrog

Using the database of automatic hydrometeorological stations, installed in the Don RIver delta and Taganrog Bay seashore, the sources of the anomalois scale water negative surge and salinization of the Azov Sea under conditions of low river flow in 2015–2016 are studied. The new schemes of stratification and advection of salty sea waters in the Don River mouth under different weather conditions, water discharge and levels are given.     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

电信网络两类逻辑图自动生成算法

为解决电信手工绘制的逻辑图不易维护的问题,研究了电信带方向逻辑图和不带方向逻辑图的自动生成算法.算法是从数据库提取缆线的路由信息,用深度优先的方式遍历缆线经过的所有坐标,采取回溯的方式,计算缆线在每个点的偏移值,调整缆线顺序,自动绘制逻辑图.产生的逻辑图没有不必要的交叉,避免了传统的手工绘制.结果表明,对于常用的电信逻辑图,自动生成算法是可以满足要求的.      

Excitation of Under-Ice Seiches of a Sea Port of the Sea of Okhotsk

Doklady Earth Sciences - Multiannual observations (2009–2017) by means of autonomous seiche recorders have been conducted for marine seiches with one second discreteness in the boot basin for...     

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法