计算矩形大定源回线瞬变电磁测深全区视电阻率

介绍了矩形大定源瞬变电磁测深全区视电阻率的计算方法。根据大定源感应电动势关于电阻率非线性变化特点,利用二分法计算感应电动势的全区视电阻率。在此过程中,计算首先从早期感应电动势开始,并将计算得到的视电阻率作为下个窗口视电阻率计算的依据。以此类推,得到全部时窗的视电阻率。理论模型和实际数据的计算结果表明文中提出方法的有效性和正确性。     

磁源瞬变电磁法中心回线装置全程视电阻率定义及计算

讨论了磁源TEM中心回线装置全程视电阻率的计算问题。利用电偶极子的瞬变场叠加来计算磁源TEM中心回线装置的瞬变场,均匀半空间的计算结果表明,该方法是可行的。用电偶极子的瞬变场定义了磁源瞬变电磁法中心回线装置全程视电阻率,通过对二、三、四层电性断面的视电阻率的计算,并与TEM晚期视电阻率进行了对比,说明了该定义的优越性。实测磁源TEM资料的计算表明,用电偶极子的瞬变场定义的磁源瞬变电磁法中心回线装置全程视电阻率较好地反映了电性断面特征,其对浅层磁源TEM探测及资料解译具有使用价值。     

全空间瞬变电磁法资料解释方法中的平移算法

全空间视电阻率公式及全区求解方法是地下瞬变电磁理论的重要组成部分。为探讨全空间全区视电阻率求解方法,笔者将半空间全区视电阻率的平移算法引入至全空间并对其进行了研究。同时,分别采用理论数据、正演模拟数据及实测数据对该方法在全空间资料处理与解释中的有效性和精确性进行了分析。结果显示,平移算法与常规的全区视电阻率解法具有相同的理论基础,但前者避免了传统全区计算方法中的多解性;此外,采用由截距参数绘制的时间剖面图和视电阻率断面图进行联合解释产生了较好的应用效果。     

TEM框内回线装置发射框边界影响及消除方法

探讨了大定源回线框内观测装置的视电阻率计算问题。中心回线晚期视电阻率计算由于受边界影响及不满足晚期条件,计算的视电阻率值有畸变,影响了煤田采空区及富水性的解释。利用电偶源时域垂直磁场求取全区、全期视电阻率,解决了发射线框的边界影响问题及早期道视电阻率计算问题。理论模型计算和实测资料的处理结果表明了该方法的有效性。      

电性源瞬变电磁全区视电阻率定义

电性源瞬变电磁法是目前研究的热点,研究该方法的全区视电阻率定义同样很重要。文中给出了一种快速有效的全区视电阻率定义方法,利用该方法可以快速计算电性源瞬变电磁多分量的全区视电阻率。利用均匀半空间磁场强度响应曲线的平移伸缩特性实现视电阻率定义,该算法无需迭代就可以直接计算,速度快,精度高。最后,对比了水平磁场分量和垂直磁场分量定义的全区视电阻率对地层的分辨能力。结果表明,两种分量定义的全区视电阻率有相同的分辨能力,给电性源瞬变电磁法的多分量综合解释技术奠定了理论基础。     

TEM全区视电阻率在煤田采空区注浆效果检测的应用

文章分析总结了前人各类求取回线源瞬变电磁全区视电阻率的方法。首先介绍一维层状介质中心回线瞬变电磁测深法的全区视电阻率定义,研究总结常用的全区视电阻率数值计算方法,介绍了隐函数的反函数求解问题,并对比两种反演方法,对典型地电断面模型进行正演模拟。最后将全区视电阻率及反演运用到煤矿采空区注浆效果验证中,取得较好效果。     

瞬变电磁全区视电阻率数值计算方法研究

基于均匀半无限介质空间中场的表达式,给出了一维层状介质中心回线装置瞬变电磁测深法的全区视电阻率定义,研究了三种全区视电阻率的数值计算方法。特别是在计算的过程中,利用分段计算,解决了全区视电阻率解的非唯一性问题。同时引入二分查找法,解决了隐函数的反函数求解问题。通过正演建模和烟圈反演,进一步对比分析了三种计算方法,这在瞬变电磁测深资料的定性解释中,正确反映地下一维电性层电性垂向的变化特征具有重要的意义。     

大定源瞬变电磁三分量全域视电阻率定义与三分量联合反演

多数瞬变电磁仪器可以同时采集三个分量的响应,但目前瞬变电磁的解释方法却主要以中心回线垂直分量为主。为了充分利用瞬变电磁三分量的信息,介绍一种大定源瞬变电磁三分量约束反演方法。基本思想是利用瞬变电磁三个分量的场分别求取全区视电阻率,并以此为基础建立初始模型,利用阻尼最小二乘法进行三分量反演。首先,将回线源看成是多个电偶极子产生的场的叠加,利用偶极子叠加的原理求取回线的三分量电磁响应,并推导视电阻率的全区与全期定义(即全域视电阻率定义);然后,利用烟圈理论求取视深度,并以此建立反演的初始地电模型。最后,利用阻尼最小二乘法实现三分量约束反演。利用文中视电阻率定义方法对已知模型数据进行处理,视电阻率曲线与模型的电性参数变化一致;以视电阻率结果为依据建立初始模型,利用三分量约束反演方法对层状模型数据进行反演,结果与模型吻合较好。     

大定圆回线频率域电磁场特征研究

大定圆回线非中心回线装置的电磁场计算公式中含有双重贝塞尔函数,因其高震荡性,使用常规方法计算误差较大。为提高计算精度,本文使用一种通过改进核函数的快速汉克尔变换算法来进行求解,计算结果表明:该算法能够很好的压制贝塞尔函数的强震荡性对计算结果的干扰,即使在非准静态近似的条件下,也能精确地计算出大圆回线源位于地表及接收点位于地面任意一点的电磁场,此外,大定源回线电磁场的分布规律不仅与接收装置位于回线内、外有关,还与圆回线所处的位置成一定的线性关系。     

矩形大定源瞬变测深数据模型约束反演

利用一维光滑模型约束反演对矩形大定源回线的数据进行反演,大定源回线瞬变电磁响应采用计算速度很快的小方形回线等效叠加方法获得.光滑模型反演分别采用基于平坦约束和光滑约束进行,通过对层状模型理论的垂直磁场感应电动势数据和实际数据进行反演,结果表明,采用光滑模型约束,能较稳定地获得地下电阻率的合理分布.     

TEM中心回线法瞬变磁场求取和全区视电阻率计算

依据TEM法中心回线装置瞬变电磁场的衰减规律,给出了由实测的dB/dt求取瞬变磁场B的简单方法.详细推导了发射电流为斜阶跃波条件下由B迭代反演全区视电阻率的方法.为提高计算精度和速度,对计算过程中的积分、求导推导了解析表达式和采用标准数据库的方法计算含有概率积分的函数.对实测资料进行了计算,通过算例进行了验证和分析.     

近区瞬变电磁测深曲线畸变段的消除方法

讨论了利用P.Raab和F.Frischknecht给出的中心回线装置瞬变衰减电压表达式,采用编程求解的方法来计算适用于全时域的视电阻率ρ_τ。这种视电阻率曲线克服了以往使用晚期计算视电阻率公式所得到的曲线,对于任何模型的左支曲线都将产生畸变的缺点,它的早期渐近线为趋于第一层电阻率值的水平直线。      

对瞬变电磁测深几个问题的思考(三)——瞬变电磁场关断效应及全期视电阻率的普适算法

首先给出瞬变电磁（TEM）场卷积计算式,并将关断电流变化率函数看成冲击函数,垂直阶跃脉冲激发的场分量看成输入信号,卷积结果就是任意形状关断电流的TEM响应。把电流激发脉冲看成是许多垂直阶跃函数的叠加,通过叠加这些阶跃函数的响应求得任意电流脉冲激发的瞬变电磁响应。当已知TEM场的实测值和关断电流形状,就可按编程的方法计算全期视电阻率。文中给出了2个斜坡关断电流激发的电偶极源TEM垂直磁场算例,一个为均匀半空间的响应及其全期视电阻率,另一个为三层K型地电断面的响应及其全期视电阻率。      

