深部咸水层CO2地质储存适宜性评价方法研究

随着温室效应的加剧,CO₂地质储存已成为减缓全球气候变暖的有效方法之一. 可用于CO₂ 地下储存的场地主要有枯竭的油气田、深部咸水层和深部不可开采的煤层等,我国深部咸水层CO₂地质储存潜力占总潜力的98%以上. 在全面分析CO₂ 地质储存适宜性影响因素的基础上,建立了适宜于沉积盆地深部咸水层CO₂ 地质储存的适宜性评价体系,主要包括地质安全性、储存规模、社会环境风险和经济适宜性4大指标层,共计28个评价指标,评价方法以层次分析法及指标叠加法为主. 以西宁盆地为研究实例,通过基于排除法的地质条件综合分析与层次分析法的定量评价相互验证,表明该指标体系和评价方法具有广泛的应用价值. 结果表明,西宁盆地一级构造单元中双树坳陷最适宜CO₂地质储存,可作为CO₂地质储存的优先选区.     

Numerical simulation of CO2 disposal by mineral trapping in deep aquifers

Carbon dioxide disposal into deep aquifers is a potential means whereby atmospheric emissions of greenhouse gases may be reduced. However, our knowledge of the geohydrology, geochemistry, geophysics, and geomechanics of CO₂ disposal must be refined if this technology is to be implemented safely, efficiently, and predictably. As a prelude to a fully coupled treatment of physical and chemical effects of CO₂ injection, the authors have analyzed the impact of CO₂ immobilization through carbonate mineral precipitation. Batch reaction modeling of the geochemical evolution of 3 different aquifer mineral compositions in the presence of CO₂ at high pressure were performed. The modeling considered the following important factors affecting CO₂ sequestration: (1) the kinetics of chemical interactions between the host rock minerals and the aqueous phase, (2) CO₂ solubility dependence on pressure, temperature and salinity of the system, and (3) redox processes that could be important in deep subsurface environments. The geochemical evolution under CO₂ injection conditions was evaluated. In addition, changes in porosity were monitored during the simulations. Results indicate that CO₂ sequestration by matrix minerals varies considerably with rock type. Under favorable conditions the amount of CO₂ that may be sequestered by precipitation of secondary carbonates is comparable with and can be larger than the effect of CO₂ dissolution in pore waters. The precipitation of ankerite and siderite is sensitive to the rate of reduction of Fe(III) mineral precursors such as goethite or glauconite. The accumulation of carbonates in the rock matrix leads to a considerable decrease in porosity. This in turn adversely affects permeability and fluid flow in the aquifer. The numerical experiments described here provide useful insight into sequestration mechanisms, and their controlling geochemical conditions and parameters.     

砂岩含水层CO2封存中考虑盐沉淀反馈作用的数值模拟:以鄂尔多斯盆地为例

盐沉淀是含水层CO2封存中需要关注的问题.当前,大多数数值模拟没有考虑盐沉淀引起的地层孔隙度和渗透率变化对流体流动的反馈作用.以鄂尔多斯盆地刘家沟组地层为例,利用TOUGH2软件建立了一个二维模型.通过修改程序源代码,使得模型能考虑盐沉淀对流体流动的反馈作用.模拟结果表明,刘家沟组地层在CO2注入20 a时,盐沉淀的反...     

深部咸水层二氧化碳封存的干化效应及选址启示

将二氧化碳注入到深部咸水层中,形成复杂的多组分、多相流系统。二氧化碳在压力梯度、浓度差作用下不断扩散,逐步带走盐溶液中的水分,导致各组分的相态变化,盐结晶析出,阻塞了咸水层孔隙通道,从而降低了二氧化碳的注入效率,研究该干化效应的影响因素并为工程选址提供依据具有重要意义。采用二维径向模型建立多相流体的流动方程,并结合相对渗透率和毛细压力方程探讨二氧化碳注入速率、咸水层含盐量、毛细压力的特征参数对干化效应的影响,干化效应可用固体饱和度值进行定量描述。结果表明：二氧化碳运移分3个区域：干涸区、气液相混合区及液相咸水区,干化效应主要发生在井周的干涸区。在毛细作用下固体饱和度随注入速率的减小而增大,随咸水层含盐量增大而增大,随毛细作用增大而增大。因此,提高二氧化碳的注入速率,向咸水层中注水稀释含盐量或选择粒径较大的均质咸水层减小毛细作用,均可降低盐结晶对孔隙通道的阻塞,提高注入效率。     

