GeoFrame软件平台在测井解释

GeoFrame软件平台在测井解释中的应用争论中石化星公司西南石油局测井公司二OO三年三月GeoFrame软件平台在测井解释中的应用争论工程名称：GeoFrame软件平台在测井解释中的应用争论托付单位：中石化星公司西南石油局科技处担当单位：中石化星公司西南石油局测井公司工程负责：徐炳高林绍文报告编写：徐炳高何传亮2023年6月西南石油局测井公司目录\l“_TOC_250018“前言 1\l“_TOC_250017“争论任务 1\l“_TOC_250016“任务完成状况 2\l“_TOC_250015“1.4主要争论成果 3\l“_TOC_250014“2、软件概述及流程建立 3\l“_TOC_250013“软件概述 3\l“_TOC_250012“进入GeoFrame软件 5流程建立 9\l“_TOC_250011“3、测井数据库治理 10\l“_TOC_250010“数据预备 10\l“_TOC_250009“数据分析 11\l“_TOC_250008“数据输入 11\l“_TOC_250007“数据输出 14\l“_TOC_250006“数据治理 15\l“_TOC_250005“4、测井资料预处理 17\l“_TOC_250004“测井曲线的预处理〔WellEdit模块〕 17\l“_TOC_250003“启动WellEdit程序 18\l“_TOC_250002“根本曲线编辑 20\l“_TOC_250001“4.4.4、曲线充填显示 23\l“_TOC_250000“地层层位比照〔WELLPIX〕 245、测井资料处理解释 28、测井资料标准化 28测井曲线交会图分析〔geoplot、utilityplot、crossplot模块〕 29岩性与含气性识别 29测井数据的多微显示 31岩相分析〔rockcell〕 33ElanPlus处理 36根本原理 36、川西地区碎屑岩EANLANPLUS方程建立 40川西场气田ELANPlus处理 41ELAN计算参数结果统计〔RESUM〕 45ELAN处理与常规处理结果比照 48第六章、测井地质建模 48剖面图显示 48平面图显示 497、结论与建议 65前言《GeoFrame2023年5月份成都龙泉培训，2023年12月美国休斯顿培训，2023年3月江苏扬州培训。通过培训和长时间对软件的摸索和学习，现在对软件有了较好的生疏和把握，本工程的争论是在对软件不断生疏和把握的根底上开展，同时，通过工程的争论，进一步加深了对软件的理解和把握。争论任务按工程设计，工程需要完成以下7项内容的争论。1、数据加载。软件平台的数据加载格式与现行使用的测井数据格式不同，因此，首先要争论解决两个系统数据格式的原理以及格式间的相互转换，最终目的是将现行使用的测井数据便利快捷地加载进GeoFrame软件系统。2(UtilityPlot&GeoPlot)牢靠的处理解释参数。3、依据测井公司资料采集设备的状况，重点争论测井资料精细评价模块〔ElanPlus〕在测井处理解释中的应用：如何准确确定地层岩性，尤其在简洁岩性地层〔低孔低渗/灰岩地层；应用已有的区域地质分析资料，如何精细计算储层参(孔，渗，饱等)；如何有效判别储层流体性质；如何实现从单井精细评价到多井区域争论的转换。等，并依据比照结果进展地质特征及储层参数统计〔ResSu。5、利用上述解释成果，结合地震解释，绘制过井油藏剖面图。 StratClass)。6J2S层位的57、将ElanPlus处理成果与常规测井资料处理成果进展比照分析。以得到准确的测井解释分析结果，有效地削减测试层数，并确保不漏失有工业价值的储层。任务完成状况本课题的任务完成状况如下：〔一、模块功能争论与流程建立要使用该软件，需要对GeoFrame软件系统主要模块的功能进展生疏把握，并学会建立合理有效的处理流程。分析争论了GeoFrame软件的各模块功能，对各模块的主要功能有了较具体的了解。依据现有测井资料建立起了适合现有测井设备系统的测井处理解释流程。〔二、数据分析与预备同的测井设备具有不同的采集数据格式。目前，西南石油局测井公司多承受国产测井设备，野外采集系统的测井数据格式与GeoFrame软件的数据格式不全都。因此，需要分析数据格式，解决数据格式间的相互转换问题。为此开发了相应的数据格式转换软件来实现数据的正确输入输出。〔三、主要处理解释程序争论GeoFrame软件系统主要模块的处理程序是整个软件系统的核心，因此对它12口井沙溪庙组地层开展争论，具体争论如下：争论数据装载程序(DataLoad&ASCIILoad)，将场沙溪庙测井数据加载进入数据库。争论测井资料的预处理程序(Preplus，WellEdit&WellCompositePlus)。对测井资料进展有效的预处理，包括环境校正，深度匹配，曲线滤波等。争论GeoFrame的交会图程序(UtilityPlot&GeoPlot)，为测井处理解释的参数确定供给牢靠的依据，包括泥质参数、骨架参数等。〔ELANPlus〕的原理及在测井处理解释中的应用。包括具体争论ELANPlus〔低孔低渗/灰岩地层〕的岩性；应用已有的区域地质分析资料，如何精细计算储层参数(孔，渗，析结果建立合理的多井解释模型，从而实现从单井精细评价扩展到多井区域的争论。进展了多井地层比照WellPi数统计ResSum。(StratLog)。本工程圆满完成了设计和合同要求的各项争论任务。在某些方面工作量超过11口。设计中模块的应用争论，实际争论中认为，这局部内容较重要，开展了较多的争论工作。主要争论成果GeoFrame软件各模块的主要功能，建立了测井处理解释流程。测井数据格式转换为GeoFrame软件认知的数据格式。现已将场地区沙溪庙十余口井的测井数据转换为GeoFrame数据格式。GeoFrame软Oracle理争论打下了根底。解释模块的分析争论。处理成果图。作图也未进展。2、软件概述及流程建立软件概述2023年从斯伦贝谢公司引进了GeoFrame程的很多学科，是多学科多专业学问穿插结合的集大成者，具有世界领先水平。该软件是基于SolarisOracle进展资料处理和解释的软件系统。整个软件系统由GeoFrame处理软件，Geonet挂靠工具，Oracle数据库，FLEXLM许可证治理器，Fortran和C/C++语言编译器等构成。我们常说的GeoFrame软件是指它的处理软件和Geonet界面，而其它局部都在后台运行。图图2-1 应用程序治理界面该软件处理解释的学科涵盖了物探、地质、测井〔岩石物理〔图2-。结合测井资料计算的地震合成记录，对地震勘探资料进展层位标定；地震资料三维可视化成图等。地质及测井成像等资料可以对某个地区的构造进展分析争论；进展区域地层比照；沉积相争论；井眼地质分析和井旁地质构造外推等。单井作出油气层精细解释、薄层分多井处理解释模型的结果，对某区域的地层和油气储层等进展分析比照及储量计算。在油藏工程方面，可以依据测井资料并结合物探、地质等方面的资料，对某个区域的油气层在横向和纵向上的展布进展数字模拟，对潜在的有利油气聚拢带作出定性及定量的分析推想，并可三维可视化成图；可对试井资料，重复地层测试器资料进展处理解释。此外，还可依据测井资料计算岩石机械特征〔动态弹性模量、静为邻区将来安全稳定地钻井供给合理的泥浆比重范围；进展裂开压力分析，为油气层的加砂压裂供给牢靠依据；进展砂岩强度分析，为油气开采供给储层在不出砂所能承受的最大压降等。GeoFrame软件功能强大，其处理手段先进，处理方法多样，操作便利灵敏。释。绘出的各种图件精巧。不过，正由于该软件是多学科学问穿插运用的结晶，其系统也格外浩大和简洁根底学问。GeoFrame软件在测井解释中的应用。进入GeoFrame软件GeoFrame软件进展测井资了的处理解释，需要经过如下的过程：〔一、进入Geonet界面：Geonet是一个应用软件挂接工具。它供给了运行地学科学应用软件的界面以及在软件之间进展数据传输的工具〔图2-。图2-2Geonet启动界面Geonet是在软件系统完成安装后，进入指定的用户,自动弹出该界面的。假设没有启动该界面，点击鼠标左键选择快捷菜单的Geonet选项就可进入到该界面。从Geonet的启动界面上可以看到，主要挂接了GeoFrame3.8.1软件和Geoshare数据传输工具。固然，还可以挂接其它软件。GeoFrame3.8.1GeoFrame应用程序管Oracle〔建议对Oracle有确定了解后再进入使用；第三项翻开GeoFrame软件系统的联机帮助文献，可以在此查到该GeoFrame软件是一的帮助信息〔2-。图图2-3联机帮助界面的主要功能是进展不同软件或不同工区工程间数据的输入输出，便于大量的数据在Oracle中的传输。进入工区〔工程〕治理器：GeoFrame软件的工区〔或工程〕治理器。要进展一个区域的争论或要进展某一个科研工程，就可以建立相应的工区。这样做的目的是，便于在Oracle数据库中对数据进展统一治理，数据调用便利牢靠，各种资料数据能有机结合〔图2-。要选用和争论区域相对应的投影坐标系统，这样作出的图件才能与实际区域的相一致。我们已建立了三个工区，登录进入场工区〔图2-。图2-4 器图2-5 登录窗口器ApplicationManage〕就可进入GeoFrame软件的各个模块了〔图2-。GeoFram〔图2-可以看到，GeoFrame软件的应用程序治理器由六大处理分析模块和三大治理工(通常把这两个模)，油藏工程模块，地震处理模块，可视化成图模块和有用工具模块。三大治理工具模块为：进成治理模块，数据治理模块和工区治理模块。2.3流程建立软件后要进展测井资料处理，首先就要建立工作流程。建立合理的工作流程，不仅使处理解释工作的思路清楚，完成工作快速准确，而且该处理流程同样可用于以后的井，它提高了工作效率和工作质量。原始测井资料的验收，原始测井数据的格式转换，现。建立合理的处理流程是一件很价值和意义的事，的功能。在计算中心可建立流程来处理常规测井资料。GeoFrame软件的进程治理器(ProcessManager)的主要功能就是建立处理解释的。现场测井数据经格式转换后，依据从测井数据输如图2-6所示。可以看到，进程治理器主窗口右边有三个小图标按钮。第一个按钮的作用是数据聚焦(DataFocus)，所谓数据聚焦，就是选取什么样的数据来进展处理，可选的数据有油田数据，井场数据，井眼数据等。其次个按钮的作用是选取处理程序或默认的处理流程链。第三个按钮的作用是用户建立自定义处理流程链；在全部程序中选取一部份组成处理流程链，并可进展保存，编辑修改或删除等；可以把选取程序形成的流程链保存下来用于以后的处理工作，以提高工作效率。断程序进程，这在程序运行出错时把握程序的运行状态是格外有用的，比直接在Solaris下中断进程更安全快速。中的挨次符合数据流向的规律，能使工作省时省力即可。图图2-6 进程治理器主窗口及处理流程链至此，我们已经建立资料处理解释的流程链。我们下面的工作，主要探讨怎样应用流程链中的程序来实现具体的处理解释。3、测井数据库治理数据预备GeoFrame软件进展区域争论工作及测井处理解释时，同样需要大量的预备工作。预备工作包括：对区域的坐标系统进展。有了这些能结合多种资料，解释结果更准确牢靠。为了更好的进展测井处理解释工作，我们除了收集常规的测井曲线资料外，还需要收集岩芯分析数据等。数据分析司的MAXIS-500测井设备承受的是DLIS格式；西方阿特拉斯公司的ECLIPS-5700测井设备承受的是XTF格式；哈里伯顿公司的EXCELL-2023测井设备承受的是LIS格式。此外，还有很多其它的测井格式。西南石油局测井公司野外测井设备多为国产设备数据格式与GeoFrame软件所支持的数据格式不同。为了能把我们现有的测井数据加载进 Oracle数据库供GeoFrame软件进展处理必需进展格式转换为此编制了测井数据格式转换程序，使所转换的数据能准确加载入库。数据格式转换程序的准确性已得到验证。假设是斯伦贝谢公司的测井数据，可以直接加载入库，而其它公司的测井数据需要转换成GeoFrame能认知的数据格式，不再赘述。数据输入Oracle数据库了。通常，用如下的两个程序将数据加载入库：〔一〕ASCII码数据加载—ASCIILoad对于ASCII数据，可以承受ASCIILoad程序将其加载进入Oracle数据库。该程序能加载多种类型的数据，加载方式灵敏多样。进入ASCIILoad程序：位坐标数据，井斜数据，地层倾角数据，岩相分析数据，岩芯分析数据和岩芯照片等数据。由于各种数据的排列方式不一样，因此加载的方式也较灵敏。在加载过程中，需要输入的主要有如下项：原始数据文件，格式把握文件，出存储坐标系统，数据采集者，记录日期格式等〔图3-1。有些项是可选项，依据需要选取。图3-1 ASCII Load程序主界面图图3-2ASCII Load把握文件关键是要选择确定数据类型并创立好相应的加载把握文件。把握文件的主要目的是描述ASCII数据的排列方式，以使系统按该排列方式正确地3-2，主要的输入项包括：加载的起始行，完毕行，数据记录长度，测井曲线的条数，深度单位及起止深度，采样间隔，输入曲线的名称及对应的单位和列数等。对于不同类型的数据，其把握文件虽然不同，但加载的方式都大同小异。〔二〕二进制数据加载—DataLoad对于测井数据为DLISDataLoad加载程序将其加载进入Oracle数据库。这种加载方式相对简洁，由于很多油田信息都在数据文件的头部赐予了定义，从而不用创立把握文件。DataLoad程序要求的输入项主要有：原始数据文件，别名文件，油气田名称，井场及井眼名称，数据采集者等。程序运行的界面如图3-3所示。图3-3 DataLoad程序它可以从本地磁盘文件，磁带或通过网络进展数据加载，我们一般选用从本地件信息是否正确。数据输出GeoFrame软件供给了一个很便利有用的程序—DataSave程序进展数据输出。运行界面如下：图3-4 数据保存Oracle可以选择多种输出数据格式，并可确定输出的深度范围和采样间隔〔图3-。数据治理在GeoFrame如此简洁浩大软件系统中，假设没有良好的数据治理功能，那么数于各种数据的治理工作。进入数据治理器Oracle数据库直接相联系，可以对Oracle数据库中的数据进展查询，编辑修改及删除等工作。使用这些数据治理工具进展数据删除时确定要留神，只有确定要删除某些数据对象前方可删除。由于一旦删除了数据，在数据库中是不行再恢复的。比方我们要场工区现有井眼的各种信息，并可进展修改〔图3-。图3-5井场和井眼数据治理器界面型的数据进展治理，这就通用数据治理器〔GeneralDataManage〔图3-。括了其它数据治理工具并有一些功能。如何进入？图图3-6GeneralDataManager程序主界面在通用数据治理器里，可以用列表或图形的方式对全部数据进展查询，还可进能：ITC数据传输和工区间数据传输。主界面右下角的两个按钮用于ITC数据传输。该功能主要用于在不同处理程序间进展数据传输，便利快捷。除了使用Geoshare进展数据传输外，应用通用数据治理器中的工区间数据传输间进展数据的输入输出，比Geoshar操作更便利〔图3-。图3-7工区间的数据传输4、测井资料预处理Edit模块〕在将数据加载进GeoFrame软件后，通常需要对数据进展预处理,包括曲线编辑、曲线校正等。该预处理工作可在WellEdit程序中完成。析和编辑(或拼接或深度匹配)测井与岩芯数据，对其进展精细加工预处理。等，可进展横向或纵向的拉伸压缩等。置某段测井曲线值为常值，进展基线移动，测井曲线滤波(去尖峰或移动平均)，曲计算等。因此，它包含了多种对应的编辑模式。、启动WellEdit程序、为每口井眼选择一个文档模板(5)、保存图形显示作为一个模板文档(9)、进展其它编辑(10)、保存编辑后的测井曲线到数据库启动WellEdit程序WellEdit中显示的一个ASCII码可调整多井眼显示的挨次，还可通过ITC进展井眼选择(图4-1)。进入WellEdit程序主窗口，我们可以看到在左边有二十多个图形小按钮。每个按钮对应不同的编辑处理模式。这包括：深度移动模式，基线移动模式，滤波器模模式，数据函数。图4-1启动窗口图4-2 WellEdit主窗口根本曲线编辑1、深度移动模式图4-3曲线移动编辑具有比较强大的深度移动功能，可以进展整体曲线的上下移动〔深度移动刻度变化能。2、基线偏移模式基线漂移主要针对自然电位曲线进展〔图4-自然电位曲线的基线指大段泥岩地层比较平直的段，以其为基准，高于该值为正偏，低于该值为负偏。在曲线的测量过程中，由于各种缘由可能导致其基线不是垂直的，需要对其进展偏移校正，使基线处于直立状态，有利于对自然电位特别层位的推断。 进展基线偏移在简洁，首先选择基线漂移的井段，然后定义不同井段的基线点，最终执行基线偏移程序将这些基线点定义为0值曲线形态得到修正完成曲线的基线偏移。3、滤波模式消退曲线的统计起伏，去掉曲线中无用的噪音信息。可以对滤波的方法作选择，比方平均滤波、加权滤波等。同时，可以选择滤波时采样点数，依据滤波幅度的要求作不同选择。图4-4 基线偏移校正模式4、插入/替换模式形成一条较好的完整的曲线。这个功能在常规曲线的DOS编辑中虽然也能实现，但不如WELLEdit简洁和易于操作。5、重采样模式0.1及0.12510个点或8设定采样率，进展曲线的重采样。6、平均模式曲线。、深度匹配〔曲线深度对齐〕(图4-5)。首先选择根底曲线，其它曲线以此为标准进展深度移动，然后选择一条需要进展深度移动的曲线以及与其深度系统全都的其他曲线。深度移动有手工和自动移动两种方式，手工移动用鼠标点击，把曲线的最相像的特征点连接起来，也可以把相运行程序后，曲线间深度自动匹配。对于曲线相关性较好的曲线可以用这种方法，但相关性差的曲线间不宜用这种方法。图4-6为以GR为根底的自动深度对齐图，其中1、2、3道为原始曲线图，横线为自动比照线，假设手工比照时，需要进展手工拾取曲线间的比照标志点，然后联线。线比照相关性判别结果，可以看出，相关性较好。图4-5 曲线深度匹配界面4.4.4、曲线充填显示

深度配窗口图4-6 深度自动对齐无限充填显示为曲线的一种显示功能，不是编辑功能。该功能挂接在WELLEDIT模块上，在这里作简要介绍。道边界间充填，充填的图案和颜色可以任意选择。曲线充填是一格外有用的工具。与密度曲线间的充填反映出中子的自然气挖掘效应〔4-。图4-7 曲线充填显示2 2 2 2 场沙溪庙组发育有5〔Js1Js2Js3-1Js3-JsWELLPIX2 2 2 2 2其中maker相当于单井分层点，zone为单井两个maker间的段，surface为多口井一样maker组成的面，layer为两个surface间的层。比方对于Js1砂体，单井的砂岩顶底叫2maker，单井的Js21砂体叫zone，区域的Js21砂体顶底面叫surface，区域上的Js21砂体。中常用的14项功能，下面简要介绍主要功能。图5-1 WELLPIX程序功能1、曲线的假彩色显示的特征，对层位比照是一格外适用的工具〔图5-。2、曲线截取比照定出比照点。3、曲线maker标注个点的层位比照。也可以在深度的状况下，输入深度数据。本次工作进展了〔Js21、Js22、Js23-1、Js23-2、Js24〕五层砂体顶底共10个maker的层位比照和标注〔图5-。图5-2 曲线的假彩色显示4、层拉平显示maker连接起来，成为剖面的比照线，可以以任一个maker为标准进展一个面的拉平，可以显示古构造特征，各井也可以以一样的深度比照，反映今构造特征。图5-4为以JS21顶为标志的层拉平显示。5、创立layerlayelaye的顶底makersurfac，确定出layer，即Js21Js2Js23-1、Js3-2Js24五层砂体。6、比照结果5-1。图图5-3maker标注图5-4 Maker的层拉平显示井号maker井深(米)井号maker井深井号maker井深(米)井号maker井深(米)js21-top2245.6js21-top2163.1js2-1bottom2270.0js2-1bottom2187.3js2-2top2324.5js23-1top2274.5js2-2bottom2344.3js23-1bottom2285.2js23-2top2380.3js23-2top2295.4js23-2bottom2398.8js23-2bottom2314.3js2-4top2436.9js2-4top2355.6js2-4bottom2454.9js2-4bottom2371.2js21-top2211.5js21-top2170.4js2-1bottom2230.4js2-1bottom2194.2js2-2top2280.3js2-2top2248.5js2-2bottom2306.4js2-2bottom2274.6js23-2top2337.4js23-1top2279.8js23-2bottom2352.0js23-1bottom2294.4js2-4top2399.5js23-2top2300.3js2-4bottom2414.4js23-2bottom2314.5js21-top2251.3js2-4top2356.5js2-1bottom2277.0js2-4bottom2378.1js2-2top2329.0js21-top2215.9js2-2bottom2356.7js2-1bottom2239.3js23-2top2392.9js2-2top2294.4js23-2bottom2403.9js2-2bottom2318.8js2-4top2451.9js23-2top2349.1js2-4bottom2459.9js23-2bottom2364.3js21-top2161.9js2-4top2410.9js2-1bottom2179.6js2-4bottom2428.0js21-top2138.2js21-top2187.9js2-1bottom2153.8js2-1bottom2212.0js23-1top2238.3js2-2top2265.4js23-1bottom2257.3js2-2bottom2295.5js23-2top2262.4js23-2top2325.8js23-2bottom2281.2js23-2bottom2338.4js2-4top2322.0js2-4top2383.6js2-4bottom2343.9js2-4bottom2399.9js21-top2130.8js21-top2191.7js2-1bottom2151.5js2-1bottom2215.5js23-1top2226.4js2-2top2264.8js23-1bottom2245.8js2-2bottom2294.5js23-2top2250.3js23-2top2319.0js23-2bottom2268.7js23-2bottom2337.8js2-4top2315.0js2-4top2385.5js2-4bottom2329.3js2-4bottom2399.35、测井资料处理解释、测井资料标准化测井曲线在测量过程中，由于井眼环境的差异，造成曲线测值存在误差，影响资料的进一步处理，因此需要对曲线进展标准化。方法是选择测量质量较好的曲线作为标准井，制作标准模板，其它井以此为依据，将曲线刻度到与标准井一样的值域区间，使之具有相近的均值和方差。标准化工作主要争对GR、CNL、DEN在场地区选择了11口井，通过比照分析觉察CX164井测井曲线质量相对较好，将其作为标准井，对其它井进展标准化，这个过程在utilityplot模块上完成。标准化工作通过挪动曲线直方图来完成。以标准井为基准，挪动目标井的直方曲线直方图可以自动进展最正确匹配。另一种方法是手工挪动，可以进展曲线直方图〔转变均值〔转变直方图方差。通过手工挪动同样可以到达曲线的最正确匹配效果。挪动直方图的同时，相对应的曲线值得到转变，到达标准化的目的。图5-1为CX164井与CX165井标准化前后的直方图比照，可以觉察，标准化前两组地层曲线进展了标准化。测井曲线交会图分析〔geoplot、utilityplot、crossplot模块〕的交会图分析功能较强大，在很多模块上都挂接了交会图分析程序。交会图可以进展岩性识别、含油气性识别等，能进展二维平面及三维立体显示。5.2.1岩性与含气性识别具有明显的“挖掘效应”。交会图与测井曲线间可以进展相互通讯〔5-Z值〔GR〕交会图上识别出气层的范围，并且将其圈定起来，在剖面图上可以明显显示出气层的井段。气层表现为低GR、低中子孔隙度、低密度特征，手工用鼠标圈定〔交会图上可以以不同的颜色同时显示出不同储层的位置，有利于含气响应模式的建立。5-3步显示出泥岩的区域。图5-2气层的交互识别图5-3泥岩层的交互识别对于简洁岩性剖面，比方同时有砂岩、灰岩、白云岩的剖面，岩性识别比较困5-2及图5-3交会图模板中有3条岩性识别及孔隙度确实定。5.2.2测井数据的多维显示geoplot模块可以进展二、三、四维显示，最多可显示四微，包括三维的空间显示和颜色显示。在一个屏幕上可以同时显示2、4、6、9张相互独立的统计图形，图5-4同时显示出了三口井〔X850、X806、X807井〕四张图，左上为中子密度的二维交会图，右下为声波、中子孔隙度二维交会图，不同颜色反映了不同的伽玛值，该图显示出3口井的中子孔隙度与目的曲线的总体特征；左下为三口井的GR频率分布直方图，不同的颜色代表不同的井，反映了三口井GR数据分布规律；图5-4右下及图5-5为三维可视化图，颜色代表不同的自然伽玛值，ZX、Y轴分别为中子孔隙度和密度。立体图形的显示格外生动，可以进展手工和自动的全3口井中，以X805井的高GR成分较多，反映泥岩较纯，曲线值反映为高中子孔隙度和中等密度。图5-4 多窗口统计图的显示图5-5 交会图三维可视化显示岩相分析〔rockcell〕岩相分析是利用神经网络对的岩相进展网络学习，建立网络模型，然后对其它井的岩相类型进展推想。区类砂泥岩剖面的岩石分析数据较少，依据试验分析结果建立岩相识别网络模型比较困难。这里选择CX164井的沙溪庙组作为网络学习样本，依据曲线特征，人工建立了砂岩、粉砂岩和泥岩的样本，利用的测井曲线有薄、GR值相对较低的层段作为粉沙岩样本，把高GR、低RT〔LLD〕层段作为泥岩样本。图5-6：网络学习样本选择1、网络参数设置在样本选择根底上，进展网络参数设置〔图5-，网络参数包括：网络层数：最多可以设置3层网络，本次设置了2层网络。网络节点数：每层网络均设置10个节点。样本筛选：样本数太多会影响运算速度，可以对样本进展均匀筛选。本次选择50%样本参与网络训练。训练样本与检验样本安排比照，以检验模型牢靠性。本次将筛选出的样本各50%分别用于样本训练和样本检验。网络学习次数：以网络学习能够收敛为止。本次网络学习次数为460次。图5-7：网络参数设置2、网络训练参数设置完成后，进展网络训练〔学习〔图5-，随着学习次数的增加，误差5-8在学习过程中，可以观看两条误差线的动态变化过程。随着学习次数的增加，误差渐渐降低，两条线几乎重合，误差线的尾部较平直，说明收敛效果较好。由此建立起网络模型。图5-8 网络训练过程及误差显示3、岩相网络计算利用建立的网络模型对区内的11口井沙溪庙组进展岩相的网络计算，分别得到各井的岩相〔性〕剖面，图5-9列出了其中5口井的剖面结果，可以觉察剖面岩性与自然伽玛所反映的岩性特征是全都的。图5-9：川西场沙溪庙组岩相判别结果ElanPlus处理5.4.1根本原理1、ELAN原理概述ELAN处理的根本原理见图5-10，其中：t：工具矢量，代表测井仪器数据，即测井曲线〔包括计算的测井曲线；v：体积矢量，代表各组分的体积；R：反映距阵，代表每条曲线的参数值，全部体积值的和为100%。对于三角形V和R计算t叫做“正问题”t=Rv，由t和R计算叫反问题。在EALN处理中，首先是解决“反问题”，既利用t和R计算v。这是ELAN计算的主要工作。R和v计算t，利用计算出的各成分体积及曲线参数值重计算出条测井曲线。将计算出的曲线与原始曲线进展比照以验证计算的体积的牢靠性。图5-10 Elanplus处理原理图解2、解释模型ELAN解释模型包括4个局部，组成成分、响应方程、参数、约束条件。组成成分指需要计算的成分的体积，包括三个组成成分：矿物、岩石和流体。可以通过CaSO4等，对砂泥岩地层而言，一般选择石英和伊利石进展计算。岩石有多种矿物组成，是多种矿物的组成的混合物，砂泥岩地层中选择砂岩和泥岩进展计算〔岩石需要计算的体积参数个数有计算方程个数确定，体积的个数不能多于方程数。洗带地层水、气。计算有六个参数体积，需要至少六个方程。响应方程响应方程是一个数学关系式。如下式：n测量值i1

ViRiVi：第i成分的体积；i成分的响应值。系。〔测井曲线〕数据。比方对于密度测井曲线〔DEN〕DEN=POR流体密度+〔1-POR〕骨架密度在表达式中全部成分体积的和为1，即：n Vi1=i1这个关系是ELAN处理中自动包含并执行了的方程。侵入模式程序中，有一个根本的假定，即钻井剖面由冲洗带〔X带〕和原状地层带〔U带〕组成，并且假定冲洗带与原状地层的的孔隙度相等。线反映原状地层信息。入因子为1.0表示该曲线完全受原状地层影响。表5-2列出了ELANPus所假定了各曲线所反映的地层深度。密度。表5-2 曲线侵入特征方程数与待求成分体积个数间的关系在联立求解EALN方程中，方程的个数不能少于所求未知数〔成分体积〕的个数。当方程数与未知数一样时，方程有唯一解，我们把这种状况叫做“完全限定”sistemtobedetermined；当方程数少于未知数时，方程有无穷解，没有意义，把欠限定〔nuderdetermine叫“超限定”〔overdetermined，从严格数学意义上讲，超限定下方程无解，是由于不同的曲线或曲线组合计算的成分体积含量有可能是不一样的。单实际状况是，在多数ELAN方程中需要计算“超限定”下的解。因此需要解决不同方程计算参数中的不协调问题，如不同方程或方程组计算的同一参数值不全都时，如何确定计算的参数值。为了更好的理解ELANPlus方程求取参数的原理，我们设想在一个直角坐标系中，一条直线与一系列点的匹配状况。如图5-11左所示，A、B两点相当于方程个数，线段AB相当于需要求取的未知数〔某一成分体积。两点可以唯一确定一条直线，这相当于“完全限定”的状况。当数据点多于2个时，一般状况下，这些点不会在一条直线上，特别是在有的曲线在测量过程中记录了一些“噪音”信息的状况下。这属于“超限定”的状况，通常用线性回归的方法来确定一条回归线〔图5-1右，这相当于所求未知数的解。简洁的回归方法，可以对不同的点赐予不同的权重。图5-12的数据与图5-11右侧图的数据是一样的，细的回归线与上图的全都，为等权重的回归。粗线表示、点的权重很大〔100，而另外两个点很小〔1〕得出得一条回归线。权重可以换算为“不确定性”〔uncertainty，如图5-12中A、D点的不确定性换算为1，B、C换算为0.01。ELANPLus程序对每个响应方程赐予了固定的“不确定性”值。当把不确定性转—权系数weightmultiplier最终的权重代入方程进展最优化运算。资料处理中可以直接引用这些隐含值，也可以修改这些值。图5-12 不同权重的线性回归5.4.2、川西地区碎屑岩EANLANPLUS方程建立GRA〔DCN〔NPHDE〔RHOR〔LL、RS〔LLS〕计算，计算参数包括石英〔QUAR、伊利石〔ILL、原状水UWA、原状气〔UGA、冲洗带水〔XGA、冲洗带气〔XGA。可以列出以下方程组：GR=V

*GR +V

*GR

+V

+V *GR

*GR +石英V *GR

石英

冲洗带水冲洗带气DT=V

冲洗带气*DT

+V *DT +V

DT +V

*DT +V

*DT +石英 石英V *DT

冲洗带水冲洗带气 冲洗带气NPHI=V *NPHI +V

*NPHI +V

NPHI +V

*NPHI +V*NPHI

石英 石英+V

冲洗带水冲洗带水RHOB=V

冲洗带气*RHOB

+V

*RHOB

+V

+V *RHOB +V石英 石英

伊利石

原状气*RHOB +V *RHOB冲洗带水LLD=V

冲洗带水*LLD

冲洗带气+V

冲洗带气*LLD +V

LLD +V

*LLD +V

*LLD石英 石英+V *LLD

冲洗带水冲洗带水

冲洗带气LLS=V

LLS +V

LLS +V

LLS +V

LLS +V

*LLS石英 石英+V *LLS

冲洗带水冲洗带水V

冲洗带气+V

冲洗带气+V +V +V =1石英 伊利石V +V

=V

+V

冲洗带气原状水

冲洗带气注：GR

：GR的石英骨架值*相应权系数石英V ：石英体积百分比石英最优解。通过以下方法来实现。方程的最优解，即得出计算结果。其中，总误差GR测量GR计算GR权系数= GR不确定性

RHOB测量RHOB计算RHOB权系数RHOB不确定性

+．．=〔测量最大值-最小值〕1.5%。当测井曲线质量较差时，其误差较大，可以给定其较小的权系数使总误差变小。5.4.3川西场气田ELANPlus处理1、参数设置首先进展参数设置，对方程、体积、和限制条件进展设置。方程包括有CUDC-IN、CXDC-IND、DT、RHOB、NPHIGR；体积有石英〔QUAR、伊利〔ILL〔XGAS〔XW〔UGAS〔UWA；限制条件设置了水基泥浆的限制条件〔图5-1。、初始化参数数，场地区设置参数分别为：RMF〔电阻率〕=0.25Ωm，MST〔温度〕=55。储层平均孔隙度=0.12。〔包括泥浆滤液电阻率等、〔粘土束缚水电导率等〔〔图5-1。这些参数在ELAN计算中是必不行少的。图图5-13场沙溪庙组参数设置图5-14 初始化计算参数响应参数冲洗带气、水参数。这些参数多数运用程序隐含值，少局部进展了修改，如长石、GR响应参数进展了修改，使计算出的岩性剖面更加合理〔5-1不确定性参数ELAN5-16为各曲线的不确定参数，均用程序的隐含值。可以依据曲线特征进展修改。图5-15 成分响应参数值图5-16 不确定性参数权系数参数权系数关系到ELAN方程中各测井曲线的权值大小，程序中设有隐含值，可作适当人工调整，对于测井质量较好的曲线可调高权系数，曲线质量较差时，可以降低权系数。此外，依据处理目的不同，可以调整权系数。图5-17为场沙溪庙组测井曲线权系数，其中，GR权系数的隐含值为0.3，利用隐含值处理的剖面不太合理，将权系数调高到1时，处理的岩性剖面比较符合客观实际。5-17 权系数参数2、EALNPlus计算与成果显示5-18显示出计算结果。图左侧显示出各条曲线与计算曲线，实现为原始曲线，虚线为计算曲线砂岩储层段原始曲线与计算〔反算曲线〕近于重合，说明计算结果牢靠；砂岩上下的5-18为ELAN计算的剖面结果，包括了各种成分的相对体积含量；图右侧为计算牢靠性的评价指标，是由各条原始曲线与反算曲线间的差异程度确定的。图5-18 ELANPlus初步处理结果3、储层参数计算ELAN隙度、渗透率、含水饱和度参数需要在ELAN处理的根底上进展计算。其中：孔隙度为原状地层水与气体积的和，或者冲洗带气与水体积的和；含水饱和度为水的体积除以总孔隙度，这个过程由FUCTION模块计算完成。图5-19显示出了储层参数特征，其中孔隙度在第三道的剖面中已经得到显示，第四道显示出孔隙度和含水饱和度值。图5-19储层参数特征ELAN计算参数结果统计〔RESUM〕ELAN计算参数结果较多，可以统计的参数很多，这个过程在RESUM模块中进展。统计的参数包括砂岩、储层、有效储层的厚度、孔隙度、渗透率、含水饱和对需要统计的参数进展选择。对于场沙溪庙组地层，主要统计储层及有效储层的厚度及孔、渗、饱参数。首先，需要对下限值进展设定。1、储层参数下限值常规测井解释的有效储层参数统计中，运用的储层参数下限值有孔隙度、渗GEOFRAMEGR值、泥质含量、孔隙度、渗透率、含水饱和度、泥岩薄夹层、储层薄夹层。依据过去的争论成果和地区阅历，场地区沙溪庙组储层参数〔包括砂岩层、储层、有效储层〕的下限值分别为GR＜70，Vsh＜0.3，PERM＜0.07，POR＜7%，Sw＜70%，泥岩薄夹层剔除＜0.1m0.1〔图5-2砂岩中夹有泥岩薄层，当薄层厚度较小予以剔除，统计中以砂岩计算；砂岩薄夹层是指在泥岩中有较薄的砂岩层，当砂岩厚度较薄时予以剔除。图5-20 储层下限值设定2、参数统计结果参数统计结果见表，其中有效储层是运用了全部储层参数下限值限制条件后和度下限值限制，砂岩层没用孔隙度、渗透率和含水饱和度限制条件。从统计结果看，有效储层孔隙度一般10%-15%，渗透率0.5-m2103和度30%-60%。表5-3 沙溪庙组储层参数统计表砂岩层 储层 有效储层井号 顶深

砂岩平均度度度CX1642245.6度度度CX1642245.62270.023.70.1121.90.121.390.4917.50.122.490.44CX1652211.52230.418.40.1418.40.140.90.3513.60.152.080.34CX1672251.32277.021.80.0817.70.10.410.687.30.113.520.56X8012161.92179.611.70.1111.60.120.050.220.0000X8042138.22153.815.50.115.50.12.870.5912.60.14.410.55X8052130.82151.519.10.1318.90.131.770.5115.40.143.510.49

孔隙渗透度 率

含水 饱和厚度 度

渗透 含水率 饱和X8062163.12187.323.60.1322.60.131.930.3818.60.143.340.34X8072170.42194.223.20.1322.90.135.580.3121.60.146.820.31X8082215.92239.316.50.1116.20.130.720.388.10.132.340.35X8512191.72215.523.30.1323.20.130.360.379.80.141.410.36CX1642324.52344.318.30.0913.60.110.570.516.10.123.050.43CX1652280.32306.419.50.1319.50.141.070.3613.10.152.730.35CX1672329.02356.724.00.120.20.121.580.4415.70.123.180.42X8072248.52274.624.20.1120.70.130.360.258.40.141.290.25X8082294.42318.89.00.098.50.120.180.371.40.130.650.37X8162265.42295.524.50.1323.40.140.370.3211.40.151.320.31X8512264.82294.522.40.1322.40.140.330.3511.90.161.330.32X8042238.32257.316.30.0811.30.091.330.753.70.113.20.61X8052226.42245.815.20.0813.00.090.070.80.50.130.820.67X8062274.52285.29.30.060.00000.0000X8072279.82294.47.80.064.10.090.030.510.0000CX1642380.32398.818.40.1218.40.120.990.4318.40.120.990.43CX1652337.42352.014.60.1214.60.120.240.414.60.120.240.4CX1672392.92403.911.00.0711.00.070.040.8911.00.070.040.89X8042262.42281.218.80.118.80.11.290.7318.80.11.290.73X8052250.32268.718.40.084.80.080.010.844.80.080.010.84X8062295.42314.318.80.0818.80.080.210.6418.80.080.210.64X8072300.32314.514.20.1114.20.110.560.4714.20.110.560.47X8082349.12364.315.20.115.20.10.110.5215.20.10.110.52X8162325.82338.412.70.112.70.10.70.5812.70.10.70.58X8512319.02337.818.80.118.80.10.110.8118.80.10.110.81CX1642436.92454.916.30.1114.80.123.050.4711.00.127.110.44CX1652399.52414.413.40.1413.40.142.020.3910.80.154.390.39CX1672451.92459.97.90.117.80.114.160.567.20.115.110.55X8042322.02343.920.60.1420.30.1410.30.4319.70.1411.660.42X8052315.02329.310.1X8062355.62371.215.20.1113.00.137.280.4810.50.1312.850.44X8072356.52378.116.90.0710.00.11.080.477.30.112.110.44X8082410.92428.014.80.1114.20.121.210.479.60.132.360.42X8162383.62399.915.50.1315.50.145.550.4313.90.1490.41X8512385.52399.313.70.1313.10.134.220.4911.90.145.870.47ELAN处理与常规处理结果比照常规处理中是指用经典的测井方法进展测井处理，对于场沙溪庙组地层，一般用GR计算泥质含量，用AC、