E-mail Alert Rss
 

物探与化探, 2019, 43(3): 584-588 doi: 10.11720/wtyht.2019.1312

方法研究·仪器研制

大回线源瞬变电磁一维自适应反演方法及应用

姚伟华

中煤科工集团 西安研究院有限公司,陕西 西安 710077

The one-dimensional adaptive inversion method for large loop source TEM and its application

YAO Wei-Hua

Xi’an Research Institute of China Coal Technology & Engineering Group,Xi’an 710077,China

责任编辑: 沈效群

收稿日期: 2018-09-3   修回日期: 2018-11-8   网络出版日期: 2019-06-20

基金资助: 国家重点基础研究发展计划(973计划)资助项目.  2017YFC0804100
中国煤炭科工集团有限公司科技创新创业资金专项.  2018MS006

Received: 2018-09-3   Revised: 2018-11-8   Online: 2019-06-20

摘要

提出了一种适用于大回线源瞬变电磁法框内任意测点垂直Z分量的一维反演方法。该算法采用CMD正则化因子优化方案、光滑度函数和自然边界条件对反演过程进行了约束,减小了初值的依赖度,保证了反演过程的自主性和反演结果的稳定性;在反演过程中,在保证模型残差沿梯度下降的同时,对模型修正量的最大值进行了限制,使反演模型不易陷入局部极小,保证了反演结果全局的收敛性;同时反演的灵敏度矩阵采用牛顿法进行更新,减少了整个反演计算的工作量。最后通过五层的双低阻层模型、双高阻层模型和实测资料证明了该算法可靠。

关键词： 瞬变电磁法 ; 自适应 ; 一维约束反演

Abstract

A one-dimensional inversion method for the vertical Z component of any measuring point in the frame of the large loop source transient electromagnetic method is proposed in this paper. The algorithm uses the CMD regularization factor optimization scheme, smoothness function and natural boundary conditions to constrain the inversion process, reduce the initial value dependence, and ensure the autonomy of the inversion process and the stability of the inversion results. In the inversion process, while ensuring that the model residual decreases along the gradient, the maximum value of the model correction is limited, so that the inversion model is not easy to fall into the local minimum, which ensures the global convergence of the inversion results. The sensitivity matrix is updated by Newton's method, thus reducing the workload of the entire inversion calculation. Finally, the reliability of the algorithm is guaranteed by five-layer double low-resistance layer model, double high-resistance layer model and measured data.

Keywords： TEM ; adaptive ; one-dimensional constraint inversion