Normalized impedance function and the straightforward inversion scheme for magnetotelluric data

This paper investigates the performance of normalized response function obtained by normalizing the Cagniard impedance function by a suitable factor and then rotating the phase by 45‡ to make it purely real for homogeneous half-space and equal to the square root of the half-space resistivity. Two apparent resistivity functions based on respectively the real and imaginary parts of this response function are proposed. The apparent resistivity function using the real part contains almost the same information as that yielded by the Cagniard expression while the one using the imaginary part qualitatively works as an indicator of the number of interfaces in the earth model. The linear straightforward inversion scheme (SIS), developed by the authors employing the concept of equal penetration layers, has been used to validate the proposed apparent resistivity functions. For this purpose, several synthetic and field models have been examined. Five synthetic models are studied to establish the veracity of the new functions and two well-studied published field data sets are inverted through SIS for comparison. We noticed that the new function and SIS compliment each other and lead to better understanding of the data information and model resolution.     

适用于袖珍计算机频率测深数字滤波方法

详细推导了电偶源频率测深单分量电场定义的视电阻率数字滤波计算公式,并且应用最小二乘法设计了只有30个系数的滤波器,编制出适用于袖珍计算机的视电阻率计算程序,使计算速度大为提高.与用均匀半空间的精确公式计算的视电阻率对比,最大相对误差<2%.     

基于核最小二乘模型的矿产靶区预测

地质统计单元的含矿性与地质找矿证据之间存在复杂的非线性关系,建立这种复杂关系的多元非线性统计模型并预测矿产靶区,对矿产勘查具有重要指导意义。以核函数为理论工具,在核最小二乘原理基础上提出了矿产靶区预测的核最小二乘模型;在GDAL数字图像输入输出函数库和CLAPACK线性代数软件包基础上,用VC++语言开发了面向栅格数据的矿产靶区预测核最小二乘模型算法程序,并把模型应用于新疆阿勒泰地区的矿产靶区预测研究。在MapInfo中生成包含100×151个网格统计单元的栅格图层,把栅格化后的15种找矿证据图层转化成100×151×15的数字图像数据立方体,用自行开发的程序计算每个网格统计单元的核最小二乘判别得分。结果表明,网格统计单元判别得分的高值区与已知矿床（点）的空间分布基本一致。     

烟圈反演和视纵向电导解释方法对比

基于浮动薄板理论,研究瞬变电磁场层状介质中视纵向电导的计算方法,利用遗传算法代替常用的阻尼最小二乘法来获取关键参数。进行“烟圈”反演计算,并分别采用3层H (K)型和4层HK型地电模型对比视纵向电导曲线和电阻率—深度曲线。结果表明,“烟圈”反演具有不需初始模型、计算简便等优点,在定性反演方面,视纵向电导解释方法可以更准确地划分电性层,且能够比较直观地体现电性变化趋势。      

Extreme Learning Machine Based Modeling of Resilient Modulus of Subgrade Soils

This paper investigates the potential of two variants of extreme learning machine based regression approaches in predicting the resilient modulus of cohesive soils. Support vector regression was used to compare the performance of the proposed extreme learning machine based regression approaches. The dataset used in this study was derived from literature and consists of 9 input parameters with a total of 891 cases. For testing, two methods i.e. train/test and tenfold cross validation was used. In case of train and test methods, a total of 594 randomly selected cases were used to train different algorithms and the remaining 297 data were used to test the created models. Correlation coefficient value of 0.991 (root mean square error = 3.47 MPa) was achieved by polynomial kernel based extreme learning machine in comparison to 0.990 and 0.990 (root mean square error = 4.790 and 4.290 MPa) by simple extreme learning machine and radial basis kernel function based support vector regression respectively with test dataset. Comparisons of results with tenfold cross validation also suggest that polynomial kernel based extreme learning machine works well in terms of root mean square error and computational cost with the used dataset. Sensitivity analysis suggests the importance of confining stress and deviator stress in predicting the resilient modulus when using with polynomial kernel based extreme learning machine modeling approach.