一种基于穿层孔测井的煤矿透明工作面构建方法与流程

本发明属于煤田地球物理测井技术领域，具体涉及一种基于煤矿井下穿层孔测井技术的透明工作面构建方法。

背景技术：

随着煤矿智能化开采技术的发展，要求构建煤矿井下透明工作面，为智能开采提供精准的地质模型。目前透明工作面的探测方法有地面三维地震、工作面电磁波透视、工作面槽波层析成像以及井下钻探等，测井探测精度高但是探测深度浅，在煤矿透明工作面构建中尚未受到重视。

技术实现要素：

针对现有技术中的缺陷和不足，本发明提供了一种基于穿层孔测井的煤矿透明工作面构建方法，通过该方法可以利用煤矿井下穿层孔的测井信息构建透明工作面，为煤矿智能开采提供高精度的地质模型。

为达到上述目的，本发明采取如下的技术方案：

一种基于穿层孔测井的煤矿透明工作面构建方法，该方法包括以下步骤：

步骤1，收集前期勘探的煤层分布和底抽巷的资料，然后设计煤矿井下底抽巷的钻孔分布，根据设计的钻孔位置打钻，现场测井；

步骤2，对获得的测井数据进行处理与分析得到分层曲线进而得到岩性曲线，将岩性曲线与钻孔轨迹联合得到目标层顶底界面数据集；

步骤3，对前期勘探得到的目标层内的构造分布区域的钻孔的测井曲线进行构造识别，在分层的基础上，提取各钻孔中测井曲线上构造的位置，将该区域内的构造进行联合，形成构造的数据集；

步骤4，根据测井数据，计算各钻孔中目标层段的孔隙度、渗透率、含水量和含气量，结合钻孔的轨迹数据得到孔隙度、渗透率、含水量和含气量的位置数据集；

步骤5，根据所述目标层顶底界面数据集、构造的数据集以及孔隙度、渗透率、含水量和含气量的位置数据集，绘制三维空间的分布。

本发明还包括如下技术特征：

具体的，所述步骤1包括：

步骤1.1，收集目标区域以及底抽巷的资料，包括目标区域的：前期勘探的煤层的分布，煤层的厚度hq、宽度wq、长度lq、陷落柱和断层的位置坐标集sg(x,y,z)，设计的底抽巷分布位置数据集sh(x,y,z)，底抽巷巷道的数量m，各底抽巷巷道的长度lhi、宽度whi、高度hhi、顶部距煤层的距离hji以及方位角azimhi；

步骤1.2，设计各个巷道内的钻孔数量、开孔点和钻孔的轨迹，测井系统的探测半径为d，一个钻孔能探测的面积为d²，根据建立目标层透明工作面的长度lq和宽度wq，目标层的面积为sq，一个巷道对应探测的目标层的长度为lqh，宽度为wqh，需要的钻孔的数量总数psum最大值：

一个钻孔群内的钻孔数量p的最大值为pmax：

钻孔群数量n的最大值nmax为：

巷道数量m的最大值mmax为：

相邻两个钻孔之间的距离dpi满足2d≤dpi≤3d，钻孔群之间的距离dzkq：

上式中，构造区域为所述步骤1.1中陷落柱和断层区域；变化相对平缓区域为没有陷落柱和断层区域；

钻孔群的数量n：

一个钻孔群内的钻孔数量p：

巷道的顶部弧形部分的对应的扇形弧度为θ，在钻孔群中第一个钻孔的钻孔倾角devip1为：

第一个钻孔的深度dp1：

上式中，hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；θ为巷道的顶部弧形部分的对应的扇形弧度角；devip1为钻孔群中第一个钻孔的钻孔倾角；

第i个钻孔的倾角devipi为：

上式中，devipi为第i个钻孔的倾角，hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；dp1为第一个钻孔的深度；devip1为钻孔群中第一个钻孔的钻孔倾角；dpi为钻孔群内相邻钻孔在穿出目标层位置的距离；θ为巷道的顶部弧形部分的对应的扇形弧度角；

第i个钻孔的深度dpi为：

上式中，dpi为第i个钻孔的深度；hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；devipi为第i个钻孔的倾角，θ为巷道的顶部弧形部分的对应的扇形弧度角；

方位角azimpi以与巷道的走向的方位角垂直作为其方位，因此azimpi为：

上式中，azimpi为第i个钻孔的方位角，azimhj是第j个巷道的方位角；

步骤1.3，现场根据设计的钻孔位置、钻孔的轨迹进行打钻，一个一个钻孔群打钻；

步骤1.4，在完成一个钻孔群的打钻之后，进行现场测井，获得测井数据。

具体的，所述步骤2包括：

步骤2.1，测井数据预处理，在完成打钻以及获得测井数据后，首先进行测井数据预处理，预处理包括对不同的测井仪器进行深度校正、无效数据剔除和数据的平滑处理；

步骤2.2，提取岩性分层特征，测井曲线经过预处理后，分析目标层与其他岩层的差异性特征，根据自然伽马测井、电阻率测井、孔隙度测井在不同岩性中表现的特征集f(gr，rt，den)；然后进行测井数据的处理，处理测井曲线，建立分层曲线；

步骤2.3，进行钻孔的分层并形成顶底界面的数据集，根据步骤2.2的分层曲线，结合测井特征值数据集f(gr，rt，den)对各个层进行测井数据的对比，并且形成岩性曲线lith，然后利用测斜测量的倾角和方位角计算钻孔三维空间轨迹trj(n，e，d)：

上式中，d、n和e分别为钻孔中第c个测点的海拔、纬度和经度；

其中：

上式中，δd、δn和δe分别为钻孔中第c个测点和c-1个测点之间的垂深增量(即海拔增量)海拔差、坐标北方向增量(纬度增量)和坐标东方向增量(经度增量)；rc是校正系数；

上式中，rc校正系数；dlac为钻孔中第c个测点的狗腿角，单位是弧度；mdc是第c个测点的孔深；

dlac＝arccos[cos(devic-1)cos(devic)+sin(devic-1)sin(devic)cos(azimc-azimc-1)]

(20)

上式中，dlac为钻孔中第c个测点的狗腿角，devic为该钻孔第c个测点的测量的倾角，azimc为该钻孔第c个测点的测量的方位角；

将岩性曲线与三维轨迹联合，提取所有钻孔的目标层的界面分层点的坐标形成底界面的数据集slb(n，e，d)和顶界面数据集slt(n，e，d)。

具体的，所述步骤3包括：

构造识别并构建构造的数据集，对于前期勘探收集得到的目标层内的构造分布区域，在这些区域的钻孔的测井曲线上进行构造识别，在分层的基础上，分析小层内的测井曲线的变化，然后进行人机交互识别构造，并提取各个钻孔中测井曲线上构造的位置，将该区域内的构造进行联合，形成构造的数据集sg(n,e,d)。

具体的，所述步骤4包括：

孔隙度、渗透率、含水量、含气量的计算并构建数据集，在前面分层、构造的建立基础上，计算目标层的孔隙度、渗透率、含水量、含气量的曲线，孔隙度ф的计算公式为：

式中：ф为孔隙度，小数；ρma＝2.65kg/m³；ρf＝1kg/m³；ρb为密度测井仪测得的密度值，kg/m³,；

渗透率perm的计算公式为：

式中perm为渗透率，md；系数a1和b1根据现场的实际情况来确定，ф为孔隙度，小数；

含水饱和度sw的计算公式为：

式中sw为含水饱和度，小数；系数a2、b2、m’、n’根据现场的实际情况来确定，ф为孔隙度，小数；rw为该地区的地层水的电阻率，rt为电阻率测井测得的视电阻率；

含水量qw的计算公式为：

式中，qw为煤层的含水量，单位m³/t；ρ为目标层煤层的密度，单位为g/m³；sw为含水饱和度，小数；

含气量qg的计算包括游离气q1和吸附气q2的两个部分，游离气的含量q1计算公式为：

式中，q1为游离气含量，单位m³/t,ρ为目标层煤层的密度，单位为g/m³；sw为含水饱和度，小数.

吸附气的含量q2计算公式为：

式中：q2—吸附气含量,m³/t；vl一干煤无灰兰氏含气量,m³/t；va一含灰量；vw—含水量；b—兰氏常数,1/pa；p—压力,pa/m²；

目标层煤层的含气量qg为：

