CN102508294A - 一种利用时移地震勘探资料进行差异avo分析的方法 - Google Patents
一种利用时移地震勘探资料进行差异avo分析的方法 Download PDFInfo
- Publication number
- CN102508294A CN102508294A CN2011103209208A CN201110320920A CN102508294A CN 102508294 A CN102508294 A CN 102508294A CN 2011103209208 A CN2011103209208 A CN 2011103209208A CN 201110320920 A CN201110320920 A CN 201110320920A CN 102508294 A CN102508294 A CN 102508294A
- Authority
- CN
- China
- Prior art keywords
- alpha
- prime
- difference
- rho
- avo
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种利用时移地震勘探资料进行差异AVO分析的方法,从Zoeppritz方程的近似弹性模量简化公式出发(Bortfeld,1961),推导出了四维差异AVO的表达式。进而解释清楚了两次时移地震AVO差异的物理意义,解决了时移地震中差异AVO分析的关键问题。利用这个表达式,可以直接对两次地震资料的振幅差异进行分析和进行储层预测。这种差异AVO技术可以用于时移地震(四维地震)资料的差异振幅分析。快速识别差异振幅所反映的流体变化。进一步可以预测流体的饱和度。
Description
技术领域:
本发明专利涉及地震分析领域,用于地震监测,具体涉及一种利用时移地震勘探资料进行差异AVO分析的方法。
背景技术:
在时移(四维)地震监测中,AVO(Amplitude versus Offset)分析技术具有识别孔隙流体变化的潜力,并且可以用于估计孔隙流体的饱和度和压力(Tura and Lumley,1999;Landro,2001;Ma and Morozov,2010)。与传统的基于各向同性介质理论的二维和三维AVO分析类似,不同的对于平面波Zoeppritz方程的近似公式(Shuey,1986;Goodway,et al,1997)可用于四维AVO分析。然而,无论对于Shuey的近似或其他近似公式,四维地震的差异AVO并没有一个简洁的、物理意义明确的表达式。就是说时移地震振幅差异究竟意味着什么呢?如果知道这个差异的表达式含义,就可以直接利用这个表达式来进行四维地震差异AVO的解释。
常规二维、三维地震的AVO技术,具有明显的表达式和物理意义,在储层预测中起了重要的作用。而四维地震资料的差异AVO,究竟具有什么样的物理含义?如何利用差异AVO进行正、反演研究是目前四维地震解释面临的问题。
参考文献
Bortfield,R.,1961,Approximations to the reflection and transmission coefficients of plane wave longitudinal and transverse waves:Geophysical Prospecting,9,485-502.
Brekhovskikh,L.M.,1960,Waves in Layered Media,Academic Press,London.
Goodway,W.,Chen,T.,and Downton,J.,1997,Improved AVO fluid detection and lithology discrimination using Lam′epetrophysical parameters;“λρ”,“μρ”,&“λ/μ fluid stack”,from P and S inversions:67th Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,183-186.
Landro,M.,2001,Discrimination between pressure and fluid saturation changes from time lapse seismic data:Geophysics,66,836-844.
Ma,J.,and Morozov,I.B.,2004,Ray-path elastic impedance:2004 CSEG National Convention.
Ma,J.,and Morozov,I.B.,2006,Property of Zoeppntz elastic impedance:2006 CSPG-CSEG-CWLS National Convention.
Ma,J.,and Morozov,I.B.,2010,AVO modeling in of Pressure-Saturation effects in Weyburn CO2 Sequestration:The Leading Edge,29,178-183
Shuey,R.T.,1985,Asimplification of the Zoeppritz equations:Geophysics,50,609-614
Tura,A.,and Lumley,D.E.,1999,Estimating pressure and saturation changes from time-lapse AVO data:69thAnn.Internat.Mtg.,Soc.Expl.,Geophys.,Expanded Abstracts,1655-1658.
Wang,Y.,1999,Approximations to the Zoeppritz equations and their use in AVO analysis:Geophysics,64,1920-1927.
Wang,Z.,Cates,M.E.&Langan,R.T.,1998,Seismic monitoring of CO2 flood in a carbonate reservoir:Arock physics study:Geophysics,63,1604-1617.
发明内容:
本发明从Zoeppritz方程的近似弹性模量简化公式出发(Bortfeld,1961),推导出了四维差异AVO的表达式。进而解释清楚了两次时移地震AVO差异的物理意义,解决了时移地震中差异AVO分析的关键问题。利用这个表达式,可以直接对两次地震资料的振幅差异进行分析和进行储层预测。
具体如下:
一种利用时移地震勘探资料进行差异AVO分析的方法,包括如下步骤:
根据Zoeppritz方程的近似弹性模量简化公式推导两次地震的反射系数差异;Zoeppritz方程的近似弹性模量简化公式为:
由此,计算两次地震的反射系数差异
其中,R’2(p)和R1(p)分别为四维地震监测后、前的反射系数;R′f2(p)和Rf1(p))分别为四维地震监测后、前的流体反射系数;α1,β1,ρ1为盖层的纵波速度、横波速度和密度;α2,β2,ρ2为储层注入流体前或者开发前的纵波速度、横波速度和密度;α’2,β’2,ρ’2为储层注入流体后或者开发后的纵波速度、横波速度和密度;
这个公式也可以变换为:
差异AVO的正演计算:
结合流体替换模型,计算流体替换或者注入前后弹性参数的变化,再根据上述公式计算差异AVO的变化。
该方法可以直接用于差异AVO的反演计算,利用公式(3)将差异AVO或者反射系数转换为Brekhovskikh定义的声波阻抗ρα/cosθ,将声波阻抗转化为流体饱和度;或者直接建立流体饱和度变化与公式(2)的关系,从差异反射系数中直接反演流体饱和度
本发明的有益效果是:
这种差异AVO技术可以用于时移地震(四维地震)资料的差异振幅分析。快速识别差异振幅所反映的流体变化。进一步可以预测流体的饱和度。
附图说明:
图1是两次地震的反射系数差异图。
图2是第一类含气砂岩的例子反射系数比较图。
图3是三类含气砂岩模型中反射系数比较图。
图4是刚性反射系数的误差比例比较图。
图5是三类含气砂岩的差异反射振幅比较图。
具体实施方式:
下面结合附图对本发明做详细描述。
Zoeppritz方程的近似弹性模量简化公式(Bortfeld,1961;Wang,1999;Maand Morozov,2004)),可以表示为流体的反射系数和刚性部分的反射系数两部分。即
若四维地震监测后、前的反射系数分别表示为R’2(p)和R1(p),那么两次地震的反射系数差异为:
其中,R’2(p)和R1(p)分别为四维地震监测后、前的反射系数;R′f2(p)和Rf1(p))分别为四维地震监测后、前的流体反射系数;α1,β1,ρ1为盖层的纵波速度、横波速度和密度;α2,β2,ρ2为储层注入流体前或者开发前的纵波速度、横波速度和密度;α’2,β’2,ρ’2为储层注入流体后或者开发后的纵波速度、横波速度和密度。
很明显,上式是一个递推公式。就是说两次地震的反射系数差异可以表示为一个递推公式。这个递推公式的形式与传统垂直入射的反射系数公式、弹性阻抗(EI)的反射系数表达式类似。那么,基于垂直入射的正、反演理论和方法,就适合于该公式。这个公式也可以写为:
其中ρα/cosθ为Brekhovskikh(1960)定义的声阻抗。ρα/cosθ和ΔRf’(p)都可以用来分析时移地震观测中流体的变化。θ2为四维地震监测前地震波以θ1入射到储层时的透射角。
这个公式也可以表示为角度的函数。这个公式的图1示说明如下,即两次时移地震的反射系数差异,相当于波在两次流体变化的分解面上的反射系数。
这里通过三类含流体砂岩模型,验证差异AVO公式的精度,如表1所示:
表1.流体替换模型(Ma and Morozov,2006).
图2.第一类含气砂岩的例子。Zoeppritz(黑色)代表从Zoeppritz方程计算的反射系数。Bortfeld(红色)代表从Bortfeld公式(1)计算的反射系数.Rf(p)(绿色)和刚性部分(蓝色)代表Bortfeld公式(1)中的流体与刚性项。
图3.三类含气砂岩模型中,Bortfeld公式(1)中流体反射系数Rf(p)(绿色)和刚性反射系数(蓝色)的比较。
图4.刚性反射系数的误差比例。其中,Rigid_Error(Oil)为(Rigid(water)-Rigid(oil))/ΔRf’(p),而ΔRf’(p)是Oil/Water边界的反射系数。Error(Gas)为(Rigid(Water)-Rigid(Oil))/ΔRf’(p),而ΔRf’(p)为Gas/Water边界反射系数。
图5.三类含气砂岩的差异反射振幅比较。Zoeppritz(Water-Oil)代表从Zoeppritz方程计算的含水砂岩反射系数减去含油砂岩反射系数。Rf(Water)-Rf(Oil)表示从流体反射系数Rf(p)计算的含水砂岩反射系数减去含油砂岩反射系数。ΔR’f(Oil/Water)表示从公式(2)计算的含水砂岩反射系数减去含油砂岩反射系数。含水与含气砂岩反射系数的差异与含水与含油砂岩反射系数差异类似,含水与含水砂岩反射系数差异为0。
公式(2)-(3)可以用于计算四维地震监测中差异AVO的振幅。这些公式既可以用于预测差异AVO的振幅,也可以用于解释差异AVO的振幅。其中,计算差异AVO的正演问题时,需要结合流体替换模型,计算流体替换或者注入前后弹性参数的变化,再计算差异AVO的变化。这可以用于预测注入流体或者油藏不同开发阶段的四维地震反射系数差异。
公式(3)可以直接用于差异AVO的反演计算。可以利用公式(3)将差异AVO(或者反射系数),转换为Brekhovskikh定义的声波阻抗(ρα/cosθ),进一步可以将声波阻抗转化为流体饱和度。这里我们基于GASSMANN方程的流体替换理论,直接建立流体饱和度变化与公式(2)的关系,从差异反射系数中直接反演流体饱和度。表1为我们计算的流体替换模型。这里我们省略了从混合流体饱和度(油、气、水、CO2)与弹性参数之间关系的计算方法。
这里公式(2)、(3)应用的前提是储层的压力处于自然状态,即流体自然替换。按照美国能源部(DOE Annual Report,2004)报告,自然状态产出的油占整个石油生产的10%。注入水、热蒸汽或二氧化碳驱油过程中,储层的压力将有所改变。弹性模量将随压力的变化,因此公式(1)中第二项会有变化。但是注入流体一段时间后,储层的压力将恢复到储层原来的压力。因此,本公式在注入流体前、后储层压力大致相同的条件下,仍然实用。
Claims (2)
1.一种利用时移地震勘探资料进行差异AVO分析的方法,其特征在于:包括如下步骤:
1)根据Zoeppritz方程的近似弹性模量简化公式推导两次地震的反射系数差异;Zoeppritz方程的近似弹性模量简化公式为:
2)由此,计算两次地震的反射系数差异
其中,R’2(p)和R1(p)分别为四维地震监测后、前的反射系数;R′f2(p)和Rf1(p))分别为四维地震监测后、前的流体反射系数;α1,β1,ρ1为盖层的纵波速度、横波速度和密度;α2,β2,ρ2为储层注入流体前或者开发前的纵波速度、横波速度和密度;α’2,β’2,ρ’2为储层注入流体后或者开发后的纵波速度、横波速度和密度;
这个公式也变换为:
3)差异AVO的正演计算:
结合流体替换模型,计算流体替换或者注入前后弹性参数的变化,再根据上述公式计算差异AVO的变化。
2.根据权利要求1所述的方法,其特征在于,该方法可以直接用于差异AVO的反演计算,利用公式(3)将差异AVO或者反射系数转换为Brekhovskikh定义的声波阻抗ρα/cosθ,将声波阻抗转化为流体饱和度;或者直接建立流体饱和度变化与公式(2)的关系,从差异反射系数中直接反演流体饱和度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103209208A CN102508294B (zh) | 2011-10-20 | 2011-10-20 | 一种利用时移地震勘探资料进行差异avo分析的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103209208A CN102508294B (zh) | 2011-10-20 | 2011-10-20 | 一种利用时移地震勘探资料进行差异avo分析的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508294A true CN102508294A (zh) | 2012-06-20 |
CN102508294B CN102508294B (zh) | 2013-11-27 |
Family
ID=46220399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103209208A Active CN102508294B (zh) | 2011-10-20 | 2011-10-20 | 一种利用时移地震勘探资料进行差异avo分析的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508294B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104656133A (zh) * | 2015-03-10 | 2015-05-27 | 安徽赛斯米克能源技术有限公司 | 一种油藏模型限定下的四维地震反演解释方法 |
CN104678438A (zh) * | 2015-03-27 | 2015-06-03 | 西北大学 | 一种co2地质封存中四维地震资料co2分布预测的方法 |
CN106383362A (zh) * | 2016-08-19 | 2017-02-08 | 中国海洋石油总公司 | 一种提高薄储层时移地震差异识别能力的方法 |
CN106646603A (zh) * | 2017-01-04 | 2017-05-10 | 中海石油(中国)有限公司 | 一种实际时移地震资料处理差异的可靠性判断方法 |
CN106842326A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107918213A (zh) * | 2016-10-07 | 2018-04-17 | 科视数字系统美国股份有限公司 | 一种用于组合光束的装置 |
CN110471106A (zh) * | 2019-09-20 | 2019-11-19 | 西南石油大学 | 一种基于滤波器设计的时移地震反演方法 |
US10534100B2 (en) | 2017-10-26 | 2020-01-14 | Chevron U.S.A. Inc. | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data |
CN111239805A (zh) * | 2020-02-13 | 2020-06-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
US10718876B2 (en) | 2017-10-26 | 2020-07-21 | Chevron U.S.A. Inc. | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions |
CN112578457A (zh) * | 2020-11-24 | 2021-03-30 | 中油奥博(成都)科技有限公司 | 一种二氧化碳驱动煤层气的光纤井中地震监测方法 |
CN112649844A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5661697A (en) * | 1995-12-18 | 1997-08-26 | Atlantic Richfield Company | Method and apparatus for detection of sand formations in amplitude-versus-offset seismic surveys |
US6058074A (en) * | 1998-07-31 | 2000-05-02 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude-versus-offset analysis with improved measurement of background statistics |
US6223126B1 (en) * | 1999-10-20 | 2001-04-24 | Phillips Petroleum Company | Multi-attribute seismic waveform classification |
US6757216B1 (en) * | 2003-05-15 | 2004-06-29 | Exxonmobil Upstream Research Company | Method for post processing compensation of amplitude for misaligned and misstacked offset seismic data |
US20040246817A1 (en) * | 2003-06-04 | 2004-12-09 | Subhashis Mallick | Method and apparatus of performing amplitude variation with offset analysis |
CN101281253A (zh) * | 2007-04-06 | 2008-10-08 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
-
2011
- 2011-10-20 CN CN2011103209208A patent/CN102508294B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5661697A (en) * | 1995-12-18 | 1997-08-26 | Atlantic Richfield Company | Method and apparatus for detection of sand formations in amplitude-versus-offset seismic surveys |
US6058074A (en) * | 1998-07-31 | 2000-05-02 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude-versus-offset analysis with improved measurement of background statistics |
US6223126B1 (en) * | 1999-10-20 | 2001-04-24 | Phillips Petroleum Company | Multi-attribute seismic waveform classification |
US6757216B1 (en) * | 2003-05-15 | 2004-06-29 | Exxonmobil Upstream Research Company | Method for post processing compensation of amplitude for misaligned and misstacked offset seismic data |
US20040246817A1 (en) * | 2003-06-04 | 2004-12-09 | Subhashis Mallick | Method and apparatus of performing amplitude variation with offset analysis |
CN101281253A (zh) * | 2007-04-06 | 2008-10-08 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104656133A (zh) * | 2015-03-10 | 2015-05-27 | 安徽赛斯米克能源技术有限公司 | 一种油藏模型限定下的四维地震反演解释方法 |
CN104678438B (zh) * | 2015-03-27 | 2017-04-05 | 西北大学 | 一种co2地质封存中四维地震资料co2分布预测的方法 |
CN104678438A (zh) * | 2015-03-27 | 2015-06-03 | 西北大学 | 一种co2地质封存中四维地震资料co2分布预测的方法 |
CN106842326B (zh) * | 2015-12-04 | 2020-10-13 | 中国石油化工股份有限公司 | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 |
CN106842326A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107479090B (zh) * | 2016-06-08 | 2019-11-12 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN106383362B (zh) * | 2016-08-19 | 2018-06-19 | 中国海洋石油集团有限公司 | 一种提高薄储层时移地震差异识别能力的方法 |
CN106383362A (zh) * | 2016-08-19 | 2017-02-08 | 中国海洋石油总公司 | 一种提高薄储层时移地震差异识别能力的方法 |
CN107918213A (zh) * | 2016-10-07 | 2018-04-17 | 科视数字系统美国股份有限公司 | 一种用于组合光束的装置 |
CN106646603A (zh) * | 2017-01-04 | 2017-05-10 | 中海石油(中国)有限公司 | 一种实际时移地震资料处理差异的可靠性判断方法 |
CN106646603B (zh) * | 2017-01-04 | 2018-11-27 | 中海石油(中国)有限公司 | 一种实际时移地震资料处理差异的可靠性判断方法 |
US10534100B2 (en) | 2017-10-26 | 2020-01-14 | Chevron U.S.A. Inc. | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data |
US10718876B2 (en) | 2017-10-26 | 2020-07-21 | Chevron U.S.A. Inc. | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions |
CN110471106A (zh) * | 2019-09-20 | 2019-11-19 | 西南石油大学 | 一种基于滤波器设计的时移地震反演方法 |
CN110471106B (zh) * | 2019-09-20 | 2022-03-11 | 西南石油大学 | 一种基于滤波器设计的时移地震反演方法 |
CN112649844A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
CN112649844B (zh) * | 2019-10-10 | 2023-08-15 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
CN111239805A (zh) * | 2020-02-13 | 2020-06-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
CN111239805B (zh) * | 2020-02-13 | 2021-02-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
CN112578457A (zh) * | 2020-11-24 | 2021-03-30 | 中油奥博(成都)科技有限公司 | 一种二氧化碳驱动煤层气的光纤井中地震监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102508294B (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508294B (zh) | 一种利用时移地震勘探资料进行差异avo分析的方法 | |
Lee | A simple method of predicting S-wave velocity | |
Johnston | Methods and applications in reservoir geophysics | |
Zhang et al. | Direct estimation of the fluid properties and brittleness via elastic impedance inversion for predicting sweet spots and the fracturing area in the unconventional reservoir | |
CN103163553B (zh) | 基于多重孔隙介质模型的地震烃类检测方法及装置 | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
CN104252007A (zh) | 一种相容性岩石物理建模方法 | |
CN106019375B (zh) | 一种页岩气地层层理地球物理评价方法 | |
CN106842326B (zh) | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 | |
Ma et al. | AVO modeling of pressure-saturation effects in Weyburn CO 2 sequestration | |
CN102156297A (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
Katterbauer et al. | Synergizing crosswell seismic and electromagnetic techniques for enhancing reservoir characterization | |
CN104570065A (zh) | 一种利用地震波阻抗定量反演孔隙度的方法 | |
Lim et al. | Inference of geomechanical properties of shales from AVO inversion based on the Zoeppritz equations | |
Lan et al. | An improved seismic fluid identification method incorporating squirt flow and frequency-dependent fluid-solid inversion | |
Li et al. | Reservoir multiparameter prediction method based on deep learning for CO2 geologic storage | |
Allawi et al. | Developing A Semi-Analytical Model for Estimating Mechanical Properties of Sandstone Reservoirs: Enhancing Applications in hydrocarbon production and Underground Gas Storage | |
White et al. | Stochastic fluid modulus inversion | |
Khalid et al. | Effect of kerogen and TOC on seismic characterization of lower cretaceous shale gas plays in lower Indus Basin, Pakistan | |
Ravazzoli et al. | AVA seismic reflectivity analysis in carbon dioxide accumulations: Sensitivity to CO2 phase and saturation | |
Wang et al. | Velocity prediction and secondary-pores quantitative inversion for complex carbonate reservoir | |
Phan et al. | Porosity estimation from seismic data at Dickman Field, Kansas for carbon sequestration | |
Nair et al. | Seismic inversion and its applications in reservoir characterization | |
Zong et al. | Robust AVO inversion for elastic modulus and its application in fluid factor calculation | |
Khatiwada et al. | A numerical sensitivity analysis to monitor CO2 sequestration in layered basalt with coda waves |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |