CN105443122B - 一种测井解释模型的处理方法及装置 - Google Patents
一种测井解释模型的处理方法及装置 Download PDFInfo
- Publication number
- CN105443122B CN105443122B CN201511000241.7A CN201511000241A CN105443122B CN 105443122 B CN105443122 B CN 105443122B CN 201511000241 A CN201511000241 A CN 201511000241A CN 105443122 B CN105443122 B CN 105443122B
- Authority
- CN
- China
- Prior art keywords
- parameter
- expression
- model
- log
- stratum
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000015572 biosynthetic process Effects 0.000 claims description 71
- 239000011435 rock Substances 0.000 claims description 14
- 230000005251 gamma ray Effects 0.000 claims description 8
- 238000005481 NMR spectroscopy Methods 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 3
- 239000011707 mineral Substances 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 230000035699 permeability Effects 0.000 claims description 2
- 239000000470 constituent Substances 0.000 claims 1
- 238000001514 detection method Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 7
- 238000011156 evaluation Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 238000010195 expression analysis Methods 0.000 abstract 1
- 238000005755 formation reaction Methods 0.000 description 59
- 235000019738 Limestone Nutrition 0.000 description 10
- 239000006028 limestone Substances 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 239000010459 dolomite Substances 0.000 description 6
- 229910000514 dolomite Inorganic materials 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 235000010755 mineral Nutrition 0.000 description 2
- 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 1
- 101100454393 Homo sapiens LCOR gene Proteins 0.000 description 1
- 102100038260 Ligand-dependent corepressor Human genes 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 235000002639 sodium chloride Nutrition 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000004457 water analysis Methods 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种测井解释模型的处理方法及装置，属于油气勘探技术领域。所述方法包括：根据预定地层分析资料确定待测地层的地层参数解释模型；根据所述地层参数解释模型配置参数表达式；将参数表达式进行后缀表达式解析；根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线。本发明通过确定的地层参数模型配置参数表达式，并将该参数表达式进行解析及求解，从而根据实际地层评价的需要，灵活构建出适应的流程框架；在处理过程中，也可根据储层类型，自由定义和添加更具有针对性的地层参数模型，储层适用范围更广，并有效提高了处理效率和计算精度。
Description
技术领域
本发明涉及一种测井解释模型的处理方法及装置，属于油气勘探技术领域。
背景技术
在现有的测井资料处理解释领域中，处理方法与程序是整个工作流程的核心。通常情况下，解释人员根据地层地质情况，选取合适的处理方法，通过输入相应的测井曲线信息计算获得地层岩性、孔隙度、饱和度等参数，以满足油气勘探开发的需要。
长期以来，测井资料处理发展形成了各种岩性识别、油气层判别以及储层物性、饱和度评价等解释模型方法，并形成了针对不同储层类型评价的处理程序，如单孔隙分析程序POR、粘土分析处理程序CLASS及复杂岩性解释程序CRA等。
从方法原理来看，这些程序在整体流程上基本一致，都是选定几条相关的测井曲线，分步计算地层矿物成分、孔隙度、饱和度等参数，不同的只是其中采用的解释模型。但在实际应用中，由于地层岩性、物性等因素影响，对不同目标储层处理时往往需要使用不同的处理程序，影响了处理解释效率。同时，在上述处理程序中，通常只包含了几种固定类型的解释模型，由于不同油田区块地质情况不一样，这些处理程序中提供的解释模型往往不通用，需要重新建立各地区针对性的解释模型和处理程序，较为费时费力。此外，在测井解释过程中，还包括了许多根据岩石物理分析资料建立的经验公式，现有处理程序欠缺对这些公式的处理能力，因而限制了新解释模型的验证和推广应用。
发明内容
本发明为解决现有的测井解释模型的处理方法存在的针对不同目标地层处理时已有的解释模型及处理程序无法通用的问题，进而提出了一种测井解释模型的处理方法及装置，具体包括如下的技术方案：
一种测井解释模型的处理方法，包括：
根据预定地层分析资料确定待测地层的地层参数解释模型；
根据所述地层参数解释模型配置参数表达式；
将所述参数表达式进行后缀表达式解析；
根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线。
一种测井解释模型的处理装置，包括：
模型确定单元，用于根据预定地层分析资料确定待测地层的地层参数解释模型；
表达式配置单元，用于根据所述地层参数解释模型配置参数表达式；
解析单元，用于将所述参数表达式进行后缀表达式解析；
处理单元，用于根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线。
本发明的有益效果是：通过确定的地层参数解释模型配置参数表达式，并将该参数表达式进行解析及求解，从而根据实际储层评价的需要，灵活构建出适应的流程框架；在处理过程中，也可根据储层类型，自由定义和添加更具有针对性的地层参数解释模型，储层适用范围更广，并有效提高了处理效率和计算精度。
附图说明
图1是以示例的方式示出了测井解释模型的处理方法的流程图。
图2是实施例一提出的测井解释模型的处理方法的流程图。
图3是实施例一提出的A井含水饱和度-电阻增大率实验测量结果示意图，其中横坐标表示含水饱和度，单位为V/V，纵坐标表示电阻增大率。
图4是实施例一提出的A井地层参数计算结果示意图。
图5是实施例二提出的测井解释模型的处理装置的结构图。
具体实施方式
结合图1所示，本实施例提出的测井解释模型的处理方法包括：
步骤11，根据预定地层分析资料确定待测地层的地层参数解释模型。
根据待测地层评价的需要，该预定地层分析资料可以包括测井资料和地层岩心实验分析数据。其中的测井资料可以包括测井仪器采集到的自然伽马测井曲线、电阻率测井曲线、三孔隙度测井曲线、元素俘获能谱测井曲线或核磁测井曲线中的至少一种；其中的实验分析数据可以包括岩心岩电分析数据、核磁分析数据、声波分析数据、孔渗分析数据或全岩矿物分析数据中的至少一种。
其中，地层参数可以包括用于储层评价所需的地层矿物组分、地层孔隙度、饱和度、渗透率等参数中的至少一种，地层参数解释模型可以通过岩心实验数据拟合或已有的测井解释方法确定。
步骤12，根据地层参数解释模型配置参数表达式。
由于步骤11中建立的地层参数解释模型通常包含的变量类型较多，例如可能既包含参与计算的测井曲线，也可能包含有各类模型参数，还可能包含有一些中间变量。另外，变量名称也包括较多种，例如同一变量名称既可能对应于模型参数，也可能为采集到的测井曲线，会对后续的参数表达式解析和求解的准确性造成较大的影响。
因此，本步骤针对地层参数解释模型的特点，将组成地层参数解释模型的元素划分为运算操作符、普通数字、模型参数符号、中间变量和测井曲线等五类，并将模型参数符号通过第一关键字表示，将测井曲线通过第二关键字表示；另外，地层参数解释模型中包含的运算操作符可为预先定义的操作符类型；中间变量则可采用与第一关键字或第二关键字类似的方法定义，通常可将不同地层参数解释模型中公用的一些的变量定义为中间参数，这样既能简化模型公式的编辑，还能提高后续计算效率。
步骤13，将参数表达式进行后缀表达式解析。
对参数表达式的解析可采用计算机领域的公式解析方法，即将参数表达式从中缀形式转换为后缀形式，并分别以堆栈的方式记录存储。
在转换的过程中，首先对参数表达式中包含的中间变量进行逐次替换，直至参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线；然后通过对参数表达式中的各组成元素依次扫描，并根据关键字类型进行匹配，识别出参数表达式中各元素类型；最后将参数表达式由中缀形式转换为后缀形式并以堆栈的方式存储。
可选的，后缀表达式堆栈中的每个元素都是一个数据结构体，该数据结构体可以包括四个字段，第一个字段为字符串型(string)，用于存放表达式元素形式；第二个字段为双精度浮点型(double)，用于存放该元素的值，数字元素的值即为其本身，模型参数和测井曲线值可通过后续计算中给定；第三个字段为整型(int)，用于标识元素类型，普通数字为0，预定义的运算操作符为1，模型参数符号为2、测井曲线为3；第四个字段为整型(int)，当元素类型为运算操作符时，该数值为预定义运算操作符类型的编号，否则，记为-1。该数据结构可表示为：
步骤14，根据经过解析的后缀表达式确定待测地层的地层参数结果曲线。
对待测地层的每个处理深度点，依次从以堆栈形式存储的后缀表达式中取出元素，当取出的元素类型为模型参数符号或测井曲线时，获取对应处理深度点的模型参数值和测井曲线值，并代入该后缀表达式计算，以确定预定深度点的处理结果。直到获得待测地层的每个处理深度点的处理结果，即可获得连续的全井段的地层参数结果曲线。根据该地层参数结果曲线即可实现对地层岩性、储层发育情况及流体性质的判别，进而指导储层综合评价。
采用本实施例提出的技术方案，通过确定的地层参数解释模型配置参数表达式，并将该参数表达式进行解析及求解，从而根据实际地层评价的需要，灵活构建出适应的流程框架；在处理过程中，也可根据储层类型，自由定义和添加更具有针对性的地层参数解释模型，储层适用范围更广，并有效提高了处理效率和计算精度。
下面通过具体的实施例对本发明提出的技术方案进行详细说明：
实施例一
结合图2所示，本实施例提出的测井解释模型的处理方法包括：
步骤21，根据预定地层分析资料确定待测地层的地层参数解释模型。
本实施例以某油田A井资料处理为例，该井4090.0-4120.0m井段位于奥陶系马家沟组，该区岩性为台地相云岩、云灰岩，夹泥岩薄层，孔隙类型以晶间孔为主。针对该区碳酸盐岩储层评价，所需基本地层参数包括泥质含量、灰岩、白云岩含量、地层孔隙度、饱和度。根据该地区经验认识，常规自然伽马曲线能够较好的反映地层泥质含量变化；光电吸收界面指数曲线对灰岩和白云岩成分反映明显；密度、中子曲线反映了地层孔隙大小；深、浅侧向电阻率曲线体现了原状地层和冲洗带部分孔隙中流体成分的变化。因此，可利用上述测井曲线建立解释模型以获取相应的地层参数。
综合岩心实验分析结果，本实施例确定的各地层参数解释模型如下：
泥质含量采用自然伽马曲线GR计算，并可通过以下公式计算获得：
其中，SH表示自然伽马相对值，无量纲；GR表示自然伽马曲线值，单位为GAPI(以标准单位API计量的伽马射线值)；GRmax、GRmin表示模型参数，分别为纯泥岩和纯岩性地层自然伽马测井值，单位为GAPI；GCUR表示与地层年代有关的经验系数，新地层例如取3.7，老地层例如取2；Vsh表示计算获得的泥质含量，单位为V/V(体积比)。
地层灰岩和白云岩含量采用光电吸收界面指数曲线PE计算，并可通过以下公式计算获得：
Vdolo＝1-Vsh-Vlime (4)
其中，Vsh表示根据式(2)计算获得的泥质含量，单位为V/V；PE表示光电吸收界面指数曲线，单位为B/E；PEmax、PEmin表示模型参数，分别为纯灰岩和纯云岩地层光电吸收界面指数响应值，例如分别取5.1B/E、3.14B/E；Vlime和Vdolo分别表示计算获得的灰岩含量和白云岩含量，单位为V/V。
地层孔隙度采用密度-中子交会方法计算，并可通过以下公式计算获得：
其中，DENma、DENf、LCOR表示模型参数，分别为岩石骨架密度值、流体骨架密度值和岩石骨架中子值，一般分别取值2.83G/CM3、1.0G/CM3和0V/V；DEN、CNL分别表示密度和中子测井值；表示密度孔隙度，单位为V/V；表示中子孔隙度，单位为V/V；表示计算获得的地层总孔隙度，单位为V/V。
含水饱和度和含油气饱和度可采用电阻率曲线进行计算。为了准确确定该井段地层岩石电性特征，本实施例对4101.3m岩心进行了岩电分析实验，图3所示的是含水饱和度(Sw)-电阻增大率(I)实验测量结果，其中的实线是采用Archie公式获得的拟合结果。结合图3所示，该井段岩石电阻增大率与含水饱和度变化较好的符合Archie关系，其计算公式为：
So＝1-Sw (11)
其中，I表示电阻增大率，无量纲；a、b、m、n、Rw分别表示模型系数，a、b、m、n根据岩电实验的拟合结果，取值分别为1、1.0052、2、1.701，无量纲；Rw表示地层水电阻率，可根据水分析资料确定，单位为OHMM；Rt表示电阻率测井值；Sw表示含水饱和度，单位为V/V；So表示含油气饱和度，单位为V/V。
步骤22，根据地层参数解释模型确定输出地层参数类型和输出曲线名称。
根据步骤21中确定的地层参数解释模型以及储层评价需求，可配置测井解释模型的实际输出结果为泥质含量、灰岩含量、白云岩含量、孔隙度、含水饱和度、含油气饱和度，输出名称分别为VSH、VLIME、VDOLO、POR、SW、SO。而相应计算式(1)、(5)、(6)等都为参数计算中的中间变量，并且由式(7)计算的地层总孔隙度同时也是式(10)计算含水饱和度Sw的中间变量。
步骤23，根据地层参数解释模型配置参数表达式。
根据步骤21中确定的地层参数解释模型，对各地层参数解释模型进行配置。在地层参数解释模型中，组成模型的元素可分为五类，即运算操作符、普通数字、模型参数符号、中间变量以及测井曲线，中间变量可由运算操作符、普通数字、模型参数符号、其它中间变量以及测井曲线组成。在进行模型公式编辑时，模型参数符号可使用第一关键字“PARA<>”包围起来；而测井曲线可使用第二关键字“CURVE<>”包围起来；运算操作符则可为预先定义的操作符类型；对于中间变量，通常可将解释模型中公用的一些变量定义为中间变量，这样既能简化模型公式的编辑，还能提高后续计算效率。
例如，利用电阻率曲线计算含水饱和度的Archie公式可记为：
SW＝(PARA<A>*PARA<B>*PARA<RW>/(CURVE<RT>*(FAI∧
PARA<M>)))∧(1/PARA<N>)
其中，A、B、M、N、RW分别为模型参数；RT为测井曲线；FAI为中间变量，通过类似方法定义，还可以在其它模型公式中重复使用。
本实施例以泥质含量、灰岩含量及白云岩含量的计算模型为例进行说明，其中，泥质含量计算模型可进行如下的配置：
SH＝(CURVE<GR>-PARA<GRMIN>)/(PARA<GRMAX>-PARA<GRMIN>) (12)
VSH＝(2∧(PARA<GCUR>)*SH-1)/(2∧PARA<GCUR>-1) (13)
其中，SH表示中间变量，GR表示自然伽马曲线值，GRMIN和GRMAX分别表示模型参数，VSH表示泥质含量，GCUR表示与地层年代有关的经验系数。
通过式(12)定义了中间变量SH，并在式(13)中调用。
另外，灰岩和白云岩的含量计算模型可分别配置为：
VLIME＝(1-VSH)*(CURVE<PE>-PARA<PEMIN>)/(PARA<PEMAX>-PARA<PEMIN>)(14)
VDOLO＝1-VSH-VLIME (15)
其中，VLIME表示灰岩含量，PE表示光电吸收界面指数曲线，PEMAX和PEMIN分别表示模型参数，VDOLO表示白云岩含量。
式(13)配置的泥质含量计算模型既可作为输出结果，同时也可作为其它参数计算模型所需的中间变量。
同理，可对式(5)-(13)进行配置，配置方法和过程同上，本实施例不再赘述。
步骤24，检查参数表达式的合法性。
本步骤是对步骤23中配置完成的各参数表达式的合法性进行检查。通过对参数表达式中的各元素进行逐次匹配，识别出参数表达式中包含的操作符、模型参数符号、中间变量、测井曲线等，同时对参数表达式中的左右括号的完整性进行检验。
以步骤23中配置的式(12)为例，若中间变量SH定义为：
SH＝(GR-PARA<GRMIN>)/(PARA<GRMAX>-PARA<GRMIN>)(16)
则在合法性检查过程中，通过关键字匹配，字符GRMIN与GRMAX会被自动判别为模型参数，而字符GR会被识别为中间变量，但同时对已定义的中间变量进行查找，未发现该中间变量名，则字符GR会被判定为非法，提示用户检查和/或修正。
步骤25，将参数表达式进行后缀表达式解析。
根据计算机领域的公式解析方法，将步骤24验证合法的参数表达式从中缀形式转换为后缀形式，并分别以堆栈的形式记录存储。在转换过程中，可首先对参数表达式中包含的中间变量进行逐次替换，直至该参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线。然后，通过对参数表达式的各组成元素依次扫描，并根据关键字类型进行匹配，识别出公式中各元素类型。最后将该参数表达式由中缀形式转换为后缀形式并放入堆栈存储。
以泥质含量计算模型的式(13)为例，首先对式(13)中包含的中间变量SH进行替换，替换后表达式为：
VSH＝(2∧(PARA<GCUR>*(CURVE<GR>-PARA<GRMIN>)/(PARA<GRMAX>-PARA<GRMIN>))-1)/2∧PARA<GCUR>-1) (17)
替换后，式(17)中不再包含有中间变量，则可对该参数表达式进行元素类型匹配并转换为后缀表达式。对该参数表达式从左至右依次扫描，其中的元素“(”、“2”、“^”、“(”分别为括号、数字、操作符、括号，直接依照中缀表达式转后缀表达式规则入栈，而元素PARA<GCUR>通过关键字匹配为模型参数，入栈后其名称为GCUR，元素类型为2(参数符号)。同理，对上述参数表达式扫描并转换完成后，后缀表达式堆栈结构可表示为：
步骤26，根据经过解析的后缀表达式确定待测地层的地层参数结果曲线。
在实际处理过程中，可对每个处理深度点获取当前深度各参数表达式中包含的测井曲线值和模型参数值。对每个参数表达式，依次从步骤25存储的堆栈中取出元素，当元素类型为模型参数符号或测井曲线时，则将具体值代入参数表达式计算，依照后缀表达式值计算方法获得该参数表达式在当前深度点的结果并输出。对全井段处理后，即可得到连续的地层参数结果曲线。
根据步骤22的配置，结合图4所示，本实施例提出的测井解释模型的处理方法的输出结果包括：左起第一部分表示岩性曲线，其中的曲线a表示GR，曲线b表示PE，曲线c表示CAL)；第二部分表示深度；第三部分表示三孔隙度曲线，其中的曲线d表示CNC，曲线e表示DEN，曲线f表示DT；第四部分表示电阻率曲线，其中的曲线g表示RLLS，曲线h表示RLLD；第五部分表示孔隙度，其中的曲线i表示POR；第六部分表示饱和度，其中的曲线j表示SO；第七部分表示岩性剖面。根据上述的输出结果可指导储层综合评价。
本实施例提出的测井解释模型的处理方法具有以下的优点：
1)首次将公式解析技术引入到测井资料处理中，成功实现了对用户自定义解释模型的添加和处理，与现有技术相比具有更强的适用性；
2)针对测井解释模型特点，设计了一套模型公式编辑、解析与存储方法，有效提升了表达式解析的准确性与计算效率；
3)提出了一套完整的基于公式解析的测井解释模型输入与处理计算的方法和流程，操作更为方便、灵活。
实施例二
本实施例提供了一种测井解释模型的处理装置，结合图5所示，包括：
模型确定单元51，用于根据预定地层分析资料确定待测地层的地层参数解释模型；
表达式配置单元52，用于根据地层参数解释模型配置参数表达式；
解析单元53，用于将参数表达式进行后缀表达式解析；
处理单元54，用于根据经过解析的后缀表达式确定待测地层的地层参数结果曲线。
可选的，在解析单元53中包括：
变量替换子单元，用于将参数表达式中包含的中间变量进行逐次替换，直至参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线；
类型匹配子单元，用于通过对参数表达式中的元素进行逐次扫描，并根据关键字类型进行匹配，以确定参数表达式中的元素类型；
堆栈存储子单元，用于将参数表达式由中缀形式的表达式转换为后缀表达式，并将后缀表达式以堆栈形式存储。
在模型确定单元51中，地层分析资料可以包括测井资料和实验分析数据。其中的测井资料可以包括自然伽马测井曲线、电阻率测井曲线、三孔隙度测井曲线、元素俘获能谱测井曲线或核磁测井曲线中的至少一种；其中的实验分析数据可以包括岩心岩电分析数据、核磁分析数据、声波分析数据、孔渗分析数据或全岩矿物分析数据中的至少一种。
由于模型确定单元51建立的地层参数解释模型通常包含的变量类型较多，例如可能既包含参与计算的测井曲线，也可能包含有各类模型参数，还可能包含有一些中间变量。另外，变量名称也变化多样，例如同一变量名称既可能对应于模型参数，也可能为采集到的测井曲线，会对后续的参数表达式解析和求解的准确性造成较大的影响。对此，针对地层参数模型的特点，可通过表达式配置单元52将组成地层参数模型的元素划分为运算操作符、普通数字、模型参数符号、中间变量和测井曲线等五类，并将参数符号通过第一关键字表示，将所述测井曲线通过第二关键字表示；另外，地层参数模型中包含的运算操作符可为预先定义的操作符类型；中间变量则可采用与第一关键字或第二关键字类似的方法定义，通常可将解释模型中公用的一些的变量定义为中间参数。
对参数表达式的解析可采用计算机领域的公式解析方法，即通过解析单元53将预定地层的参数表达式从中缀形式转换为后缀形式，并分别以堆栈的方式记录存储。在转换的过程中，首先通过变量替换子单元对参数表达式中包含的中间变量进行逐次替换，直至参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线；然后通过类型匹配子单元对参数表达式中的各组成元素依次扫描，并根据关键字类型进行匹配，识别出参数表达式中各元素类型；最后通过堆栈存储子单元将参数表达式由中缀形式转换为后缀形式并以堆栈形式存储。
对待测地层的每个处理深度点，可通过处理单元54依次从解析单元53建立的堆栈存储中取出元素，当取出的元素类型为模型参数符号或测井曲线时，获取对应处理深度点的模型参数或测井曲线值并代入该后缀表达式计算，以确定预定深度点的处理结果。直到获得待测地层的每个处理深度点的处理结果，即可获得连续的全井段的地层参数结果曲线。根据该地层参数结果曲线即可实现对地层岩性、储层发育情况及流体性质判别，进而指导储层综合评价。
采用本实施例提出的技术方案，通过确定的地层参数解释模型配置参数表达式，并将该参数表达式进行解析及求解，从而根据实际地层评价的需要，灵活构建出适应的流程框架；在处理过程中，也可根据储层类型，自由定义和添加更具有针对性的地层参数模型，储层适用范围更广，并有效提高了处理效率和计算精度。
本具体实施方式是对本发明的技术方案进行清楚、完整地描述，其中的实施例仅仅是本发明的一部分实施例，而并不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有经过创造性劳动的前提下所获得的所有其它实施方式都属于本发明的保护范围。
Claims (7)
1.一种测井解释模型的处理方法，其特征在于，包括：
根据预定地层分析资料确定待测地层的地层参数解释模型；
根据所述地层参数解释模型配置参数表达式；
将所述参数表达式进行后缀表达式解析；
根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线；
所述配置参数表达式包括：
将组成所述地层参数解释模型的元素划分为运算操作符、普通数字、模型参数符号、中间变量和测井曲线；
将所述模型参数符号通过第一关键字表示，将所述测井曲线通过第二关键字表示；
将所述参数表达式进行后缀表达式解析包括：
将所述参数表达式中包含的中间变量进行逐次替换，直至所述参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线；
通过对所述参数表达式中的元素进行逐次扫描，并根据关键字类型进行匹配，以确定所述参数表达式中的元素类型；
将所述参数表达式由中缀形式的表达式转换为后缀表达式，并将所述后缀表达式以堆栈形式存储。
2.如权利要求1所述的方法，其特征在于，所述方法还包括：检查所述参数表达式的合法性，并当检查出不合法的参数表达式时提示用户检测和/或修正。
3.如权利要求2所述的方法，其特征在于，检查所述参数表达式的合法性包括：
通过对所述参数表达式包含的运算操作符、普通数字、模型参数符号、中间变量以及测井曲线进行分类匹配，并检查所述参数表达式的配置是否完整、操作符是否为预定义符号之内、模型参数符号是否使用所述第一关键字表示、测井曲线是否使用所述第二关键字表示以及中间变量是否已经定义；
若所述检查的结果中至少有一项的为否，则确定所述参数表达式不合法。
4.如权利要求1所述的方法，其特征在于，根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线包括：
依次从以所述堆栈形式存储的后缀表达式中取出元素，当所述元素类型为模型参数符号或测井曲线时获取预定深度点的对应模型参数或测井曲线值，并根据所述对应模型参数或测井曲线值确定预定深度点的处理结果；
根据所述待测地层的全部深度点的处理结果建立地层参数结果曲线。
5.如权利要求1所述的方法，其特征在于，所述地层分析资料包括测井资料和实验分析数据；所述测井资料包括测井仪器采集到的自然伽马测井曲线、电阻率测井曲线、三孔隙度测井曲线、元素俘获能谱测井曲线或核磁测井曲线中的至少一种；所述实验分析数据包括岩心岩电分析数据、核磁分析数据、声波分析数据、孔渗分析数据或全岩矿物分析数据中的至少一种。
6.如权利要求1所述的方法，其特征在于，所述地层参数包括用于储层评价所需的地层矿物组分、地层孔隙度、饱和度或渗透率中的至少一种。
7.一种测井解释模型的处理装置，其特征在于，包括：
模型确定单元，用于根据预定地层分析资料确定待测地层的地层参数解释模型；
表达式配置单元，用于根据所述地层参数解释模型配置参数表达式；
解析单元，用于将所述参数表达式进行后缀表达式解析；
处理单元，用于根据经过所述解析的后缀表达式确定所述待测地层的地层参数结果曲线；
表达式配置单元，具体用于将组成所述地层参数解释模型的元素划分为运算操作符、普通数字、模型参数符号、中间变量和测井曲线；将所述模型参数符号通过第一关键字表示，将所述测井曲线通过第二关键字表示；
在所述解析单元中包括：
变量替换子单元，用于将所述参数表达式中包含的中间变量进行逐次替换，直至所述参数表达式中只包括运算操作符、普通数字、模型参数符号和测井曲线；
类型匹配子单元，用于通过对所述参数表达式中的元素进行逐次扫描，并根据关键字类型进行匹配，以确定所述参数表达式中的元素类型；
堆栈存储子单元，用于将所述参数表达式由中缀形式的表达式转换为后缀表达式，并将所述后缀表达式以堆栈形式存储。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511000241.7A CN105443122B (zh) | 2015-12-28 | 2015-12-28 | 一种测井解释模型的处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511000241.7A CN105443122B (zh) | 2015-12-28 | 2015-12-28 | 一种测井解释模型的处理方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105443122A CN105443122A (zh) | 2016-03-30 |
| CN105443122B true CN105443122B (zh) | 2019-01-18 |
Family
ID=55553704
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511000241.7A Active CN105443122B (zh) | 2015-12-28 | 2015-12-28 | 一种测井解释模型的处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105443122B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108074282B (zh) * | 2016-11-08 | 2021-02-23 | 中国石油化工股份有限公司 | 基于三维地质建模软件的物性参数的厚度校正方法 |
| CN107092719B (zh) * | 2017-03-17 | 2020-09-04 | 中国石油天然气股份有限公司 | 注水优势通道识别和微球封堵粒径选择的方法及装置 |
| CN111562630B (zh) * | 2019-02-13 | 2023-06-20 | 中国石油化工集团有限公司 | 基于网格划分的储层参数测井评价方法 |
| CN113153285B (zh) * | 2021-05-13 | 2023-11-24 | 陕西正通煤业有限责任公司 | 一种煤层围岩富水性电性评价标准计算方法 |
| CN113374472B (zh) * | 2021-06-11 | 2024-03-19 | 中国石油天然气股份有限公司 | 一种应用元素录井判别火成岩岩性方法 |
| CN115434698B (zh) * | 2022-09-01 | 2024-05-28 | 西南石油大学 | 基于光电吸收截面指数的地层钙质胶结物含量计算方法 |
| CN117649529B (zh) * | 2024-01-30 | 2024-07-09 | 中国科学技术大学 | 一种基于多维信号分析神经网络的测井数据解释方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
| CN104675392A (zh) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | 一种基于叠前多参数降维的储层岩性识别方法 |
| CN104881435A (zh) * | 2015-05-05 | 2015-09-02 | 中国海洋石油总公司 | 一种基于数据挖掘的研究流程自动化测井评价专家系统 |
| CN104880737A (zh) * | 2015-04-27 | 2015-09-02 | 北京师范大学 | 测井资料识别地下流体类型的多元Logistic方法 |
| CN105089661A (zh) * | 2015-07-16 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种连续判断快慢地层的方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103914027B (zh) * | 2013-01-07 | 2016-07-06 | 上海宝信软件股份有限公司 | 利用表达式解析实现工业现场底层设备监测管理的方法 |
| CN204047597U (zh) * | 2013-10-24 | 2014-12-31 | 林梅叶 | 一种新型自动节能的水产养殖给药装置 |
| CN104866636B (zh) * | 2014-02-24 | 2018-05-01 | 中国石油化工集团公司 | 一种随钻测井数据实时处理方法 |
-
2015
- 2015-12-28 CN CN201511000241.7A patent/CN105443122B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1040245A (zh) * | 1988-07-25 | 1990-03-07 | 石油工业部石油勘探开发科学研究院石油地质研究所 | 一种测井信息的解释处理方法 |
| CN104675392A (zh) * | 2013-12-02 | 2015-06-03 | 中国石油化工股份有限公司 | 一种基于叠前多参数降维的储层岩性识别方法 |
| CN104880737A (zh) * | 2015-04-27 | 2015-09-02 | 北京师范大学 | 测井资料识别地下流体类型的多元Logistic方法 |
| CN104881435A (zh) * | 2015-05-05 | 2015-09-02 | 中国海洋石油总公司 | 一种基于数据挖掘的研究流程自动化测井评价专家系统 |
| CN105089661A (zh) * | 2015-07-16 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种连续判断快慢地层的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 分形测井解释理论方法研究;文环明;《中国优秀博硕士学位论文全文数据库（博） 基础科学辑》;20060115(第01期);A011-1 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105443122A (zh) | 2016-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105443122B (zh) | 一种测井解释模型的处理方法及装置 | |
| CN104879103B (zh) | 一种分层注水效果分析方法 | |
| CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
| CN103603659B (zh) | 利用常规测井资料识别复杂孔隙结构储层流体类型的方法 | |
| CN112133377B (zh) | 一种海域天然气水合物赋存状态的判别方法及系统 | |
| CN110318745B (zh) | 一种沉积微相约束下粒径岩性测井评价方法 | |
| CN105488248B (zh) | 一种深部矿层和岩层的判定方法 | |
| CN105445441B (zh) | 储层含气饱和度确定方法及装置 | |
| CN108897066A (zh) | 碳酸盐岩裂缝密度定量预测方法及装置 | |
| CN106355571B (zh) | 一种白云岩储层质量的确定方法及装置 | |
| CN106353820A (zh) | 一种基于知识模型的测井解释方法 | |
| CN104678432A (zh) | 砂砾岩裂缝识别的方法 | |
| CN105114067A (zh) | 一种岩性测井相的方法 | |
| Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN107762483A (zh) | 一种基于测井曲线的相关系数及包络面积的流体识别方法 | |
| CN107121448A (zh) | 碳酸盐岩储层沥青含量计算方法 | |
| CN112990320A (zh) | 一种岩性的分类方法、装置、电子设备及存储介质 | |
| CN110688781B (zh) | 一种低渗非均质气藏储层测井解释方法 | |
| CN109738955B (zh) | 一种基于成分-结构分类下的变质岩岩性综合判别方法 | |
| CN111236934B (zh) | 水淹级别确定方法和装置 | |
| CN112069444B (zh) | 一种利用测井资料计算储集层试井渗透率方法及计算机 | |
| CN113945992B (zh) | 泥岩、油页岩识别方法、装置、电子设备及介质 | |
| CN116930023A (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 |