CN117236232B - 天然气水合物与浅层气、深层气联合开采模拟方法和系统 - Google Patents
天然气水合物与浅层气、深层气联合开采模拟方法和系统 Download PDFInfo
- Publication number
- CN117236232B CN117236232B CN202311515925.5A CN202311515925A CN117236232B CN 117236232 B CN117236232 B CN 117236232B CN 202311515925 A CN202311515925 A CN 202311515925A CN 117236232 B CN117236232 B CN 117236232B
- Authority
- CN
- China
- Prior art keywords
- gas
- mass
- hydrate
- layer
- component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000007789 gas Substances 0.000 title claims abstract description 301
- 238000004088 simulation Methods 0.000 title claims abstract description 51
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000004519 manufacturing process Methods 0.000 claims abstract description 71
- 230000035699 permeability Effects 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 83
- 239000012071 phase Substances 0.000 claims description 70
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 47
- 239000008346 aqueous phase Substances 0.000 claims description 21
- 238000000354 decomposition reaction Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 10
- 238000004134 energy conservation Methods 0.000 claims description 9
- 238000005065 mining Methods 0.000 claims description 9
- 230000009969 flowable effect Effects 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 abstract description 3
- 238000010168 coupling process Methods 0.000 abstract description 3
- 238000005859 coupling reaction Methods 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000012512 characterization method Methods 0.000 abstract description 2
- 230000001186 cumulative effect Effects 0.000 description 6
- 230000007423 decrease Effects 0.000 description 6
- 239000003345 natural gas Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- -1 Natural gas hydrates Chemical class 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 230000036962 time dependent Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 150000004677 hydrates Chemical class 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Landscapes
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明公开了天然气水合物与浅层气、深层气联合开采模拟方法和系统，涉及天然气水合物开采技术领域，包括：构建浅层气层与水合物层联合开采模拟模型，包括设定浅层气层与水合物层并网格划分，通过设置地层参数、生产参数构建联合开采模拟模型，求解模型获取产能数据；构建浅层气层、深层气层与水合物层联合开采模拟模型，包括设定浅层气层、深层气层与水合物层并网格划分，通过设置地层参数、生产参数构建联合开采模拟模型，求解模型获取产能数据。本发明基于渗流场、温度场、力学场相互作用机制进行三场耦合流动模拟，建立天然气水合物储层孔隙度、渗透率、饱和度、毛管力等物性参数的表征方法，分别探明了两气合采和三气合采的产能变化规律。
Description
技术领域
本发明涉及天然气水合物开采技术领域，尤其涉及一种天然气水合物与浅层气、深层气联合开采模拟方法和系统。
背景技术
随着经济的发展和社会的进步，能源需求逐渐增加。美国能源信息署在2019年国际能源展望中报告称，2018年至2050年间，世界能源消耗量将增加近50％。随着常规化石燃料储量的减少，探索和开发非常规油气资源势在必行。天然气水合物作为非常规能源的重要组成部分，受到了全世界的关注和研究。与其他传统化石能源相比，它是一种新型清洁燃烧化石能源，产生的CO2更少。除此之外，天然气水合物中的总碳量是传统化石能源的两倍。
天然气水合物是在低温高压条件下由天然气和水发生化学反应而形成的类冰状的结晶物质，俗称可燃冰，是一种非常规的、清洁的天然气资源，其分布面广、资源量大。标准状况下，1m3天然气水合物分解产出164m3天然气和0.8m3水。
我国天然气水合物储量丰富，但目前仍处于实验室模拟、地质调查等阶段，水合物储层的数值模拟技术是一种有效的评价方法。鉴于目前有关开采水合物模拟的研究还不够全面，且缺少对相关开采动力学规律以及多气合采相关开采方案的系统性研究，需要对各种水合物藏的开采方法进行系统的研究，并对其开采规律进行分析，厘清多气合采的特征规律，并确定最佳的开采参数，从而为多气合采提供科学的依据。
数值模拟技术是研究天然气水合物储层的关键技术之一，通过数值模拟技术可以对各种开采方式的采收率进行一定的预测，并对其进行合理的开发和设计；同时，通过数值模拟，可以对水合物的各个参数进行分析，从而对水合物的动力学特性进行再认识，从而为今后的开发利用提供依据。
因此，亟需研究一种天然气水合物与浅层气、深层气联合开采模拟方法和系统。
发明内容
为解决上述技术问题，本发明公开了一种天然气水合物与浅层气、深层气联合开采模拟方法和系统，能够厘清多气合采的特征，解决试采中存在的日产气量低，稳产时间短的问题。本发明基于渗流场、温度场、力学场相互作用机制进行三场耦合流动模拟，建立天然气水合物储层孔隙度、渗透率、饱和度、毛管力等物性参数的表征方法，可以有效阐明多气合采储层物性参数时变规律及演化特征，明确不同开采阶段储层物性参数变化规律。
为实现上述目的，本发明采用下述技术方案：
本发明第一方面提出了天然气水合物与浅层气、深层气联合开采模拟方法，包括以下步骤：
S1、构建浅层气层与水合物层联合开采模拟模型，包括设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数等，构建联合开采模拟模型，求解模型获取产能数据；
S2、在步骤S1的基础上，在水合物层下方增加深层气层，并进一步对深层气层进行网格划分，设定地质参数、生产参数等，之后构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据。
可选地，步骤S1中，水合物分解动力学方程为：
式中，Nh为水合物结合的水分子数；
水合物层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，t为时间，s；κ为组分标识，本方程中表示水合物组分h、甲烷组分m、水组分w或能量e；Mκ为κ的所有组分的和，kg/m3或J/m3；Fκ为κ的可流动的组分，kg/(m2·s)；qκ为κ的源汇，kg/(m3·s)或J/(m3·s)；
水合物组分质量守恒方程为：
Mh＝φSHρH；
式中，Mh为水合物组分的质量之和，kg/m3；SH为水合物相的饱和度；ρH为水合物相的密度，kg/m3；
甲烷组分质量守恒方程为：
式中，Mm为甲烷组分的质量之和，kg/m3；φ为储层孔隙度；SA为水相的饱和度；SG为气相的饱和度；ρA为水相的密度，kg/m3；ρG为气相的密度，kg/m3；为甲烷组分的质量与水相的质量之比；/>为甲烷组分的质量与气相的质量之比；/>为甲烷组分的质量与水合物相的质量之比；FA为水相的在单位截面积通过的质量流量，kg/(m2·s)；FG为气相在单位截面积通过的质量流量，kg/(m2·s)；Fm为甲烷组分在单位截面积通过的质量流量，kg/(m2·s)；qm为甲烷组分的源汇，kg/(m3·s)；qA为水相的源汇，kg/(m3·s)；qG为气相的源汇，kg/(m3·s)；/>为甲烷组分的质量与水相源汇的质量之比；/>为甲烷组分的质量与气相源汇的质量之比；
水组分质量守恒方程为：
式中，Mw为水组分的质量之和，kg/m3；为水组分的质量与水相的质量之比；/>为水组分的质量与气相的质量之比；/>为水组分的质量与水合物相的质量之比；Fw为水组分在单位截面积通过的质量流量，kg/(m2·s)；qw为水组分的源汇，kg/(m3·s)；/>为水组分的质量与水相源汇的质量之比；/>为水组分的质量与气相源汇中质量之比；
能量守恒方程为：
式中，Me为能量之和，J/m3；ρR为岩石的密度，kg/m3；HR为岩石的焓，J/kg；sβ为β相的饱和度；Hβ为β相的焓，J/kg；ΔsH为当前时间步内水合物饱和度的变化值；ΔH0为水合物的分解/形成焓，J/kg；Fe为能量流速，J/(m2·s)；Kc为体系的综合导热系数，W/(m·K)；Fβ为β相在单位截面积通过的质量流量，kg/(m2·s)；qe为能量的源汇，J/(m3·s)。
可选地，步骤S2中，浅层气层、深层气层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，x为气组分g或能量e；Mx为x的所有组分的和，kg/m3或J/m3；Fx为x的可流动的组分，kg/(m2·s)；qx为x的源汇，kg/(m3·s)或J/(m3·s)；
气组分质量守恒时为：
Mg＝φSGρGXG；
Fg＝XGFG；
qg＝Xq,GqG；
式中，Mg为气组分的质量之和，kg/m3；XG为气组分的质量与气相的质量之比；Fg为气组分在单位截面积通过的质量流量，kg/(m2·s)；qg为气组分的源汇，kg/(m3·s)；Xq,G为气组分的质量与气相源汇的质量之比；
能量守恒时为：
Me＝(1-φ)ρRHR+φsGρGHG；
qe＝HGqG；
式中，sG为气相的饱和度；HG为气相的焓，J/kg；FG为气相在单位截面积通过的质量流量，kg/(m2·s)。
可选地，步骤S2中，深水气层温度较低，实际气体压缩性与理想气体存在较大差异，因此在气体状态方程中引入偏差因子：
pV＝nZRT；
式中，p为气体的绝对压力，Mpa；T为气体的绝对温度，K；V为气体的体积，m3；n为气体的物质的量，mol；R为气体常数，8.314×10-3Mpa/(mol·K)；Z为气体的偏差因子，无量纲；
气体渗流与液体渗流相类似，当气体处于层流状态时，采用达西渗流定律描述其流动状态，对于均质地层中，在三维渗流空间中，广义达西定律为：
气体流动在三维空间中的流动速度表示为：
式中，v为气体渗流速度，m/s；K为地层渗透率，D；μ为气体粘度，mPa·s；g为重力加速度，g/cm3；x、y、z为空间坐标轴；
当气体的渗流速度增加到一定程度之后，紊流和惯性的影响将会愈发显著，其渗流速度与压力梯度呈非线性关系，不满足达西渗流定律；在水平方向上，气体渗流过程中有紊流和惯性阻力存在时，满足其动力学规律的非线性二次项式运动方程为：
式中，ζ为影响紊流和惯性阻力的孔隙结构特征参数；
上式右端第一项为黏滞阻力，与渗流速度成正比，第二项为惯性阻力，与渗流速度平方成正比，所以当渗流速度较小之时，可以不考虑惯性阻力的影响；当渗流速度较大时，紊流和惯性的影响逐渐显著，当渗流速度增加到偏离达西定律时，惯性将会起主导作用；因此，上式为广义的运动方程，将达西定律变为：
其中，
式中，δ为紊流修正系数；
水合物与浅层气、深层气联合开采过程中，水相相对渗透率krw、气相相对渗透率krg和毛细管压力Pc的计算式分别为：
其中，
式中，Swr为束缚水饱和度，Sgr为残余气饱和度，krw0和krg0为渗透率端点值，Pc0是毛管压力端点值，m为van Genuchten参数。
本发明第二方面提出了天然气水合物与浅层气、深层气联合开采模拟系统，包括：
天然气水合物与浅层气联合开采模拟模块，用于设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数，构建联合开采模拟模型，求解模型获取产能数据；
天然气水合物与浅层气、深层气联合开采模块，在天然气水合物与浅层气联合开采模拟模块的基础上，水合物层下方增加深层气层，并对深层气层进行网格划分，设定地质参数、生产参数，构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据。
本发明第三方面提出了一种计算机设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述任一项所述方法的步骤。
本发明第四方面提出了一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现上述任一项所述方法的步骤。
本发明的有益效果是，
本发明提供的天然气水合物与浅层气、深层气联合开采模拟方法，基于渗流场、温度场、力学场相互作用机制进行三场耦合流动模拟，分别探明了两气合采和三气合采的产能变化规律。构建两气合采与三气合采数值模拟模型，可以有效阐明多气合采储层物性参数时变规律及演化特征，明确不同开采阶段储层物性参数变化规律，为实现海洋天然气水合物多气合采提供理论与技术支撑。
附图说明
图1是本发明实施例提供的天然气水合物与浅层气、深层气联合开采模拟方法和系统流程图；
图2是本发明实例构建的浅层气层与水合物层联合开采模拟模型示意图；
图3是本发明实例井底流压为3.75MPa时浅层气与水合物联合开采的累产气与产气速率随时间变化的规律图；
图4是本发明实例构建的浅层气层、深层气层与水合物层联合开采模拟模型示意图；
图5是本发明实例井底流压为3.75MPa时浅层气、深层气与水合物联合开采的累产气与产气速率随时间变化的规律图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
天然气水合物与浅层气、深层气联合开采模拟方法，如图1所示，包括以下步骤：
S1、构建浅层气层与水合物层联合开采模拟模型，包括设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数等，构建联合开采模拟模型，求解模型获取产能数据；
水合物分解动力学方程为：
式中，Nh为水合物结合的水分子数；
水合物层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，t为时间，s；κ为组分标识，本方程中表示水合物组分h、甲烷组分m、水组分w或能量e；Mκ为κ的所有组分的和，kg/m3或J/m3；Fκ为κ的可流动的组分，kg/(m2·s)；qκ为κ的源汇，kg/(m3·s)或J/(m3·s)；
水合物组分质量守恒方程为：
Mh＝φSHρH；
式中，Mh为水合物组分的质量之和，kg/m3；SH为水合物相的饱和度；ρH为水合物相的密度，kg/m3；
甲烷组分质量守恒方程为：
式中，Mm为甲烷组分的质量之和，kg/m3；φ为储层孔隙度；SA为水相的饱和度；SG为气相的饱和度；ρA为水相的密度，kg/m3；ρG为气相的密度，kg/m3；为甲烷组分的质量与水相的质量之比；/>为甲烷组分的质量与气相的质量之比；/>为甲烷组分的质量与水合物相的质量之比；FA为水相的在单位截面积通过的质量流量，kg/(m2·s)；FG为气相在单位截面积通过的质量流量，kg/(m2·s)；Fm为甲烷组分在单位截面积通过的质量流量，kg/(m2·s)；qm为甲烷组分的源汇，kg/(m3·s)；qA为水相的源汇，kg/(m3·s)；qG为气相的源汇，kg/(m3·s)；/>为甲烷组分的质量与水相源汇的质量之比；/>为甲烷组分的质量与气相源汇的质量之比；
水组分质量守恒方程为：
式中，Mw为水组分的质量之和，kg/m3；为水组分的质量与水相的质量之比；/>为水组分的质量与气相的质量之比；/>为水组分的质量与水合物相的质量之比；Fw为水组分在单位截面积通过的质量流量，kg/(m2·s)；qw为水组分的源汇，kg/(m3·s)；/>为水组分的质量与水相源汇的质量之比；/>为水组分的质量与气相源汇中质量之比；
能量守恒方程为：
式中，Me为能量之和，J/m3；ρR为岩石的密度，kg/m3；HR为岩石的焓，J/kg；sβ为β相的饱和度；Hβ为β相的焓，J/kg；ΔsH为当前时间步内水合物饱和度的变化值；ΔH0为水合物的分解/形成焓，J/kg；Fe为能量流速，J/(m2·s)；Kc为体系的综合导热系数，W/(m·K)；Fβ为β相在单位截面积通过的质量流量，kg/(m2·s)；qe为能量的源汇，J/(m3·s)。
S2、在步骤S1的基础上，在水合物层下方增加深层气层，并进一步对深层气层进行网格划分，设定地质参数、生产参数等，之后构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据。
浅层气层、深层气层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，x为气组分g或能量e；Mx为x的所有组分的和，kg/m3或J/m3；Fx为x的可流动的组分，kg/(m2·s)；qx为x的源汇，kg/(m3·s)或J/(m3·s)；
气组分质量守恒时为：
Mg＝φSGρGXG；
Fg＝XGFG；
qg＝Xq,GqG；
式中，Mg为气组分的质量之和，kg/m3；XG为气组分的质量与气相的质量之比；Fg为气组分在单位截面积通过的质量流量，kg/(m2·s)；qg为气组分的源汇，kg/(m3·s)；Xq,G为气组分的质量与气相源汇的质量之比；
能量守恒时为：
Me＝(1-φ)ρRHR+φsGρGHG；
qe＝HGqG；
式中，sG为气相的饱和度；HG为气相的焓，J/kg；FG为气相在单位截面积通过的质量流量，kg/(m2·s)。
深水气层温度较低，实际气体压缩性与理想气体存在较大差异，因此在气体状态方程中引入偏差因子：
pV＝nZRT；
式中，p为气体的绝对压力，Mpa；T为气体的绝对温度，K；V为气体的体积，m3；n为气体的物质的量，mol；R为气体常数，8.314×10-3Mpa/(mol·K)；Z为气体的偏差因子，无量纲；
气体渗流与液体渗流相类似，当气体处于层流状态时，采用达西渗流定律描述其流动状态，对于均质地层中，在三维渗流空间中，广义达西定律为：
气体流动在三维空间中的流动速度表示为：
式中，v为气体渗流速度，m/s；K为地层渗透率，D；μ为气体粘度，mPa·s；g为重力加速度，g/cm3；x、y、z为空间坐标轴；
当气体的渗流速度增加到一定程度之后，紊流和惯性的影响将会愈发显著，其渗流速度与压力梯度呈非线性关系，不满足达西渗流定律；在水平方向上，气体渗流过程中有紊流和惯性阻力存在时，满足其动力学规律的非线性二次项式运动方程为：
式中，ζ为影响紊流和惯性阻力的孔隙结构特征参数；
上式右端第一项为黏滞阻力，与渗流速度成正比，第二项为惯性阻力，与渗流速度平方成正比，所以当渗流速度较小之时，可以不考虑惯性阻力的影响；当渗流速度较大时，紊流和惯性的影响逐渐显著，当渗流速度增加到偏离达西定律时，惯性将会起主导作用；因此，上式为广义的运动方程，将达西定律变为：
其中，
式中，δ为紊流修正系数；
水合物与浅层气、深层气联合开采过程中，水相相对渗透率krw、气相相对渗透率krg和毛细管压力Pc的计算式分别为：
其中，
式中，Swr为束缚水饱和度，Sgr为残余气饱和度，krw0和krg0为渗透率端点值，Pc0是毛管压力端点值，m为van Genuchten参数。
本发明实施例还提供了一种天然气水合物与浅层气、深层气联合开采模拟系统，包括：
天然气水合物与浅层气联合开采模拟模块，用于设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数，构建联合开采模拟模型，求解模型获取产能数据；
天然气水合物与浅层气、深层气联合开采模块，在天然气水合物与浅层气联合开采模拟模块的基础上，水合物层下方增加深层气层，并对深层气层进行网格划分，设定地质参数、生产参数，构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据。
本发明实施例还提出了一种计算机设备，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述方法的步骤。
本发明实施例还提出了一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现上述方法的步骤。本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本发明所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory，ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random AccessMemory，SRAM)或动态随机存取存储器(Dynamic Random Access Memory，DRAM)等。
下面结合实施例对本发明的技术方案作进一步的描述。
如图2所示，该发明实例仅考虑浅层气层和天然气水合物层，模型的上、下盖层是不可渗透的。海底水平面距离海面水平面的深度为1200m，海底处压力取11760kPa，海底温度取3.8℃。海底沉积物密度取2600kg/m3，沉积物比热容为1000J/(kg·℃)。同时，将天然气水合物层距离海面的深度设置为1530m，将浅层气层距离海面深度设置为1700m。
模型大小为560m×560m×285m，长宽都为560m，厚度为285m。XY方向上划分为56个网格，纵向上划分为35个网格。模型共有4层，分别为上下盖层，水合物层和浅层气层。上下盖层厚度均为40m，划分为4个网格，水合物层厚度为15m，划分为5个网格，浅层气层厚度为35m，划分为5个网格。天然气水合物层与浅层气层之间为170m的其他各类地质类型，将其设置为无效网格，忽略其地层参数对结果的影响。模型中央一口直井进行降压生产，从浅层气层顶部钻至水合物层底部。按地层厚度划分的地质模型如图2所示，图中水合物层和浅层气层之间的区域为无效网格区。
水合物层厚度为15m，包括水，天然气水合物和甲烷气，有效孔隙度为33％，水合物饱和度为31％，渗透率为1.5mD。浅层气层厚度为35m，孔隙度取25％，渗透率取6mD，储层中包含甲烷相与气相。
图3展示了井底流压为3.75MPa时浅层气与天然气水合物联合开采的累产气与产气速率随时间的变化规律图，从图中可以看出产气速率为一条先上升后快速下降，之后平缓下降的曲线，而累产气量为一条先快速上升后平缓上升的曲线，因此生产过程可分为两个阶段：①自由气大量产出阶段。投产最初一段时间获得极高的产气速率后快速下降，最大产气速率达到了1.83×106m3/d。在此阶段内，浅层气层内存在大量自由气，故初始时刻的产气速率与产气量突增。之后由于储层的自由气的大量产出，储层压力降低，水合物开始分解，产出的气体包括自由气和水合物分解产生的气体两部分。该阶段的产出气以自由气为主，开采初期能获得很大的产气速率，而后续由于自由气大量减少，产气速率迅速下降，累产气量开始逐渐变得平缓；②水合物分解产气阶段。水合物开始分解产生分解气，而产气速率主要受储层中天然气水合物的分解速率影响。由于浅层气的大量采出导致气体压力快速下降，导致井底压力与地层压力之差快速增大，促进水合物的分解，使水合物分解速度加快。然而，由于水合物分解吸热导致了水合物层温度降低，在没有外界热量的补充的情况下，水合物分解速率逐渐下降，产气速率也一直呈下降趋势，但下降速率缓慢，基本保持在稳定值3.03×104m3/d左右，有较好的商业开采价值。
如图4所示，该发明实例考虑浅层气层和天然气水合物层和深层气层，模型的上、下盖层是含水饱和层。将天然气水合物层距离海面的深度设置为1530m，将浅层气距离海面的深度设置为1700m，将深层气距离海面的深度设置为4500m。
模型大小为560m×560m×3077m，长宽都为560m，厚度为3077m。XY方向上划分为56个网格，纵向上划分为131个网格。模型共有5层，分别为上下盖层，浅层气层，水合物层和深层气层。上下盖层厚度均为40m，划分为4个网格，水合物层厚度为15m，划分为5个网格，浅层气层厚度为35m，划分为5个网格，深层气层为27m，划分为3个网格。浅层气层与天然气水合物层之间为155m以及天然气水合物层与深层气层之间2765m的其他各类地质类型，将其设置为无效网格，忽略其地层参数对结果的影响。模型中央一口直井进行降压生产，从浅层气层顶部钻至深层气层底部。按地层厚度划分的地质模型如图4所示，图中水合物层和浅层气层之间以及浅层气层和深层气层之间的区域为无效网格区。
水合物层厚度为15m，包括水和天然气水合物和甲烷气，有效孔隙度为33％，水合物饱和度为31％，渗透率为1.5mD。浅层气层厚度为35m，孔隙度取25％，渗透率取6mD，储层中包含甲烷气。深层气层厚度为27m，包含甲烷气，有效孔隙度32％，渗透率为7.4mD。
图5展示了井底流压为3.75MPa时多气合采模型累产气与产气速率随时间的变化规律图，从图中可以看出产气速率为一条先上升后快速下降，之后平缓下降的曲线，而累产气量为一条先快速上升后平缓上升的曲线，因此生产过程可分为两个阶段：①自由气大量产出阶段。生产最初一段时间获得极高的产气速率后快速下降，最大产气速率达到了2.82×107m3/d。这一阶段，浅层气层和深层气层中自由气含量较高，因此初始时刻产气速率和产气量都发生了突增。开采气主要来自于储层中的自由气以及由水合物转化而来的甲烷气。这一阶段产出气主要为自由气，在采出程度前期可取得较大产气速率，随后因自由气的大幅降低，产气速率快速降低，累产气量也开始渐趋平缓。②水合物分解产气阶段。水合物开始分解生成甲烷气，产气速率则主要受到储层内天然气水合物分解速率的影响。由于浅层气和深层气大量开采造成气体压力迅速降低，致使井底压力和地层压力之差迅速增加，促使水合物分解并加速其分解。因此，为了研究水合物分解对气井产能的影响，需要准确测定不同条件下产气量随时间的变化规律。而水合物分解吸热造成水合物层内温度下降，无外界热量补充时水合物分解速率会逐渐减小，产气速率始终呈现减小趋势，只是减小速率较慢，基本保持在稳定值1.81×104m3/d左右。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (4)
1.天然气水合物与浅层气、深层气联合开采模拟方法，其特征在于，包括以下步骤：
S1、构建浅层气层与水合物层联合开采模拟模型，包括设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数，构建联合开采模拟模型，求解模型获取产能数据；
S2、在步骤S1的基础上，水合物层下方增加深层气层，并对深层气层进行网格划分，设定地质参数、生产参数，构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据；
步骤S1中，水合物分解动力学方程为：
式中，Nh为水合物结合的水分子数；
水合物层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，t为时间，s；κ为组分标识，本方程中表示水合物组分h、甲烷组分m、水组分w或能量e；Mκ为κ的所有组分的和，kg/m3或J/m3；Fκ为κ的可流动的组分，kg/(m2·s)；qκ为κ的源汇，kg/(m3·s)或J/(m3·s)；
水合物组分质量守恒方程为：
Mh＝φSHρH；
式中，Mh为水合物组分的质量之和，kg/m3；SH为水合物相的饱和度；ρH为水合物相的密度，kg/m3；
甲烷组分质量守恒方程为：
式中，Mm为甲烷组分的质量之和，kg/m3；φ为储层孔隙度；SA为水相的饱和度；SG为气相的饱和度；ρA为水相的密度，kg/m3；ρG为气相的密度，kg/m3；为甲烷组分的质量与水相的质量之比；/>为甲烷组分的质量与气相的质量之比；/>为甲烷组分的质量与水合物相的质量之比；FA为水相的在单位截面积通过的质量流量，kg/(m2·s)；FG为气相在单位截面积通过的质量流量，kg/(m2·s)；Fm为甲烷组分在单位截面积通过的质量流量，kg/(m2·s)；qm为甲烷组分的源汇，kg/(m3·s)；qA为水相的源汇，kg/(m3·s)；qG为气相的源汇，kg/(m3·s)；/>为甲烷组分的质量与水相源汇的质量之比；/>为甲烷组分的质量与气相源汇的质量之比；
水组分质量守恒方程为：
式中，Mw为水组分的质量之和，kg/m3；为水组分的质量与水相的质量之比；/>为水组分的质量与气相的质量之比；/>为水组分的质量与水合物相的质量之比；Fw为水组分在单位截面积通过的质量流量，kg/(m2·s)；qw为水组分的源汇，kg/(m3·s)；/>为水组分的质量与水相源汇的质量之比；/>为水组分的质量与气相源汇中质量之比；
能量守恒方程为：
式中，Me为能量之和，J/m3；ρR为岩石的密度，kg/m3；HR为岩石的焓，J/kg；sβ为β相的饱和度；Hβ为β相的焓，J/kg；ΔsH为当前时间步内水合物饱和度的变化值；ΔH0为水合物的分解/形成焓，J/kg；Fe为能量流速，J/(m2·s)；Kc为体系的综合导热系数，W/(m·K)；Fβ为β相在单位截面积通过的质量流量，kg/(m2·s)；qe为能量的源汇，J/(m3·s)；
步骤S2中，浅层气层、深层气层遵循质量守恒方程和能量守恒方程，体系守恒关系表示为：
式中，x为气组分g或能量e；Mx为x的所有组分的和，kg/m3或J/m3；Fx为x的可流动的组分，kg/(m2·s)；qx为x的源汇，kg/(m3·s)或J/(m3·s)；
气组分质量守恒时为：
Mg＝φSGρGXG；
Fg＝XGFG；
qg＝Xq,GqG；
式中，Mg为气组分的质量之和，kg/m3；XG为气组分的质量与气相的质量之比；Fg为气组分在单位截面积通过的质量流量，kg/(m2·s)；qg为气组分的源汇，kg/(m3·s)；Xq,G为气组分的质量与气相源汇的质量之比；
能量守恒时为：
Me＝(1-φ)ρRHR+φsGρGHG；
Fe＝-Kc▽T+HGFG；
qe＝HGqG；
式中，sG为气相的饱和度；HG为气相的焓，J/kg；FG为气相在单位截面积通过的质量流量，kg/(m2·s)；
步骤S2中，气体渗流规律如下：
在气体状态方程中引入偏差因子：
pV＝nZRT；
式中，p为气体的绝对压力，Mpa；T为气体的绝对温度，K；V为气体的体积，m3；n为气体的物质的量，mol；R为气体常数，8.314×10-3Mpa/(mol·K)；Z为气体的偏差因子，无量纲；
气体渗流与液体渗流相类似，当气体处于层流状态时，采用达西渗流定律描述流动状态，对于均质地层中，在三维渗流空间中，广义达西定律为：
气体流动在三维空间中的流动速度表示为：
式中，v为气体渗流速度，m/s；K为地层渗透率，D；μ为气体粘度，mPa·s；g为重力加速度，g/cm3；x、y、z为空间坐标轴；
当气体的渗流速度增加到渗流速度与压力梯度呈非线性关系，不满足达西渗流定律；在水平方向上，气体渗流过程中有紊流和惯性阻力存在时，满足其动力学规律的非线性二次项式运动方程为：
式中，ζ为影响紊流和惯性阻力的孔隙结构特征参数；
当渗流速度增加到偏离达西定律时，达西定律变为：
其中，
式中，δ为紊流修正系数；
水合物与浅层气、深层气联合开采过程中，水相相对渗透率krw、气相相对渗透率krg和毛细管压力Pc的计算式分别为：
其中，
式中，Swr为束缚水饱和度，Sgr为残余气饱和度，krw0和krg0为渗透率端点值，Pc0是毛管压力端点值，m为van Genuchten参数。
2.采用权利要求1所述方法的天然气水合物与浅层气、深层气联合开采模拟系统，其特征在于，包括：
天然气水合物与浅层气联合开采模拟模块，用于设定浅层气层与水合物层，并进行网格划分，设定地质参数、生产参数及井控参数，构建联合开采模拟模型，求解模型获取产能数据；
天然气水合物与浅层气、深层气联合开采模块，在天然气水合物与浅层气联合开采模拟模块的基础上，水合物层下方增加深层气层，并对深层气层进行网格划分，设定地质参数、生产参数，构建浅层气层、深层气层与水合物层联合开采模拟模型，求解模型获取产能数据。
3.一种计算机设备，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现权利要求1所述方法的步骤。
4.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述指令被执行时实现权利要求1所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311515925.5A CN117236232B (zh) | 2023-11-15 | 2023-11-15 | 天然气水合物与浅层气、深层气联合开采模拟方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311515925.5A CN117236232B (zh) | 2023-11-15 | 2023-11-15 | 天然气水合物与浅层气、深层气联合开采模拟方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117236232A CN117236232A (zh) | 2023-12-15 |
| CN117236232B true CN117236232B (zh) | 2024-02-20 |
Family
ID=89089784
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311515925.5A Active CN117236232B (zh) | 2023-11-15 | 2023-11-15 | 天然气水合物与浅层气、深层气联合开采模拟方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117236232B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108241777A (zh) * | 2017-12-27 | 2018-07-03 | 青岛海洋地质研究所 | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 |
| CN111236908A (zh) * | 2020-01-09 | 2020-06-05 | 西南石油大学 | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 |
| CN111353242A (zh) * | 2020-04-26 | 2020-06-30 | 青岛杰瑞工控技术有限公司 | 一种可燃冰开采时甲烷监测传感器的布置方法 |
| CN111411946A (zh) * | 2020-05-12 | 2020-07-14 | 中国石油大学(北京) | 一种致密气藏气井开采方式的确定方法及装置 |
| CN111665546A (zh) * | 2019-03-07 | 2020-09-15 | 中普宝信(北京)科技有限公司 | 用于可燃冰探测的声学参数获取方法 |
| CN112855129A (zh) * | 2021-01-15 | 2021-05-28 | 青岛海洋地质研究所 | 海洋天然气水合物系统多气合采模拟装置及方法 |
| CN113327647A (zh) * | 2021-04-13 | 2021-08-31 | 中国石油大学(华东) | 一种模拟水合物藏中不同来源甲烷运移轨迹的方法 |
| CN114692472A (zh) * | 2022-06-01 | 2022-07-01 | 中国石油大学(华东) | 储层改造辅助降压开发天然气水合物藏的数值模拟方法 |
| CN114776264A (zh) * | 2022-03-28 | 2022-07-22 | 中国科学院广州能源研究所 | 一种天然气水合物开采过程中的固相控制方法 |
| CN115408889A (zh) * | 2022-11-01 | 2022-11-29 | 中国石油大学(华东) | 一种注热-压裂-降压联合开采海域天然气水合物的方法 |
| CN116167302A (zh) * | 2023-04-18 | 2023-05-26 | 中国石油大学(华东) | 一种天然气水合物增产模拟中人工复杂裂缝的描述方法 |
| CN116291333A (zh) * | 2023-02-15 | 2023-06-23 | 西南石油大学 | 一种海洋天然气水合物、浅层气、深部气合采系统及方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0412060D0 (en) * | 2004-05-28 | 2004-06-30 | Univ Newcastle | Process for stimulating production of methane from petroleum in subterranean formations |
| CN114706118A (zh) * | 2022-04-02 | 2022-07-05 | 中国石油大学(华东) | 海域天然气水合物体系地震响应特征的正演物理模拟方法 |
-
2023
- 2023-11-15 CN CN202311515925.5A patent/CN117236232B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108241777A (zh) * | 2017-12-27 | 2018-07-03 | 青岛海洋地质研究所 | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 |
| CN111665546A (zh) * | 2019-03-07 | 2020-09-15 | 中普宝信(北京)科技有限公司 | 用于可燃冰探测的声学参数获取方法 |
| CN111236908A (zh) * | 2020-01-09 | 2020-06-05 | 西南石油大学 | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 |
| CN111353242A (zh) * | 2020-04-26 | 2020-06-30 | 青岛杰瑞工控技术有限公司 | 一种可燃冰开采时甲烷监测传感器的布置方法 |
| CN111411946A (zh) * | 2020-05-12 | 2020-07-14 | 中国石油大学(北京) | 一种致密气藏气井开采方式的确定方法及装置 |
| CN112855129A (zh) * | 2021-01-15 | 2021-05-28 | 青岛海洋地质研究所 | 海洋天然气水合物系统多气合采模拟装置及方法 |
| CN113327647A (zh) * | 2021-04-13 | 2021-08-31 | 中国石油大学(华东) | 一种模拟水合物藏中不同来源甲烷运移轨迹的方法 |
| CN114776264A (zh) * | 2022-03-28 | 2022-07-22 | 中国科学院广州能源研究所 | 一种天然气水合物开采过程中的固相控制方法 |
| CN114692472A (zh) * | 2022-06-01 | 2022-07-01 | 中国石油大学(华东) | 储层改造辅助降压开发天然气水合物藏的数值模拟方法 |
| CN115408889A (zh) * | 2022-11-01 | 2022-11-29 | 中国石油大学(华东) | 一种注热-压裂-降压联合开采海域天然气水合物的方法 |
| CN116291333A (zh) * | 2023-02-15 | 2023-06-23 | 西南石油大学 | 一种海洋天然气水合物、浅层气、深部气合采系统及方法 |
| CN116167302A (zh) * | 2023-04-18 | 2023-05-26 | 中国石油大学(华东) | 一种天然气水合物增产模拟中人工复杂裂缝的描述方法 |
Non-Patent Citations (5)
| Title |
|---|
| 基于数值模拟的天然气水合物和浅层气合采产能影响因素研究;温慧芸;天然气技术与经济;第16卷(第5期);第1-3节 * |
| 杨宝善.凝析气藏开发工程.石油工业出版社,1995,第6-7页. * |
| 王启立 等.煤基多孔介质的表征、流动及应用.中国矿业大学出版社,2020,第107页. * |
| 碳达峰目标下中国深海天然气水合物开发战略;魏纳 等;天然气工业;第42卷(第2期);第3.4节，表4，图6 * |
| 考虑井筒与地层耦合的低渗气藏气井动态模拟;彭希宁;中国优秀硕士学位论文全文数据库 工程科技I辑(第06期);第2.4.2节、3.3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117236232A (zh) | 2023-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| You et al. | Mechanisms of methane hydrate formation in geological systems | |
| Janicki et al. | Simulation of Methane Recovery from Gas Hydrates Combined with Storing Carbon Dioxide as Hydrates. | |
| Burwicz et al. | Estimation of the global amount of submarine gas hydrates formed via microbial methane formation based on numerical reaction-transport modeling and a novel parameterization of Holocene sedimentation | |
| Bhatnagar et al. | Generalization of gas hydrate distribution and saturation in marine sediments by scaling of thermodynamic and transport processes | |
| CN107133452B (zh) | 油藏渗流数值模拟方法及装置 | |
| CN108133080B (zh) | 一种考虑非稳态吸附的非均质裂缝页岩气藏数值模拟方法 | |
| Liang et al. | Effects of vertical center well and side well on hydrate exploitation by depressurization and combination method with wellbore heating | |
| CN113327647B (zh) | 一种模拟水合物藏中不同来源甲烷运移轨迹的方法 | |
| CN110952976B (zh) | 一种气藏开发模式下单井开采稳产潜力评价方法 | |
| Luboń et al. | The influence of the first filling period length and reservoir level depth on the operation of underground hydrogen storage in a deep aquifer | |
| CN115034489B (zh) | 一种考虑溶解作用的气藏co2埋存潜力预测方法 | |
| Gu et al. | Numerical analysis on gas production from silty hydrate reservoirs in the South China sea by depressurizing: The effect of permeability reduction caused by pore compression | |
| CN105464634A (zh) | 一种利用埋存二氧化碳开采甲烷水合物的方法 | |
| CN112507279A (zh) | 一种探究浅层咸水越流补给深层淡水的方法 | |
| Vishal et al. | Sensitivity analysis of methane hydrate bearing Class 3 reservoirs during thermal injection | |
| Gambelli et al. | Kinetic considerations and formation rate for carbon dioxide hydrate, formed in presence of a natural silica-based porous medium: How initial thermodynamic conditions may modify the process kinetic | |
| CN117236232B (zh) | 天然气水合物与浅层气、深层气联合开采模拟方法和系统 | |
| Xiao et al. | Numerical analysis of production behaviors and permeability characteristics on the second gas hydrate production test in the South China sea | |
| Zhao et al. | Mathematical model and simulation of gas hydrate reservoir decomposition by depressurization | |
| Nikolaev et al. | Investigation of underground bio-methanation using bio-reactive transport modeling | |
| Kawada et al. | Formation of a hydrothermal reservoir due to anhydrite precipitation in an arc volcano hydrothermal system | |
| CN113010568B (zh) | 页岩吸附气产出比例的确定方法、装置、存储介质及设备 | |
| CN115273994A (zh) | 一种天然气水合物岩芯分解和输运物性的预测方法 | |
| Guo et al. | Numerical simulation of depressurization production of natural gas hydrate in different well types | |
| CN111101930B (zh) | 一种气藏开发模式下单井开采提产潜力评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |