本文關(guān)鍵詞： 非結(jié)構(gòu)網(wǎng)格 電阻率 三維正演 有限單元法 巷道超前探測 蒙特卡羅方法 三維反演 地形　出處：《中國科學(xué)技術(shù)大學(xué)》2016年博士論文　論文類型：學(xué)位論文

【摘要】：地表起伏地形在野外礦產(chǎn)資源勘察中不可避免,其對直流電法勘探影響巨大。我們通過在電阻率三維正演中采用非結(jié)構(gòu)網(wǎng)格,能夠進(jìn)行任意復(fù)雜地形和幾何模型的數(shù)值模擬。與結(jié)構(gòu)化網(wǎng)格相比,非結(jié)構(gòu)網(wǎng)格具有可局部加密的特點(diǎn),只需要在重點(diǎn)觀測計(jì)算區(qū)域內(nèi)添加有限數(shù)量的節(jié)點(diǎn),即可充分發(fā)揮有限單元法的優(yōu)勢,在減少計(jì)算量的同時(shí),大幅度提高求解精度。球體和無限大板體模型的電阻率三維非結(jié)構(gòu)有限單元數(shù)值模擬結(jié)果表明,誤差均在0.9%以內(nèi)。不同起伏地形模型的電阻率三維非結(jié)構(gòu)有限單元數(shù)值模擬結(jié)果表明,地形影響巨大,甚至掩蓋地下真實(shí)異常體的響應(yīng)。快速、精確的電阻率三維非結(jié)構(gòu)有限單元數(shù)值模擬為其在巷道超前探測中的應(yīng)用及帶地形電阻率三維反演奠定堅(jiān)實(shí)基礎(chǔ)。直流電阻率法通過觀測單極-偶極裝置(三極裝置)視電阻率曲線最小值位置估計(jì)掌子面前方低阻異常體距離,在巷道超前探測及礦井水害預(yù)測與防治工作中有廣泛應(yīng)用。由于地下巷道本身即是三維模型,必須應(yīng)用電阻率三維數(shù)值模擬才能獲得地下復(fù)雜結(jié)構(gòu)的全空間三維地電場的分布特征。我們基于電阻率三維非結(jié)構(gòu)有限單元數(shù)值模擬,給出低阻異常體超前探測的線性預(yù)測模型。然而,實(shí)際井下地質(zhì)結(jié)構(gòu)復(fù)雜,異常大小不定且任意分布,而且可能存在有電阻率各向異性,使得目前由簡單模型實(shí)驗(yàn)或數(shù)值模擬獲得的預(yù)測公式不確定性較大,而且可靠性難于評價(jià)。蒙特卡羅方法用隨機(jī)化思想解決非確定性問題,我們將該方法與并行算法結(jié)合,對大量不確定的隨機(jī)模型進(jìn)行三維數(shù)值模擬,以驗(yàn)證各預(yù)測模型的準(zhǔn)確度及可靠性。結(jié)果表明,本文提出的超前探線性預(yù)測模型較其他預(yù)測模型更為準(zhǔn)確、可靠。同時(shí),模擬結(jié)果表明介質(zhì)各向異性的影響很大,尤其當(dāng)圍巖電阻率具有各向異性時(shí),所有預(yù)測模型的準(zhǔn)確度及可靠性較差,給安全生產(chǎn)帶來很大隱患。我們首次將并行蒙特卡羅方法應(yīng)用于直流電阻率巷道超前探測研究,為巷道超前探測準(zhǔn)確性和可靠性的定量化描述提供了新思路。現(xiàn)階段,面向?qū)嶋H應(yīng)用的起伏地形下電阻率三維反演依然困難。考慮到應(yīng)用GPS/GNSS定位技術(shù)時(shí),區(qū)域地球物理調(diào)查中可非規(guī)則布設(shè)測網(wǎng)的實(shí)際特點(diǎn),我們在基于非結(jié)構(gòu)網(wǎng)格的電阻率三維正演數(shù)值模擬基礎(chǔ)上,首次實(shí)現(xiàn)了任意偶極-偶極觀測數(shù)據(jù)的電阻率、激發(fā)極化三維帶地形反演。任意偶極-偶極觀測方式有效實(shí)現(xiàn)大規(guī)模、隨機(jī)化、高密度的三維觀測數(shù)據(jù)采集,同時(shí),非結(jié)構(gòu)網(wǎng)格為地形建模采集地表高程數(shù)據(jù)及其靈活、高效的野外電阻率三維數(shù)據(jù)解釋創(chuàng)造了有利條件。理論模型合成數(shù)據(jù)的帶地形電阻率、激發(fā)極化三維反演結(jié)果表明,無論常規(guī)的中間梯度數(shù)據(jù)、還是任意偶極-偶極數(shù)據(jù)均取得良好反演效果,起伏地形下的不同異常體電性反映清晰、定位準(zhǔn)確。通過不斷調(diào)整合成觀測數(shù)據(jù)中高斯噪聲大小,展示了本文的反演方法對不同質(zhì)量觀測數(shù)據(jù)的反演效果,顯示了對噪聲的良好兼容性,是電阻率三維反演應(yīng)用于實(shí)際勘探的關(guān)鍵。我們將本文的反演方法應(yīng)用于甘肅某地的多金屬礦三維勘察,對實(shí)際觀測的電阻率三維數(shù)據(jù)進(jìn)行反演解釋,給出了礦區(qū)的地下電性三維分布,其中一個(gè)低阻異常體與正在開采的金礦體吻合,其余三個(gè)低阻異常推斷為隱伏礦體,以待其它勘探手段查證。電阻率三維反演解釋成功應(yīng)用于實(shí)際礦產(chǎn)資源勘探的結(jié)果國內(nèi)、外尚不多見,是實(shí)際電法勘探三維反演解釋向前推進(jìn)的重要一步。
[Abstract]:The surface topography in the field of mineral resources exploration is inevitable, the electrical prospecting has great impact. We passed in 3-D resistivity forward using unstructured grid, numerical simulation can be arbitrary complex terrain and geometric model. Compared with the structured meshes, unstructured grid has the characteristics of partial encryption, only need to calculate the node area in addition a limited number of observations in the key, can give full play to the advantages of the finite element method, to reduce the amount of computation at the same time, greatly improve the accuracy of solution. The sphere and the infinite plate model of the electrical resistivity of three-dimensional unstructured finite element numerical simulation results show that the errors are within 0.9%. The resistivity of three-dimensional unstructured finite element numerical different terrain model simulation results show that the terrain impact response and even obscure the real underground abnormal body. Fast, accurate resistivity Numerical simulation of three-dimensional unstructured finite element for its application in tunnel in advance in the detection of terrain and 3D resistivity inversion and lay a solid foundation. The DC resistivity observation by pole dipole device (transistor device) apparent resistivity curve minimum position estimation in front of the tunnel face low resistivity body distance in tunnel in advance are widely used in detection and prediction mine water disaster and prevention work. Because the underground roadway itself is a three-dimensional model must be applied to get 3D resistivity 3-D numerical simulation to the whole space of underground complex structure of electric field distribution. Our numerical simulation of 3-D resistivity structure based on linear non finite element, given the low resistance abnormal body advanced detection model. However, the actual underground the geological structure is complex, the abnormal size is uncertain and arbitrary distribution, but there may be a resistivity anisotropy, making the current by Jane Single model test or numerical simulation prediction formula for the uncertainty and reliability evaluation. The Monte Carlo method to solving the uncertain problems with random thoughts, we will combine with the method of parallel algorithm, a large number of uncertain stochastic model of three-dimensional numerical simulation, to verify the accuracy and reliability of the prediction model. The results show that in this paper, the advance on linear prediction model than other prediction model is more accurate and reliable. At the same time, the simulation results show that the influence of anisotropy, especially when the rock resistivity is anisotropic, accuracy and poor reliability of all forecasting models, bring great hidden trouble to the safe production. We will be the first application of parallel Monte Carlo method in the DC resistivity in advanced detection research, provides a quantitative detection for the roadway ahead of the accuracy and reliability of New ideas. At this stage, 3-D resistivity inversion topography and the practical application is still difficult. Considering the application of GPS/GNSS positioning technology, the actual characteristics of irregular measuring network layout can be the regional geophysical survey, we in resistivity based on unstructured grids of 3D numerical modeling based on first realized arbitrary dipole resistivity dipole induced polarization observation data, three-dimensional topography inversion. Arbitrary dipole dipole observation way of effective implementation of large-scale, randomized, 3-D high density data acquisition, at the same time, the unstructured grid for terrain modeling and data acquisition of surface elevation and flexible 3D field resistivity data interpretation, and create favorable conditions. The theoretical model of synthetic data the terrain with 3D resistivity, induced polarization inversion results show that both the conventional intermediate gradient data, or arbitrary dipole dipole number According to the good inversion effect, undulating terrain under different electric anomalies reflect clear, accurate positioning. By continuously adjusting the Gauss noise synthesis data size, showing the effect of the inversion inversion method for different quality observation data, showing good compatibility to the noise, the key is applied to 3-D resistivity inversion we will survey and exploration. The 3D multi metal mine inversion method in this paper was used to in Gansu, inversion interpretation of 3D resistivity data of actual observations, given the mine underground electric distribution in which a low resistivity body and gold mining is consistent, the remaining three low resistivity anomaly for for concealed orebodies, other exploration methods verified. 3D resistivity inversion applied to explain the results of domestic mineral resources exploration, is rare, is the actual power The three dimensional inversion of the method of exploration is an important step forward.

