伽马测井值，单位为GAPI;GCUR表示与地层年代有关的经验系数，新地层例如取3.7,老地层 例如取2; Vsh表示计算获得的泥质含量，单位为V/V(体积比）。
[0048] 地层灰岩和白云岩含量采用光电吸收界面指数曲线PE计算，并可通过W下公式计 算获得：
[0051]其中，Vsh表示根据式(2)计算获得的泥质含量，单位为V/V;PE表示光电吸收界面指 数曲线，单位为B/E ; PEmax、PEmin表示模型参数，分别为纯灰岩和纯云岩地层光电吸收界面指 数响应值，例如分别取5.1B/E、3.14B/E;Vlime和Vdolo分别表示计算获得的灰岩含量和白 云岩含量，单位为V/V。
[0052]地层孔隙度采用密度-中子交会方法计算，并可通过W下公式计算获得：
[0056]其中，DENma、DE化、LCOR表示模型参数，分别为岩石骨架密度值、流体骨架密度值和 岩石骨架中子值，一般分别取值2.83G/CM3、1.0G/CM3和0V/V;DEN、C化分别表示密度和中子 测井值；0.舶《表示密度孔隙度，单位为V/V; 0?表示中子孔隙度，单位为V/V; 0表示计算 获得的地层总孔隙度，单位为v/v。
[0057]含水饱和度和含油气饱和度可采用电阻率曲线进行计算。为了准确确定该井段地 层岩石电性特征，本实施例对4101.3m岩屯、进行了岩电分析实验，图3所示的是含水饱和度 (Sw)-电阻增大率（I)实验测量结果，其中的实线是采用ArcMe公式获得的拟合结果。结合 图3所示，该井段岩石电阻增大率与含水饱和度变化较好的符合ArcMe关系，其计算公式 为：
[0062] 其中，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。
[0063]步骤22,根据地层参数解释模型确定输出地层参数类型和输出曲线名称。
[0064]根据步骤21中确定的地层参数解释模型W及储层评价需求，可配置测井解释模型 的实际输出结果为泥质含量、灰岩含量、白云岩含量、孔隙度、含水饱和度、含油气饱和度， 输出名称分别为VSH、化IME、VD化0、？03、5胖、50。而相应计算式（1)、巧）、（6)等都为参数计算 中的中间变量，并且由式(7)计算的地层总孔隙度0同时也是式（10)计算含水饱和度Sw的 中间变量。
[0065]步骤23,根据地层参数解释模型配置参数表达式。
[0066]根据步骤21中确定的地层参数解释模型，对各地层参数解释模型进行配置。在地 层参数解释模型中，组成模型的元素可分为五类，即运算操作符、普通数字、模型参数符号、 中间变量W及测井曲线，中间变量可由运算操作符、普通数字、模型参数符号、其它中间变 量W及测井曲线组成。在进行模型公式编辑时，模型参数符号可使用第一关键字"PARAO" 包围起来;而测井曲线可使用第二关键字乂URVEO"包围起来;运算操作符则可为预先定义 的操作符类型;对于中间变量，通常可将解释模型中公用的一些变量定义为中间变量，运样 既能简化模型公式的编辑，还能提高后续计算效率。
[0067] 例如，利用电阻率曲线计算含水饱和度的Archie公式可记为：
[006引SW=(PARA<A〉*PARA<B〉*PARA<RW〉/(CURVE<RT〉*(FAI八PARA<M>)))八（1/PARA<N ?
[0069]其中，A、B、M、N、RW分别为模型参数;RT为测井曲线;FAI为中间变量，通过类似方法 定义，还可W在其它模型公式中重复使用。
[0070]本实施例W泥质含量、灰岩含量及白云岩含量的计算模型为例进行说明，其中，泥 质含量计算模型可进行如下的配置：
[0071]甜=(CURVE<GR〉-PARA<GRMIN〉）/(PARA<GRMAX〉-PARA<GRMIN〉）(12)
[0072] V甜=(2 八（PARA<GCUR>)巧H-1)/(2APARA<GCUR〉-1) (13)
[0073] 其中，SH表示中间变量，GR表示自然伽马曲线值，GRMIN和GRMAX分别表示模型参 数，VSH表示泥质含量，GCUR表示与地层年代有关的经验系数。
[0074]通过式(12)定义了中间变量甜，并在式(13)中调用。
[007引另外，灰岩和白云岩的含量计算模型可分别配置为：
[0076] VLIME=(1-VW)*(CURVE<PE〉-PARA<PEMIN〉）/(PARA<PEMAX〉-PARA<PEMIN〉）（14)
[0077]Vdolo= I-VSH-VLIME (15)
[007引其中，VLIME表示灰岩含量，阳表示光电吸收界面指数曲线，阳MAX和阳MIN分别表 示模型参数，VDOLO表示白云岩含量。
[0079]式（13)配置的泥质含量计算模型既可作为输出结果，同时也可作为其它参数计算 模型所需的中间变量。
[0080]同理，可对式(5)-(13)进行配置，配置方法和过程同上，本实施例不再寶述。
[0081] 步骤24，检查参数表达式的合法性。
[0082]本步骤是对步骤23中配置完成的各参数表达式的合法性进行检查。通过对参数表 达式中的各元素进行逐次匹配，识别出参数表达式中包含的操作符、模型参数符号、中间变 量、测井曲线等，同时对参数表达式中的左右括号的完整性进行检验。
[0083] W步骤23中配置的式（12)为例，若中间变量甜定义为：
[0084]甜=(GR-PARA<GRMIN〉）/(PARA<GRMAX〉-PARA<GRMIN〉） (16)
[008引则在合法性检查过程中，通过关键字匹配，字符GRMIN与GRMAX会被自动判别为模 型参数，而字符GR会被识别为中间变量，但同时对已定义的中间变量进行查找，未发现该中 间变量名，则字符GR会被判定为非法，提示用户检查和/或修正。
[0086]步骤25,将参数表达式进行后缀表达式解析。
[0087] 根据计算机领域的公式解析方法，将步骤24验证合法的参数表达式从中缀形式转 换为后缀形式，并分别W堆找的形式记录存储。在转换过程中，可首先对参数表达式中包含 的中间变量进行逐次替换，直至该参数表达式中只包括运算操作符、普通数字、模型参数符 号和测井曲线。然后，通过对参数表达式的各组成元素依次扫描，并根据关键字类型进行匹 配，识别出公式中各元素类型。最后将该参数表达式由中缀形式转换为后缀形式并放入堆 找存储。
[0088] W泥质含量计算模型的式（13)为例，首先对式（13)中包含的中间变量SH进行替 换，替换后表达式为：
[0089] V甜=(2A(PARA<GCUR〉*(CURVE<GR〉-PARA<GRMIN〉）/(PARA<GRMAX〉-PARA<GRMIN 〉））-1)/2八PARA<GCUR〉-1) (17)
[0090]替换后，式（17)中不再包含有中间变量，则可对该参数表达式进行元素类型匹配 并转换为后缀表达式。对该参数表达式从左至右依次扫描，其中的元素"("、"2"、("分 别为括号、数字、操作符、括号，直接依照中缀表达式转后缀表达式规则入找，而元素PARA< GCUR〉通过关键字匹配为模型参数，入找后其名称为GCUR，元素类型为2(参数符号）。同理， 对上述参数表达式扫描并转换完成后，后缀表达式堆找结构可表示为：
[0093] 步骤26,根据经过解析的后缀表达式确定待测地层的地层参数结果曲线。
[0094] 在实际处理过程中，可对每个处理深度点获取当前深度各参数表达式中包含的测 井曲线值和模型参数值。对每个参数表达式，依次从步骤25存储的堆找中取出元素，当元素 类型为模型参数符号或测井曲线时，则将具体值代入参数表达式计算，依照后缀表达式值 计算方法获得该参数表达式在当前深度点的结果并输出。对全井段处理后，即可得到连续 的地层参数结果曲线。
[0095] 根据步骤22的配置，结合图4所示，本实施例提出的测井解释模型的处理方法的输 出结果包括：左起第一部分表示岩性曲线，其中的曲线a表示GR，曲线b表示PE曲线C表示 CAL);第二部分表示深度;第S部分表示S孔隙度曲线，其中的曲线d表示CNC，曲线e表示 DEN，曲线f表示DT;第四部分表示电阻率曲线，其中的曲线g表示化LS，曲线h表示化LD;第五 部分表示孔隙度，其中的曲线i表示P0R;第六部分表示饱和度，其中的曲线j表示SO;第屯部 分表示岩性剖面。根据上述的输出结果可指导储层综合评价。
[0096] 本实施例提出的测井解释模型的处理方法具有W下的优点：
[0097] 1)首次将公式解析技术引入到测井资料处理中，成功实现了对用户自