首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The ultra-low-permeability shale gas reservoir has a lot of well-developed natural fractures. It has been proven that hydraulic fracture growth pattern is usually a complex network fracture rather than conventional single planar fractures by micro-seismic monitoring, which can be explained as the shear and tensile failure of natural fractures or creation of new cracks due to the increase in reservoir pore pressure caused by fluid injection during the process of hydraulic fracturing. In order to simulate the network fracture growth, a mathematical model was established based on full tensor permeability, continuum method and fluid mass conservation equation. Firstly, the governing equation of fluid diffusivity based on permeability tensor was solved to obtain the reservoir pressure distribution. Then Mohr–Coulomb shear failure criterion and tensile failure criterion were used to decide whether the rock failed or not in any block on the basis of the calculated reservoir pressure. The grid-block permeability was modified according to the change of fracture aperture once any type of rock failure criterion was met within a grid block. Finally, the stimulated reservoir volume (SRV) zone was represented by an enhancement permeability zone. After calibrating the numerical solution of the model with the field micro-seismic information, a sensitivity study was performed to analyze the effects of some factors including initial reservoir pressure, injection fluid volume, natural fracture azimuth angle and horizontal stress difference on the SRV (shape, size, bandwidth and length). The results show that the SRV size increases with the increasing initial pore reservoir and injection fluid volume, but decreases with the increase in the horizontal principal stress difference and natural fracture azimuth angle. The SRV shape is always similar for different initial pore reservoir and injection fluid volume. The SRV is observed to become shorter in length and wider in bandwidth with the decrease in natural fracture azimuth angle and horizontal principal stress difference.  相似文献   

2.
随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝(射孔孔眼),选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示:增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。  相似文献   

3.
Multistage fracturing of the horizontal well is recognized as the main stimulation technology for shale gas development. The hydraulic fracture geometry and stimulated reservoir volume (SRV) is interpreted by using the microseismic mapping technology. In this paper, we used a computerized tomography (CT) scanning technique to reveal the fracture geometry created in natural bedding-developed shale (cubic block of 30 cm × 30 cm × 30 cm) by laboratory fracturing. Experimental results show that partially opened bedding planes are helpful in increasing fracture complexity in shale. However, they tend to dominate fracture patterns for vertical stress difference Δσ v  ≤ 6 MPa, which decreases the vertical fracture number, resulting in the minimum SRV. A uniformly distributed complex fracture network requires the induced hydraulic fractures that can connect the pre-existing fractures as well as pulverize the continuum rock mass. In typical shale with a narrow (<0.05 mm) and closed natural fracture system, it is likely to create complex fracture for horizontal stress difference Δσ h  ≤ 6 MPa and simple transverse fracture for Δσ h  ≥ 9 MPa. However, high naturally fractured shale with a wide open natural fracture system (>0.1 mm) does not agree with the rule that low Δσ h is favorable for uniformly creating a complex fracture network in zone. In such case, a moderate Δσ h from 3 to 6 MPa is favorable for both the growth of new hydraulic fractures and the activation of a natural fracture system. Shale bedding, natural fracture, and geostress are objective formation conditions that we cannot change; we can only maximize the fracture complexity by controlling the engineering design for fluid viscosity, flow rate, and well completion type. Variable flow rate fracturing with low-viscosity slickwater fluid of 2.5 mPa s was proved to be an effective treatment to improve the connectivity of induced hydraulic fracture with pre-existing fractures. Moreover, the simultaneous fracturing can effectively reduce the stress difference and increase the fracture number, making it possible to generate a large-scale complex fracture network, even for high Δσ h from 6 MPa to 12 MPa.  相似文献   

4.
The hydraulic fracturing technique has been widely applied in many fields, such as the enhanced geothermal systems (EGS), the improvement of injection rates for geologic sequestration of CO2, and for the stimulations of oil and gas reservoirs. The key points for the success of hydraulic fracturing operations in unconventional resources are to accurately estimate the redistribution of pore pressure and stresses around the induced fracture and predict the reactivations of preexisting natural fractures. The pore pressure and stress regime around hydraulic fracture are affected by poroelastic and thermoelastic phenomena as well as by fracture opening compression. In this work, a comprehensive semi-analytical model is used to estimate the stress and pore pressure distribution around an injection-induced fracture from a single well in an infinite reservoir. The model allows the leak-off distribution in the formation to be three-dimensional with the pressure transient moving ellipsoidically outward into the reservoir from the fracture surface. The pore pressure and the stress changes in three dimensions at any point around the fracture caused by poroelasticity, thermoelasticity, and fracture compression are investigated. With Mohr-Coulomb failure criterion, we calculate the natural fracture reactivations in the reservoir. Then, two case studies of constant water injection into a hydraulic fracture are presented. This work is of interest in the interpretation of microseismicity in hydraulic fracturing and in the estimation of the fracture spacing for hydraulic fracturing operations. In addition, the results from this study can be very helpful for the selection of stimulated wells and further design of the refracturing operations.  相似文献   

5.
Sedimentary porous rocks can be used for long-term subsurface containment of CO2. Before injecting CO2 to sedimentary reservoirs, it is necessary to perform stability analysis of the reservoir and to estimate the maximum sustainable pore fluid pressures. In order to avoid the reservoir damage during the CO2 injection, the effective stresses in the reservoir should be evaluated. In this paper, numerical modeling techniques are used for the evaluation of stresses and deformations in a naturally fractured carbonate sedimentary reservoir. The developed numerical modeling scheme couples the behavior of the CO2 injection and the change in the geomechanical behavior of the sedimentary carbonate reservoir for a case study from Saudi Arabia. The present investigation extends the previous studies by considering the sorption-based deformation during the injection of the compressed CO2 fluid into the Arab-D naturally fractured carbonate reservoir. The change in permeability during the injection of CO2 is evaluated. The adopted CO2 injection scenario was shown to be within the safe maximum occupancy, and that the increase in the pore pressure does not result in the reservoir failure.  相似文献   

6.
This paper presents the first published 3D geomechanical modelling study of the CO2CRC Otway Project, located in the state of Victoria, Australia. The results of this work contribute to one of the main objectives of the CO2CRC, which is to demonstrate the feasibility of CO2 storage in a depleted gas reservoir. With this aim in mind, a one-way coupled flow and geomechanics model is presented, with the capability of predicting changes to the in situ stress field caused by changes in reservoir pressure owing to CO2 production and injection. A parametric study investigating the pore pressures required to reactivate key, reservoir-bounding faults has been conducted, and the results from the numerical simulation and analytical analysis are compared. The numerical simulation indicates that the critical pore fluid pressure to cause fault reactivation is 1.15 times the original pressure as opposed to 1.5 times for the comparable analytical model. Possible reasons for the differences between the numerical and analytical models can be ascribed to the higher degree of complexity incorporated in the numerical model. Heterogeneity in terms of lateral variations of hydrological and mechanical parameters, effect of topography, presence of faults and interaction between cells are considered to be the main sources for the different estimation of critical pore pressure. The numerical model, which incorporates this greater complexity, is able then to better describe the state of stress that acts in the subsurface compared with a simple 1D analytical model. Moreover, the reactivation pressures depend mainly on the state of stress described; therefore we suggest that numerical models be performed when possible.  相似文献   

7.
The multi-stage fracturing in horizontal well is a common technique for shale gas reservoir exploitation, in which cluster spacing governs the fracturing performance. Undersized cluster spacing might make the stimulated reservoir volume (SRV), activated by the respective hydraulic fracture, excessively overlap with each other, while oversized cluster spacing might leave a large unstimulated volume between neighboring hydraulic fractures; in either case, fracturing would be inefficient. Previous design of cluster spacing has failed to maximize the SRV due to the absence of a dynamic SRV evaluation model. A numerical model of SRV evaluation in shale reservoir was established by integrating four main modules, including fracture propagation, reservoir pressure distribution, formation stress distribution, and natural fracture failure criterion. Then, a method to optimize cluster spacing was developed with the goal of maximizing SRV. In order to validate this method, it was applied in Fuling shale gas reservoir in Southwest China to determine the optimal cluster spacing. The sensitivity of key parameters on the optimal cluster spacing has been analyzed. This research proposed a compelling cluster spacing optimization method, which could reduce the uncertainty in cluster spacing design, and provides some new insights on the optimal design of multi-stage fracturing in horizontal shale gas well.  相似文献   

8.
水力压裂会对煤储层造成水锁伤害,在压裂液中加入表面活性剂是减缓水锁伤害的有效途径。采用1.5% KCl溶液为基液,在其中加入不同浓度的阴离子型表面活性剂AS和非离子型表面活性剂NS,配置出8种压裂液,分别对河东煤田柳林沙曲矿的焦煤样和太原西山屯兰矿的瘦煤样进行静置沉降实验、毛细管压力测试和离心分离实验,最终优选出煤储层水基压裂液用最佳表面活性剂为0.05% AN复配溶液(AS:NS=9:1),由此构成了表面活性剂压裂液(1.5% KCl+0.05% AN)。研究结果表明,在水基压裂液中加入0.05% AN复配溶液(AS:NS=9:1),可以大幅度降低压裂液的表面张力,改变压裂液的界面状态,从而增加煤表面的亲水性,降低煤孔隙的毛细管压力,使得压裂液的可排性增强,进而能够有效控制储层的水锁伤害,实现增产。   相似文献   

9.
综合压裂技术在低压致密气藏储层中的研究与应用   总被引:1,自引:0,他引:1  
大牛地气田储层属于低压致密气层,非均质性强、孔喉直径小、毛管压力高、多层叠置且层间跨度不等导致储层保护和加砂压裂难度大。通过研究提出了采用液氮增能低伤害压裂液体系,降低对储层的伤害,有效地保护了储层;采用大型加砂压裂技术,增加了裂缝长度,扩大了泄气面积,提高了单层产量,气产量由200 m3/d提高到7.23×104 m3/d;采用不动管柱连续分层压裂工艺技术实现了多层叠置气层均衡改造,压裂投产作业时间缩短了近67%,费用降低了20.69%。  相似文献   

10.
渤海湾盆地渤海D区块沙河街组发育低渗-特低渗透油藏,地层具有异常高温高压特征。通过分析储层填隙物特征和孔隙结构特征,认为渤海D区块沙河街组储层属于中、低孔-低渗型。储层填隙物中易于膨胀及分散运移的伊利石含量较高,储层孔隙空间主要为残余原生孔及溶蚀孔,孔隙和喉道大小为细孔喉级别,孔喉连通性差,非均质性强;通过研究压裂液对储层的敏感性伤害,认为渤海D区块沙河街组储层属中等偏强水敏、弱速敏、中等偏弱碱敏、易水锁伤害。为降低压裂液对储层的敏感伤害、压裂液残留物伤害及与地层环境不配伍造成的储层伤害,渤海D区块高温低渗储层的压裂液须具有耐高温、弱碱性特性,且配方中宜添加防膨剂、稳定剂、缓冲剂和助排剂。  相似文献   

11.
Carbon dioxide (CO2) sequestration in depleted sandstone hydrocarbon reservoirs could be complicated by a number of geomechanical problems associated with well drilling, completions, and CO2 injection. The initial production of hydrocarbons (gas or oil) and the resulting pressure depletion as well as associated reduction in horizontal stresses (e.g., fracture gradient) narrow the operational drilling mud weight window, which could exacerbate wellbore instabilities while infill drilling. Well completions (casing, liners, etc.) may experience solids flowback to the injector wells when injection is interrupted due to CO2 supply or during required system maintenance. CO2 injection alters the pressure and temperature in the near wellbore region, which could cause fault reactivation or thermal fracturing. In addition, the injection pressure may exceed the maximum sustainable storage pressure, and cause fracturing and fault reactivation within the reservoirs or bounding formations. A systematic approach has been developed for geomechanical assessments for CO2 storage in depleted reservoirs. The approach requires a robust field geomechanical model with its components derived from drilling and production data as well as from wireline logs of historical wells. This approach is described in detail in this paper together with a recent study on a depleted gas field in the North Sea considered for CO2 sequestration. The particular case study shows that there is a limitation on maximum allowable well inclinations, 45° if aligning with the maximum horizontal stress direction and 65° if aligning with the minimum horizontal stress direction, beyond which wellbore failure would become critical while drilling. Evaluation of sanding risks indicates no sand control installations would be needed for injector wells. Fracturing and faulting assessments confirm that the fracturing pressure of caprock is significantly higher than the planned CO2 injection and storage pressures for an ideal case, in which the total field horizontal stresses increase with the reservoir re-pressurization in a manner opposite to their reduction with the reservoir depletion. However, as the most pessimistic case of assuming the total horizontal stresses staying the same over the CO2 injection, faulting could be reactivated on a fault with the least favorable geometry once the reservoir pressure reaches approximately 7.7 MPa. In addition, the initial CO2 injection could lead to a high risk that a fault with a cohesion of less than 5.1 MPa could be activated due to the significant effect of reduced temperature on the field stresses around the injection site.  相似文献   

