CN107831540A - 储层物性参数直接提取新方法 - Google Patents
储层物性参数直接提取新方法 Download PDFInfo
- Publication number
- CN107831540A CN107831540A CN201710753641.8A CN201710753641A CN107831540A CN 107831540 A CN107831540 A CN 107831540A CN 201710753641 A CN201710753641 A CN 201710753641A CN 107831540 A CN107831540 A CN 107831540A
- Authority
- CN
- China
- Prior art keywords
- parameters
- reservoir
- physical property
- physical
- porosity
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000000704 physical effect Effects 0.000 claims abstract description 64
- 238000000605 extraction Methods 0.000 claims abstract description 19
- 238000012417 linear regression Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000012706 support-vector machine Methods 0.000 claims description 7
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 4
- 238000000611 regression analysis Methods 0.000 claims description 3
- 239000011435 rock Substances 0.000 abstract description 11
- 238000011160 research Methods 0.000 description 10
- 239000012530 fluid Substances 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6224—Density
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种储层物性参数直接提取新方法。
背景技术
储层物性反演一直是储层地球物理研究的难点和热点,目前常用的方法有两种,一是地震多属性物性反演(吴媚等,2008;吴秋波等,2011),其思路是建立物性和地震属性之间定量的关系,进而反演物性参数,二是地质统计学物性反演(孙思敏等,2007;刘百红等,2009),主要是通过地质统计学随机模拟的思路实现物性参数的反演。它们的不足之处在于都需要大量的测井信息,岩石物理的机理也不清楚,且只能实现单一物性参数的反演。
针对这些不足,国内外的学者们开展了基于岩石物理的储层物性反演研究。斯坦福大学的Mukerji等(2001)和Eidsvik等(2004)等人率先将岩石物理和叠前反演相结合,估计储层物性参数,Bachrach(2006)进一步地将岩石物理模型用于储层孔隙度和饱和度反演,Spikes等(2008)等人提出了用井资料和地震资料对储层物性进行估计的概率地震反演方法,Larsen等(2006)等人则提出了一种基于马尔科夫链的储层流体反演方法,随后Grana等(2010)等人发展了这些前人的研究成果,将岩石物理和贝叶斯反演相结合,实现了物性的反演,取得了较好的效果。国内在基于岩石物理的储层物性反演方面的研究较少,郑继新等(2009)首次将统计岩石物理用于含气储层饱和度与孔隙度联合反演,取得了一定的成果,但是其对岩石物理模型和实际应用的研究不够深入,属于一种前瞻性的研究,胡华锋等(2012)进一步发展了这种方法,提出了基于贝叶斯分类的储层物性参数联合反演方法,也取得了较好的应用效果,但是方法研究仍需要进一步深入。杨培杰等(2017)在前人研究的基础上,针对砂泥岩储层,基于岩石物理理论、贝叶斯理论、蒙特卡罗和遗传算法,提出了一种孔隙度和含水饱和度同步反演方法。
如何从地震数据中有效地提取孔隙度、泥质含量、含水饱和度等储层物性参数是储层地球物理研究的热点和难点,目前的物性提取方法主要是通过两步法实现的,即首先反演弹性参数,然后通过弹性参数再次计算物性参数,该方法的不足是效率低,并且准确度不高。为提高物性提取的效率和准确度,提出了一种从叠前数据直接提取物性参数的新方法。
发明内容
本发明的目的是提供一种可以同时反演储层孔隙度、泥质含量和密度等参数,提高储层物性估计结果的客观性和准确性的储层物性参数直接提取新方法。
本发明的目的可通过如下技术措施来实现:储层物性参数直接提取新方法,该储层物性参数直接提取新方法包括:步骤1,输入地震数据、地震子波、测井数据以及地质层位;步骤2,以统计岩石物理分析为基础,从测井数据出发,建立起储层物性参数和弹性参数之间的定量关系;步骤3,推导物性提取的近似公式,里面包含有孔隙度、泥质含量、密度三个参数;步骤4,以叠前分析为主要手段,结合物性提取的近似公式和点对点方程,从叠前数据中直接提取物性参数。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,输入的地震数据为叠前偏移道集或叠前角道集,测井数据包括纵波速度、横波速度、密度、孔隙度和泥质含量这几种曲线。
在步骤1中,通过测井二次解释获得孔隙度,泥质含量以及饱和度这些物性参数,测井二次解释是指在地区井资料丰富的情况下,利用各种信息,来建立更精确的孔隙度模型、渗透率模型和饱和度模型。
在步骤1中,在测井二次测井解释中,孔隙度曲线为通过怀利Wyllie、雷默Raymer这些估算方法,通过声波、密度、中子这些数据计算;泥质含量曲线为通过伽玛或自然电位曲线,通过克拉维尔Clavier、斯泰贝Steiber、拉里奥诺夫Larionov这些估算方法进行计算;饱和度是在原始状态下储积层中石油体积占有效孔隙体积的百分数,利用电阻率曲线,通过阿尔奇Archies、费尔和汉迈克ert&Hammack这些公式计算饱和度。
在步骤2中,利用测井数据,通过回归分析得到统计的岩石物理模型,首先通过线性回归初步确定回归的系数,然后通过支持向量机SVM的思路,对这些回归系数进行微调,以提高物性参数和弹性参数间的定量关系的准确度。
在步骤2中,物性参数和弹性参数间的定量关系的表达式为:
其中,M表示弹性参数,包括纵波速度、横波速度;N表示物性参数,包括孔隙度、泥质含量;C和D表示回归过程所得到的不同的系数。
在步骤3中,将待提取的物性参数和弹性参数相结合,从叠前理论出发,基于阿基方程,得到包含孔隙度、泥质含量、密度的近似方程,简记为PCD公式,该方程为物性参数直接反演的基础。
在步骤3中,所构建的物性参数近似方程为:
其中,Δ表示取数据的差值,Aφ表示与孔隙度反射有关的常数,φ表示孔隙度;AN表示与泥质含量反射有关的常数,N表示泥质含量;Aρ表示与密度反射有关的常数,ρ表示密度,Rpp表示地震数据反射系数,θ地震波入射角。
在步骤4中,将物性参数近似方程和叠前反演有效结合,开发了一种点对点的叠前反演策略,通过点对点方程实现约束。
在步骤4中,建立的点对点方程为:
其中,ξ表示待提取的物性参数,n是约束点的数量,是变化的,m≥n>1,m是待反演的地震道的总点数,r表示地层反射系数。
本发明中的储层物性参数直接提取新方法,主要用于石油地球物理勘探的储层地球物理领域,首先,以统计岩石物理分析为基础,从测井数据出发,建立起储层物性参数和弹性参数之间的定量关系,其次,推导物性提取的PCD近似公式,里面包含有孔隙度、泥质含量、密度三个参数,最后,以叠前分析为主要手段,结合PCD公式和点对点方程,即可从叠前数据中直接提取物性参数。整个提取过程简洁高效,在实际生产中具有很好的应用效果。该方法将岩石物理和地质统计分析相结合,结合PCD公式和点对点方程,应用叠前道集、地震子波、控制层位,可以同时反演储层孔隙度、泥质含量和密度等参数,提高储层物性估计结果的客观性和准确性,降低投资风险,为地质人员估算油气储量、确定井位提供更为可靠的依据。
附图说明
图1为本发明的储层物性参数直接提取新方法的一具体实施例的流程图;
图2为本发明的一具体实施例中点对点约束模式的示意图;
图3为本发明的一具体实施例中直接提取物性参数中的孔隙度的结果的示意图;
图4为本发明的一具体实施例中直接提取物性参数中的泥质含量的结果的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的储层物性参数直接提取新方法的流程图。
步骤101、数据输入
主要包括地震数据、地震子波、测井数据以及地质层位,地震数据需要是叠前偏移道集,也可以是叠前角道集,子波是在物性参数提取前就在软件中提取出来的,测井数据主要包括纵波速度、横波速度、密度、孔隙度和泥质含量等几种曲线,地质层位是由地质人员解释出来的。
常规的测井无法直接测量储层的孔隙度和泥质含量等物性参数,需要通过二次解释来获得。测井二次解释是指在地区井资料比较丰富的情况下进行的,充分利用各种信息,如:试油结果,来建立更精确的孔隙度模型、渗透率模型和饱和度模型,目的是要弄清地区的各项特征,精细描述。是理论与实际相结合的过程,是不断认识的过程,使解释更接近生产状态,为以后解释建立相关模型。
测井二次测井解释的关键是测井曲线标准化处理、“四性”关系分析和解释模版确定。并结合岩心分析资料、试油试采资料.研究出适合该区的测井解释模型和油水层识别标准.提高了解释精度。具体地说,这个工作其实就是进行第二次测井解释。通常,是对开发井进行的。油田打开发井时,会有一个随钻测井解释,投产后,会进行储量复算。复算时,测井人员会把全油田的所有开发井,系统考虑,重新解释处理,包括定量处理、定性结论的复检等。这个过程,也就是大家遇见最多的测井二次解释了,其实最常用的叫法是测井复算。另外,还有些零碎的,从定义上讲,应该也算二次解释。比如,打完探井,进行储量评价时,会需要对所有探井、评价井重新做统一处理,部分开发井随钻并投产后,证实原先解释成果有误,需要及时纠正等等。
在测井二次测井解释中,孔隙度曲线可通过Wyllie、Raymer等估算方法,通过声波、密度、中子等数据计算;泥质含量曲线可通过伽玛或自然电位曲线,通过Clavier、Steiber、Larionov、Bateman等估算方法进行计算;饱和度是在原始状态下储积层中石油体积占有效孔隙体积的百分数,它是测井解释中的一个重要的参数,利用电阻率曲线,通过Archies、fert&Hammack等公式计算饱和度。
步骤102、岩震分析
将这些解释出的测井曲线(孔隙度和泥质含量曲线)和地震弹性参数相结合,通过岩石物理理论建立物性参数和弹性参数间的定量关系。
表征岩石物理学特征的地震参数主要有岩石的弹性模量、纵波速度、横波速度、密度和衰减等等,它是识别岩性及油气的重要参数,也是联系储集层特征的参数,进行定量地震油藏描述的桥梁,例如,孔隙度φ、含水饱和度Sw、渗透率K和地层压力P等。岩石物理学用于烃类直接检测的主要问题是:当孔隙流体成分改变时,速度、密度和衰减是如何变化的。
为了实现这个目的,有必要估算不同饱和流体状态之间的岩石特性的差异,因而需要建立一些基本的岩石物理关系。这些关系包括岩石特性和弹性参数之间的经验公式、统计关系和理论关系。另一方面,了解地震波特性与岩石、流体性质的关系,有助于模拟地震波在复杂介质中的传播规律。这就要求对反映岩石物理学特征的地震参数与岩性和烃类的关系有深刻的理解,通过岩石物理正演,可以由储层岩石、流体参数得到储层弹性参数,反之,通过岩石物理反演,可以由储层弹性参数反演储层岩石、流体参数。
岩石物理就是建立物性参数和弹性参数间的定量关系。一般来说,可以通过理论、经验以及统计关系来描述弹性参数和物性参数之间的关系。然而,对于理论岩石物理模型,总是缺少应用的条件,对于经验公式,使用范围会受到限制。因此,在实际应用中,往往是利用测井数据,通过回归分析得到统计的岩石物理模型。
线性回归一般是通过最小二乘的方法建立不同参数间线性关系,线性关系虽然简单,但是鲁棒性很好,特别是在做地震储层预测的时候。但是在实际应用中,仅仅通过线性回归方法,准确度不够,因此将线性回归和SVM相结合,首先是通过线性回归初步确定回归的系数,然后通过SVM的思路,对这些回归系数进行微调,可进一步提高物性参数和弹性参数间的定量关系的准确度。物性参数和弹性参数间的定量关系的一般表达式为:
其中,M表示弹性参数,比如纵波速度、横波速度等;N表示物性参数,如孔隙度、泥质含量等;C和D表示回归过程所得到的不同的系数。
步骤103、建立近似方程
将待提取的物性参数和弹性参数相结合,从叠前理论出发,基于Aki方程,得到包含孔隙度、泥质含量、密度的近似方程,简记为PCD公式,该方程即为物性参数直接反演的基础。所构建的物性参数近似方程如公式(2)所示。
其中,Δ表示取数据的差值,Aφ表示与孔隙度反射有关的常数,φ表示孔隙度;AN表示与泥质含量反射有关的常数,N表示泥质含量;Aρ表示与密度反射有关的常数,ρ表示密度,Rpp表示地震数据反射系数,θ地震波入射角。
该公式是实现物性参数直接提取的关键所在,然后设定近似公式中的各项参数,其中的地震波入射角是个非常关键的问题,可用现有的软件进行处理。
步骤104、物性参数直接提取
储层物性参数直接提取的两个理论基础是物性参数近似方程和叠前反演的方法,对于叠前反演方法来说,目前已经比较成熟,不是本专利探究的重点,但是,如何将物性参数近似方程和叠前反演有效结合,是提高反演结果分辨率、特别是稳定性非常关键的问题,为了实现这个目的,开发了一种点对点的叠前反演策略,通过点对点方程实现约束,点对点方程如公式(3)所示。
其中,ξ表示待提取的物性参数,n是约束点的数量,是变化的,m≥n>1,m是待反演的地震道的总点数,r表示地层反射系数。通过点对点的叠前反演策略,可以大大提高物性提取过程的效率和提取结果的准确性。
图2为点对点约束模式的示意图,通过将待反演数据的头、尾进行点对点约束,从而将反演结果约束在真实值的附近,并且提高了横向物性取结果的稳定性。
以叠前分析为主要手段,结合PCD公式和点对点方程,应用叠前道集、地震子波、地质层位,即可实现孔隙度与泥质含量的直接提取,整个提取过程简介高效,具有很好的推广应用价值。对于直接提取的物性参数结果,可输入专业软件进行专业的绘制,然后进行数据的分析和解释,最终为预测有利储层,确定有利勘探目标提供数据的基础。
图3为直接提取的孔隙度示意图,图4为直接提取的泥质含量示意图。一般来说,对于砂泥岩储层,孔隙度高、泥质含量低的地层为好储层,有利于流体的运移或是储存,通过物性参数直接提取的方法,有效地提高了这两个参数提取结果的准确度和稳定性,同时,通过这两个物性参数的联合应用,可以有效地减小地震反演多解性,提高储层预测准确性。
本发明的储层物性参数直接提取新方法,通过岩震关系分析,从测井数据出发,采用线性回归和SVM相结合的方法,建立储层物性和弹性参数之间的定量关系,所建立起来的方程更加准确,更符合研究区实际情况。从叠前理论出发,推导得到包含孔隙度、泥质含量、密度等参数的近似方程,简记为PCD方程,该方程即为物性参数直接反演的基础。以叠前分析为主要手段,结合PCD公式和点对点方程,应用叠前道集、地震子波、控制层位,即可实现储层物性参数的直接提取,过程简单明了,可大大提高物性提取的效率和可靠性。
Claims (10)
1.储层物性参数直接提取新方法,其特征在于,该储层物性参数直接提取新方法包括:
步骤1,输入地震数据、地震子波、测井数据以及地质层位;
步骤2,以统计岩石物理分析为基础,从测井数据出发,建立起储层物性参数和弹性参数之间的定量关系;
步骤3,推导物性提取的近似公式,里面包含孔隙度、泥质含量、密度三个参数;
步骤4,以叠前分析为主要手段,结合物性提取的近似公式和点对点方程,从叠前数据中直接提取物性参数。
2.根据权利要求1所述的储层物性参数直接提取新方法,其特征在于,在步骤1中,输入的地震数据为叠前偏移道集或叠前角道集,测井数据包括纵波速度、横波速度、密度、孔隙度和泥质含量这几种曲线。
3.根据权利要求2所述的储层物性参数直接提取新方法,其特征在于,在步骤1中,通过测井二次解释获得孔隙度,泥质含量以及饱和度这些物性参数,测井二次解释是指在地区井资料丰富的情况下,利用各种信息,来建立更精确的孔隙度模型、渗透率模型和饱和度模型。
4.根据权利要求3所述的储层物性参数直接提取新方法,其特征在于,在步骤1中,在测井二次测井解释中,孔隙度曲线为通过怀利Wyllie、雷默Raymer这些估算方法,通过声波、密度、中子这些数据计算;泥质含量曲线为通过伽玛或自然电位曲线,通过克拉维尔Clavier、斯泰贝Steiber、拉里奥诺夫Larionov这些估算方法进行计算;饱和度是在原始状态下储积层中石油体积占有效孔隙体积的百分数,利用电阻率曲线,通过阿尔奇Archies、费尔和汉迈克ert&Hammack这些公式计算饱和度。
5.根据权利要求1所述的储层物性参数直接提取新方法,其特征在于,在步骤2中,利用测井数据,通过回归分析得到统计的岩石物理模型,首先通过线性回归初步确定回归的系数,然后通过支持向量机SVM的思路,对这些回归系数进行微调,以提高物性参数和弹性参数间的定量关系的准确度。
6.根据权利要求5所述的储层物性参数直接提取新方法,其特征在于,在步骤2中,物性参数和弹性参数间的定量关系的表达式为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,M表示弹性参数,包括纵波速度、横波速度;N表示物性参数,包括孔隙度、泥质含量;C和D表示回归过程所得到的不同的系数。
7.根据权利要求1所述的储层物性参数直接提取新方法,其特征在于,在步骤3中,将待提取的物性参数和弹性参数相结合,从叠前理论出发,基于阿基方程,得到包含孔隙度、泥质含量、密度的近似方程,简记为PCD公式,该方程为物性参数直接反演的基础。
8.根据权利要求7所述的储层物性参数直接提取新方法,其特征在于,在步骤3中,所构建的物性参数近似方程为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<msub>
<mi>A</mi>
<mi>&phi;</mi>
</msub>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&phi;</mi>
</mrow>
<mi>&phi;</mi>
</mfrac>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>N</mi>
</msub>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>N</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>&rho;</mi>
</msub>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>&rho;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Δ表示取数据的差值,Aφ表示与孔隙度反射有关的常数,φ表示孔隙度;AN表示与泥质含量反射有关的常数,N表示泥质含量;Aρ表示与密度反射有关的常数,ρ表示密度,Rpp表示地震数据反射系数,θ地震波入射角。
9.根据权利要求1所述的储层物性参数直接提取新方法,其特征在于,在步骤4中,将物性参数近似方程和叠前反演有效结合,开发了一种点对点的叠前反演策略,通过点对点方程实现约束。
10.根据权利要求9所述的储层物性参数直接提取新方法,其特征在于,在步骤4中,建立的点对点方程为:
其中,ξ表示待提取的物性参数,n是约束点的数量,是变化的,m≥n>1,m是待反演的地震道的总点数,r表示地层反射系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710753641.8A CN107831540A (zh) | 2017-08-28 | 2017-08-28 | 储层物性参数直接提取新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710753641.8A CN107831540A (zh) | 2017-08-28 | 2017-08-28 | 储层物性参数直接提取新方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107831540A true CN107831540A (zh) | 2018-03-23 |
Family
ID=61643219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710753641.8A Pending CN107831540A (zh) | 2017-08-28 | 2017-08-28 | 储层物性参数直接提取新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831540A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109993219A (zh) * | 2019-03-21 | 2019-07-09 | 长江大学 | 基于支持向量机的辫状河致密砂岩储渗单元划分方法 |
CN112415616A (zh) * | 2019-08-22 | 2021-02-26 | 中国石油天然气集团有限公司 | 深埋储层孔隙度反演方法及装置 |
CN112965103A (zh) * | 2021-02-09 | 2021-06-15 | 中国石油大学(华东) | 一种多重孔隙储层叠前地震概率化多道反演方法 |
CN113933910A (zh) * | 2020-07-13 | 2022-01-14 | 中国石油化工股份有限公司 | 物性参数同时计算的方法、系统、存储介质与电子设备 |
CN114428355A (zh) * | 2020-09-23 | 2022-05-03 | 中国石油化工股份有限公司 | 储层物性参数直接反演方法、装置、电子设备及介质 |
CN114428294A (zh) * | 2020-09-23 | 2022-05-03 | 中国石油化工股份有限公司 | 储层参数统计反演方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104950331A (zh) * | 2015-06-30 | 2015-09-30 | 中国石油天然气集团公司 | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 |
CN106556867A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油天然气股份有限公司 | 基于贝叶斯分类的相控孔隙度反演方法 |
-
2017
- 2017-08-28 CN CN201710753641.8A patent/CN107831540A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104950331A (zh) * | 2015-06-30 | 2015-09-30 | 中国石油天然气集团公司 | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 |
CN106556867A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油天然气股份有限公司 | 基于贝叶斯分类的相控孔隙度反演方法 |
Non-Patent Citations (5)
Title |
---|
DONG LI ET AL.: "Direct estimation of petrophysical properties based on AVO inversion", 《2015 SEG NEW ORLEANS ANNUAL MEETING》 * |
刘宝和: "《中国石油勘探开发百科全书 勘探卷》", 30 November 2008, 北京:石油工业出版社 * |
宋延杰等: "《地球物理测井》", 31 May 2011, 北京:石油工业出版社 * |
崔维等: "基于叠前地震道集实现储层物性参数直接反演", 《中国地球科学联合学术年会2014》 * |
杨培杰等: "敏感流体因子定量分析与直接提取", 《石油地球物理勘探》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109993219A (zh) * | 2019-03-21 | 2019-07-09 | 长江大学 | 基于支持向量机的辫状河致密砂岩储渗单元划分方法 |
CN112415616A (zh) * | 2019-08-22 | 2021-02-26 | 中国石油天然气集团有限公司 | 深埋储层孔隙度反演方法及装置 |
CN113933910A (zh) * | 2020-07-13 | 2022-01-14 | 中国石油化工股份有限公司 | 物性参数同时计算的方法、系统、存储介质与电子设备 |
CN113933910B (zh) * | 2020-07-13 | 2024-04-09 | 中国石油化工股份有限公司 | 物性参数同时计算的方法、系统、存储介质与电子设备 |
CN114428355A (zh) * | 2020-09-23 | 2022-05-03 | 中国石油化工股份有限公司 | 储层物性参数直接反演方法、装置、电子设备及介质 |
CN114428294A (zh) * | 2020-09-23 | 2022-05-03 | 中国石油化工股份有限公司 | 储层参数统计反演方法、装置、电子设备及存储介质 |
CN114428294B (zh) * | 2020-09-23 | 2024-08-06 | 中国石油化工股份有限公司 | 储层参数统计反演方法、装置、电子设备及存储介质 |
CN112965103A (zh) * | 2021-02-09 | 2021-06-15 | 中国石油大学(华东) | 一种多重孔隙储层叠前地震概率化多道反演方法 |
CN112965103B (zh) * | 2021-02-09 | 2022-08-23 | 中国石油大学(华东) | 一种多重孔隙储层叠前地震概率化多道反演方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106368691B (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
AU2010332157B2 (en) | Workflow for petrophysical and geophysical formation evaluation of wireline and LWD log data | |
EP2616979B1 (en) | Production estimation in subterranean formations | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
RU2661489C1 (ru) | Способ комплексирования исходных данных для уточнения фильтрационного строения неоднородных карбонатных коллекторов | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN106842301B (zh) | 一种凝灰质砂岩有利储层的定量识别与预测方法 | |
OA12151A (en) | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data. | |
CN114746774A (zh) | 预测油气和运移路径的综合地质力学模型 | |
CN110646850B (zh) | 隔夹层地震预测方法及装置 | |
CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
CN106842326B (zh) | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 | |
CN112034521B (zh) | 一种欠压实与生烃混合成因地层超压的预测方法 | |
CN105572747A (zh) | 一种识别高盐富钾区注水开发油藏碎屑岩岩性的方法 | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
CN116047602B (zh) | 基于生烃数值模拟的ii型水合物饱和度预测方法 | |
CN115586572B (zh) | 一种孔隙参数与储层参数的地震岩石物理解析反演方法 | |
CN113534263B (zh) | 一种不依赖测井资料的含油气饱和度预测方法 | |
Davies | Permeability Modelling of a Sandstone Reservoir in Parts of the Niger Delta | |
CN113376692B (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
Worthington et al. | Optimizing the value of reservoir simulation through quality-assured initialization | |
CN108875109B (zh) | 一种异常地层压力的预测方法及系统 | |
CN111594156A (zh) | 一种天然气水合物饱和度计算方法及系统 | |
Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180323 |
|
WD01 | Invention patent application deemed withdrawn after publication |