PDF (1371KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

姚伟华. 大回线源瞬变电磁一维自适应反演方法及应用. 物探与化探[J], 2019, 43(3): 584-588 doi:10.11720/wtyht.2019.1312

YAO Wei-Hua. The one-dimensional adaptive inversion method for large loop source TEM and its application. Geophysical and Geochemical Exploration[J], 2019, 43(3): 584-588 doi:10.11720/wtyht.2019.1312

0 前言

瞬变电磁测深是一种应用较广泛的时间域电磁探测方法,具有灵敏度高、分辨率强、探测深度大、施工效率高等诸多优点,近年来被广泛应用在金属矿产勘探、煤矿地质灾害调查、隧道超前预报、工程地质勘察等领域[1,2,3]。随着计算机技术与仪器硬件水平的发展,该方法得到了快速发展,但是由于瞬变电磁数值模拟理论的复杂性,二、三维的反演方法还未达到应用程度[2]。大回线源是实际生产中应用最多的瞬变电磁工作装置之一,目前主要的数据处理技术是基于视电阻率的定性解释和一维反演的半定量解释[4,5]。一维反演技术方法较多,在实际生产中应用较广泛的有:共轭梯度法(ICCG)、高斯-牛顿法、阻尼最小二乘法、OCCAM反演法等[6,7,8,9,10,11,12,13,14],它们均为线性化反演方法,优点是收敛速度快、计算量相对较小,但也存在反演过程对初值依赖度高、反演结果容易陷入局部极小、灵敏度矩阵计算量大等问题[9]。针对上述问题,广大地球物理学者开展了大量的研究工作:翁爱华[10]将Occam反演方法应用到中心回线观测装置瞬变电磁测深中,该方法在反演过程中考虑到了实际地层电性纵向上的连续性,减弱了反演结果对初始模型的依赖度,缺点是拉格朗日因子计算量大;张维[8]采用阻尼最小二乘法实现了大定源回线瞬变电磁一维反演,并引入迭代控制上界、下界以及阻尼因子调节的衡量因子实现了阻尼因子的自动调节,提升了反演迭代的速度,缺点是迭代控制上下界人为性强,在实际中不易取舍;徐玉聪,赵宁,秦策等[9]等采用自适应正则化反演法实现了大定源瞬变电磁一维反演,反演过程中正则化因子通过数据目标函数和模型目标函数自适应得到,并采用模型对深度的二阶导数极小的模型约束,使反演收敛速度快、稳定性好,但是计算雅克比矩阵需要花费大量的时间,反演效率不高;戴锐、张达、冀虎[7]采用高斯-牛顿法实现了大定源瞬变电磁法一维反演,也取得了良好的效果,但是该算法对反演初始值依赖度高,需要人为不断的尝试。针对瞬变电磁法反演技术现状,本文根据大回线源瞬变电磁场特点,采用CMD正则化因子优化方案[15]、光滑度函数、自然边界条件和模型修正量的最大值对大回线源瞬变电磁梯度反演过程进行了优化,灵敏度矩阵采用牛顿法[16]进行更新,减少了反演计算的工作量,最后理论模型和实测资料结果表明该方法有效。

1 大回线源瞬变电磁一维反演方法

1.1 一维正演方法

在求解大回线源瞬变电磁响应时,先求解频率域响应,再通过傅里叶变换求解时间域响应。首先采用剖分的思想,把回线源分成若干个电偶极子,然后对所剖分的各个电偶极子产生的频率域响应进行矢量叠加、求和,就可得层状大地大回线源在地表产生的频率域响应,表达式为

Bz(x,y,ω)=∑j=1mPEj2πμ0sinφj∫0∞λ2λ+u1/R1J1(λrj)dλR1=cthu1h1+arcthu1u2cthu2h2+…+artchun-1un,

其中:ui= λ2+ki2, ki2=-iωμσi,σi是第i层地层的电导率,hi为第i层地层的厚度,φj为第j个电偶极子与测点的夹角,rj为测点与第j个电偶极子的距离, PEj为第j个电偶极子的磁距。最后对式(1)进行余旋变换就可求得大回线源地表任意测点的时间域响应。

1.2 一维反演方法

在进行大回线源瞬变电磁法一维反演时,在使观测数据与模型数据的拟合误差达到极小的条件下,同时也要使反演模型的粗糙度达到极小,反演的目标函数为:

U=R+μ-1{‖W·d-W·F[m]‖2-X∞2}。

式中:μ是正则化因子,d是大回线源瞬变电磁框内任意测点的实测数据,F是大回线源瞬变电磁一维正演算子,m为反演模型参数向量, X∞2是反演所要求达到的拟合差,W=diag{1/ε1,1/ε2,…,1/εj,…,1/εM},εj是第j个数据的标准差,是反演模型的粗糙度函数,将其写成矩阵的形式为

R=‖P·m‖2

其中P为粗糙度矩阵,表示为

P=00-11-11……0-11

对公式(2)中的模型m求偏导数,并取 ▽Um=0,就可以求得Um在极小值情况下对应的模型参数修正量Δm的迭代方程

[((WJ)TWJ+μPTP]Δmk=(WJ)TWΔdk。

式中Δmk是第k次反演迭代的模型修改量,Δdk是第k次迭代模型响应与实测数据的残差,J是雅可比矩阵,其元素就为Jij=∂ Fi[m]/∂ mj。

对于公式(4)中的正则化因子选取文献[10]提出的CMD自适应调节方案,其第k次迭代的公式为:

μk=‖PTPmk‖2‖PTPmk‖2+‖WTWΔdk‖2。

在瞬变电磁一维反演过程中,对地层厚度按照等对数间隔进行离散,来反演每一层的电阻率。为了保证反演的电阻率符合自然规律,将每一层的电阻率限制在[0.01,10 000]Ω·m之间,为了不使反演结果快速的陷入局部极小,对模型修正量的最大值限制在50 Ω·m以内。

在反演迭代中每次都必须更新雅克比矩阵J,其计算量是巨大的,在此采用牛顿法进行灵敏度矩阵的更新:

Jk+1=Jk+(Δdk-JkΔmk)ΔmTkΔmTkΔmk。

在反演过程中还要对矩阵方程的条件数进行考虑,如果反演系数矩阵的条件数较大,那么就必须增大正则化因子,以改善其病态性。

2 理论模型计算

为了验证大回线源瞬变电磁一维反演算法的有效性,选取五层的双低阻层和双高阻层模型进行理论模型计算。正演模型的发射外框设置为360 m×360 m,采样时间为0.1~10 ms,等对数间隔40道。反演过程中最小反演深度设置为50 m,最大反演深度设置为600 m,将地层等对数剖分为30层,反演计算这30层大地的电阻率。反演最大迭代次数设置为20次。

双低阻层模型从浅至深各层的电阻率设置为100、10、150、10、100 Ω·m,各层厚度分别为100、30、100、60 m;双高阻层模型从浅至深各层电阻率分别为100、300、80、500、100 Ω·m,厚度为80、100、100、150 m。反演初始模型均设置为100 Ω·m的均匀半空间。双低阻层模型不同反演次数的反演结果及反演模型和理论模型响应的晚期视电阻率拟合结果见图1,双高阻层模型的计算、拟合结果见图2。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   双低阻层模型反演结果(a)及晚期视电阻率拟合结果(b)

Fig.1   Inversion results(a) and fitting diagram of late apparent resistivity(b) of double low resistance layer model


图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   双高阻层模型反演结果(a)及晚期视电阻率拟合结果(b)

Fig.2   Inversion results(a) and fitting diagram of late apparent resistivity(b) of double high resistance layer model


分析图1、图2可见,不管是双低阻层模型还是双高阻层模型,随着反演迭代次数的增加,反演模型结果快速趋近于理论模型;反演模型响应的晚期视电阻率和理论模型响应的晚期视电阻率拟合较好,第20次反演模型结果与理论模型电性变化规律一致,高阻极值与低阻极值的深度与理论模型深度基本一致。但是,综合分析两种模型的反演结果发现:对于低阻层,大回线瞬变电磁一维反演的结果更贴合理论模型;对于高阻层反演电阻率相对于理论模型偏小,这主要是因为瞬变电磁场对低阻体敏感、对高阻体不敏感造成的。总体来说,两种模型的反演结果均能反映低阻层和高阻层的相对电性差异,证明本文提出的大回线瞬变电磁一维反演方法对多层地电模型的反演结果是可靠的。

3 实测数据反演

测区位于山西北部的宁武煤田,地形复杂,地表有基岩出露,大部分被新生界黄土层覆盖,同时沟谷发育,局部冲沟两侧偶有分支沟谷,沟谷底大部分比较开阔。

浅部第四系以亚黏土、亚砂土为主,向下至二叠系上石盒子组及下石盒子组上部,岩性以砂岩、砂泥岩互层为主,呈低阻电性特征;至下石盒子组下部,以粗粒厚层状砂岩、砂质泥岩、砂岩互层为主,呈相对高阻电性特征;再向下至二叠系山西组,岩性为泥岩、粉砂岩,电性相对偏低;继续向下至本区主采煤层石炭系太原组,岩性为灰岩、砂岩、砂泥岩及本区主采煤层,与上述地层电阻率相比,其整体表现为较高阻特征。总体来说,测区地层电性整体呈现为“低阻—高阻—低阻—高阻—低阻”的变化特征。区内小煤窑分布较多,主要采用房柱式开采2号煤层,原生煤层被采动破坏范围较大,且局部位置可能积水,为矿方安全生产埋下了安全隐患,因此,采用地面瞬变电磁法对测区进行了采空区范围勘查。本次施工区域控制面积约0.49 km2,采用发射线框边长为240 m×240 m、发射频率为25 Hz、点距20 m、线距40 m,共布设测线32条。

笔者以28线为例来验证反演算法的有效性。28线工作布置如图3所示,测线长度1 500 m,地形大体呈东北高、西南低趋势,在水平距离600~700 m间分布有相对较深的沟谷,煤层亦表现为东北高、西南低的变化特征。28线反演电阻率断面如图4所示。在反演过程中反演层数剖分为40层,初始电阻率模型均为100 Ω·m,在联想T440P PC机上,76个测点的反演计算时间约为20 min,计算速度较快。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   大回线源瞬变电磁28线工作布置

Fig.3   Large loop source transient electromagnetic 28 line work layout


由图4可见,在探测深度范围内,地层电阻率由浅到深总体上呈“低—高—低—高”的变化特征,与实际地层的电性变化规律一致;横向上反演电阻率等值线走向与含煤地层基本一致,较好地反映了沉积地层的电性分布特征。在2号煤层附近发现1处低阻异常区,其中低阻异常区位于水平距离100~600 m之间,电阻率等值线扭曲,变形明显,呈低阻异常特征,且异常中间有阻断或不连续,推断该低阻异常可能与2号煤层采空破坏区积水有关。后期通过收集该矿的地质资料,发现该区域以前为小煤窑主采区,且采煤历时较长,主要采用房柱式开采2号煤层,采动破坏范围较大;由于该采空区埋深较浅,易于受大气降水和上部含水层中水的补给,所以形成了充水采空区。该反演结果较好地反映了地下采空区的实际位置,证明了该反演算法的有效性。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   大回线瞬变电磁28线反演电阻率断面

Fig.4   Large loop transient electromagnetic 28 line inversion section


4 结论

本文实现了大回线源瞬变电磁法一维自适应反演,通过理论模型和实测资料验证,得出如下结论:

1) 该算法考虑到沉积地层电性纵向的连续性,采用CMD方案、光滑度函数、自然边界条件和模型修正量的最大值对大回线源瞬变电磁的梯度反演过程进行了优化和约束,在保证反演过程自主和稳定的前提下,使反演结果对初值要求不高,同时使反演结果不易陷入局部极小,最后采用牛顿法更新灵敏度矩阵提高了反演效率。

2) 在两层低阻层或者高阻层的反演中,大回线瞬变电磁一维反演算法反演的低阻层结果更贴合实际模型,而对于高阻层反演电阻率相对于实际模型偏小,这主要是由瞬变电磁场对低阻敏感、对高阻不敏感的特性造成的,是瞬变电磁场反演的固有特点。

总体来说,理论模型与实测数据的反演结果均能客观反映地层的相对电性差异,说明该算法对于复杂地电结构具有较好的有效性和适应性,证明本文提出的大回线瞬变电磁一维反演方法是可靠的。

The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。

参考文献

薛国强, 于景邨 .

瞬变电磁法在煤炭领域的研究与应用新进展

[J]. 地球物理学进展, 2017,32(1):319-326.

[本文引用: 1]

Xue G Q, Yu J C .

New development of TEM research and application in coal mine exploration

[J]. Progress in Geophysics, 2017,32(1):0319-0326.

[本文引用: 1]

程建远, 石显新 .

中国煤炭物探技术的现状与发展

[J]. 地球物理学进展, 2013,28(4):2024-2032.

DOI:10.6038/pg20130446      Magsci     [本文引用: 2]

<p>总结了近年来我国煤田物探与矿井物探技术的6大突出进展,即三维地震应用领域不断扩大、地震资料精细处理与动态解释得到推广、地震勘探正在实现从构造勘探向岩性勘探的跨越、地面瞬变电磁技术得到广泛应用、矿井物探技术重新得到重视、以&quot;嵌入式&quot;服务为标志的煤矿安全生产地质保障模式出现了可喜创新等;通过与国内外同类技术发展水平的对比分析,找出了我们存在的差距,提出了加快我国煤炭物探技术发展的对策与建议,并对今后煤炭物探技术与装备的发展前景进行了展望.</p>

Chen J Y, Shi X X .

Current status and development of coal geophysical technology in China

[J]. Progress in Geophys, 2013,28(4):2024-2032.

Magsci     [本文引用: 2]

陈卫营, 薛国强 .

瞬变电磁法多装置探测技术在煤矿采空区调查中的应用

[J]. 地球物理学进展, 2013,28(5):2709-2717.

DOI:10.6038/pg20130554      Magsci     [本文引用: 1]

瞬变电磁法在采空区调查中得到广泛应用. 本文对比了回线源和电性源两种装置在响应特征、探测能力等方面的不同特点,认为回线源瞬变电磁信号对地层灵敏度较高,电性源瞬变电磁探测深度较大.尝试将两种装置联合应用于某变电站选址地层稳定性及周边采空区调查,利用回线源数据可获得浅部地层信息,利用电性源装置探测数据可获得深部地层信息,达到了很好的效果.同时,研究结果认为电性源瞬变电磁法具有探测深度大、施工方便、工作效率高等优点,是一种可用于深部探测的值得推广的电磁方法.

Chen W Y, Xue G Q .

Application on coal-mine voids detection with multi-device TEM technology

[J]. Progress in Geophys, 2013,28(5):2709-2717.

Magsci     [本文引用: 1]

赵越, 李貅, 王祎鹏 .