12.
在低渗、低孔隙率页岩气储层压裂改造过程中,压裂液体系的pH值直接影响地层中的造缝、稳缝效果。pH值对阴离子表面活性剂压裂液的粘度、表面张力和页岩膨胀抑制性影响很大。通过改变阴离子表面活性剂的浓度,测定了在不同pH值条件下压裂液体系基本性能的变化。研究表明:SDBS压裂液体系的粘度随着pH值的增加先增大后减小,在pH值=9时达到最大;体系在pH值=5时降低体系表面张力的能力最好;页岩的膨胀量随着pH值的增大逐渐增大,在pH=7时膨胀量最大,且随着pH值的进一步增大,页岩膨胀量先减小后增大,在pH值=9时膨胀量最小;当体系的pH值一定时,SDBS在临界胶束浓度状态能够最大程度抑制页岩的膨胀。从化学微观角度上分析了pH值对SDBS压裂液性能的影响,有利于SDBS压裂液体系配方优化。  相似文献   

13.
Quick-look assessments to identify optimal CO2 EOR storage sites   总被引:1,自引:0,他引:1  
A newly developed, multistage quick-look methodology allows for the efficient screening of an unmanageably large number of reservoirs to generate a workable set of sites that closely match the requirements for optimal CO2 enhanced oil recovery (EOR) storage. The objective of the study is to quickly identify miscible CO2 EOR candidates in areas that contain thousands of reservoirs and to estimate additional oil recovery and sequestration capacities of selected top options through dimensionless modeling and reservoir characterization. Quick-look assessments indicate that the CO2 EOR resource potential along the US Gulf Coast is 4.7 billion barrels, and CO2 sequestration capacity is 2.6 billion metric tons. In the first stage, oil reservoirs are screened and ranked in terms of technical and practical feasibility for miscible CO2 EOR. The second stage provides quick estimates of CO2 EOR potential and sequestration capacities. In the third stage, a dimensionless group model is applied to a selected set of sites to improve the estimates of oil recovery and storage potential using appropriate inputs for rock and fluid properties, disregarding reservoir architecture and sweep design. The fourth stage validates and refines the results by simulating flow in a model that describes the internal architecture and fluid distribution in the reservoir. The stated approach both saves time and allows more resources to be applied to the best candidate sites.  相似文献   

14.
The problem of predicting the geometric structure of induced fractures is highly complex and significant in the fracturing stimulation of rock reservoirs. In the traditional continuous fracturing models, the mechanical properties of reservoir rock are input as macroscopic quantities. These models neglect the microcracks and discontinuous characteristics of rock, which are important factors influencing the geometric structure of the induced fractures. In this paper, we simulate supercritical CO2 fracturing based on the bonded particle model to investigate the effect of original natural microcracks on the induced‐fracture network distribution. The microcracks are simulated explicitly as broken bonds that form and coalesce into macroscopic fractures in the supercritical CO2 fracturing process. A calculation method for the distribution uniformity index (DUI) is proposed. The influence of the total number and DUI of initial microcracks on the mechanical properties of the rock sample is studied. The DUI of the induced fractures of supercritical CO2 fracturing and hydraulic fracturing for different DUIs of initial microcracks are compared, holding other conditions constant. The sensitivity of the DUI of the induced fractures to that of initial natural microcracks under different horizontal stress ratios is also probed. The numerical results indicate that the distribution of induced fractures of supercritical CO2 fracturing is more uniform than that of common hydraulic fracturing when the horizontal stress ratio is small.  相似文献   

15.
以进一步探索表面活性剂压裂液(1.5% KCl+0.05% AN)的防水锁增产机理为目的,采集河东煤田柳林地区沙曲矿的焦煤作为研究对象,通过表面张力和接触角测试、渗透率伤害实验、变压解吸实验,探讨煤储层水锁伤害的控制方法,揭示表面活性剂压裂液防水锁增产作用机理。研究结果表明:表面活性剂压裂液具有极低的表面张力和界面接触角,具备降低煤储层孔喉毛细管压力、增强压裂液可排性的能力,使得其水锁伤害率比活性水压裂液降低约40%,能够有效抑制煤储层的水锁伤害;这种防水锁作用不仅能够提高煤储层束缚水状态下的气相渗透率,实现增透增产,而且还能降低煤层气产出的临界孔径以促进煤层气解吸,实现增解增产。   相似文献   

