CN106199695B - 基于空变目标函数的叠前三参数反演方法 - Google Patents
基于空变目标函数的叠前三参数反演方法 Download PDFInfo
- Publication number
- CN106199695B CN106199695B CN201610489375.8A CN201610489375A CN106199695B CN 106199695 B CN106199695 B CN 106199695B CN 201610489375 A CN201610489375 A CN 201610489375A CN 106199695 B CN106199695 B CN 106199695B
- Authority
- CN
- China
- Prior art keywords
- mrow
- space
- object function
- mfrac
- prestack
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 5
- 230000007812 deficiency Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于空变目标函数的叠前三参数反演方法,该基于空变目标函数的叠前三参数反演方法包括:步骤1,基于贝叶斯理论,构建贝叶斯目标函数;步骤2,采用用于叠前三参数反演的鲁棒初始模型进行三参数空间交叉建模;步骤3,求取空变因子,修正观测地震数据不足以修正模型的问题;步骤4,采用三参数空间交叉约束建模得到的初始模型约束下的叠前三参数反演,而目标函数则是加入空变因子的空变目标函数。该基于空变目标函数的叠前三参数反演方法,在反演应用过程中,通过三参数模型空间交叉约束的方法提高目标函数的收敛性,适用于井网分布不均匀或目的层横向快速变化地区。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种基于空变目标函数的叠前三参数反演方法。
背景技术
叠前地震反演技术因其对叠前地震信息的挖掘和丰富的弹性结果而受到重视,在复杂隐蔽性油气勘探和开发过程中发挥了重要作用。但是地震资料反映的是地下反射界面的信息,对地下剩余油气的分布和丰度的预测存在多解性;测井资料能够更精确地反映井点处地下岩层信息,但空间外推性较差。地震资料和测井资料的二者结合虽然能够提高反演技术的空间预测能力,但是现有方法没有针对性地对两种资料的空间匹配问题开展研究。
由于地震资料精度有限,对于反映尺度较小或变化较快的地质体时多解性较大,也会导致反演结果不理想。尽管随着测井资料等附加信息的引入在一定程度上提高了反演精度,但是实际应用过程中,尤其是在开展叠前地震反演时同样存在问题。一方面,工区内井位分布不均匀,初始模型的精度在空间上差异较大;另一方面,受到横波测井资料匮乏的影响,基于估算横波构建的初始模型容易引起反演过程中的收敛问题。为此我们发明了一种新的基于空变目标函数的叠前三参数反演方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种基于空变目标函数的叠前三参数反演方法,针对地震资料和测井资料的空间匹配问题提出一种新型的地震反演空变目标函数,利用空变权值约束残差项来改善反演结果的局部质量,同时为反演过程提供了修正低频模型的方法。
本发明的目的可通过如下技术措施来实现:基于空变目标函数的叠前三参数反演方法,该基于空变目标函数的叠前三参数反演方法包括:步骤1,基于贝叶斯理论,构建贝叶斯目标函数;步骤2,采用用于叠前三参数反演的鲁棒初始模型进行三参数空间交叉建模;步骤3,求取空变因子,修正观测地震数据不足以修正模型的问题;步骤4,采用三参数空间交叉约束建模得到的初始模型约束下的叠前三参数反演,而目标函数则是加入空变因子的空变目标函数。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,基于贝叶斯理论,将待反演地层参数的后验概率分解为其先验分布和实际地震资料与合成地震记录的似然函数。
在步骤1中,考虑地震褶积模型:
d=Gr+n
其中d=[d1,d2,...,dN]T是观测到的地震数据,r=[r1,r2,...,rM]T是反射系数序列,G是N×M维子波褶积矩阵,n=[n1,n2,...,nN]T表示观测噪声,由贝叶斯公式可得以下近似:
p(r|d)∝p(r)p(d|r)
其中p(r|d)表示反射系数的后验概率,p(r)表示反射系数先验信息,p(d|r)表示似然函数;
已知先验信息p(r)和似然函数,根据贝叶斯公式得到待反演参数的后验分布,将后验分布取对数,通过约束观测地震数据和模型数据的差值,使之达到最小,构建目标函数如下:
多项式中,第一项为主要优化项,第二项为反射系数系数约束项,第三项为模型数据约束项,其中G是子波褶积矩阵,r是反射系数序列,K是模型矩阵,ξ是扰动项,α、β和μ分别是三项的权值,在反演过程中人为给定。
在步骤3中,利用提取子波与空间交叉建模得到的三个模型进行褶积运算,获得合成地震记录,并进一步获得反演过程的初始残差,针对权值采取空变策略,通过引入空变因子C,修正观测地震数据不足以修正模型的问题,新构建的空变目标函数为:
其中,x和y为空间坐标,θ为入射角度。
在步骤2中,首先把工区内测井资料中三参数完备的井与不完备的井分为两组;然后再按照不完备组所缺少的测井资料类型再进行细分,包括缺密度资料组、缺横波速度资料组;最后,降低完备组井参与建模的资料数,尽可能保留附近不完备组所缺资料,同时去掉一种附近存在较多的一项资料类型。
本发明中的基于空变目标函数的叠前三参数反演方法,在反演应用过程中,通过三参数模型空间交叉约束的方法提高目标函数的收敛性,使得基于新目标函数的叠前反演方法适用于井网分布不均匀或目的层横向快速变化地区。本发明形成了以三参数空间交叉建模和贝叶斯空变目标函数为特色的新型叠前三参数地震反演技术方法,其主要优点是:①按照贝叶斯思想将先验模型信息单独表达出来,用于构建空变目标函数;②不减少工区内可用约束井位的数量,通过建立优势参数的空间分配策略避免了由于估算横波等步骤可能带来的线性相关问题;③通过模型残差构建空变因子来进一步提高新模型约束能力的空间平均化,同时避免由减少单个模型参与井位带来的反演精度损失,形成优势互补。因此,对于目标地层横向变化较快,且原始资料不足时开展这种基于空变目标函数的叠前三参数反演不失为一种更为稳定可靠的方法。
附图说明
图1为本发明的基于空变目标函数的叠前三参数反演方法的一具体实施例的流程图;
图2为本发明的一具体实施例中采用常规反演方法出现局部极值扩散的反演剖面的示意图;
图3为本发明的一具体实施例中采用新方法得到的正确反演剖面的示意图;
图4为本发明的一具体实施例中采用常规方法预测的目标井位A钻遇的储层平面图;
图5为本发明的一具体实施例中采用新方法预测的目标井位A钻遇的储层平面图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的基于空变目标函数的叠前三参数反演方法的流程图。
在步骤101,构建贝叶斯目标函数。
步骤101实现空变目标函数的构建。基于贝叶斯理论,将待反演地层参数的后验概率分解为其先验分布和实际地震资料与合成地震记录的似然函数。
考虑地震褶积模型:
d=Gr+n
其中d=[d1,d2,...,dN]T是观测到的地震数据,r=[r1,r2,...,rM]T是反射系数序列,G是N×M维子波褶积矩阵,n=[n1,n2,...,nN]T表示观测噪声。由贝叶斯公式可得以下近似:
p(r|d)∝p(r)p(d|r)
其中p(r|d)表示反射系数的后验概率,p(r)表示反射系数先验信息,p(d|r)表示似然函数。
已知先验信息p(r)和似然函数,根据贝叶斯公式就可以得到待反演参数的后验分布。将后验分布取对数。地震反演目标函数的基本构建思路是通过约束观测地震数据和模型数据的差值,使之达到最小。构建目标函数如下:
多项式中,第一项为主要优化项,第二项为反射系数系数约束项,第三项为模型数据约束项。其中G是子波褶积矩阵,r是反射系数序列,K是模型矩阵,ξ是扰动项,α、β和μ分别是三项的权值,在反演过程中人为给定。流程进入步骤102。
在步骤102,三参数空间交叉建模。
步骤102中采用了用于叠前三参数反演的鲁棒初始模型构建策略。对于测井资料不足的工区反演,在构建叠前三参数反演所需要的三个初始模型时可用的约束井位数量不一,通常是缺少横波测井井位。此时,三个模型的空间约束力不均性就会带来如局部极值扩散等问题。
为了充分保障横向内各个位置附近的井位可以参与初始反射系数计算,本方案中采用一种新的三参数建模策略。首先把工区内测井资料中三参数完备的井与不完备的井分为两组;然后再按照不完备组所缺少的测井资料类型再进行细分,如缺密度资料组、缺横波速度资料组等;最后,降低完备组井参与建模的资料数,具体做法如下:尽可能保留附近不完备组所缺资料,同时去掉一种附近存在较多的一项资料类型。流程进入到步骤103。
在步骤103,求取空变因子。
利用提取子波与空间交叉建模得到的三个模型进行褶积运算,获得合成地震记录,并进一步获得反演过程的初始残差。本发明主要针对权值采取空变策略,通过引入空变因子C,修正观测地震数据不足以修正模型的问题。新构建的空变目标函数为:
其中,x和y为空间坐标,θ为入射角度。流程进入步骤104。
在步骤104,叠前空变反演实现。
步骤104即采用三参数空间交叉约束建模得到的初始模型约束下的叠前三参数反演,而目标函数则是加入空变因子的空变目标函数。
图2为本发明的一具体实施例中反演实验剖面上出现异常值扩散的问题,其主要原因是空间内三参数约束不均匀导致反演迭代过程中局部极值不能在空间内得到合理的纠正,而影响了空间地层反演结果。利用本发明方法对反演目标函数进行改进,则有效地消除了局部极值问题,如图3所示。图4和图5则分别给出了改进前后两个反演结果在目标井位A附近的平面图。图4所示的常规反演方法在刻画A井钻遇的薄互地层时,边界不够清楚,且与周边多口开发井位钻遇情况不符;图5所示结果则清晰地揭示了薄互地层的边界,及与周边其他储层发育位置间的关系,符合实际钻井情况。
Claims (3)
1.基于空变目标函数的叠前三参数反演方法,其特征在于,包括以下步骤:
步骤1,基于贝叶斯理论,构建贝叶斯目标函数;
步骤2,采用三参数空间交叉建模,获得用于叠前三参数反演的初始模型;
步骤3,求取空变因子,实现迭代过程中模型的空间差异化修正;
步骤4,联合步骤1与步骤3获得空变目标函数,实现叠前三参数反演方法;
在步骤1中,考虑地震褶积模型:
d=Gr+n
其中d=[d1,d2,...,dN]T是观测到的地震数据,r=[r1,r2,...,rM]T是反射系数序列,G是N×M维子波褶积矩阵,n=[n1,n2,...,nN]T表示观测噪声;
已知先验信息p(r)和似然函数,根据贝叶斯公式得到待反演参数的后验分布,将后验分布取对数,通过约束观测地震数据和模型数据的差值,使之达到最小,构建目标函数如下:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&alpha;</mi>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mn>2</mn>
</mfrac>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>r</mi>
<mi>T</mi>
</msup>
<mi>r</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mi>&mu;</mi>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mi>r</mi>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mi>r</mi>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
</mrow>
多项式中,第一项为主要优化项,第二项为反射系数系数约束项,第三项为模型数据约束项,其中K是模型矩阵,ξ是扰动项,α、β和μ分别是三项的权值。
2.根据权利要求1所述的基于空变目标函数的叠前三参数反演方法,其特征在于,在步骤3中,利用提取子波与三参数空间交叉建模得到的三个模型进行褶积运算,获得合成地震记录,并进一步获得反演过程的初始残差,针对权值采取空变策略;通过引入空变因子,实现迭代过程中模型的空间差异化修正;新构建的空变目标函数为:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mn>2</mn>
</mfrac>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>r</mi>
<mi>T</mi>
</msup>
<mi>r</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mi>&mu;</mi>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mi>r</mi>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mi>r</mi>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,x和y为空间坐标,θ为入射角度。
3.根据权利要求1所述的基于空变目标函数的叠前三参数反演方法,其特征在于,在步骤2中,首先把工区内测井资料中三参数完备的井与不完备的井分为两组;然后再按照不完备组所缺少的测井资料类型再进行细分,包括缺密度资料组、缺横波速度资料组;最后,降低完备组井参与建模的资料数,尽可能保留附近不完备组所缺资料,同时去掉一种附近存在较多的一项资料类型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610489375.8A CN106199695B (zh) | 2016-06-29 | 2016-06-29 | 基于空变目标函数的叠前三参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610489375.8A CN106199695B (zh) | 2016-06-29 | 2016-06-29 | 基于空变目标函数的叠前三参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106199695A CN106199695A (zh) | 2016-12-07 |
CN106199695B true CN106199695B (zh) | 2018-03-09 |
Family
ID=57462158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610489375.8A Active CN106199695B (zh) | 2016-06-29 | 2016-06-29 | 基于空变目标函数的叠前三参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199695B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932819B (zh) * | 2017-02-23 | 2019-01-01 | 河海大学 | 基于各向异性马尔科夫随机域的叠前地震参数反演方法 |
CN110858005B (zh) * | 2018-08-24 | 2021-11-05 | 中国石油化工股份有限公司 | 一种基于基追踪横向多道约束的各向异性参数反演方法 |
CN113608263B (zh) * | 2021-09-02 | 2023-06-20 | 中海石油(中国)有限公司深圳分公司 | 一种叠前反演的可行性评价方法、装置、设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8422335B2 (en) * | 2009-06-17 | 2013-04-16 | Seoul National University R&Db Foundation | Apparatus and method for imaging subsurface structure of target area by using waveform inversion |
CN104007467A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 一种基于混合范数正则化的叠前三参数反演实现的储层与流体预测方法 |
CN103245970B (zh) * | 2012-02-08 | 2015-05-27 | 中国石油化工股份有限公司 | 叠前地震宽角度反演方法 |
CN104808243A (zh) * | 2015-05-08 | 2015-07-29 | 中国石油大学(华东) | 一种叠前地震贝叶斯反演方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EA017177B1 (ru) * | 2008-03-21 | 2012-10-30 | Эксонмобил Апстрим Рисерч Компани | Эффективный способ инверсии геофизических данных |
-
2016
- 2016-06-29 CN CN201610489375.8A patent/CN106199695B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8422335B2 (en) * | 2009-06-17 | 2013-04-16 | Seoul National University R&Db Foundation | Apparatus and method for imaging subsurface structure of target area by using waveform inversion |
CN103245970B (zh) * | 2012-02-08 | 2015-05-27 | 中国石油化工股份有限公司 | 叠前地震宽角度反演方法 |
CN104007467A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 一种基于混合范数正则化的叠前三参数反演实现的储层与流体预测方法 |
CN104808243A (zh) * | 2015-05-08 | 2015-07-29 | 中国石油大学(华东) | 一种叠前地震贝叶斯反演方法和装置 |
Non-Patent Citations (2)
Title |
---|
基于空变阻抗约束的贝叶斯稀疏脉冲反演;肖开攀 等;《CT理论与应用研究》;20151130;第24卷(第6期);第809-817页 * |
基于贝叶斯理论的薄层反演方法;黄捍东 等;《石油地球物理勘探》;20111231;第46卷(第6期);第919-924页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106199695A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105954804B (zh) | 页岩气储层脆性地震预测方法及装置 | |
US10795053B2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
CN104866653B (zh) | 一种获取地下三维密度结构的方法 | |
CN103105624B (zh) | 基于数据库技术的纵横波时差定位方法 | |
CN106597533A (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
CN107817535A (zh) | 薄夹层的确定方法和装置 | |
US10466388B2 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
US20060052938A1 (en) | Method for creating facies probability cubes based upon geologic interpretation | |
CN106199695B (zh) | 基于空变目标函数的叠前三参数反演方法 | |
CN106716181A (zh) | 使用多点统计模拟调节基于对象或事件的储层模型 | |
CN104297785A (zh) | 岩相约束储层物性参数反演方法及装置 | |
CN102222365A (zh) | 复杂空间的曲面重构方法 | |
CN104504754B (zh) | 一种油气储层多点统计建模的方法及装置 | |
CN103886216B (zh) | 一种基于地质矢量信息的多点地质统计方法 | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
US9411915B2 (en) | Method of selecting positions of wells to be drilled for petroleum reservoir development | |
CN103105622B (zh) | 基于数据库技术的同型波时差定位方法 | |
CN111815773A (zh) | 适用于机器学习算法的三维复杂地质模型标签制作方法 | |
Alhuthali et al. | Optimal rate control under geologic uncertainty | |
CN104297792B (zh) | 一种扇上叠置水道储层的相控反演方法 | |
CN113933899B (zh) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 | |
CN107939372A (zh) | 针对小断块油藏的最优井位部署方法和装置 | |
CN107884831A (zh) | 一种基于近似支持向量机的avo类型判别方法 | |
CN111768503B (zh) | 一种基于三维地质模型的海砂资源量估算方法 | |
CN106338764A (zh) | 生物启发计算地层圈闭油气藏超剥线识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |