欧式距离法在电阻率测深法找水有效性研究.pdf

2023年第8期西部探矿工程*收稿日期：2022-09-19修回日期：2022-09-29第一作者简介：王超（1983-），男（汉族），河南长垣人，工程师，现从事地球物理勘探工作。欧式距离法在电阻率测深法找水有效性研究王超*1,2，李科强1,2，高晶3（1.河南省地质矿产勘查开发局测绘地理信息院，河南 郑州 450006；2.河南省地质勘查信息化工程技术研究中心，河南 郑州 450006；3.河南农业大学，河南 郑州 450002）摘要：在电阻率测深法找水资料处理解释时，通过多个测点电阻率曲线类型图与已知井位相对比，只能定性分析相似程度，在多条电阻率曲线类型和已知井位相似时，不好判断最佳井位。为了能从定量分析的角度解决这一问题，模糊数学中的欧式距离公式提供一种数学方法，以已知井位测量数据为标准样本，利用相对欧式距离法公式计算工作区测点测量数据与已知井位的距离值，距离值越小，测点越可能找到水。通过实测数据，利用欧式距离法进行计算处理，为打井优劣位置的选取提供了物探依据。通过井位抽水试验，达到了预定指标要求。结果表明，欧式距离法在电阻率测深法找水资料处理解释的有效性，证明欧式距离法在电阻率测深法找水确定最佳井位是行之有效的。关键词：模糊数学；欧式距离法；电阻率测深法；找水；有效性中图分类号：P641 文献标识码：A 文章编号：1004-5716(2023)08-0151-03地球物理勘探方法在找水有其独特优势，电阻率测深法因研究地下介质的垂向变化有其独特优势，应用最广泛。电阻率测深法测量数据成果资料处理解释对于找水井位确定至关重要。对于电阻率测深法找水资料处理解释，物探工作者做了大量研究工作1-5。常规方法是结合水文地质情况，用测点电阻率曲线类型图去定性判断含水层的位置，成果解释从一维发展到二维6。对于多测点测量，通过和已知井位电阻率曲线类型图进行对比，进行定性分析判断解释，当多条电阻率曲线类型相类似时，确定最佳井位难度较大。许多物探工作者为了能对测量数据从定量分析角度处理解释，确保确定井位的出水量能够到达预定指标，做了大量研究工作。模糊数学7中贴近度8-9、灰色关联度分析法10、模糊相似优先比11-12、欧式距离法13、电反射系数（K）法14、梯度K法15、衍生指数16等数学方法都应用到电测深数据处理解释。通过对比分析，欧式距离法具有数据运算简单、实用，且与实际情况符合较好的优点。本文通过运用欧式距离法对野外测点实测数据进行处理解释，验证其应用效果。1方法原理模糊数学是用数学方法研究和处理具有“模糊性”现象的数学，其模糊性是指客观事物差异间的过渡界限“不分明性”。有模糊熵、距离测度和相似性测度对模糊集进行计算处理，计算结果可以比较两个模糊集的近似程度。在实际数据计算中，模糊熵应用不太方便和贴近度计算过程相对复杂，应用不太方便，利用闵氏距离公式的特殊形式欧氏距离法，计算相对方便、简洁13。在泛函分析17中，对于距离空间的定义，设X是非空集合，对于X中任意两个元素x和y，按某一法则都对应唯一的实数d(x，y)，且满足：（1）非负性：d(x，y)0；d(x，y)=0，当且仅当x=y；（2）对称性：d(x，y)=d(x，y)；（3）三角不等性：对于任意的x，y，zX，d(x，y)d(x，z)+d(y，z)，则称d(x，y)为x与y的距离。设A、B为论域U中的两个模糊子集，当论域U为有限集时，可用闵可夫斯基公式计算A、B的距离，闵可夫斯基距离公式又称为闵氏距离公式，闵氏公式是一组距离的定义。模糊子集A、B的距离公式为：d(A,B)=i=1n|UA(xi)-UB(xi)q1q（1）式中：n论域U中的个数；1512023年第8期西部探矿工程q指数系数，且为正整数。当q=2时，得到绝对欧氏距离公式，变换如下：d(A,B)=i=1nUA(xi)-UB(xi)2（2）由文献13可知，相对欧氏距离公式为：(A,B)=1ni=1nUA(xi)-UB(xi)2（3）根据模糊数学运算要求，需通过隶属函数对模糊子集映射，使U(xi)的取值在（0,1）之间，由于隶属函数的选取对数据处理结果影响加大。为了真实反映测点地下介质的垂向变化，尽量得到可靠有用信息，用U(xi)取电测深测量数据进行计算。2电阻率测深法测量数据的计算根据物探方法“由已知到未知”的解释原则，在使用闵氏距离公式对电阻率测深法测量数据的计算时，具体测量数据计算步骤：（1）首先在出水量满足找水预定指标的井位进行测量试验，确定合适的电极距，测量结果作为UA(xi)；（2）在预定区域找水点位，保持与已知井位相同的电极距，测量结果为Ui(xi)（i=1,2，n)；（3）利用公式（3）计算相对欧氏距离；（4）对所有计算结果进行有小到大排序。一般情况下，工程物探找水的区域离已知井位范围较近，水文地质情况差别不是很大。值越小，表示其与已知井位的距离越小，其水文地质情况与已知井位越接近，其测点处找到水可能性越大，可以将较小作为预设井位。3应用实例为尽快解决豫西地区人畜饮水困难，河南省政府积极组织动员各方力量找水打井，我院作为技术支撑单位，积极参与本次抗旱找水任务，并开展水文地质和物探勘查工作，力争以最快的速度和最好的找水效果解决当地村庄居民和牲畜生活饮用水问题。依据收集当地水文地质和矿产地质资料及前人工作资料，查得1 20万水文地质图，地下水类型为浅层孔隙水和基岩裂隙水。地下水径流自北部山区流向南部平原，补给平原区松散岩类孔隙水。区域地层岩性为：松散层岩性为粉质粘土、粘土夹钙质结核层和砂砾石、中细砂等；基岩为砂岩、砂砾岩和泥岩。主要的含水层为松散岩类地层的钙质结核层和砂层，预测下部的基岩裂隙水水量较小。根据水文地质情况，在任务区范围内开展对称四极电阻率测深工作，在工作区有一眼机井，单井涌水量712m3/h，其用水量基本到达找水指标。将该井电阻率测深法数据作为已知井位，测量结果作为UA(xi)。在找水任务区范围内测点 11 个，依次作为 Ui(xi)（i=1,2，11)。根据找水的有效数据信息，去掉了前两个电极距数据，电阻率测深法原始数据详见表1。A(试验点）及B1（最终确定出水量最佳井位）电测深曲线图分别如图1、图2所示。利用公式（3），计算待定井位测点与已知井位间的相对欧氏距离，并对由小到大依次排序，详见表2。由表2统计可知，相对欧氏距离最小值为B11测点，后面依次为B1、B9、B5、B10、B3、B2、B6、B7、B4、B8。物探组根据计算结果向工程组提交了5个打井位置。为表1电阻率测深原始数据s值(单位：m)点号AB1B2B3B4B5B6B7B8B9B10B11AB/21010.623.511.710.916.618.322.812.113.412.623.013.9159.918.010.210.512.814.119.010.910.911.919.711.5259.118.810.410.712.112.014.610.511.911.213.39.4359.116.010.211.411.312.012.210.312.69.99.78.3509.114.112.111.511.612.010.310.515.28.98.47.87510.311.810.214.115.012.511.812.921.49.89.98.810011.59.619.618.322.012.314.816.626.711.511.610.515015.613.829.428.234.012.722.524.137.616.316.314.720018.316.838.438.142.713.229.731.846.421.921.619.325022.523.545.545.253.713.538.238.752.927.026.824.430027.529.346.442.566.715.044.945.159.032.431.929.935026.527.444.038.548.516.351.252.764.637.336.234.240024.425.539.637.046.317.856.062.276.542.741.338.91522023年第8期西部探矿工程了能够真正满足实际用水需求，现场定打井位置时，征求了当地村委会的意见，B11、B9、B10不作为打井位置，B5与B6、B7、B8处于同一个村庄，当地村委会急需在B4位置打井，最后确定打井位置B1、B5、B3、B2、B4。根据物探测量的建议，工程组在上述测点位置进行了打钻施工。经过成井、洗井后，进行了抽水试验，B1测点单井出水量24m3/h，B5测点单井出水量18m3/h，B3测点单井出水量14m3/h，B2测点单井出水量12m3/h，B4测点单井出水量5m3/h。B1、B5、B3、B2到达指标要求，B4没有达到指标要求。从抽水试验出水量分析，井位出水量和相对闵氏距离成正比关系。通过对水质分析，各种常规元素含量、微量元素含量测定，属于合格饮用水，同时满足农业灌溉水质标准。4结论通过实测数据的计算分析，欧氏距离法对于电阻率测深法找水快速确定井位是切实可行的，其可以对电阻率测深法曲线进行定量解释，避免了单纯从定性资料解释曲线类型相似无法确定最佳井位的缺点，为电阻率测深法找水提供了理论依据。通过该方法实践，其理论基础可靠，方法易于掌握，计算过程简单，并且计算结果与实际情况符合较好，充分显示了欧氏距离法在电测深找水中的有效性和优越性。但在使用该方法需注意：真实掌握找水区域水文地质情况，大致掌握含水层的深度；确定该区域出水量已知井位，因预设井位测量数据以该井测量数据为标准，至关重要；保持预设井位测量电极距和已知井位一致，便于资料处理解释。总之，欧氏距离法对于电阻率测深法测量资料解释，是行之有效的一种方法，对于野外快速确定最佳井位有很强的指导意义。参考文献：1刘天佑.应用地球物理数据采集与处理M.武汉:中国地质大学出版社,2009:170-182.2王聿军.贫水山区电法找水的应用实例J.地质与勘探，2001,37(3):73-74.3陈灿华,陈绍求.电测深法在岩溶探测中的应用J.中南工业大学学报,2002,31(1):9-2.4何公民.谈电测深法野外工作J.煤炭技术,2009,28(7):157-158.5刘福臣,黄怀峰,马祥配.电法勘探找水野外常见干扰问题分析与处理J.水利水电科技进展,2008,28(5):74-77.6裴正文.试论地球物理与现代水文、工程、环境地质J.水文地质工程地质,1999,26(3):58-60.7贺仲雄.模糊数学及其应用M.天津:天津科学技术出版社,1983:20-30.8才庆喜.模糊数学方法在电测深曲线对比中的应用J.河北地质学院学报,1987,11(2):69-76.9才庆喜.Fuzzy贴近度在电测深找水数据解释中的应用J.中国煤炭地质,1991,3(1):67-70.10 庞绪贵.应用灰色关联度分析对电测深找水资料进行数据处理J.工程勘察,1996(6):67-68.11 李明山.模糊相似比在电测深找水中的应用J.EngineeringSurvey,1996(6):69,36.12 常聚友.模糊相似比在第四纪电探找水中的应用J.物探化探计算技术,1996,18(3):254-259.13 丁华,张志勇,潘绍明.欧氏距离法在电测深找水中应用的可行性探讨J.工程勘察,2001,25(4):285-289.14 李保国.电反射系数（K）法在解释地质效果中的应用J.水文地质工程地质,2001,28(6):66-6.15 孙金龙.梯度K法在电测深找水中的应用J.水文地质工程地质,2002(6):58-60.16 黄真萍,褚高强,俞侃.衍生指数在电测深找水中的研究与应用J.工程地质学报,2010,18(S1):244-248.17 许天周.应用泛函分析M.北京:科技出版社,2002:47.图1A电测深曲线图图2B1电测深曲线图表2相对欧式距离值点号B11B1B91.2391.4621.636点号B5B10B31.6431.8942.929点号B2B6B73.2813.6653.912点号B4B84.8636.313153