地质封存二氧化碳沿断层泄漏数值模拟研究

安全性和封存效果是二氧化碳（CO2）地质封存（GCS）最受关注的问题,其中以断层/裂缝引起的场地不确定性因素较为复杂。利用储层多相流模拟软件TOUGH2/ECO2N,研究了不同情境下注入深部咸水层中的CO2沿断层发生泄漏的时间、速率、泄漏量等特征。通过15个场地尺度的模型结果分析表明,CO2注入方案、断层性质（发育位置、产状、几何形态、内部结构）、系统内岩层的组合形态对CO2泄漏均有影响。相同储存环境和断层发育条件下,CO2注入速率从1.59 kg/s增加至6.34 kg/s,CO2泄漏时间提前3 706 d,泄漏量在20 a后增加至注入总量的32.43%。断层发育位置对CO2封存影响极为显著,本次研究,距注入井100 m位置的断层可造成CO2的泄漏量在20 a后高达总注入量的63.39%。相同条件下,倾斜/窄小的断层比垂直/宽厚的断层对CO2及储层咸水的泄漏影响更小。断层渗透率增加一倍,可导致CO2泄漏量增加2.11%,泄漏速率约提高0.006~0.01 kg/s。     

煤层气羽状水平井开采数值模拟研究

建立了煤层气定向羽状水平井开采的数学模型,该模型考虑了地层非均质和各向异性的影响、渗透率的压力敏感性和主支、分支井筒内压降影响;本文给出了数值模拟方法,实现了数值模拟计算;并分别对理论模型和樊庄区块煤层气试验区所设计的羽状水平井进行产能预测,证明了模型及算法的正确性。     

规模化深部咸水含水层CO2地质储存选址方法研究

本文依据中国沉积盆地CO2地质储存潜力评价结果,认为深部咸水含水层是实现规模化CO2地质储存的主体,进而对适宜CO2地质储存的深部咸水含水层属性进行了界定。提出了深部咸水含水层CO2地质储存选址原则,合理划分了选址工作阶段。建立了选址技术指标、安全性评价指标、经济适宜性和地面地质-社会环境选址指标4个指标层,60余个指标的选址指标体系,提出了基于层次分析(AHP)的多因子排序选址评价方法。本文的研究成果对中国深部咸水含水层CO2地质储存场地选址具有一定的指导意义。     

Post-injection spreading and trapping of CO2 in saline aquifers: impact of the plume shape at the end of injection

We use an analytical model for the post-injection spreading of a plume of CO₂ in a saline aquifer under the action of buoyancy and capillary trapping to show that the spreading behavior is at all times strongly influenced by the shape of the plume at the end of the injection period. We solve the spreading equation numerically and confirm that, at late times, the volume of mobile CO₂ is given by existing asymptotic analytical solutions. The key parameters governing plume spreading are the mobility ratio, M, and the capillary trapping number, Γ—the former sets the shape of the plume at the end of the injection period, and the latter sets the amount of trapping. As a quantitative measure of the dependence of the spreading behavior on the initial shape, we use a volume ratio. That is, we evolve the plume from a true end-of-injection initial shape and also from an idealized “step” initial shape, and we take the ratio of these mobile plume volumes in the asymptotic regime. We find that this volume ratio is a power-law in M, where the exponent is governed exclusively by Γ. For conditions that are representative of geologic CO₂ sequestration, the ratio of mobile volumes between “true” and “step” initial plume shapes can be 50% or higher.     

基于数值模拟探讨提高咸水层CO2 封存注入率的途径

咸水层CO2地质封存是减少大气中CO2排放量的有效途径。CO2注入率是衡量咸水层中CO2注入能力的有效因素,因此,研究注入速率的变化规律及提高的措施是很有工程价值的。在很多区域,地层的低渗透性限制了CO2的注入率。针对鄂尔多斯盆地的水文地质条件,通过数值模拟,探讨在低渗透性咸水层中提高CO2注入率的途径,包括改变储层中的盐度、采用水平井注入、增加注入井段的长度以及采取水力压裂等工程措施。其中改变储层中的盐度可通过在注入CO2前向储层中注入一定量的水来实现。模拟结果表明,这些方式可以有效地提高CO2注入率,其中水平井改造方式和水力压裂工程措施效果显著,盐度改造措施在地层初始含盐度较高时,会有更好的效果。研究结果可为鄂尔多斯盆地和类似地区的咸水层CO2地质封存项目提供参考。     

