三维地形频率域井筒电磁场区域积分方程法模拟
|
摘要
井筒电磁法作为一种高效的地球物理勘探技术特别适合我国地形复杂地区(沙漠、高山等)的油气资源勘探.地形起伏区域对井筒电磁响应的观测具有严重影响,但到目前为止人们对三维井筒电磁地形效应特征的研究十分有限.本文提出基于区域划分的积分方程法模拟带地形频率域井筒电磁系统响应,与基于偏微分方程的有限差分、有限单元法相比,该方法能更高效地模拟地形响应.首先根据地形起伏情况定义感应数,将地形条件下目标体的井筒电磁场模拟区域划分为参考模型、背景介质及目标体介质分布子区域,针对各子区域的模拟计算特点,配置Anderson算法、稳定型双共轭梯度-快速傅里叶积分方程算法,从而获得三维地形频率域井筒电磁场响应.通过将计算结果与半空间模型的Anderson算法解析解、带山谷地形模型的其他已发表的三维边界积分方程结果进行对比,检验了本文算法的精度及高效性.最后,系统分析了山谷地形对井筒电磁地井观测系统电磁场响应的影响特征.本文研究结果对三维井筒电磁地形效应的识别和校正具有指导意义.
As an efficient geophysical exploration technology, well electromagnetic method is particularly applicable to oil and gas exploration in China's complex terrain areas(deserts, mountains, etc.). A serious influence of topographic relief area on the electromagnetic response of well is inevitable but challenging. To the best of our knowledge, there is no literature on modeling the electromagnetic response of three-dimensional(3D)topography with well electromagnetic method. Based on the domain decomposition, an integral equation method is presented to simulate the electromagnetic response of 3D topography in frequency domain via the well electromagnetic method. Compared with the finite difference and finite element method based on partial differential equation, this method is very efficient in simulating topographic response without huge computation or truncation boundary error accumulation or special boundary condition requirements. Firstly, an induction coefficient is defined according to the topographic relief situation. Then the computational domain consisting of the target body, background medium and 3D topography is divided into reference model, background medium and the distribution of target body medium area. According to the characteristics of each sub-region, Anderson algorithm is an analytic solution based on Gaussian filtering, which is used to provide the primary field from the excited sources in surface. And then, the stable double conjugate gradient-fast Fourier transform is incorporated into integral equation algorithm to obtain the fast 3D terrain shaft frequency domain electromagnetic responses. By comparing the calculation results using the new algorithm presented in this paper with the analytical solutions of Anderson algorithm for half-space model with surface electromagnetic method,the precision and the efficiency of this new algorithm are demonstrated. And the ability to model the electromagnetic responses of 3D topography is shown by comparing with the published results of 3D boundary integral equation. Thus, the high accuracy and high efficiency of the new algorithm presented in this paper are validated. Finally, the influence of 3D valley topography on electromagnetic field response of surface to borehole electromagnetic(SBEM) observation system is presented and analyzed. It is observed that the response of SBEM is seriously disturbed by the field of 3D valley topography which is necessarily removed. The research results presented in this paper are of significance for guiding the identification and correction of electromagnetic topographic effect from 3D SBEM.
As an efficient geophysical exploration technology, well electromagnetic method is particularly applicable to oil and gas exploration in China's complex terrain areas(deserts, mountains, etc.). A serious influence of topographic relief area on the electromagnetic response of well is inevitable but challenging. To the best of our knowledge, there is no literature on modeling the electromagnetic response of three-dimensional(3D)topography with well electromagnetic method. Based on the domain decomposition, an integral equation method is presented to simulate the electromagnetic response of 3D topography in frequency domain via the well electromagnetic method. Compared with the finite difference and finite element method based on partial differential equation, this method is very efficient in simulating topographic response without huge computation or truncation boundary error accumulation or special boundary condition requirements. Firstly, an induction coefficient is defined according to the topographic relief situation. Then the computational domain consisting of the target body, background medium and 3D topography is divided into reference model, background medium and the distribution of target body medium area. According to the characteristics of each sub-region, Anderson algorithm is an analytic solution based on Gaussian filtering, which is used to provide the primary field from the excited sources in surface. And then, the stable double conjugate gradient-fast Fourier transform is incorporated into integral equation algorithm to obtain the fast 3D terrain shaft frequency domain electromagnetic responses. By comparing the calculation results using the new algorithm presented in this paper with the analytical solutions of Anderson algorithm for half-space model with surface electromagnetic method,the precision and the efficiency of this new algorithm are demonstrated. And the ability to model the electromagnetic responses of 3D topography is shown by comparing with the published results of 3D boundary integral equation. Thus, the high accuracy and high efficiency of the new algorithm presented in this paper are validated. Finally, the influence of 3D valley topography on electromagnetic field response of surface to borehole electromagnetic(SBEM) observation system is presented and analyzed. It is observed that the response of SBEM is seriously disturbed by the field of 3D valley topography which is necessarily removed. The research results presented in this paper are of significance for guiding the identification and correction of electromagnetic topographic effect from 3D SBEM.
引文
[1]Augustin A,Kennedy W,Morrison H 1989 Geophysics 54 90
[2]Li J H,He Z X,Lii Y Z 2011 Chin.J.Eng.Geophys.8 303(in Chinese)[李静和,何展翔,吕玉增2011工程地球物理学报8 303]
[3]Jahandari H,Ansari S,Farquharson C 2017 J.Appl.Geophys.138 185
[4]Tang W W,Liu J X,Ye Y X 2018 Oil Geophys.Prosp.53617(in Chinese)[汤文武,柳建新,叶益信2018石油地球物理勘探53 617]
[5]Nornikman H,Pee N,Ahmad B 2018 JTEC10 35
[6]Zhang Y,Lin L,Chen G B 2018 Chinese J.Geophys.61 1639(in Chinese)[张烨,林蔺,陈桂波2018地球物理学报61 1639]
[7]Wei B J,Chen T,Hou X L 2014 J.China Uni.Petrol. 38 57(in Chinese)[魏宝君,陈涛,侯学理2014中国石油大学学报38 57]
[8]Guo C B,Xiao C H,Liu D M 2008 Acta Phys.Sin.57 4182(in Chinese)[郭成豹,肖昌汉,刘大明2008物理学报57 4182]
[9]Chen G B,Wang H N,Yao J J,Han Z Y 2009 Acta Phys.Sin.58 3848(in Chinese)[陈桂波,汪宏年,姚敬金,韩子夜2009物理学报58 3848]
[10]Li J H,He Z X 2012 Oil Geophys.Prosp.47 653(in Chinese)[李静和,何展翔2012石油地球物理勘探47 653]
[11]Li J H,He Z X,Xu Y X 2017 Appl. Geophys.14 559
[12]Kruglyakov M,Kuvshinov A 2018 Geophys.J.Int.213 1387
[13]Tiberi G,Monorchio A,Manara G,Mittra R A 2006 IEEE T.Anteen.Propag.54 2508
[14]Nie X C,Yuan N,Liu C R 2010 IEEE Geosci.Remote.S.4872
[15]Chobanyan E,Notaros B M,Ilic M M 2014 APSURSI213 4
[16]Chen G B,Bi J,Wang J B,Chen X Y,Sun G C,Lu J 2011Acta Phys.Sin.60 094102(in Chinese)[陈桂波,毕娟,汪剑波,陈新邑,孙贯成,卢俊2011物理学报60 094102]
[17]Zhdanov M S,Wan L,Gribenko A,Martin C,Key K,Constable S 2011 Geophysics 76
[18]Yin C C,Zhang B,Liu Y H,Cai J 2015 Chine.se J.Geophys.58 1411(in Chinese)[殷长春,张博,刘云鹤,蔡晶2015地球物理学报58 1411]
[19]Li X J,Lei L,Chen Y P,Jiang M,Rong Z,Hu J 2019 The Chinese J of Rad.Sci.11(in Chinese)[李先进,雷霖,陈涌频,江明,荣志,胡俊2019电波科学学报1 1]
[20]Li J H,Song L P,Liu Q H 2016 Pure Appl.Geophys.173 607
[21]Anderson W L 1979 Geophysics 44 1287
[22]Ruan B Y,Wang Y X 2005 Chinese J.Geophys.48 1197(in Chinese)[阮百尧,王有学2005地球物理学报48 1197]
[2]Li J H,He Z X,Lii Y Z 2011 Chin.J.Eng.Geophys.8 303(in Chinese)[李静和,何展翔,吕玉增2011工程地球物理学报8 303]
[3]Jahandari H,Ansari S,Farquharson C 2017 J.Appl.Geophys.138 185
[4]Tang W W,Liu J X,Ye Y X 2018 Oil Geophys.Prosp.53617(in Chinese)[汤文武,柳建新,叶益信2018石油地球物理勘探53 617]
[5]Nornikman H,Pee N,Ahmad B 2018 JTEC10 35
[6]Zhang Y,Lin L,Chen G B 2018 Chinese J.Geophys.61 1639(in Chinese)[张烨,林蔺,陈桂波2018地球物理学报61 1639]
[7]Wei B J,Chen T,Hou X L 2014 J.China Uni.Petrol. 38 57(in Chinese)[魏宝君,陈涛,侯学理2014中国石油大学学报38 57]
[8]Guo C B,Xiao C H,Liu D M 2008 Acta Phys.Sin.57 4182(in Chinese)[郭成豹,肖昌汉,刘大明2008物理学报57 4182]
[9]Chen G B,Wang H N,Yao J J,Han Z Y 2009 Acta Phys.Sin.58 3848(in Chinese)[陈桂波,汪宏年,姚敬金,韩子夜2009物理学报58 3848]
[10]Li J H,He Z X 2012 Oil Geophys.Prosp.47 653(in Chinese)[李静和,何展翔2012石油地球物理勘探47 653]
[11]Li J H,He Z X,Xu Y X 2017 Appl. Geophys.14 559
[12]Kruglyakov M,Kuvshinov A 2018 Geophys.J.Int.213 1387
[13]Tiberi G,Monorchio A,Manara G,Mittra R A 2006 IEEE T.Anteen.Propag.54 2508
[14]Nie X C,Yuan N,Liu C R 2010 IEEE Geosci.Remote.S.4872
[15]Chobanyan E,Notaros B M,Ilic M M 2014 APSURSI213 4
[16]Chen G B,Bi J,Wang J B,Chen X Y,Sun G C,Lu J 2011Acta Phys.Sin.60 094102(in Chinese)[陈桂波,毕娟,汪剑波,陈新邑,孙贯成,卢俊2011物理学报60 094102]
[17]Zhdanov M S,Wan L,Gribenko A,Martin C,Key K,Constable S 2011 Geophysics 76
[18]Yin C C,Zhang B,Liu Y H,Cai J 2015 Chine.se J.Geophys.58 1411(in Chinese)[殷长春,张博,刘云鹤,蔡晶2015地球物理学报58 1411]
[19]Li X J,Lei L,Chen Y P,Jiang M,Rong Z,Hu J 2019 The Chinese J of Rad.Sci.11(in Chinese)[李先进,雷霖,陈涌频,江明,荣志,胡俊2019电波科学学报1 1]
[20]Li J H,Song L P,Liu Q H 2016 Pure Appl.Geophys.173 607
[21]Anderson W L 1979 Geophysics 44 1287
[22]Ruan B Y,Wang Y X 2005 Chinese J.Geophys.48 1197(in Chinese)[阮百尧,王有学2005地球物理学报48 1197]