物探与化探
物探與化探
물탐여화탐
Geophysical and Geochemical Exploration
2015年
5期
1059-1062
,共4页
航空瞬变电磁%2.5维正演%有限单元法%异常场法%响应特征
航空瞬變電磁%2.5維正縯%有限單元法%異常場法%響應特徵
항공순변전자%2.5유정연%유한단원법%이상장법%향응특정
ATEM%2.5-D forward modeling%finite element method%anomalous field algorithm%response characteristic
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。
採用三角網格剖分的有限元法,研究瞭2.5維航空瞬變電磁法正縯模擬問題。利用時頻變換數值方法將時間域電磁場轉換到拉氏域,再利用傅裏葉變換將三維問題降維變為2.5維問題,然後由有限元法求解得到拉氏域二維電磁場,逆拉氏變換後得到時間域航空瞬變響應。為瞭迴避正縯模擬中總感應磁場在場源處的奇異性問題,採用異常場算法,場源響應通過在微分方程中施加揹景電磁場實現。由于瞬變電磁信號具有較大的動態範圍,而且需要經過兩次正、逆拉氏變換和傅裏葉變換,每箇環節的計算精度和速度要嚴格控製在較高的水平上,否則積纍誤差會非常大。模型計算錶明均勻大地和層狀大地模型解析解與數值解吻閤很好。這證明該算法是正確可行的,可作為研究二維複雜地質體的方法手段。
채용삼각망격부분적유한원법,연구료2.5유항공순변전자법정연모의문제。이용시빈변환수치방법장시간역전자장전환도랍씨역,재이용부리협변환장삼유문제강유변위2.5유문제,연후유유한원법구해득도랍씨역이유전자장,역랍씨변환후득도시간역항공순변향응。위료회피정연모의중총감응자장재장원처적기이성문제,채용이상장산법,장원향응통과재미분방정중시가배경전자장실현。유우순변전자신호구유교대적동태범위,이차수요경과량차정、역랍씨변환화부리협변환,매개배절적계산정도화속도요엄격공제재교고적수평상,부칙적루오차회비상대。모형계산표명균균대지화층상대지모형해석해여수치해문합흔호。저증명해산법시정학가행적,가작위연구이유복잡지질체적방법수단。
Based on the triangular mesh dissection of the finite element method, this paper presents the 2.5?D airborne transient elec?tromagnetic method forward modeling. Firstly the time domain electromagnetic partial equations are converted into laplace domain with the numerical time?frequency transform algorithm, then the 2.5?D expressions are obtained from 3?D case by adopting the fourier trans?form. The airborne data are acquired from the results of the digital inverse laplace transform of the 2?D electric and magnetic compo?nents in laplacian domain, which is computed using finite element method. Instead of the conventional algorithm, the anomalous field method is employed throughout the finite element process, so that the source singularity of total field is avoided. the source impact is embodied in the differencial equation by appling the background electromagnetic field item. The computing accuracy and efficiency would be managed rigorously due to the wide dynamic range of the transient singnal, as well as the two?time inverse laplace and fourier transforms. Otherwise the accumulative error will soar acceptably. The synthetic model experiments shows that the numeric modeling so?lutions match the analytical results of homogeneous and layered erath responses well, which also proves the effectiveness of the algo?rithm.