Abstract：Multiplicity of solutions to inversion is a inherent problem of 3D electrical resistivity detection，which always causes errors and even mistakes in the geological interpretation. A research thought is presented to solve the problem，in which the structural and morphology information of anomalous body gained by other geophysical methods，e.g. seismic prospecting，ground penetrating radar，is used as prior constraints in resistivity inversion. First，construction method of a novel constraint named spatial structural constraint is proposed. The known anomalous regions obtained with other geophysical methods are mapped into 3D inversion model；and the difference between grids in corresponding regions is restricted to minimum. Thus，spatial structural constraint matrix is established to characterize typical geological structures，especially the more complicated structures；and it is relatively easy to be expressed mathematically and constructed. Based on above research，the objective function and inversion tomography equation of 3D resistivity inversion with spatial structural constraint are constructed. Then，numerical tests and engineering application are implemented；and the spatial structural constraint gained by ground penetrating radar(GPR) is applied to 3D resistivity inversion. It is found that the prior information increases greatly and false anomaly in inversion tomography is eliminated effectively. And the location precision and interface identification-effect are improved significantly compared with traditional method，making the inversion result be consistent with geological model or actual situation and the multiplicity of solutions be pressed obviously. And it is proved that the 3D resistivity inversion tomography with spatial structural constraint is a feasible way to reduce multiplicity of solutions and improve location precision and interface identification-effect.
底青云，王妙月. 积分法三维电阻率成像[J]. 地球物理学报，2001，44(6)：843-852.(DI Qingyun，WANG Miaoyue. 3D resistivity tomography by integral method[J]. Chinese Journal of Geophysics，2001，44(6)：843-852.(in Chinese))
吴小平，徐果明. 利用共轭度法的电阻率三维反演研究[J]. 地球物理学报，2000，43(3)：420-426.(WU Xiaoping，XU Guoming. Study on 3D resistivity inversion using conjugate gradient method[J]. Chinese Journal of Geophysics，2000，43(3)：420-426.(in Chinese))
[12]
MUSIL M，MAURER H R，GREEN A G. Discrete tomography and joint inversion for loosely connected or unconnected physical properties：application to crosshole seismic and georadar data sets[J]. Geophysical Journal International，2003，153(2)：389-402.
<body>Your browse does not support frame!</body> " target="_blank">
刘 斌，李术才，聂利超，等. 基于自适应加权光滑约束与PCG算法的三维电阻率探测反演成像[J]. 岩土工程学报，2012，34(9)：1 646-1 653.(LIU Bin，LI Shucai，NIE Lichao，et al . Study on inversion method of 3D resistivity detection using adaptive-weighted smooth constraint and PCG algorithm[J]. Chinese Journal of Geotechnical Engineering，2012，34(9)：1 646-1 653.(in Chinese))
<body>Your browse does not support frame!</body> " target="_blank">
宛新林，席道瑛，高尔根，等. 用改进的光滑约束最小二乘正交分解法实现电阻率三维反演[J]. 地球物理学报，2005，48(1)：439-444.(WAN Xinlin，XI Daoying，GAO Ergen，et al. 3-D resistivity inversion by the least-squares QR factorization method under improved smoothness constraint condition[J]. Chinese Journal of Geophysics，2005，48(1)：439-444.(in Chinese))
<body>Your browse does not support frame!</body> " target="_blank">
[13]
刘 斌，李术才，李树忱，等. 基于不等式约束的最小二乘法三维电阻率反演及其算法优化[J]. 地球物理学报，2012，55(1)：260-268.(LIU Bin，LI Shucai，LI Shuchen，et al. 3D electrical resistivity inversion with least-squares method based on inequality constraint and its computation efficiency optimization[J]. Chinese Journal of Geophysics，2012，55(1)：260-268.(in Chinese))
阮百尧，村上峪，徐世浙. 电阻率/激发极化率数据的二维反演程序[J]. 物探化探计算技术，1999，21(2)：116-125.(RUAN Baiyao，MURAKAMI Y，XU Shizhe. 2D inversion program of induced polarization data[J]. Computing Techniques for Geophysical and Geochemical Exploration，1999，21(2)：116-125.(in Chinese))
<body>Your browse does not support frame!</body> " target="_blank">
2012, Vol. 31 