大回线源瞬变电磁全域视电阻率定义

[J]. 地球物理学进展, 2015,30(4):1856-1863.

[本文引用: 1]

Zhao Y, Li X, Wang Y P .

Full-domain apparent resistivity definition for large-loop TEM

[J]. Progress in Geophys, 2015,30(4):1856-1863.

[本文引用: 1]

薛国强, 李貅, 底青云 .

瞬变电磁法正反演问题研究进展

[J]. 地球物理学进展, 2008,23(4):1165-1172.

DOI:      Magsci     [本文引用: 1]

<FONT face=Verdana>对瞬变电磁法的方法发展概况和仪器研制状况做出了综述性评价.对瞬变电磁法正反演问题的研究成果进行了系统总结.目前的数值模拟正演方法主要有一维滤波系数法，三维积分方程法，二维，三维有限差分法，2.5维有限元法等，主要的反演方法有：一维浮动薄板解释法，人机对话自动反演法，烟圈理论解释法，神经网络反演法，成像类反演等，论述了瞬变电磁法各种计算方法的特点.瞬变电磁法的正反演发展趋势主要是研究三维正反演的计算方法和目标体成像系统.</FONT>

Xue G Q, Li X, Di Q Y .

Research progress in TEM forward modeling and inversion calculation

[J]. Progress in Geophys, 2008,23(4):1165-1172.

Magsci     [本文引用: 1]

刘晓 .

瞬变电磁法三维正演及非线性共轭梯度反演研究

[D]. 北京:中国地质大学( 北京), 2016.

[本文引用: 1]

Li X .

Study of Three-dimensional TEM Forward Modeling and Non-linear Conjugate Gradients Inversion

[D]. Beijing: China University of Geosciences( Beijing), 2016.

[本文引用: 1]

戴锐, 张达, 冀虎 .

大定源瞬变电磁法反演效果分析

[J]. 有色金属:矿山部分, 2017,69(3):1-4.

[本文引用: 2]

Dai R, Zhang D, Ji H .

Analysis of Large loop transient electromagnetic method inversion effect

[J]. Nonferrous Metals:Mining Section, 2017,69(3):1-4.

[本文引用: 2]

张维 .

大定源瞬变电磁法一维正反演研究

[D]. 长沙:中南大学, 2013.

[本文引用: 2]

Zhang W .

The Study about 1D Forward and Inversion of large loop TEM

[D]. Changsha:Central South University, 2013.

[本文引用: 2]

徐玉聪, 赵宁, 秦策 , 等.

大定源瞬变电磁一维自适应正则化反演

[J]. 地质与勘探, 2015,51(2):360-365.

[本文引用: 3]

Xu Y C, Zhao N, Qin C , et al.

One-dimensional adaptive regularization inversion of transient electromagnetic sounding with a large fixed source

[J]. Geology and Exploration, 2015,51(2):360-365.

[本文引用: 3]

翁爱华 .

Occam反演及其在瞬变电磁测深申的应用

[J]. 地质与勘探, 2007,43(5):74-76.

[本文引用: 3]

Weng A H .

OCCAM’s inversion and its application to transient electromagnetic method

[J]. Geology and Prospecting, 2007,43(5):74-76.