16.
砂砾岩水力压裂裂缝扩展规律的数值模拟分析   总被引:1,自引:0,他引:1  
李连崇  李根  孟庆民  王昊  王振 《岩土力学》2013,34(5):1501-1507
砂砾岩储层一般具有岩性和渗透性变化大、孔隙度低、连通性差、孔隙结构复杂和非均匀性严重等特点,因此,在水力压裂过程中,裂缝扩展形态难以控制,大规模改造难度大。针对国内某典型砂砾岩油藏特征,采用数值计算方法对砂砾岩压裂裂缝的扩展规律进行了研究,包括地应力场、砾石含量和粒径等对裂缝扩展形态及压裂压力的影响。研究表明,砾石的存在增加了压裂裂纹扩展的复杂性,裂纹主要有止裂、偏转、穿透和吸附4种表现模式,但主应力差严格控制着裂纹的走向,随着主应力差的增大,裂纹由总体绕砾扩展转变为总体穿砾扩展,失稳压力随着主应力差的增大而明显减小;砾石含量的多少体现了砂砾岩试样宏观的非均匀性,含量越高均匀性越差,随着砾石含量的提高,裂纹与砾石的相互作用占据主导地位,失稳压力随砾石含量的增加而增大;当砾石体积含量一定时,砾石粒径对压裂压力的影响主要取决于砾石排列的随机性,失稳压力随砾石粒径的增大而略有增大。  相似文献   

17.
The reaction kinetics and fluid expulsion during the decarbonation reaction of calcite+quartz=wollastonite+CO2 in water-absent conditions were experimentally investigated using a Paterson-type gas apparatus. Starting materials consisted of synthetic calcite/quartz rock powders with variable fractions of quartz (10, 20, and 30 wt%) and grain sizes of 10 µm (calcite) and 10 and 30 µm (quartz). Prior to reaction, samples were HIPed at 700 °C and 300 MPa confining pressure and varying pore pressures. Initial porosity was low at 2.7–6.3%, depending on pore pressure during HIP and the amount and grain size of quartz particles. Samples were annealed at reaction temperatures of 900 and 950 °C at 150 and 300 MPa confining pressures, well within the wollastonite stability field. Run durations were between 10 min and 20 h. SEM micrographs of quenched samples show growth of wollastonite rims on quartz grains and CO2-filled pores between rims and calcite grains and along calcite grain boundaries. Measured widths of wollastonite rims vs. time indicate a parabolic growth law. The reaction is diffusion-controlled and reaction progress and CO2 production are continuous. Porosity increases rapidly at initial stages of the reaction and attains about 10–12% after a few hours. Permeability at high reaction temperatures is below the detection limit of 10–21 m2 and not affected by increased porosity. This makes persistent pore connectivity improbable, in agreement with observed fluid inclusion trails in form of unconnected pores in SEM micrographs. Release of CO2 from the sample was measured in a downstream reservoir. The most striking observation is that fluid release is not continuous but occurs episodic and in pulses. Ongoing continuous reaction produces increase in pore pressure, which is, once having attained a critical value (Pcrit), spontaneously released. Connectivity of the pore space is short-lived and transient. The resulting cycle includes pore pressure build-up, formation of a local crack network, pore pressure release and crack closure. Using existing models for plastic stretching and decrepitation of pores along with critical stress intensity factors for the calcite matrix and measured pore widths, it results that Pcrit is about 20 MPa. Patterns of fluid flow based on mineralogical and stable isotope evidence are commonly predicted using the simplifying assumption of a continuous and constant porosity and permeability during decarbonation of the rock. However, simple flow models, which assume constant pore pressure, constant fluid filled porosity, and constant permeability may not commonly apply. Properties are often transient and it is most likely that fluid flow in a specific reacting rock volume is a short-lived episodic process.Editorial responsibility: J. Hoefs  相似文献   