【學(xué)位授予單位】：中國科學(xué)技術(shù)大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2016
【分類號】：P631.322

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 黃俊革;阮百堯;王家林;;坑道內(nèi)赤道偶極電阻率測深的可行性[J];中南大學(xué)學(xué)報(bào)(自然科學(xué)版);2006年04期

2 陳仲候 ,何昌禮;電阻率測深的數(shù)字解釋[J];物探與化探;1980年02期

3 李光明;電阻率測深曲線接頭開口的引起原因及處理方法的探討[J];物探與化探;1987年03期

4 劉煜洲;曹靜平;陳福集;;微分電阻率測深原理及應(yīng)用[J];國外地質(zhì)勘探技術(shù);1991年10期

5 劉象勤,李純志;二維直立巖層電阻率測深異常曲線的某些特征[J];礦產(chǎn)與地質(zhì);1991年03期

6 莫承彬，覃佑邦，劉康華;二階差分法在電阻率測深中的應(yīng)用[J];物探與化探;1995年01期

7 阮百堯，，徐世浙;二維直流電阻率測深曲線的快速反演[J];物探與化探;1996年06期

8 阮百堯,熊彬;電導(dǎo)率連續(xù)變化的三維電阻率測深有限元模擬[J];地球物理學(xué)報(bào);2002年01期

