频率域CSEM全区视电阻率数值计算方法.pdfVIP

中国地球物理2011 ·251· (7)电磁方法研究与应用 频率域CSEM全区视电阻率的数值计算方法 NumericalCalculationMethodsforFull—wave ApparentResistivity of DomainCSEM Frequency 王顺国1+熊彬2李帝铨1 bin2Li WangShungu01’XiongDiquanl 1．中南大学 地球科学与信息物理学院 长沙410083；2．桂林理工大学地球科学学院 桂林541004 在电法勘探中，通过引入适当的视电阻率定义使得视电阻率响应曲线比观测的场值更好地反映地下电 性界面的存在，这是地球物理勘探中电法较之于其他方法的独特之处。频率域CSEM在准静态近似条件下， 能得到场的解析表达式，然表达式中场值与电阻率之间存在着复杂的隐函数关系，难以用解析方法直接求 出电阻率与场之间的显示反函数，因此视电阻率定义通常采用各种近似方法。然各种近似定义方法在前提 条件满足时能很好地反映地下电性变化。条件不满足时，则不能真正反映地下电性分布规律。且在实际工 作中很难准确判断这些前提条件满足与否，这就给数据解释工作造成困难。国内外学者对上述问题进行深 入研究，确定全区视电阻率定义比近似定义更好地反映了地下介质的电性变化规律，故采用全区视电阻率 可有效解决上述问题。 目前，频率域CSEM全区视电阻率常用迭代法和逆插值法求解。为了对比二者的优劣，对电偶源激励 下层，H，，H：，z删，Re(乙)所定义的全区视电阻率用这两种数值方法分别进行了计算，并对计算结 果做了较深入的分析。 结果表明，迭代法受多种因素的影响，主要是迭代初值的选取、迭代期望精度￡的设置、迭代次数Ⅳ的控 制以及迭代公式的收敛性。迭代初值选取的越接近真值，迭代收敛速度越快，采用远区近似表达式的计算结果 作为迭代初值，就是比较理想的选择。迭代进行时，当视电阻率的相对误差I角一pol／／90小于迭代期望精度￡(￡ 为一个较小的正数)时，迭代结束，计算中取￡一0．5×10。2即能保证精度要求，如若￡过小，迭代可能因达到 次数上限而终止。迭代次数Ⅳ的选取也不宜过大，迭代结果在迭代进行若干次后将不再明显的提高精度。而迭 代公式的收敛性从根本上决定迭代是否成功，如果直接迭代发散，可以修正迭代步长Ap(如使△n一△风／2， 这里△朋即为新的步长，在计算H，定义的全区视电阻率时，就用到了该方法)，也可修正迭代公式，以使迭代 收敛。由于收敛性问题，使得迭代算法不稳定，其计算精度和时间会随收敛性而发生变化。 逆插值法主要受以下因素的影响，即插值目标函数的单调性、插值算法的选取和插值区间的确定。如 果使用逆插值法而不能在插值区间上保证插值目标函数严格单调，一个场值(如H，)将对应多个n值， 此时逆插值结果就会出现错误，令得到的全区视电阻率曲线出现畸变。所以在具体计算时，首先要对插值 目标函数的单调性进行验证。对于非单调目标函数f(p)，可以引入转折点∥，使得在转折点划分的各插 值区域中f(p)是单调的，然后逐段实施逆插值，即能正确得出各场值对应的全区视电阻率值。但当函数 为了能确保计算结果满足精度要求，选择精度较高的三次样条插值，这样，插值算法选取问题就得到了妥 单调性问题，逆插值法将不利于计算单调性复杂的场值表达式所导出的全区视电阻率。 总体比较，两种方法所得合理结果都能正确反映地下介质的电性变化规律。逆插值的影响因素相对较 少，计算过程较稳定，然迭代法不需要验证和处理目标函数单调性问题，具体实施相对简单。在实际应用 中，可根据目标函数的特点，择优选取方法。 参考文献 【11何继善．广域电磁法和伪随机信号电法【M】．北京：高等教育出版社，2010：72～119． 【2】熊彬．大回线瞬变电磁法全区视电阻率的逆样条插值计算【J】．吉林大学学报(地球科学版)，2005。35(4)：515—519 【3】汤井田，何继普．可控源音频大地电磁法及其应用【M】．长沙；中南大学出版社，2005