18.
CO2 pilot injection studies, with site-specific geologic assessment and engineering reservoir design, can be instrumental for demonstrating both incremental enhanced oil recovery and permanent geologic storage of greenhouse gases. The purpose of this paper is to present the geologic and reservoir analyses in support of a field pilot test that will evaluate the technical and economic feasibility of commercial-scale CO2-enhanced oil recovery to increase oil recovery and extend the productive life of the Citronelle Oil Field, the largest conventional oil field in Alabama (SE USA). Screening of reservoir depth, oil gravity, reservoir pressure, reservoir temperature, and oil composition indicates that the Cretaceous-age Donovan sand, which has produced more than 169 × 106 bbl in Citronelle Oil Field, is amenable to miscible CO2 flooding. The project team has selected an 81 ha (200 ac) 5-spot test site with one central gas injector, two producers, and two initially temporarily abandoned production wells that are now in production. Injection is planned in two separate phases, each consisting of 6,804 t (7,500 short tons) of food-grade CO2. The Citronelle Unit B-19-10 #2 well (Permit No. 3232) is the CO2 injector for the first injection test. The 14-1 and 16-2 sands of the upper Donovan are the target zones. These sandstone units consist of fine to medium-grained sandstone that is enveloped by variegated mudstone. Both of these sandstone units were selected based on the distribution of perforated zones in the test pattern, production history, and the ability to correlate individual sandstone units in geophysical well logs. The pilot injections will evaluate the applicability of tertiary oil recovery to Citronelle Field and will provide a large volume of information on the pressure response of the reservoirs, the mobility of fluids, time to breakthrough, and CO2 sweep efficiency. The results of the pilot injections will aid in the formulation of commercial-scale reservoir management strategies that can be applied to Citronelle Field and other geologically heterogeneous oil fields and the design of similar pilot injection projects.  相似文献   

19.
页岩气的工业化开发要求对致密储层进行体积压裂改造,水平井多级分段压裂技术是成功开发页岩气藏的关键技术之一。页岩储层在地应力与致裂压力联合作用下突破页岩基质,沟通天然裂隙等弱结构面,形成高导流能力缝网系统,页岩储层中赋存的游离气和吸附气得以释放,提高了页岩气井的初始产量和最终采收率。根据渝东地区页岩储层工程地质特征,结合渝页A-2HF井、渝页B-2HF井矿场压裂施工数据,探索适用于渝东地区的页岩储层压裂工艺技术,结果表明,采用前置盐酸处理储层保持近井带导流能力,能够保障后续大排量施工压力窗口;低砂比粉砂段塞多级降滤实现“控近扩远”,增加储层改造体积;压力异常层段胶液前置充分造缝,可为后续滑溜水携带支撑剂进入地层创造有利条件。  相似文献   

20.
Very limited investigations have been done on the numerical simulation of carbon dioxide (CO2) migration in sandstone aquifers taking consideration of the interactions between fluid flow and rock stress. Based on the poroelasticity theory and multiphase flow theory, this study establishes a mathematical model to describe CO2 migration, coupling the flow and stress fields. Both finite difference method (FDM) and finite element method (FEM) were used to discretize the mathematical model and generate a numerical model. A case study was carried out using the numerical model on the Jiangling sandstone aquifer in the Jianghan basin, China. The rock mechanics parameters of reservoir and overlying strata of Jiangling depression were obtained by triaxial tests. A two-dimensional model was then built to simulate carbon dioxide migration in the sandstone aquifer. The numerical simulation analyzes the carbon dioxide migration distribution rule with and without considering capillary pressure. Time-dependent migration of CO2 in the sandstone aquifer was analyzed, and the result from the coupled model was compared with that from a traditional non-coupled model. The calculation result indicates a good consistency between the coupled model and the non-coupled model. At the injection point, the CO2 saturation given by the coupled model is 15.39 % higher than that given by the non-coupled model; while the pore pressure given by the coupled model is 4.8 % lower than that given by the non-coupled model. Therefore, it is necessary to consider the coupling of flow and stress fields while simulating CO2 migration for CO2 disposal in sandstone aquifers. The result from the coupled model was also sensitized to several parameters including reservoir permeability, porosity, and CO2 injection rate. Sensitivity analyses show that CO2 saturation is increased non-linearly with CO2 injection rate and decreased non-linearly with reservoir porosity. Pore pressure is decreased non-linearly with reservoir porosity and permeability, and increased non-linearly with CO2 injection rate. When the capillary pressure was considered, the computed gas saturation of carbon dioxide was increased by 10.75 % and the pore pressure was reduced by 0.615 %.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号