9 溫永輝,何門貴;電阻率測深的數(shù)值差分法在工程上的應(yīng)用[J];工程勘察;2002年02期

10 黃俊革,阮百堯;直流電阻率測深中二維與三維反演結(jié)果的對比與分析[J];物探與化探;2004年05期

相關(guān)會議論文 前9條

1 阮百堯;徐世浙;;二維直流電阻率測深曲線的快速反演[A];1995年中國地球物理學(xué)會第十一屆學(xué)術(shù)年會論文集[C];1995年

2 黃俊革;阮百堯;;坑道內(nèi)金屬管道對電阻率測深影響的數(shù)值模擬[A];地球科學(xué)與資源環(huán)境——華南青年地學(xué)學(xué)術(shù)研討會論文集[C];2003年

3 黃俊革;阮百堯;;坑道內(nèi)金屬管道對電阻率測深影響的數(shù)值模擬[A];中國地球物理.2003——中國地球物理學(xué)會第十九屆年會論文集[C];2003年

4 申時(shí)新;雷揚(yáng)敬;譚昌勇;;利用電阻率測深劃分咸淡水[A];1997年中國地球物理學(xué)會第十三屆學(xué)術(shù)年會論文集[C];1997年

5 張致付;程志平;阮百堯;李澤林;;Zohdy方法在三維電阻率測深數(shù)據(jù)反演中的應(yīng)用[A];資源、環(huán)境與工程電法勘探方法技術(shù)交流會論文集[C];2003年

6 唐君輝;劉國興;;電阻率測深在漠河地區(qū)探測永久凍土帶分布的效果[A];中國地球物理第二十一屆年會論文集[C];2005年

7 張致付;程志平;阮百堯;;三維直流電阻率測深曲線的近似反演[A];中國地球物理學(xué)會年刊2002——中國地球物理學(xué)會第十八屆年會論文集[C];2002年

8 楊長龍;張林;何莉;;視電阻率測深在鋁土礦勘查中的應(yīng)用[A];河南省地質(zhì)調(diào)查與研究通報(bào)2007年卷（上冊）[C];2007年

9 劉斌;程久龍;徐世浙;;電阻率測深二維有限元模擬在煤田勘探中的應(yīng)用[A];1996年中國地球物理學(xué)會第十二屆學(xué)術(shù)年會論文集[C];1996年

相關(guān)博士學(xué)位論文 前3條

1 劉洋;基于非結(jié)構(gòu)網(wǎng)格的電阻率三維正反演及其應(yīng)用研究[D];中國科學(xué)技術(shù)大學(xué);2016年

2 趙廣茂;帶地形的復(fù)電阻率2.5維電磁場正反演研究[D];吉林大學(xué);2009年

3 徐凱軍;2.5維復(fù)電阻率電磁場正反演研究[D];吉林大學(xué);2007年

相關(guān)碩士學(xué)位論文 前1條

1 周峰;復(fù)電阻率1D/2.5D正反演研究[D];東華理工大學(xué);2015年

本文編號：1465187