CN115032337B - 一种盐水预抽量和回注量的计算方法及场地试验验证方法 - Google Patents
一种盐水预抽量和回注量的计算方法及场地试验验证方法 Download PDFInfo
- Publication number
- CN115032337B CN115032337B CN202210561742.6A CN202210561742A CN115032337B CN 115032337 B CN115032337 B CN 115032337B CN 202210561742 A CN202210561742 A CN 202210561742A CN 115032337 B CN115032337 B CN 115032337B
- Authority
- CN
- China
- Prior art keywords
- brine
- layer
- reinjection
- pressure
- saline
- 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
- 239000012267 brine Substances 0.000 title claims abstract description 343
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 title claims abstract description 343
- 238000005086 pumping Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 87
- 238000012360 testing method Methods 0.000 title claims abstract description 70
- 238000004364 calculation method Methods 0.000 title claims abstract description 23
- 238000012795 verification Methods 0.000 title claims abstract description 21
- 238000002347 injection Methods 0.000 claims abstract description 101
- 239000007924 injection Substances 0.000 claims abstract description 101
- 238000007789 sealing Methods 0.000 claims abstract description 62
- 238000000605 extraction Methods 0.000 claims abstract description 57
- 230000008569 process Effects 0.000 claims abstract description 27
- 230000009919 sequestration Effects 0.000 claims abstract description 18
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 143
- 239000011780 sodium chloride Substances 0.000 claims description 143
- 239000011148 porous material Substances 0.000 claims description 90
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 85
- 239000012530 fluid Substances 0.000 claims description 34
- 230000015572 biosynthetic process Effects 0.000 claims description 28
- 238000001764 infiltration Methods 0.000 claims description 16
- 230000008595 infiltration Effects 0.000 claims description 16
- 238000012544 monitoring process Methods 0.000 claims description 15
- 229910000831 Steel Inorganic materials 0.000 claims description 14
- 239000010959 steel Substances 0.000 claims description 14
- 239000007788 liquid Substances 0.000 claims description 11
- 230000008901 benefit Effects 0.000 claims description 8
- 238000002474 experimental method Methods 0.000 claims description 8
- 150000003839 salts Chemical class 0.000 claims description 8
- 230000008961 swelling Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 229910001566 austenite Inorganic materials 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 101150054854 POU1F1 gene Proteins 0.000 description 7
- 230000035699 permeability Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 239000004568 cement Substances 0.000 description 2
- 239000004570 mortar (masonry) Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/02—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B9/00—Simulators for teaching or training purposes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N2013/003—Diffusion; diffusivity between liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N2015/0007—Investigating dispersion of gas
- G01N2015/0015—Investigating dispersion of gas in solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N2015/0813—Measuring intrusion, e.g. of mercury
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/70—Combining sequestration of CO2 and exploitation of hydrocarbons by injecting CO2 or carbonated water in oil wells
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Combustion & Propulsion (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明属于CO2地质封存领域,涉及一种最大化CO2封存效率的盐水预抽量和回注量计算方法及场地试验验证方法,并提供了一种提高CO2封存效率的预抽‑回注方法,分别将计算得到的盐水预抽量和回注量作为CO2封存过程中实际盐水抽取和回注的控制条件,同时建立了CO2注入速度与盐水抽取速度关系、盐水回注速度与CO2注入速度关系。通过本发明设计的一种最大化CO2封存效率的场地试验验证方法和装置,与场地试验对比验证了所述计算方法精度超过95%;通过计算盐水预抽量和回注量实施盐水预抽和回注策略,与常规封存试验对比封存能力提高25%、封存安全提高20%。
Description
技术领域
本发明属于CO2地质封存领域,涉及一种最大化CO2封存效率的盐水预抽量和回注量计算方法及场地试验验证方法。
背景技术
CO2地质封存已成为当前研究热点,特别是最大化CO2地质封存效率已成为重中之重。CO2封存效率密切相关于封存能力和封存安全。为了提高CO2封存能力,可以在CO2注入盐水层前预先抽取盐水增大CO2储存空间。然而,预抽盐水将降低地层孔隙压力,极易诱发地层坍塌。另外,在CO2注入后,地层孔隙压力迅速抬升从低压(盐水预抽引起)转为高压状态,又将引起地层膨胀甚至破裂。为了应对地层破裂问题,可以将抽取的盐水回注到上部毗邻的盐水层,缓解下部地层孔隙压力抬升引起的地层膨胀变形。同时,预抽盐水的再次回注还可以省去复杂的盐水处理流程、节省高额的盐水处理费用。综上分析,盐水预抽和回注可以提高CO2地质封存能力和封存安全,然而实施这一策略需要计算盐水预抽量和回注量。盐水预抽量由地层坍塌条件控制,而盐水回注量由地层破裂条件控制。鉴于此,迫切需要提出一种最大化CO2封存效率的盐水预抽量和回注量计算方法,并通过场地试验验证计算方法在实际应用中的可靠性。
目前,已有学者研发了少数实验装置模拟CO2地质封存,如一种实验模拟CO2地质封存可行性的***(CN202010239445.0)、一种CO2地层封存***(CN202020387776.4)、一种CO2地层封存方法及***(CN202010215234.3)。然而,上述实验装置主要模拟CO2在盐水层内的运移规律,难以评价盐水预抽、盐水回注对CO2封存效率的影响,更无法确定具体的盐水预抽量和盐水回注量。鉴于此,提出一种最大化CO2封存效率的盐水预抽量和回注量计算方法及场地试验验证已势在必行。
发明内容
本发明提供了一种最大化CO2封存效率的盐水预抽量和回注量计算方法及场地试验验证方法,以地层坍塌和地层破裂为条件计算盐水预抽量和盐水回注量。
为实现上述目的,本发明一方面提供一种最大化CO2封存效率的盐水预抽量和回注量计算方法,其计算步骤是:
(1)计算底部盐水层盐水预抽量V预抽
盐水预抽量指在CO2注入前的盐水抽取量;CO2注入后,盐水抽取仍会持续一段时间,但该段时间内的盐水抽取量不计入盐水预抽量内;
式中,Pw c、为底部盐水层坍塌压力临界值和许用值,MPa;σH、σh分别为最大、最小水平主应力,由施加的水平围压获得,MPa;Pp为盐水层初始地层孔隙压力,由所施加地层孔隙压力获得,MPa;/>为内摩擦角,通过多组三轴压缩实验的莫尔包络线获得;c为内聚力,通过多组三轴压缩实验的莫尔包络线获得,MPa;α为比奥系数,一般为0.85;n1为安全系数,无量纲,结合工程经验设定;
式中,ΔPp预抽为盐水预抽引起的底部盐水层孔隙压力变化,MPa;H预抽为盐水预抽引起的水位高度变化,m;γ盐水为盐水重度,KN·m-3;V预抽、V渗分别为盐水预抽量、实际条件下盐水预抽引起的地下流体从无穷远处向预抽井渗流的补给量,m3;S为底部盐水层水平面积,m2;φ下为底部盐水层孔隙度,无量纲;
(2)计算顶部盐水层盐水回注量V回注
②确定由盐水回注引起的顶部盐水层孔隙压力增加ΔPp回注,以增大后的孔隙压力值Pp+ΔPp回注不高于破裂压力许用值[Pw f]为临界条件计算盐水回注量V回注;
式中,ΔPp回注为盐水回注引起的顶部盐水层孔隙压力变化,MPa;H回注为盐水回注引起的水位高度变化,m;V回注、V′渗分别为盐水回注量、实际条件下盐水回注引起的地下流体从回注井向无穷远处渗流的溢出量,m3;S′为顶部盐水层水平面积,m2;φ上为顶部盐水层孔隙度,无量纲;
(3)底部盐水层CO2注入与盐水抽取同步发生时,建立CO2注入速度v注与盐水抽取速度v抽取的关系:
式中,ΔPp注为CO2注入引起底部盐水层孔隙压力变化MPa;v注为CO2注入速度,m3·s-1;v抽取为盐水抽取速度,m3·s-1;v″渗-v渗为实际条件下底部盐水层地下流体从无穷远处向预抽井及注入井渗流的补给速度,m3·s-1;
(4)底部盐水层CO2注入与顶部盐水层盐水回注同步发生时,建立CO2注入速度v注与盐水回注速度v回注的关系:
式中,v回注为盐水回注速度,m3·s-1;v注为CO2注入速度,m3·s-1。
本发明另一方面提供一种提高CO2封存效率的预抽-回注方法,采用上述计算方法,具体包括以下步骤:
首先对底部盐水层进行盐水预抽;当盐水预抽量达到上述计算得到的盐水预抽量V预抽时,开始向底部盐水层注入CO2,同时由上述建立的CO2注入速度与盐水抽取速度关系控制CO2注入速度;当CO2检测器检测到CO2时,停止盐水抽取;当底部盐水层孔隙压力恢复到初始孔隙压力时,开始向顶部盐水层回注盐水,同时由上述建立的盐水回注速度与CO2注入速度关系控制盐水回注速度;当盐水回注量达到上述计算得到的盐水回注量V回注时,停止盐水回注;当底部盐水层孔隙压力达到破裂压力许用值时,停止CO2注入。
本发明还提供一种最大化CO2封存效率的场地试验验证方法,验证盐水预抽量和回注量计算方法,其操作步骤是:
步骤1:布设场地试验
布设CO2封存的场地试验装置;
步骤2:开展场地试验
(1)施加压力
通过液压支架及配重块给盐水层施加三向围压,通过液体增压器及承压钢管给盐水层施加初始地层孔隙压力;
(2)预抽盐水
在CO2注入前,通过预抽井从底部盐水层中抽取盐水并暂存于储存箱;在CO2注入后,当盐水预抽井井底处CO2检测器检测到CO2时,停止盐水抽取;
(3)注入CO2
当流体压力传感器监测到的底部盐水层孔隙压力达到坍塌压力许用值时,开始通过注入井向底部盐水层中注入CO2;当流体压力传感器监测到的底部盐水层孔隙压力达到破裂压力许用值时,停止CO2注入;
(4)回注盐水
当流体压力传感器监测到的底部盐水层孔隙压力值恢复到初始孔隙压力时,开始通过回注井将从底部盐水层中抽取的盐水回注到顶部盐水层中;当流体压力传感器监测到的顶部盐水层孔隙压力达到顶部盐水层破裂压力许用值时,停止盐水回注;
(5)结束试验
停止CO2注入,试验结束;记录CO2注入量、底部盐水层***变形值、底部盐水层孔隙压力监测值达到坍塌压力许用值时的盐水实际抽取量、顶部盐水层孔隙压力监测值达到破裂压力许用值时的盐水实际回注量;
步骤3:验证计算方法
(1)盐水预抽量验证
获取底部盐水层孔隙压力监测值达到坍塌压力许用值时的盐水实际抽取量,以盐水实际抽取量的等效体积与盐水预抽量理论值对比,若误差在5%以内,表明盐水预抽量计算方法可以满足现场工程;
(2)盐水回注量验证
获取顶部盐水层孔隙压力监测值达到破裂压力许用值时的盐水实际回注量,以盐水实际回注量的等效体积与盐水回注量理论值对比,若误差在5%以内,表明盐水回注量计算方法可以满足现场工程;
步骤4:对比封存效率
(1)通过步骤2中的(1)、(2)、(3)和(4)开展最大化CO2封存效率场地试验,通过步骤2中的(5)获取CO2注入量和底部盐水层最大***变形;
(2)通过步骤2中的(1)和(3)开展常规CO2封存场地试验,通过步骤2中的(5)获取CO2注入量和底部盐水层最大***变形;
(3)对比上述(1)和(2)中的CO2注入量,确认最大化CO2封存效率场地试验具有封存能力优势;对比(1)和(2)中的底部盐水层最大***变形,确认最大化CO2封存效率场地试验具有封存安全优势。
本发明中,CO2封存的场地试验装置中设置有盖层(盖层为CO2封存选址的必备条件),盖层上方的盐水层在本发明中称为顶部盐水层,盖层下方的盐水层在本发明中称为底部盐水层。
更进一步的,本发明中,所述CO2封存的场地试验装置,包括基坑、底部盐水层、盖层、顶部盐水层、位移传感器、流体压力传感器、防渗箱体、“回”字型基坑、液压支架、配重块、液体增压器、承压钢管、预抽井、储存箱、CO2检测器、定压泵、定压水箱、超临界CO2气罐、注入井、回注井;底部盐水层、盖层、顶部盐水层构成盐水层,盐水层外部包裹有防渗箱体,盐水层和防渗箱体均布置在基坑内,防渗箱体与基坑壁面间用土回填密实;在底部盐水层上表面布置位移传感器,监测盐水层变形;在基坑外四周开挖“回”字型基坑,在“回”字型基坑内放置液压支架,施加水平围压;在盐水层顶部布置配重块,施加垂向围压;液体增压器与盐水层内的承压钢管连接,负责将增压后的盐水注入盐水层,实现初始地层孔隙压力施加;预抽井设置在盐水层内,用于从底部盐水层中抽取盐水并暂存于储存箱;CO2检测器布置在预抽井井底处,用于判断CO2是否达到预抽井井底;注入井设置在盐水层中心处,并连接超临界CO2气罐,实现向底部盐水层注入CO2;回注井设置在顶部盐水层,负责将储存箱中存储的盐水回注到顶部盐水层;定压泵、定压水箱通过承压钢管与盐水层连通,由定压泵维持定压水箱内压力恒定且等于初始地层孔隙压力,模拟预抽盐水、CO2注入、回注盐水引起的实际条件下地下流体渗流过程(即无限大地层内孔隙压力的自调节功能);在预抽井、注入井及回注井的井底位置处布置流体压力传感器,监测盐水层地层孔隙压力。
本发明中,以地层坍塌和破裂为条件计算最大化CO2封存效率的盐水预抽量和回注量,以计算得到的盐水预抽量V预抽为实际盐水抽取过程中的控制条件,防止底部盐水层在预抽盐水过程中引起地层坍塌;以计算得到的回注量V回注为实际盐水回注过程中的控制条件,防止顶部盐水层在回注盐水过程中引起地层破裂;同时建立了CO2注入速度与盐水抽取速度关系,盐水回注速度与CO2注入速度关系。并且,本发明设计了一种最大化CO2封存效率的场地试验验证方法和装置,结构紧凑合理,原理可靠,可行性高。
与现有技术相比,本发明具有如下的有益效果:
1、本发明提出了一种最大化CO2封存效率的盐水预抽量和回注量计算方法,以地层坍塌和破裂为条件计算盐水预抽量和回注量,通过与场地试验对比验证计算精度超过95%;
2、本发明设计了一种最大化CO2封存效率的场地试验验证方法和装置,通过计算盐水预抽量和回注量实施盐水预抽和回注策略,与常规封存试验对比封存能力提高25%、封存安全提高20%。
附图说明
图1为本发明CO2封存的场地试验装置的结构示意图;
图2为本发明CO2封存的场地试验装置的俯视图;
图3为本发明盐水预抽量计算流程图;
图4为本发明盐水回注量计算流程图;
图5为本发明计算方法验证对比图;
图6为本发明实施例1结果对比图;
图中:1、基坑;2、底部盐水层;3、盖层;4、顶部盐水层;5、位移传感器;6、流体压力传感器;7、防渗箱体;8、“回”字型基坑;9、液压支架;10、配重块;11、液体增压器;12、承压钢管;13、预抽井;14、储存箱;15、CO2检测器;16、定压泵;17、定压水箱;18、超临界CO2气罐;19、注入井;20、回注井。
具体实施方式
结合附图和本发明具体实施方式的描述,能够更加清楚地了解本发明的细节。但是,在此描述的本发明的具体实施方式,仅用于解释本发明的目的,而不能以任何方式理解成是对本发明的限制。在本发明的教导下,技术人员可以构想基于本发明的任意可能的变形,这些都应被视为属于本发明的范围。
本实施方式一方面提供一种最大化CO2封存效率的盐水预抽量和回注量计算方法,其计算步骤是:
(1)计算底部盐水层2盐水预抽量V预抽
盐水预抽量指在CO2注入前的盐水抽取量;CO2注入后,盐水抽取仍会持续一段时间,但该段时间内的盐水抽取量不计入盐水预抽量内;
①根据弹性力学确定底部盐水层2预抽井13井壁围岩应力分布(见式(1)),并将径向有效应力σ′r、环向有效应力σ′θ代入Mohr-Coulomb准则确定底部盐水层2坍塌压力临界值,见式(2):
考虑安全要求,实际工程中的坍塌压力临界值通常用许用值代替,见式(3):
式中,Pw为底部盐水层2预抽井13井眼压力,MPa;Pw c、为底部盐水层2坍塌压力临界值和许用值,MPa;σH、σh分别为最大、最小水平主应力,由施加的水平围压获得,MPa;Pp为盐水层初始地层孔隙压力,由所施加地层孔隙压力获得,MPa;/>为内摩擦角,通过多组三轴压缩实验的莫尔包络线获得;c为内聚力,通过多组三轴压缩实验的莫尔包络线获得,MPa;α为比奥系数,一般为0.85;n1为安全系数,无量纲,结合工程经验设定;
根据水头压力计算公式,底部盐水层2孔隙压力下降ΔPp预抽表示为式(4):
式中,ΔPp预抽为盐水预抽引起的底部盐水层2孔隙压力变化,MPa;H预抽为盐水预抽引起的水位高度变化,m;γ盐水为盐水重度,KN·m-3;V预抽、V渗分别为盐水预抽量、实际条件下盐水预抽引起的地下流体从无穷远处向预抽井13渗流的补给量,m3;S为底部盐水层2水平面积,m2;φ下为底部盐水层2孔隙度,无量纲;
以V预抽为实际盐水抽取过程中的控制条件,防止底部盐水层2在预抽盐水过程中引起地层坍塌;
(2)计算顶部盐水层4盐水回注量V回注
①将式(1)中的径向有效应力σ′r、环向有效应力σ′θ代入最大拉应力准则确定地层破裂压力,见下式(6):
Pw f=3σh-σH-αPp+σt (6)
考虑安全要求,实际工程中的破裂压力临界值通常用许用值代替,见式(7):
根据水头压力计算公式,顶部盐水层4孔隙压力增加ΔPp回注表示为式(8):
式中,ΔPp回注为盐水回注引起的顶部盐水层4孔隙压力变化,MPa;H回注为盐水回注引起的水位高度变化,m;V回注、V′渗分别为盐水回注量、实际条件下盐水回注引起的地下流体从回注井20向无穷远处渗流的溢出量,m3;S′为顶部盐水层4水平面积,m2;φ上为顶部盐水层4孔隙度,无量纲;
以V回注为实际盐水回注过程中的控制条件,防止顶部盐水层4在回注盐水过程中引起地层破裂;
(3)建立CO2注入速度与盐水抽取速度关系
在抽取盐水与注入CO2同步发生时,以盐水抽取引起的底部盐水层2孔隙压力下降等于CO2注入引起的底部盐水层2孔隙压力增加为控制条件,建立CO2注入速度与盐水抽取速度的关系,见下式(10):
式中,ΔPp注为CO2注入引起底部盐水层2孔隙压力变化MPa;v注为CO2注入速度,m3·s-1;v抽取为盐水抽取速度,m3·s-1;v″渗-v渗为实际条件下底部盐水层2地下流体从无穷远处向预抽井13及注入井19渗流的补给速度,通过地层孔隙压力与定压水箱17的压力差及达西定律求得,见下式(11):
式中,k为地层渗透系数,通过渗透率实验获得,m·s-1;A为垂直于渗流方向的地层横截面积,m2;L为渗径长度,m;
(4)建立盐水回注速度与CO2注入速度关系
在注入CO2与回注盐水同步发生时,以CO2注入引起的底部盐水层2孔隙压力增加等于盐水回注引起的顶部盐水层4孔隙压力增加为控制条件,建立盐水回注速度与CO2注入速度关系,见下式(12):
式中,v回注为盐水回注速度,m3·s-1;v注为CO2注入速度,m3·s-1;v′渗为实际条件下地下流体从回注井20向无穷远处渗流的外溢速度,m3·s-1;v″渗为实际条件下地下流体从注入井19向无穷远处渗流的外溢速度,m3·s-1;
特别地,由“CO2注入引起的底部盐水层2孔隙压力增加等于盐水回注引起的顶部盐水层4孔隙压力增加”这一条件可以确定底部盐水层2、顶部盐水层4与定压水箱17之间的压差相等,进一步结合达西定律求得v′渗=v″渗,为此式(12)进一步简化为式(13):
本发明另一方面提供一种最大化CO2封存效率的场地试验,验证盐水预抽量和回注量计算方法,其操作步骤是:
步骤1:布设场地试验
如图1-2所示,所述一种最大化CO2封存效率的场地试验装置包括基坑1、底部盐水层2、盖层3、顶部盐水层4、位移传感器5、流体压力传感器6、防渗箱体7、“回”字型基坑8、液压支架9、配重块10、液体增压器11、承压钢管12、预抽井13、储存箱14、CO2检测器15、定压泵16、定压水箱17、超临界CO2气罐18、注入井19、回注井20;
开挖基坑1,体积为长3km×宽3km×深100m,在基坑1内部自下而上依次填充高度分别为50m、10m、40m的底部盐水层2、盖层3、顶部盐水层4,共同构成盐水层;底部盐水层2、顶部盐水层4的砂浆质量配比为水泥:砂子:水等于1:0.33:0.42,盖层3砂浆质量配比为水泥:砂子:水等于1:0.11:0.42;盐水层外部包裹有防渗箱体7(通过薄钢板制造),避免盐水在高压条件下从盐水层中渗出,防渗箱体7与基坑1壁面间用土回填密实;在底部盐水层2上表面每隔600m横纵布置位移传感器5,监测盐水层变形;
在基坑1外四周间隔50m开挖宽2m高100m的“回”字型基坑8,并在“回”字型基坑8内放置液压支架9,液压支架9底部固定在“回”字型基坑8外壁面(即远离基坑1的壁面),液压支架9顶部固定在“回”字型基坑8内壁面(即靠近基坑1的壁面),施加水平围压;在盐水层顶部每隔300m横纵布置配重块10,施加垂向围压;液体增压器11与盐水层内的承压钢管12连接,负责将增压后的盐水注入盐水层,实现初始地层孔隙压力施加;
注入井19一端连接超临界CO2气罐18,另一端从盐水层顶部中心处垂直深入底部盐水层2内20m,实现向底部盐水层2注入CO2;预抽井13设置在注入井19左右两侧各1km处且深入底部盐水层2内20m,用于从底部盐水层2中抽取盐水并暂存于储存箱14;CO2检测器15布置在预抽井13井底处,用于判断CO2是否达到预抽井13井底;回注井20设置在顶部盐水层4中心的左右两侧1km处且深入顶部盐水层4内15m,负责将储存箱14中存储的盐水回注到顶部盐水层4;定压泵16、定压水箱17通过承压钢管12与盐水层连通,由定压泵16维持定压水箱17内压力恒定且等于初始地层孔隙压力,用于模拟预抽盐水、CO2注入、回注盐水引起的实际条件下地下流体渗流过程(即无限大地层内孔隙压力的自调节功能);在预抽井13、注入井19及回注井20的井底位置处布置流体压力传感器6,监测盐水层地层孔隙压力;
步骤2:开展场地试验
(1)施加压力
通过液压支架9及配重块10给盐水层施加三向围压;通过液体增压器11给盐水增压,然后通过承压钢管12将增压后的盐水注入盐水层,实现初始地层孔隙压力施加;
(2)预抽盐水
通过预抽井13从底部盐水层中抽取盐水并暂存于储存箱14,根据工程经验设定盐水抽取速度v抽取;同时通过定压水箱17及承压钢管12模拟实际条件下底部盐水层由预抽盐水引起的地下流体从无穷远处向预抽井渗流的补给过程;
(3)注入CO2
当预抽井13井底位置处的流体压力传感器6监测到的底部盐水层2孔隙压力达到坍塌压力许用值时,开始通过注入井19向底部盐水层2中注入CO2,同时根据式(10)计算CO2注入速度v注;
(4)停止盐水抽取
当盐水预抽井13井底处设置的CO2检测器15检测到CO2时,停止盐水抽取;
(5)回注盐水
当注入井19井底位置处的流体压力传感器6监测到的底部盐水层2孔隙压力恢复到初始地层孔隙压力时,开始通过回注井20将从底部盐水层2中抽取的盐水回注到顶部盐水层4中,同时根据式(13)计算盐水回注速度v回注;
(6)停止盐水回注
当流体压力传感器6(回注井20井底位置处)监测到的顶部盐水层4孔隙压力达到破裂压力许用值时,停止盐水回注;
(7)停止CO2注入
当流体压力传感器6(注入井19井底位置处)监测到的底部盐水层2孔隙压力达到破裂压力许用值时,停止CO2注入;
(8)结束试验
停止CO2注入,试验结束;记录CO2注入量、底部盐水层2***变形值、底部盐水层2孔隙压力监测值达到坍塌压力许用值时的盐水实际抽取量V预抽-试验、顶部盐水层4孔隙压力监测值达到破裂压力许用值时的盐水实际回注量V回注-试验;
步骤3:验证计算方法
特别地,场地试验中盐水层的初始孔隙压力是通过液体增压器11对盐水增压实现的,而非实际条件下盐水依靠自身重力达到的初始孔隙压力,因此场地试验中注入盐水层的盐水体积(记为盐水增压体积V盐水-增压,由场地试验直接获得)远小于实际条件下产生相同孔隙压力的盐水体积(记为盐水原始体积,表示为或/>);鉴于此,应对场地试验中盐水的增压体积和实际条件下盐水的原始体积进行等效转换,以便与理论计算方法得到的盐水预抽量和盐水回注量进行比较;
(1)盐水预抽量验证
(2)盐水回注量验证
步骤4:对比封存效率
(1)通过步骤2中的(1)-(7)开展最大化CO2封存效率场地试验,通过步骤2中的(8)获取CO2注入量V1和底部盐水层2最大***变形u1;
(2)通过步骤2中的(1)、(3)和(7)开展常规CO2封存场地试验,通过步骤2中的(8)获取CO2注入量V2和底部盐水层2最大***变形u2;
(3)对比上述(1)和(2)中的CO2注入量V1、V2,确认最大化CO2封存效率场地试验具有封存能力优势;对比(1)和(2)中的底部盐水层2最大***变形u1、u2,确认最大化CO2封存效率场地试验具有封存安全优势。
实施例1:
根据步骤1布设场地试验;通过岩石强度实验获取底部盐水层2、顶部盐水层4抗拉强度σt为4MPa、内聚力c为9MPa、内摩擦角为40°;通过渗透率实验获取底部盐水层2、顶部盐水层4渗透率k为1×10-10m/s、孔隙度φ为0.15;根据步骤2中的(1)设定水平最大、水平最小、垂向围压为20MPa、16MPa、10MPa,设定初始地层孔隙压力为12MPa;根据工程经验设定安全系数n1、n2分别为1.25、1.35,根据式(3)、式(7)计算坍塌压力许用值及破裂压力许用值分别为10.9MPa、14.8MPa;通过资料获取盐水及超临界CO2重度分别为9.8KN·m-3、4.48KN·m-3;根据步骤2中的(2)-(7)开展最大化CO2封存效率场地试验,设定盐水抽取速度为0.5kg/s,根据式(10)及式(11)计算CO2注入速度为1.13kg/s,然后根据式(13)计算盐水回注速度为0.52kg/s;根据步骤2中的(3)、(7)开展常规CO2封存场地试验;根据步骤3中的(1)和(2)分别获取最大化CO2封存效率场地试验下盐水预抽量(等效体积)/>并分别与盐水预抽量理论值V预抽-理论、盐水回注量理论值V回注-理论对比,如图5所示;根据步骤4中的(1)和(2)分别获取最大化CO2封存效率场地试验、常规CO2封存场地试验下CO2注入量V1、V2、底部盐水层2最大***变形u1、u2;根据步骤4中的(3)对比两种条件下CO2注入量V1、V2、底部盐水层2最大***变形u1、u2,如图6所示。
由图5可知,本发明所述最大化CO2封存效率的盐水预抽量和回注量理论值与场地试验的实际抽取量、实际回注量相差仅为3.1%、4.3%,表明一种最大化CO2封存效率的盐水预抽量和回注量计算方法可以满足现场工程需要。
由图6可知,本发明所述最大化CO2封存效率场地试验比常规CO2封存场地试验具有封存能力优势,CO2注入量提高25%;最大化CO2封存效率场地试验比常规CO2封存场地试验具有封存安全优势,底部盐水层2最大***变形降低20%。
Claims (5)
1.一种最大化CO2封存效率的盐水预抽量和回注量计算方法,其特征在于,具体计算步骤包括:
(1)计算底部盐水层盐水预抽量V预抽
盐水预抽量指在CO2注入前的盐水抽取量;CO2注入后,盐水抽取仍会持续一段时间,但该段时间内的盐水抽取量不计入盐水预抽量内;
式中,Pw c、为底部盐水层坍塌压力临界值和许用值,MPa;σH、σh分别为最大、最小水平主应力,由施加的水平围压获得,MPa;Pp为盐水层初始地层孔隙压力,由所施加地层孔隙压力获得,MPa;/>为内摩擦角,通过多组三轴压缩实验的莫尔包络线获得;c为内聚力,通过多组三轴压缩实验的莫尔包络线获得,MPa;α为比奥系数,一般为0.85;n1为安全系数,无量纲,结合工程经验设定;
②确定由盐水预抽引起的底部盐水层孔隙压力下降ΔPp预抽,以下降后的孔隙压力值Pp-ΔPp预抽不低于坍塌压力许用值[Pw c]为临界条件计算盐水预抽量V预抽;
式中,ΔPp预抽为盐水预抽引起的底部盐水层孔隙压力变化,MPa;H预抽为盐水预抽引起的水位高度变化,m;γ盐水为盐水重度,KN·m-3;V预抽、V渗分别为盐水预抽量、实际条件下盐水预抽引起的地下流体从无穷远处向预抽井渗流的补给量,m3;S为底部盐水层水平面积,m2;φ下为底部盐水层孔隙度,无量纲;
(2)计算顶部盐水层盐水回注量V回注
式中,ΔPp回注为盐水回注引起的顶部盐水层孔隙压力变化,MPa;H回注为盐水回注引起的水位高度变化,m;V回注、V′渗分别为盐水回注量、实际条件下盐水回注引起的地下流体从回注井向无穷远处渗流的溢出量,m3;S′为顶部盐水层水平面积,m2;φ上为顶部盐水层孔隙度,无量纲;
(3)底部盐水层CO2注入与盐水抽取同步发生时,建立CO2注入速度v注与盐水抽取速度v抽取的关系:
式中,ΔPp注为CO2注入引起底部盐水层孔隙压力变化MPa;v注为CO2注入速度,m3·s-1;v抽取为盐水抽取速度,m3·s-1;v″渗-v渗为实际条件下底部盐水层地下流体从无穷远处向预抽井及注入井渗流的补给速度,m3·s-1;
(4)底部盐水层CO2注入与顶部盐水层盐水回注同步发生时,建立CO2注入速度v注与盐水回注速度v回注的关系:
式中,v回注为盐水回注速度,m3·s-1;v注为CO2注入速度,m3·s-1。
2.一种提高CO2封存效率的预抽-回注方法,其特征在于,具体包括以下步骤:
首先对底部盐水层进行盐水预抽;当盐水预抽量达到权利要求1中计算得到的盐水预抽量V预抽时,开始向底部盐水层注入CO2,同时由权利要求1中所建立的CO2注入速度与盐水抽取速度关系控制CO2注入速度;当CO2检测器检测到CO2时,停止盐水抽取;当底部盐水层孔隙压力恢复到初始孔隙压力时,开始向顶部盐水层回注盐水,同时由权利要求书1中所建立的盐水回注速度与CO2注入速度关系控制盐水回注速度;当盐水回注量达到权利要求书1中计算得到的盐水回注量V回注时,停止盐水回注;当底部盐水层孔隙压力达到破裂压力许用值时,停止CO2注入。
3.一种最大化CO2封存效率的场地试验验证方法,验证权利要求1所述的最大化CO2封存效率的盐水预抽量和回注量计算方法,其特征在于,具体操作步骤包括:
步骤1:布设场地试验
布设CO2封存的场地试验装置;
步骤2:开展场地试验
(1)施加压力
通过液压支架及配重块给盐水层施加三向围压,通过液体增压器及承压钢管给盐水层施加初始地层孔隙压力;
(2)预抽盐水
在CO2注入前,通过预抽井从底部盐水层中抽取盐水并暂存于储存箱;在CO2注入后,当盐水预抽井井底处CO2检测器检测到CO2时,停止盐水抽取;
(3)注入CO2
当流体压力传感器监测到的底部盐水层孔隙压力达到坍塌压力许用值时,开始通过注入井向底部盐水层中注入CO2;当流体压力传感器监测到的底部盐水层孔隙压力达到破裂压力许用值时,停止CO2注入;
(4)回注盐水
当流体压力传感器监测到的底部盐水层孔隙压力值恢复到初始孔隙压力时,开始通过回注井将从底部盐水层中抽取的盐水回注到顶部盐水层中;当流体压力传感器监测到的顶部盐水层孔隙压力达到顶部盐水层破裂压力许用值时,停止盐水回注;
(5)结束试验
停止CO2注入,试验结束;记录CO2注入量、底部盐水层***变形值、底部盐水层孔隙压力监测值达到坍塌压力许用值时的盐水实际抽取量、顶部盐水层孔隙压力监测值达到破裂压力许用值时的盐水实际回注量;
步骤3:验证计算方法
(1)盐水预抽量验证
获取底部盐水层孔隙压力监测值达到坍塌压力许用值时的盐水实际抽取量,以盐水实际抽取量的等效体积与盐水预抽量理论值对比,若误差在5%以内,表明盐水预抽量计算方法可以满足现场工程;所述盐水预抽量理论值是指权利要求1的步骤(1)所述的底部盐水层盐水预抽量V预抽;
(2)盐水回注量验证
获取顶部盐水层孔隙压力监测值达到破裂压力许用值时的盐水实际回注量,以盐水实际回注量的等效体积与盐水回注量理论值对比,若误差在5%以内,表明盐水回注量计算方法可以满足现场工程;所述盐水回注量理论值是指权利要求1的步骤(2)所述的顶部盐水层盐水回注量V回注;
步骤4:对比封存效率
(1)通过步骤2中的(1)、(2)、(3)和(4)开展最大化CO2封存效率场地试验,通过步骤2中的(5)获取CO2注入量和底部盐水层最大***变形;
(2)通过步骤2中的(1)和(3)开展常规CO2封存场地试验,通过步骤2中的(5)获取CO2注入量和底部盐水层最大***变形;
(3)对比上述(1)和(2)中的CO2注入量,确认最大化CO2封存效率场地试验具有封存能力优势;对比(1)和(2)中的底部盐水层最大***变形,确认最大化CO2封存效率场地试验具有封存安全优势。
4.根据权利要求3所述的一种最大化CO2封存效率的场地试验验证方法,其特征在于,所述CO2封存的场地试验装置,包括基坑、底部盐水层、盖层、顶部盐水层、位移传感器、流体压力传感器、防渗箱体、“回”字型基坑、液压支架、配重块、液体增压器、承压钢管、预抽井、储存箱、CO2检测器、定压泵、定压水箱、超临界CO2气罐、注入井、回注井。
5.根据权利要求4所述的一种最大化CO2封存效率的场地试验验证方法,其特征在于,底部盐水层、盖层、顶部盐水层构成盐水层,盐水层外部包裹有防渗箱体,盐水层和防渗箱体均布置在基坑内,防渗箱体与基坑壁面间用土回填密实;在底部盐水层上表面布置位移传感器,监测盐水层变形;在基坑外四周开挖“回”字型基坑,在“回”字型基坑内放置液压支架,施加水平围压;在盐水层顶部布置配重块,施加垂向围压;液体增压器与盐水层内的承压钢管连接,负责将增压后的盐水注入盐水层,实现初始地层孔隙压力施加;注入井一端连接超临界CO2气罐,另一端从盐水层顶部中心处垂直深入底部盐水层内,实现向底部盐水层注入CO2;预抽井设置在盐水层内,用于从底部盐水层中抽取盐水并暂存于储存箱;CO2检测器布置在预抽井井底处,用于判断CO2是否达到预抽井井底;回注井设置在顶部盐水层,负责将储存箱中存储的盐水回注到顶部盐水层;定压泵、定压水箱通过承压钢管与盐水层连通,由定压泵维持定压水箱内压力恒定且等于初始地层孔隙压力,模拟预抽盐水、CO2注入、回注盐水引起的实际条件下地下流体渗流过程,即无限大地层内孔隙压力的自调节功能;在预抽井、注入井及回注井的井底位置处布置流体压力传感器,监测盐水层地层孔隙压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210561742.6A CN115032337B (zh) | 2022-05-23 | 2022-05-23 | 一种盐水预抽量和回注量的计算方法及场地试验验证方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210561742.6A CN115032337B (zh) | 2022-05-23 | 2022-05-23 | 一种盐水预抽量和回注量的计算方法及场地试验验证方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115032337A CN115032337A (zh) | 2022-09-09 |
CN115032337B true CN115032337B (zh) | 2023-05-23 |
Family
ID=83121108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210561742.6A Active CN115032337B (zh) | 2022-05-23 | 2022-05-23 | 一种盐水预抽量和回注量的计算方法及场地试验验证方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115032337B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002340881A (ja) * | 2001-05-21 | 2002-11-27 | Tokai Univ | 大気ー海洋間における二酸化炭素交換の動的変動の計算方法 |
AU2010223803A1 (en) * | 2009-03-11 | 2011-08-25 | Roman Bilak | Process for sequestration of fluids in geological formations |
CA2751186C (en) * | 2010-12-23 | 2013-06-18 | Paxton Corporation | Zero emission steam generation process |
AU2012272496A1 (en) * | 2011-06-24 | 2014-01-16 | Bilak, Roman | Sequestration of greenhouse gasses by generating an unstable gas/saline front within a formation |
US11353621B2 (en) * | 2019-03-04 | 2022-06-07 | King Fahd University Of Petroleum And Minerals | Method and alarming system for CO2 sequestration |
CN111285374A (zh) * | 2020-03-24 | 2020-06-16 | 兰州理工大学 | 一种二氧化碳地层封存方法及*** |
CN111289698A (zh) * | 2020-03-30 | 2020-06-16 | 中国华能集团清洁能源技术研究院有限公司 | 一种实验模拟co2地质封存可行性的*** |
CN113266333B (zh) * | 2021-06-29 | 2023-04-18 | 西北大学 | 一种通过挤入饱和co2盐水改善油砂储层渗透率的方法 |
-
2022
- 2022-05-23 CN CN202210561742.6A patent/CN115032337B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115032337A (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103510501B (zh) | 一种防止回灌井中堵塞涌水的方法 | |
Zhang et al. | Model tests on tilt adjustment techniques for a mooring dolphin platform with three suction caisson foundations in clay | |
CN109024710A (zh) | 一种深基坑工程承压水管涌的封堵方法 | |
CN102011591A (zh) | 深井工作面注浆工艺 | |
CN106499384A (zh) | 煤层气定向井注入/压降试井测试装置及其方法 | |
CN107178077A (zh) | 一种岩溶地区溶洞上覆土层注浆加固结构及其设计方法 | |
CN115032337B (zh) | 一种盐水预抽量和回注量的计算方法及场地试验验证方法 | |
CN110714742B (zh) | 一种提高底水凝析气藏采收率的方法 | |
CN109931038B (zh) | 一种缝洞型油藏注氮气设计方法 | |
CN104213556A (zh) | 一种端承载抗浮锚杆组件及抗浮施工方法 | |
CN218861536U (zh) | 锚定筒锚固结构及其施工所用的装置 | |
CN104141294B (zh) | 基坑支护止水体系缺陷引发土体渗透破坏的判断方法 | |
CN114482021B (zh) | 预设反向摩阻力联合桩底土固结的方法及形成的灌注桩 | |
CN217681857U (zh) | 一种用于煤层顶板岩层变形监测的压水装置 | |
CN114293549A (zh) | 一种灰岩地区溶洞抗浮锚杆的施工方法 | |
Zhu et al. | Field test on the mechanism of composite bucket foundation penetrating sandy silt overlying clay | |
WO2016085627A1 (en) | Storage systems for storing and extracting energy | |
CN103526780B (zh) | 一种回灌井中堵水防渗方法 | |
Busari et al. | A finite element model for optimum design of plain concrete pressure tunnels under high internal pressure | |
CN210322666U (zh) | 一种注浆浆液析水实验装置 | |
CN112049105A (zh) | 富水砂层矿山法隧道地连墙封闭止水的施工方法 | |
CN113898379B (zh) | 一种实测注浆固结地层压力及参数优化的装置及方法 | |
CN114856507A (zh) | 盐穴储气库的注采气装置和注采气方法 | |
CN209779615U (zh) | 树根桩主动循踪分层注浆*** | |
Shen | Numerical Simulation Analysis of Confined Water |
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 |