Factors affecting CO2 storage capacity and efficiency with water withdrawal in shallow saline aquifers

Carbon sequestration in shallow aquifers can be facilitated by water withdrawal. The factors that optimize the injection/withdrawal balance to minimize potential environmental impacts have been studied, including reservoir size, well pattern, injection rate, reservoir heterogeneity, anisotropy ratio, and permeability sequence. The effects of these factors on CO₂ storage capacity and efficiency were studied using a compositional simulator Computer Modeling Group-General Equation of State Model, which modeled features including residual gas trapping, CO₂ solubility, and mineralization reactions. Two terms, storage efficiency and CO₂ relative breakthrough time, were introduced to better describe the problem. The simulation results show that simultaneous water withdrawal during CO₂ injection greatly improves CO₂ storage capacity and efficiency. A certain degree of heterogeneity or anisotropy benefits CO₂ storage. A high injection rate favors storage capacity, but reduces the storage efficiency and CO₂ breakthrough time, which in turn limits the total amount of CO₂ injected.     

南海西部盐水层CO2埋存潜力评估

南海西部东方1-1、崖城13-1和乐东等气田具有CO2含量高的特点,每年从生产出的天然气中分离出大量的CO2。为了减少CO2在大气中的排放,考虑在莺歌海地区选择合适的盐水层埋存体,拟进行CO2地质埋存示范工程。根据CO2在盐水层中的各种埋存机理,并考虑盐水层构造特征对CO2运移分布的影响,提出了一种较为准确的CO2埋存潜力评估方法,并利用此方法对筛选出的5个备选盐水层进行了评估。评估结果表明,各盐水层埋存潜力巨大,都远远大于示范工程期限内预计的CO2埋存总量,并可将这些盐水层作为将来海南省及临近广东省人为CO2的埋存场所。南海西部东方1-1、崖城13-1和乐东等气田具有CO2含量高的特点,每年从生产出的天然气中分离出大量的CO2。为了减少CO2在大气中的排放,考虑在莺歌海地区选择合适的盐水层埋存体,拟进行CO2地质埋存示范工程。根据CO2在盐水层中的各种埋存机理,并考虑盐水层构造特征对CO2运移分布的影响,提出了一种较为准确的CO2埋存潜力评估方法,并利用此方法对筛选出的5个备选盐水层进行了评估。评估结果表明,各盐水层埋存潜力巨大,都远远大于示范工程期限内预计的CO2埋存总量,并可将这些盐水层作为将来海南省及临近广东省人为CO2的埋存场所。     

Experimental and numerical simulation study of the mineral sequestration mechanism of the Shiqianfeng saline aquifers in the Ordos Basin,Northwest China

This study focused on typical injection layers of deep saline aquifers in the Shiqianfeng Formation used in the Carbon Capture and Sequestration Demonstration Projects in the Ordos Basin, Northwest China. The study employed experiments and numerical simulations to investigate the mechanism of CO₂ mineral sequestration in these deep saline aquifers. The experimental results showed that the dissolved minerals are plagioclase, hematite, illite–smectite mixed layer clay and illite, whereas the precipitated minerals are quartz (at 55, and 70 °C) and kaolinite (at 70 °C). There are rare carbonate mineral precipitations at the experimental time scale, while the precipitation of quartz as a product of the dissolution of silicate minerals and some intermediate minerals rich in K and Mg that transform to clay minerals, reveals the possibility of carbonate precipitation at the longer time scale. These results are consistent with some results previously reported in the literature. We calibrated the kinetic parameters of mineral dissolution and precipitation by these experimental results and then simulated the CO₂ mineral sequestration under deep saline aquifer conditions. The simulation results showed that the dissolved minerals are albite, anorthite and minor hematite, whereas the precipitated minerals are calcite, kaolinite and smectite at 55 and 70 °C. The geochemical reaction of illite is more complex. At 55 °C, illite is dissolved at the relatively lag time and transformed to dawsonite; at 70 °C, illite is precipitated in the early reaction period and then transformed to kaolinite. Based on this research, sequestrated CO₂ minerals, which are mainly related to the temperature of deep saline aquifers in Shiqianfeng Fm., are calcite and dawsonite at lower temperature, and calcite at higher temperature. The simulation results also establish that calcite could precipitate over a time scale of thousands of years, and the higher the temperature the sooner such a process would occur due to increased reaction rates. These characteristics are conducive, not only to the earlier occurrence of mineral sequestration, but also increase the sequestration capacity of the same mineral components. For a sequestration period of 10,000 years, we determined that the mineral sequestration capacity is 0.786 kg/m³ at 55 °C, and 2.180 kg/m³ at 70 °C. Furthermore, the occurrence of mineral sequestration indirectly increases the solubility of CO₂ in the early reaction period, but this decreases with the increase in temperature.     

Numerical simulation of fracture path and nonlinear closure for simultaneous and sequential fracturing in a horizontal well

Reservoir stimulation requires a model to evaluate the fracture path and closure for the simultaneous or sequential propagation of the hydraulic fracture (HF). This paper presents a fluid-solid coupled model to simulate multi-stage HF propagation. A non-linear joint model is proposed to evaluate the fracture closure when the created fractures are elastically propped. HF closure continues until the balance of external stress matches the proppant's resistance. The reservoir along the horizontal wellbore is not stimulated equally by the multi-stage fracturing. The HFs in the subsequent stage are ‘repelled’ and restrained by the HFs in the previous stage.     

鄂尔多斯深部咸水层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地质储存为例

在地下流动系统问题的研究中,热-水动力-力学(THM)耦合过程是研究的热点问题。在地下多相非等温数值模拟软件TOUGH2的框架内,基于Biot固结理论和摩尔-库仑破坏判定准则,建立了THM耦合模型;采用积分有限差和有限元联合的空间离散方法,开发了THM模拟器TOUGH2Biot。该模拟器中热和水动力过程是全耦合,力学过程是部分耦合。通过与解析解的对比,验证了其正确性。基于鄂尔多斯盆地CCS示范工程,采用TOUGH2Biot研究了CO2注入地层后的THM响应。结果显示CO2的注入引起流体压力急剧增加,地层有效应力减小,地表隆起,隆起大小在几十个厘米,同时孔渗增加,利于CO2注入引起的压力上升向外消散。CO2注入最有可能导致剪切破坏的位置位于最大速率注入点上部盖层,其次为靠近地表的位置。     

离散裂隙网络对页岩气井产能影响的数值模拟

页岩作为典型的非常规储层,基质孔隙小,渗透率极低,水平井多级水力压裂为其商业开发的主要手段。准确模拟页岩气产能,应同时考虑水力裂隙和天然裂隙的渗流。基于离散裂隙模型和等效连续模型建立页岩气渗流数学模型,利用有限元分析方法进行数值求解,研究不同走向裂隙组对页岩气井产能的影响。研究认为,页岩基质为气体的生产提供了主要气源,天然裂隙作为渗流的主要通道,将气体输送到水力裂缝,进而到达井筒。模拟结果表征,离散裂隙的渗流特征对于页岩气井的产能有重要影响。根据页岩储层的天然裂隙走向,可以优化相应的水平井方位。对于二维离散裂隙网络模型,水平井沿着2个裂隙组夹角的平分线更有利于生产。     

Geochemical effects of CO2 sequestration in sandstones under simulated in situ conditions of deep saline aquifers

The geochemical effects of brine and supercritical CO₂ (SCCO₂) on reservoir rocks from deep (1500–2000 m) saline aquifers were examined via experimental simulation at in situ conditions. Dry sandstone samples were mounted in a triaxial cell and autoclave system, evacuated, and saturated with 1 M NaCl solution. The brine-rock system was allowed to react at 30 MPa confining pressure, 15 MPa pore fluid pressure, and 60 °C while SCCO₂ was injected at a pressure gradient of 1–2 MPa. The experiment was conducted for a period of 1496 h, during which fluids were periodically sampled and analyzed. The pH measured in partially degassed fluid samples at 25 °C decreased from a starting value of 7.0–4.3 (9 days) and finally 5.1 after saturation with SCCO₂.