地下铁质管线埋深磁偶极子构造法优化反演制造技术

本发明专利技术涉及地下铁质管线埋深磁偶极子构造法优化反演方法，具体涉及一种基于磁偶极子构造法正演的地下铁质管线埋深优化反演计算方法。建立磁异常探测数据预处理、地磁场定位查询及反演搜索模型参数替代，将十二个正演模型参数简化为埋深和磁化率两个未知模型参数。埋深和管道磁化率两个未知模型参数作为全局优化反演求解参数，使用初次反演得到埋深和管道磁化率作为第二次反演参数搜索范围参考点，在新搜索范围内再次进行多次反演搜索，取埋深反演平均值作为地下铁质管线埋深优化反演结果。本发明专利技术地下铁质管线埋深磁偶极子构造法优化反演兼顾正演计算速率和精度，反演求解快速且准确度高。

Underground iron pipeline; buried depth; magnetic dipole structure method; optimal inversion

The invention relates to an optimized inversion method for the buried depth of an underground iron pipeline and a magnetic dipole structure method, in particular to an optimized inversion calculation method for buried iron pipelines forward based on the magnetic dipole structure method. The establishment of magnetic anomaly detection data preprocessing, geomagnetic field location query and inversion search model parameter substitution, twelve forward model parameters are simplified to two unknown model parameters of buried depth and magnetic susceptibility. The depth and pipe susceptibility of two unknown parameters as global optimization inversion parameters, using the first inversion depth and pipeline susceptibility as the second parameter inversion search reference point, multiple inversion search again in the new search range, the depth of the inversion average value as the underground iron pipeline buried depth optimization inversion results. The underground iron pipeline buried depth is optimized by the magnetic dipole structure method, and the forward calculation speed and the accuracy are taken into account, and the inversion solution is fast and accurate.

【技术实现步骤摘要】
地下铁质管线埋深磁偶极子构造法优化反演
本专利技术涉及一种使用探测磁异常信号反演求解埋地铁质管线埋深的方法，具体涉及地下铁质管线埋深磁偶极子构造法优化反演。
技术介绍
地下埋设铁质管线主要有污水、雨水、给水、电缆、燃气、电信、电力和石油管线，为国家工业生产和城镇居民生活中提供基础保障作用。社会经济的发展，地下铁质管线数量的增多，其建设、维护和管理的成本越来越高，因此地下管线的探测变得越来越重要。随着技术的不断发展和完善，对各种管线都有自身适应的探测方法，其勘测施工难度、探测深度范围、经济投入及准确性各不相同。磁异常铁质管线探测作为一种被动源管线探测方法，使用地球自身的磁场作为激励磁源，不需要人工提供磁源发生装置，使得该方法相对其他管线探测方法而言具有硬件设备简单，操作容易，投入成本低等，探测位置准确优点。目前，使用磁异常探测地下铁质管线的商业仪器是基于探测总磁异常分布来判断地下管线的有无及大致走向，而不能计算地下铁质管线的埋深。专门针对地下铁质管线磁异常探测，利用数据处理方式来计算管道埋深的方法报道较少，已有技术报道使用磁异常化极估算的方法计算管线埋深精度较低，且容易受噪声影响。另有论文报道使用差分能量谱估计法估算出埋深在印度某小镇附近的一段供水管的深度为3.88±0.10m，而使用磁感应强度大小与磁异常梯度张量的解析信号的比值的方法来估计这段管道的埋深时，其估计埋深为4.09m，且两种方法都未给出准确性分析。另有文献报道使用磁异常梯度的方法探测估计某水管的埋深时，深度估计误差能够小于5％。其它使用磁数据处理估算磁源模型深度的方法主要有Werner反褶积法，解析信号，Euler反褶积法，局部波数法，及源参数成像法等，这些方法尚未在地下铁质管线探测埋深计算方面应用，无法知道其准确性如何。基于磁异常反演求解地下铁质管线埋深的方法较少，本专利技术前期研究报道使用遗传算法基于四个反演参数进行过埋深反演，但该方法对于埋深较浅的管线需要进行分块单元划分，反演花费时间较长。其他相关磁异常反演方法尚未应用于地下铁质管线探测埋深反演，无法判断其有效性。
技术实现思路
本专利技术的目的是针对地下铁质管线磁异常埋深求解现有的问题，提供了一种地下铁质管线埋深磁偶极子构造法优化反演求解方法。本专利技术是通过以下技术方案实现地下铁质管线埋深磁偶极子构造法优化反演，具体步骤如下：(1)建立正演模型：以地理北为X轴，地理东为Y轴，垂直向下为Z轴，探测地平面为XY平面，建立计算坐标系；坐标系中，设所致探测磁异常有效管道长度为Lp，管道外径为φ，管道厚度为δ，管道磁化率为χm；设管道方位角为θ(0≤θ≤π)，管道轴线埋深为h1；地磁背景场由总强度F，磁倾角Ic及磁偏角D三个独立的地磁要素表征；磁异常正演时，管道轴线位于坐标系原点正下方，管道路径垂直于管道轴线，管道路径的中点在Z轴上，高度为h2，长度为Ld，管道路径方位角为η；将管道路径上总磁异常大小或总磁异常分量Bd表示为十二个参数表征的函数：Bd＝f(Lp,φ,δ,θ,h1,F,Ic,D,χm,h2,Ld,η)①(2)测点磁异常计算：将管道分割成有限数量的单元块，将管道分割所得单元视作坐标为单元几何中心的磁偶极子，使用磁偶极子磁场计算公式计算每一个单元在测点处的磁感应强度，使用矢量和计算所有单元在测点处的叠加磁感应强度，该叠加磁感应强度加上该点的地磁场强度即为测点磁异常；设单元块Ei几何中心坐标为(Eix,Eiy,Eiz)，单元体积为Vi，测点Pj的坐标为(Pjx,Pjy,Pjz)；则有磁偶极子磁场计算公式可得单元块Ei在测点Pj处的磁感应强度Bij三分量为：上式②中μ0为真空磁导率；当Vi为分块单元时当Vi为分块单元时Vi＝πφδ(φ-δ)；上式②中，Hx＝Fcos(Ic)cos(D)，Hy＝Fcos(Ic)sin(D)，Hz＝Fsin(Ic)；则根据磁场叠加原理，地下铁质管线在探测面Pj点的磁感应强度三分量为：上式③中I为构造地下铁质管道的单元总数；当管道模型参数和单元划分方式确定后，单元块Ei几何中心坐标即随之确定；当探测路径及测点间隔确定后，测点Pj坐标也为随之确定；当使用探测总磁异常大小作为观测值时，则使用模拟总磁异常作为匹配值；当使用探测总磁异常轴向分量作为观测值时，则使用③式中对应的轴向分量作为匹配值；测线上所有测点的总磁异常的大小或总磁异常分量组成的向量为Bsim；(3)探测计算管道方位角：初步选择几条探测直线路径探测地下铁质管线所致总磁异常分布，确定管道大致走向后，在垂直管道方向上选取两条相隔一定距离的探测路径测量总磁异常，使用曲线拟合的方式找到两条探测曲线的峰值点，使用两个峰值点对应的水平坐标(B#1x,B#1y)、(B#2x,B#2y)计算管道方位角为：当θ确认后，测线垂直于管道轴线方位角η等于90°；使用两个峰值点连线初步确定管道的水平位置；(4)磁异常测量：垂直管道轴线方向布设探测路径，探测路径长度Ld范围在16m到30m之间；探测路径中点设于初步确定的管道的正上方；每隔10cm左右测量记录一次总磁异常或总磁异常分离，得到探测磁异常向量Bobs；(5)建立反演模型：使用GPS定位得到探测区域坐标，查询国际地磁参考场得到地磁总强度F，磁倾角Ic及磁偏角D；反演时使用有限长管道作为反演管道长度，有效管道长度Lp在18m至30m范围选取；测线高度h2在探测磁异常时确定；反演时，使用替代管道模型厚度δ与管道外径φ作为确定的参数，管道厚度δ选取范围为5mm至15mm，管道外径φ选取范围为20cm至80cm；此时建立参数反演模型：(h1,χm)＝f-1(Bobs)⑤(6)反演数据预处理及反演目标函数的确定：使用插值方法找到模拟计算磁异常向量Bsim的最大值或者最小值，使用拟合方式找到探测磁异常向量Bobs拟合曲线上的最大值或者最小值，以上述两个最大值点或两个最小值点作为反演数据定位点，按一定测点间隔选取测线长度为16m至20m的探测磁异常向量Bobs的拟合计算得到目标磁异常向量并选取对应坐标点由计算磁异常Bsim插值计算得到匹配磁异常向量更新⑤式反演模拟为建立反演目标函数：(7)管道埋深参数优化反演：使用全局优化算法搜索找到式⑥中使得目标函数尽量小的埋深h1和磁化率χm，优化反演得到的埋深h1即为地下铁质管线的轴线埋深。对上述内容的进一步补充，使用探测总磁异常的数值大小作为目标求解磁异常向量时，计算总磁异常的数值大小作为匹配磁异常使用探测磁场某一坐标轴向分量，如探测磁异常垂直分量作为目标求解磁异常向量时，计算磁异常某一坐标轴向分量，如计算磁异常垂直分量作为匹配磁异常向量对上述内容的进一步补充，地下铁质管线埋深磁偶极子构造法优化反演，所述单元划分，将探测路径中点正下方对应的管道轴线上的点作为单元划分的开始端点，沿管道轴线方向向管道两端分割单元；首先将管道从起始点开始分割成长度为管道直径的节单元，当节单元的几何中心到测线中点的距离大于6.5倍的管道直径时，保留该节单元，继续分割节单元，直至长度为探测路径长度Ld的管道单元划分完毕；当所分节单元的几何中心到测线中点的距离小于6.5倍的管道直径时，取消该节单元划分，沿着轴向分割为环，环宽度为δ，将每一个环沿圆周分割为块单元，每一个管道本文档来自技高网...

【技术保护点】
地下铁质管线埋深磁偶极子构造法优化反演，其特征在于，具体步骤如下：(1)建立正演模型：以地理北为X轴，地理东为Y轴，垂直向下为Z轴，探测地平面为XY平面，建立计算坐标系；坐标系中，设所致探测磁异常有效管道长度为L

【技术特征摘要】
1.地下铁质管线埋深磁偶极子构造法优化反演，其特征在于，具体步骤如下：(1)建立正演模型：以地理北为X轴，地理东为Y轴，垂直向下为Z轴，探测地平面为XY平面，建立计算坐标系；坐标系中，设所致探测磁异常有效管道长度为Lp，管道外径为φ，管道厚度为δ，管道磁化率为χm；设管道方位角为θ(0≤θ≤π)，管道轴线埋深为h1；地磁背景场由总强度F、磁倾角Ic及磁偏角D三个独立的地磁要素表征；磁异常正演时，管道轴线位于坐标系原点正下方，探测路径垂直于管道轴线，探测路径的中点在Z轴上，高度为h2，长度为Ld，探测路径方位角为η；将探测路径上总磁异常大小或总磁异常分量Bd表示为十二个参数表征的函数：Bd＝f(Lp,φ,δ,θ,h1,F,Ic,D,χm,h2,Ld,η)；(2)测点磁异常计算：将管道分割成有限数量的单元块，将所得单元视作坐标为单元几何中心的磁偶极子，使用磁偶极子磁场计算公式计算每一个单元所在测点处的磁感应强度，使用矢量和计算所有单元在测点处的叠加磁感应强度，该叠加磁感应强度加上该点的地磁场强度即为测点磁异常；设测线上所有测点的总磁异常的大小或总磁异常分量组成的向量为计算磁异常向量Bsim；(3)探测计算管道方位角：初步选择几条探测直线路径探测地下铁质管线所致总磁异常分布，确定管道大致走向后，在垂直管道方向上选取两条相隔一定距离的探测路径测量总磁异常，使用曲线拟合的方式找到两条探测曲线的峰值点，使用两个峰值点坐标计算管道方位角θ，测线垂直于管道轴线方位角η等于90°,初步确定管道的水平位置；(4)磁异常测量：垂直管道轴线方向布设探测路径，探测路径长度Ld范围在16m到30m之间，探测路径中点设于初步确定的管道的正上方，每隔10cm左右测量记录一次总磁异常或总磁异常分离，得到探测磁异常向量Bobs；(5)建立反演模型：使用GPS定位得到探测区域坐标，查询国际地磁参考场得到地磁总强度F，磁倾角Ic及磁偏角D；反演时使用有限长管道作为反演管道长度，有效管道长度Lp在18m至30m范围选取；测线高度h2在探测磁异常时确定；反演时，使用替代管道厚度δ与管道外径φ作为确定的参数，管道厚度δ选取范围为5mm至15mm，管道外径φ选取范围为20cm至80cm；此时建立参数反演模型：(h1,χm)＝f-1(Bobs)；(6)反演数据预处理及反演目标函数的确定：使用插值方法找到模拟计算磁异常向量Bsim的最大值或者最小值，使用拟合方式找到探测磁异常向量Bobs拟合曲线上的最大值或者最小值，以上述两个最大值点或两个最小值点作为反演数据定位点，按一定...

【专利技术属性】
技术研发人员：郭智勇，刘先平，杜坚，梁海波，张禾，赵丹丹，刘得军，
申请(专利权)人：西南石油大学，中国石油大学北京，
类型：发明
国别省市：四川,51