qg＝q1+q2(27)

将各个孔计算的孔隙度、渗透率、含水量、含气量与该孔的三维轨迹坐标联系起来，形成孔隙度sф、渗透率sperm、含水量sqw、含气量sqg的数据集。

具体的，所述步骤5包括：

三维成图，构建透明工作面；将上述获得的各个数据集输入三维绘图的软件中，先加载各个数据集的坐标位置，然后利用克里金插值，实现钻孔区域间空白部分的充填。

具体的，所述步骤2.1中，测井数据预处理包括：

深度校正：各测井仪器串联进入孔内进行测量，第一个测井仪器的深度测量点与测井系统的深度基准点重合，第二个测井仪器深度点与测井系统的深度基准点距离d2，第三个测井仪器深度点与测井系统的深度基准点距离d3，第四个测井仪器深度点与测井系统的深度基准点距离d4；测井系统的深度基准点的深度为deppi，各个测井仪器的测井深度分别记为dep1pi、dep2pi、dep3pi、dep4pi，因此各个测井仪器的深度校正如下：

无效数据剔除：在将各个测井仪器测量数据的深度进行归位之后，进行各个测井仪器的数据的检查，根据各个测井仪器的无效数据标记，剔除无效数据；

数据的平滑处理：放射性测井原始数据由于统计的涨落会有毛刺，测井仪器还有自身的记录数据的波动性会产生毛刺，因此需要进行平滑处理，消除这些影响；采用滑动窗的平滑处理办法，假设对窗长(2t+1)个相邻点进行平滑，测井的数据点为tk，相邻的点为tk-t、tk-t+1…tk+t,平滑公式为：

具体的，所述步骤2.2中，分层曲线layerlog的建立方法为：各个测井曲线共有n个采样点，x1a、x2a分别为这两层第a点的测井值，若分层界面在b，b+1采样点间，则两层的层内差方和s为：

式中，

然后再求取s与q的极值，s'＝0，q'(b)＝0时的a值为分层点，形成分层曲线layerlog，并且对于各个层的测井值求其层内的平均值作为该层段的测井值。

本发明与现有技术相比，有益的技术效果是：

本发明测井探测精度高，在底抽巷或顶抽巷内，布局多个钻孔群，在高瓦斯矿井抽采瓦斯的同时，可利用测井精细探测煤层界面、煤层内部的构造(陷落柱、断层等)，亦可探测煤层的孔隙度、渗透率、瓦斯含量以及水的含量等，构建透明工作面，为智能开采提供精细探测的数据。

附图说明

图1是本发明基于穿层孔测井的透明工作面构建的流程图；

图2是本发明煤矿井下底抽巷道与目标层及设计钻孔的切面示意图；

图3是本发明其中一个巷道方向的目标层与穿层孔的关系俯视图；

图4是本发明目标层与其中一个巷道及钻孔群的侧面示意图；

图5是本发明各测井仪器与深度测量标准点的关系示意图；

图6是本发明基于穿层孔构建的透明工作面的示意图；

附图标号含义：

图2中，2-1-目标层(煤层)；2.1-该钻孔群内第一个钻孔，2.2-第二个钻孔，2.3-第三个钻孔，2.i-第i个钻孔，2.p-2-第p-2个钻孔，2.p-1-第p-1个钻孔，2.p-第p个钻孔，p-该钻孔群的钻孔总数；2-3-钻孔群内相邻钻孔在穿出目标层位置的距离dpi；2-4-目标层的厚度hq；2-5-各个巷道顶部距离目标层的距离hj(j为第j个巷道)；2-6-巷道的高度hhj(h巷道，j是第j个巷道)；2-7-巷道的宽度whj；2-8-巷道顶端弧形部分的弧度角θj；2-9-巷道；2-10-钻孔的倾角αi；2-11-钻孔的深度dpi。

图3中，3-1-目标层(煤层)；3-2-目标层长度lq；3-3-目标层的宽度wq；3-4-该巷道内的第1个钻孔群；3-5-该巷道内第2个钻孔群；3-6-该巷道内的第i’个钻孔群；3-7-该巷道内的第n-1个钻孔群；3-8-相邻钻孔群之间的距离dzkq；3-9-该巷道内的第n个钻孔群；3-10-钻孔群内相邻钻孔在穿出目标层位置的距离dpi；3-11-钻孔群内钻孔示意。

图4中，4-1-目标层(煤层)；4-2-煤层厚度hq；4-3-巷道顶部与目标层之间的距离hj；4-4-巷道高度hhj；4-5-巷道(侧面)；4-6-巷道长度lhj；4-7-钻孔群示意。

图5中，5-1-深度测量基准点；5-2-位于测井系统前端的第一个测井仪(测斜仪)；5-3-位于测井系统的第二个测井仪(自然伽马测井仪)；5-4-位于测井系统的第三个测井仪(电阻率测井仪)；5-5-位于测井系统的第四个测井仪(孔隙度测井仪)；5-6-第四个测井仪与深度测量基准点的距离d4；5-7-第三个测井仪与深度测量基准点的距离d3；5-8-第二个测井仪与深度测量基准点的距离d2。

图6中，6-1-目标层的顶界面示意图，6-2-目标层底界面示意图，6-3-钻孔示意图，6-4-巷道示意图，6-5-绘图的框架示意；6-6-坐标角度示意。

以下结合说明书附图和具体实施方式对本发明做具体说明。

具体实施方式

本发明所设计的基于穿层孔测井的透明工作面构建方法，包括测井钻孔布置设计、基于单孔的煤层顶底板测井解释方法、基于多孔测井的断层解释方法以及基于多孔测井的透明工作面构建方法。

测井钻孔布置是指煤矿井下巷道内钻孔布置设计，即煤矿井下巷道内钻孔布置的设计。

煤矿井下巷道内钻孔布置设计是指煤层底板中底抽巷内的钻孔设计，底抽巷内钻孔主要是高瓦斯矿井在煤层开采前，沿煤层底板施工岩巷并在巷道内打孔，将其上部煤层中的瓦斯安全抽出。本发明利用底抽巷内钻孔进行测井，获取地层信息，构建透明工作面，底抽巷内钻孔设计一是根据前期瓦斯分布的资料，结合巷道的位置设计钻孔群的位置，二是根据测井构建透明工作面的要求设计钻孔群内钻孔的开孔位置及轨迹。

单孔测井解释方法为井下钻孔测井之后进行的处理解释方法，处理解释的主要目标是在沿钻孔轨迹方向区分钻孔钻遇地层的岩性，将岩性进行层位归一，为透明工作面的构建做准备。

基于测井的断层解释方法为根据相邻钻孔测井数据的差异，判断断层的位置，根据多个钻孔岩性的变化，确定断层的大小及分布特征。

基于测井的透明工作面的构建方法为根据穿层孔的测井解释结果、断层的分布位置等，将各个钻孔解释的结果通过特定的方法，构建一个透明的三维地质体。

一种基于测井的透明工作面构建方法，包括如下的步骤：

设计煤矿井下底抽巷的钻孔分布，根据前期收集的资料，针对底抽巷巷道的分布、底抽巷与目标区域的相对位置关系，设计钻孔的分布。预计构建透明工作面区域的长lq，宽wq，高hq，底抽巷距离目标的煤层距离为hj，巷道的长度为lh，高度为hh，巷道宽度为wh，在底抽巷内有m个巷道，在各个巷道内设计n个钻孔群，各个钻孔群内设计p个钻孔，所采用的测井系统的探测半径为d，钻孔群的距离dzkq根据前期勘探资料，在有构造的区域进行密集的设计，在相对较为平缓的区域，进行稀疏设计，dzkq满足如下的公式：

上式中，构造区域为所述步骤1.1中陷落柱和断层区域；变化相对平缓区域为没有陷落柱和断层区域；

式中，d为测井系统的探测半径，m。

各个钻孔群内设计p个钻孔，钻孔的开孔位置为巷道顶部的弧形区域，如果在巷道的侧帮位置打穿层孔，孔深较大长不利于打钻，p个钻孔相互之间的距离dpi在目标层位的穿层区域点的分布要满足2d≤dpi≤3d，根据各个钻孔与顶部的目标层位的关系，设计各个钻孔群里各钻孔的轨迹(包括钻孔的倾角devii、钻孔的方位角azimi、钻孔的长度dzi)。

对设计的钻孔在井下巷道进行标记，然后根据设计的钻孔位置以及钻孔的位置进行打钻，在每一个钻孔群打钻完成后，立即进行测井，测量钻孔的实际钻孔轨迹(钻孔的实际倾角、实际的方位角)，采集各个测井方法的探测数据log1、log2……logs。

根据获得的测井数据，进行测井数据的处理与分析。首先进行数据的预处理，包括深度校正(各个测井仪器的深度差校正)、测井数据的无效数据的剔除、测井数据的平滑处理。

然后进行数据的分析，提取不同岩性的测井参数，分析不同岩性的分层特征，然后根据测井参数，分析各个钻孔的数据中不同岩性的分界点，分析构造(断层、陷落柱)的位置。结合钻孔的轨迹数据，进行各个特征的定位，将各个钻孔的穿层位置坐标提取出来lb(n,e,d)和lt(n,e,d)，形成目标层底界面的数据集和顶界面的数据集，以及根据测井系列提取出来区域内的构造的位置数据集(包括断层、陷落柱)。

根据测井数据，计算各个钻孔中目标层段的孔隙度、渗透率、含水量、含气量的大小，结合钻孔的轨迹数据进行位置的定位，形成目标区域内的孔隙度、渗透率、含水量和含气量的分布，形成孔隙度、渗透率、含水量和含气量的位置数据集。

根据形成的数据集，绘制三维空间的分布。

需要注意的是：煤矿井下智能化开采透明工作面的构建，其首要任务是明确目的层(煤层)高精度的顶底板界面、精确的构造位置，其次是煤层的含水和含气量的分布，以及煤层孔隙度、渗透率的大小等，因此测井仪器系统要求必须能测量钻孔的轨迹、自然伽马、其次是电阻率，有条件可以进行孔隙度测井(密度测井仪或中子测井仪)。

遵从上述技术方案，以下给出本发明的具体实施例，需要说明的是本发明并不局限于以下具体实施例，凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。下面结合实施例对本发明做进一步详细说明。

实施例1：

如图1至图6所示，本实施例给出、一种基于穿层孔测井的煤矿透明工作面构建方法，其特征在于，包括以下步骤：

步骤1，收集目标区域以及底抽巷的资料，包括目标区域(煤层)的前期勘探的煤层的分布，煤层的厚度hq、煤层的宽度wq、煤层的长度lq、煤层中的陷落柱和断层的位置坐标集sg(x,y,z)，设计的底抽巷分布位置数据集sh(x,y,z)，底抽巷巷道的数量m、各个底抽巷巷道的长度lhi、各个底抽巷巷道的宽度whi、各个底抽巷巷道的高度hhi、各个底抽巷巷道的顶部距煤层的距离hji、各个底抽巷巷道的方位角azimhi；如图2、图3、图4所示的关于目标层和巷道标注信息；

步骤2，设计各个巷道内的钻孔数量和钻孔的开孔点、钻孔的轨迹；在底抽巷设计钻孔，要求钻孔能将要求探测构建透明工作面的目标区域能通过孔中物探测井方法实现精细探测，因此底抽巷的巷道数量m、钻孔群数量n、一个钻孔群中钻孔的数量p和钻孔的分布与目标层的大小以及底抽巷的分布位置紧密联系在一起；测井系统的探测半径为d，一个钻孔可以探测的面积为d²，根据建立目标层透明工作面的长度lq和宽度wq，目标层的面积为sq，一个巷道对应探测的目标层的长度为lqh，宽度为wqh，需要的钻孔的数量总数psum最大值：

一个钻孔群内的钻孔数量p的最大值为pmax：

钻孔群数量n的最大值nmax为：

巷道数量m的最大值mmax为：

根据现场的实际情况，考虑打钻经济投入的问题，对于钻孔群数量、钻孔数量的设计需要根据实际的地质情况考虑，在地质变化相对平缓的区域钻孔的距离可以相对放宽，如图3所示钻孔与目标区域的关系示意图中的钻孔及钻孔群之间的距离、钻孔群内钻孔之间的距离，

钻孔相互之间的距离dpi满足2d≤dpi≤3d即可，钻孔群的距离也可以相对放宽，钻孔群之间的距离dzkq：

上式中，构造区域为所述步骤1.1中陷落柱和断层区域；变化相对平缓区域为没有陷落柱和断层区域；

钻孔群的数量n：

一个钻孔群内的钻孔数量p：

各个钻孔的轨迹包括钻孔的深度、钻孔的倾角、方位角与钻孔所在巷道与目标区域的相对位置关系紧密相关，如图2所示的示意图，巷道的顶部弧形部分的对应的扇形弧度为θ，在钻孔群中第一个钻孔的钻孔倾角devip1为：

第一个钻孔的深度dp1：

上式中，hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；θ为巷道的顶部弧形部分的对应的扇形弧度角；devip1为钻孔群中第一个钻孔的钻孔倾角；

第i个钻孔的倾角devipi为：

上式中，devipi为第i个钻孔的倾角，hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；dp1为第一个钻孔的深度；devip1为钻孔群中第一个钻孔的钻孔倾角；dpi为钻孔群内相邻钻孔在穿出目标层位置的距离；θ为巷道的顶部弧形部分的对应的扇形弧度角；

第i个钻孔的深度dpi为：

上式中，dpi为第i个钻孔的深度；hq为目标层的厚度；hj为第j个巷道的巷道顶部距离目标层的距离；hhj为巷道的高度；whj为巷道的宽度；devipi为第i个钻孔的倾角，θ为巷道的顶部弧形部分的对应的扇形弧度角；

各个钻孔的方位角azimpi以与巷道的走向的方位角垂直作为其方位，因此azimpi为：

上式中，azimpi为第i个钻孔的方位角，azimhj是第j个巷道的方位角；

步骤3，现场根据设计的钻孔位置、钻孔的轨迹进行打钻，一个一个钻孔群打钻；

步骤4，在完成一个钻孔群的打钻之后，进行现场测井，获得测井数据；

步骤5，测井数据预处理，在完成打钻以及获得测井数据后，首先进行测井数据预处理，预处理包括对不同的测井仪器进行深度校正、无效数据的剔除、数据的平滑处理；

深度校正：如图5所示，各个测井仪器串联进入孔内进行测量，与测井系统记录的深度基准点有一定的距离差，第一个测井仪器的深度测量点与测井系统的深度测量点重合，第二个测井仪器深度点与测井系统的深度测量点距离d2，第三个测井仪器深度点与测井系统的深度测量点距离d3，第四个测井仪器深度点与测井系统的深度测量点距离d4，测井系统的深度测量记录的深度为deppi，各个测井仪器的测井深度分别记为dep1pi、dep2pi、dep3pi、dep4pi，因此各个测井仪器的深度校正如下：

无效数据剔除：在测井过程中，可能会有一些无效数据需要剔除，在将各个测井仪器测量数据的深度进行归位之后，进行各个测井仪器的数据的检查，根据各个测井仪器的无效数据标记，剔除无效数据；

测井数据的平滑处理，放射性测井(自然伽马测井、孔隙度测井(密度测井))原始数据由于统计的涨落会有毛刺，测井仪器还有自身的记录数据的波动性会产生毛刺，因此需要进行平滑处理，消除这些影响；采用滑动窗的平滑处理办法，假设对窗长(2t+1)个相邻点进行平滑，测井的数据点为tk，相邻的点为tk-t、tk-t+1…tk+t,平滑公式为：

步骤6，提取岩性分层特征，在测井曲线经过预处理之后，分析目标层(煤层)与其他岩层的差异性特征，根据自然伽马测井、电阻率测井、孔隙度测井在不同岩性中表现的特征集f(gr，rt，den)，并且将其加入地区数据库中，以便以后建立其他的透明工作面参考使用；然后进行测井数据的处理，处理测井曲线，建立分层曲线；分层曲线layerlog的建立方法为：各个测井曲线共有n个采样点，x1a、x2a分别为这两层第a点的测井值，若分层界面在b，b+1采样点间，则两层的层内差方和s为：

式中，

然后再求取s与q的极值，s'＝0，q'(b)＝0时的a值为分层点，形成分层曲线layerlog，并且对于各个层的测井值求其层内的平均值作为该层段的测井值；

步骤7，进行钻孔的分层并形成顶底界面的数据集，根据s6的分层曲线，结合s6的测井特征值数据集f(gr，rt，den)对各个层进行测井数据的对比，并且形成岩性的曲线lith，然后利用测斜测量的钻孔的倾角和方位计算钻孔三维空间轨迹trj(n,e,d)：

上式中，d、n和e分别为钻孔中第c个测点的海拔、纬度和经度；

其中：

上式中，δd、δn和δe分别为钻孔中第c个测点和c-1个测点之间的垂深增量(即海拔增量)海拔差、坐标北方向增量(纬度增量)和坐标东方向增量(经度增量)；rc是校正系数；

上式中，rc校正系数；dlac为钻孔中第c个测点的狗腿角，单位是弧度；mdc是第c个测点的孔深；

dlac＝arccos[cos(devic-1)cos(devic)+sin(devic-1)sin(devic)cos(azimc-azimc-1)]

(20)

上式中，dlac为钻孔中第c个测点的狗腿角，devic为该钻孔第c个测点的测量的倾角，azimc为该钻孔第c个测点的测量的方位角；

将岩性曲线与三维轨迹联合，提取所有钻孔的目标层的界面分层点的坐标形成底界面的数据集slb(n，e，d)和顶界面分层点的坐标集slt(n，e，d)；

步骤8，构造识别并构建构造的数据集，对于前期勘探收集的到的目标层内的构造分布区域，进行详细的解释，在这些区域的钻孔的测井曲线上进行构造识别，在分层的基础上，分析小层内的测井曲线的变化，然后进行人机交互识别构造，并提取各个钻孔中测井曲线上构造的位置，将该区域内的构造进行联合，形成构造的数据集sg(n，e，d)；

步骤9，孔隙度、渗透率、含水量、含气量的计算并构建数据集，在前面分层、构造的建立基础上，计算目标层的孔隙度、渗透率、含水量、含气量的曲线，孔隙度ф的计算公式为：

式中：ф为孔隙度，小数；ρma＝2.65kg/m³；ρf＝1kg/m³；ρb为密度测井仪测得的密度值，kg/m³，；

渗透率perm的计算公式为：

式中perm为渗透率，md；系数a1和b1根据现场的实际情况来确定，ф为孔隙度，小数；

含水饱和度sw的计算公式为：

式中sw为含水饱和度，小数；系数a2、b2、m’、n’根据现场的实际情况来确定，ф为孔隙度，小数；rw为该地区的地层水的电阻率，rt为电阻率测井测得的视电阻率；

含水量qw的计算公式为：

式中，qw为煤层的含水量，单位m³/t；ρ为目标层煤层的密度，单位为g/m³；sw为含水饱和度，小数；

含气量qg的计算包括游离气q1和吸附气q2的两个部分，游离气的含量q1计算公式为：

式中，q1为游离气含量，单位m³/t,ρ为目标层煤层的密度，单位为g/m³；sw为含水饱和度，小数.

吸附气的含量q2计算公式为：

式中：q2—吸附气含量,m³/t；vl一干煤无灰兰氏含气量,m³/t；va一含灰量；vw—含水量；b—兰氏常数,1/pa；p—压力,pa/m²；

目标层煤层的含气量qg为：

qg＝q1+q2(27)

将各个孔计算的孔隙度、渗透率、含水量、含气量与该孔的三维轨迹坐标联系起来，形成孔隙度sф、渗透率sperm、含水量sqw、含气量sqg的数据集；

步骤10，三维成图，构建透明工作面，如图6所示；将上述获得的各个数据集输入三维绘图的软件中，先加载各个数据集的坐标位置，然后利用克里金插值，实现钻孔区域间空白部分的充填。

需要注意的是：1)对于煤矿井下的开采煤层的透明工作面的构建首要任务是明确目的层(煤层)的精细的顶底界面、精确的构造位置，其次是煤层的含水和含气量的分布，以及煤层孔隙度、渗透率的大小，因此测井仪器系统要求必须能测量钻孔的轨迹、自然伽马测井、其次是电阻率测井，有条件可以进行孔隙度测井(密度测井仪或中子测井仪)。2)根据现在的钻孔间距设置，钻孔之间尚存在一定的盲区，对于精度要求更高的地方，钻孔间距、钻孔群的间距可以适当缩小，以实现更精细的探测。