[本文引用: 3]

周道卿, 谭捍东, 王卫平 .

频率域航空电磁资料Occam反演研究

[J]. 物探与化探, 2006,30(2):71-74.

[本文引用: 1]

Zhou D Q, Tan H D, Wang W P .

The OCCAM inversion in faem data processing

[J]. Geophysical and Geochemical Exploration, 2006,30(2):71-74.

[本文引用: 1]

李明星, 肖林通, 张倚瑞 , 等.

瞬变电磁粒子群优化反演研究

[J]. 煤炭技术, 2014,33(9):302-304.

[本文引用: 1]

Li M X, Xiao L T, Zhang Y R , et al.

Research on particle swarm optimization inversion of transient electromagnetic method

[J]. Coal Technology, 2014,33(9):302-304.

[本文引用: 1]

Constable S C, Parker R L, Constable C G .

Occam’s inversion: A practical algorithm for generating smooth model from electromagnetic sounding data

[J]. Geophysics, 1987,52(2):289.

DOI:10.1190/1.1442303      URL     [本文引用: 1]

毛立峰, 王绪本, 陈斌 .

直升机航空瞬变电磁自适应正则化一维反演方法研究

[J]. 地球物理学进展, 2011,26(1):300-305.

DOI:10.3969/j.issn.1004-2903.2011.01.035      Magsci     [本文引用: 1]

本文以中心回线式直升机航空瞬变电磁法为例对时间域航空电磁资料的一维反演算法进行理论研究.由于时间域航空电磁反演问题是典型的不适定问题,多解性严重且稳定性能较差.为此,论文使用Tikhonov正则化反演方法进行一维解释,可以得到用最光滑模型、最平坦模型和最小构造模型约束条件下的反演结果,而正则化因子是根据各次迭代的数据目标函数与约束模型目标函数的关系,以自适应的方式给出,避免了反复地搜索正则因子,使反演能够稳定快速进行.论文介绍了自适应正则化反演方法之后,用多层理论模型的垂直感生电动势响应数据进行反演试算,对比分析不同约束模型下的反演结果.结果表明,三种模型约束下,算法对初始模型要求不高,可以使用均匀半空间模型作为反演初始模型,并且均能够在不超过十次的反演迭代计算就能较好地重构理论模型,体现了良好的稳定性和可靠性.

Mao L F, Wang X B, Chen B .

Study on an adaptive regularized 1D inversion method of helicopter TEM data

[J]. Progress in geophys., 2011,26(1):300-305.

Magsci     [本文引用: 1]

陈小斌, 赵国泽, 汤吉 , 等.

大地电磁自适应正则化反演算法

[J]. 地球物理学报, 2005,48(4):937-946.

Magsci     [本文引用: 1]

针对大地电磁正则化反演中正则化因子的选取困难问题提出了自适应正则化反演算法（Adaptive Regularized Inversion Algorithm, ARIA）. 在该算法中， ①提出了一种新的数据方差处理方法：数据方差规范化，使得数据方差的大小只对数据的拟合发生影响，不对数据目标函数和模型约束目标函数的权重产生影响，从而减少了正则化因子取值的影响因素；②提出了粗糙度核矩阵的概念，并给出了由基本结构插值基函数计算粗糙度核矩阵的公式，使得模型目标函数的构建更为简便、直接；③根据数据目标函数、模型约束目标函数和正则化因子之间的关系，提出了两种正则化因子自适应调节方法. 本文详细阐述了最平缓模型约束下的大地电磁一维连续介质反演的ARIA实现，以几个算例的分析比较来说明ARIA的有效性.

Chen X B, Zhao G Z, Tang J , et al.

An adaptive regularized inversion algorithm for magnetotelluric data. Chinese

[J]. Geophys., 2005,48(4):937-946.

Magsci     [本文引用: 1]

李刚, 潘和平, 王智 , 等.

回线源瞬变电磁法一维反演算法

[J]. 煤田地质与勘探, 2017,45(5):161-166.

[本文引用: 1]

Li G, Pan H P, Wang Z , et al.

One-dimensional inversion for loop source transient electromagnetic method

[J]. Coal Geology & Exploration, 2017,45(5):161-166.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com