多阵列感应测井+毕设文献翻译

1、复杂地底环境下任意指向的电磁线圈的磁场摘要：本文提出了一种基于有限差分的仿真复杂地质环境下任意方向的线圈产生的磁场的方法。地底传播媒介在水平和垂直两个方向上都有多个层次，并且曾现明显的各向异性。本改进的有限差分方法用一个交错的网格近似得到分散的电磁场中的向量方程。然后用广义最小裕量法迭代地求解得到的线性稀疏矩阵，同时用不完全逻辑单元预处理方法提升此线性方程的收敛性，加速求解过程。此算法已经被实验多次证实，并且应用于常见的三维感应工具中，用于电测井中各向异性的测量。1、 简介复杂地形环境下线圈产生的磁场的研究对于地球物理勘探而言十分重要。众所周知，油气储层相比起非储层岩石来有更高的电阻率(或更低

2、的电导率)。因此，从地层的电阻率特征中，地质学家可以判别岩石中的油气储层及碳氢化合物的含量。地层的电阻率特征通常通过地层对电测井工具的反应磁场获得。近年来，厚的油气储层逐渐被开采完毕，寻找和开采浅的油气储层变得越来越重要。浅的砂页岩储层通常表现出电气特性上的各向异性，突出的特点是水平方向和垂直方向上的电阻率各向异性。地质环境电气各向异性的检测是一个地质学家们研究了近70年的问题。最近开发的三维感应工具就能够检测地层媒介的各向异性。在三维感应测井工具中发射端和接收端是两个指向相互垂直的小感应线圈。各向异性媒介对三维感应测井收发装置的不同部分产生反应。因此，此类工具能够测量各向异性的地层电阻率。为

3、了能在复杂地质环境下有效地接受获取到的数据，精确地预测电磁感应线圈在三维地理条件下在水平方向和垂直方向的多个地层的电磁反应尤为重要。在这样的三维各向异性媒介中如何高效地求解麦克斯韦方程是一个具有挑战性但又吸引了极大关注的难题。已知的用于测井的三维电磁建模技术包括有限元技术、最大相关有限体积法、积分方程法、有限差分法等等。每种方法都有利有弊。在这些方法中有限元技术在复杂地理环境下的建模中灵活度最高。然而，他们不像有限差分法那样直白。积分方程法在相对简单的模型中或许更高效，但是它们计算复杂并且随着未知元的增加计算所需的内存容量快速增长。最近研究出的快速算法能够加速求解过程并且减少存储空间需求，但是

4、这个快速算法只能用于各向同性介质。有限差分法在概念和实现上都比较简单，并且已经有许多相关的技术应用于这种方法使之更灵活、高效。有限差分法已经成为地球物理电磁仿真中最流行的技术。本文中，我们提出一种基于有限元的仿真方法，用于计算嵌入在三维完全各向异性介质中的感应源产生的电磁场。介质的各向异性用一个对称的3X3维张量来描述。本算法中，我们用一个基于边缘中心的、交错排列的(staggered-grid)的网格的有限差分方法求解准静态频域里的麦克斯韦方程。给每个节点一个平均电导率张量，目的是得到相对较粗的网格划分。然后，用通用最小裕量法(GMRES)求解得到的矩阵方程，并应用不完全逻辑单元(ILU)先

5、决条件以提高方程的收敛性。开发的代码非常高效，它能确保同时进行三维各向异性介质中多空间和多频率的磁场响应的计算。本文的组织结构如下：首先，我们得到描述各向异性介质中电磁感应的偏微分方程。然后，用基于交错网格的有限差分法求解这个偏微分方程。我们提供一个完全各向异性介质中的详细的有限差分方程。在本论文的数字结果部分，我们提过一些例子以证明本方法的高效性与高能性。最后，将这种方法应用于三维各向异性介质中的三轴感应工具的电磁响应的研究。2、正文2.1 控制方程推导图1所示为一个多层的向异性介质中无限小的任意指向的感应线圈。其电磁感应的控制方程由法拉第定律和安培定律可得：图1 多层各向异性地层中的感应线

6、圈 (1a) (1b)以上两式中为源电流密度，为感应电流密度。上式还考虑了由引起的时域谐波并将其抑制。感应电流密度与电场强度有如下关系式： (2)总电场强度E可以表示为嵌入在地底参考媒介中的源电流引起的主磁场，和由于电导率和介电常数差异引起的分散磁场。在地球物理应用中，我们通常更倾向于一个由分散磁场表示的公式而不是总磁场表示的。因为前者的计算更为健壮和准确，特别是当特别接近信号源进行测量时。如果采用总磁场这一表示方法，通常就需要很好的网格，这会导致计算资源的需求量大增。除此之外，采用总磁场表示时根本不可能获取精确的同步响应。这是因为直接耦合磁场在总磁场中占据主要地位。设公式(1)中E=E，并由

7、公式(1)和(2)推出一个以分散的电磁场表示的一元、二阶偏微分方程(PDE)：(3)上式中是受分散磁场影响的源电流密度，其表达式如下： (4)其中，I是一个3X3维密度矩阵。需要强调的是式中的电导率和介电常数都与位置有关并且具备完全各向异性。比如：(5)2.2 有限差分解接下来，我们用基于如图2所示网格的有限差分方程去解(3)的偏微分方程。解空间在笛卡尔坐标系中被离散化，并且电磁场的组成单元，定义在坐标轴上。磁场的组成单元在y轴与z轴上是错开的，在x轴与z轴上是错开的，在x轴和y轴上是错开的。如图2所示：图2 有限差分方程所用网格结合公式(3)与(4)并扩展旋度操作得到下面所示的关于分散电场的

8、表达式：本文中，一个以13个点为中心的有限差分方程网被用于近似代替(3)中的旋度-旋度操作符。这样，上面的差分方程被转化为如下所示线性系统中的等式(7)其中，矩阵K是维度为的系统矩阵。E是一个长度为的向量，它包含了所有节点的次级电场值，。S(长度为)由(6)式中右边部分推出的次源(secondary-source)向量。系统矩阵K是一个每行有13个非零项的稀疏矩阵。这些非零项与网格间距和介质的与频率相关的属性有关。在线性方程的推导中，采用了电导率平均方法以获取边界节点(定义电磁场的地方)的电导率。经过此处理后，边界上的节点的电导率被表示为与其相邻的四个节点的电导率的加权和。在网格的最边界上对分

9、散磁场应用狄利克雷条件。详细的矩阵表达式可以参见附录。需要指出的是原来的系统方程是非对称的。通过将(A1)乘以，(A2)乘以，和(A3)乘以，我们得到一个形式上对称的系统方程。其中，是网格节点i，j,k的长度；，是编号为i+1与I，j+1与j，k+1与k的节点的中心距离。式(7)中的线性系统能够用最小裕量法(GMRES)高效地求解，并且用不完全逻辑单元预估(ILU)来提高矩阵方程的收敛性。一旦由公式(7)求解出电场，就可以用法拉第定律推导出磁场的值，如公式(8)所示：在实现过程中，(8)中的偏导函数用相应的差分方程表示，这些差分方程可以通过将电场中插入到观测点或几个相邻点得到。等式(7)通常缓

10、慢收敛，特别是当介质中的导电率较大的时候。为了克服这个缺点，我们用一个不完全逻辑单元预估(ILU)来提高矩阵方程的收敛性。2.3 三轴感应工具的电磁响应解公式(3)以得到分散电场。主电场加上分散电场即得到总磁场。公式(3)适用于发射线圈指向任意方向时的情况。对于一个如图3(a)所示的三轴感应线圈，通常包含三个正交的发射线圈和三个正交的接受线圈，如图3(b)所示。因为发射线圈无限小，故可以将它当作一个指向线圈表面法线方向的磁偶极子。磁偶极子源M=(可以由源电流I()表征如下：(9)其中，是发射线圈的位置，a=(，)是投射到(y，z)，(x，z)和(x，y)平面的载流回路区域。发射线圈产生的主磁场

11、通过同质背景介质来计算，这种介质的导电率和介电常数在源点处给出。倾斜线圈的总辐射场是指向x轴，y轴，z轴的源线圈产生的磁场的总和。同质背景模型中，关于频域中电磁场的明确的表达式可以推导。图3 三轴感应工具及其通用线圈结构。 (a)在各向异性地层的三轴感应工具，（b）一个感应线圈结构的通用工具线圈表面的法线方向。磁偶极子源 M = (M x ,M y ,M z )的特征是由电流源决定 I()M = jµ 0 I()a(r r s ) 其中rs是发射线圈的位置, a = (a x ,a y ,a z )是载流回路分别投射到（y，z），（x，z）和（x，y) 区域。由发射机线圈产生的主要区

12、域是由均匀在原点给出电导率和介电常数的背景介质计算得到的。倾斜线圈的总的辐射场是由x方向，z方向和y方向源的总和。对均匀背景模型，对于电磁场在频率域的显式表达式可以推导出来。2.4平均电导率张量由于我们认为地下地层通常是非常复杂的，包括钻孔，倾斜床，地层和各向异性介质，细网格对在不同媒体之间形成复杂结构和接口来说是必要的。然而，细网格需要大量的计算机资源。为了缓解这一困难，可行的办法是采用合理的粗网格形成几何形状，使用平均电导率为每个单元塑造介质的电性能。这是一个很好的在精度和计算复杂性之间的折衷。所以在大多数情况下，用于有限差分法的网格对于介质电性能是独立的，不同的介质可以包含在一个使用的平

13、均电导率单一的矩形单元中。在本文中，一项类似在 10中描述的技术被用来测量平均电导率张量。图4. 划分一个单元计算平均电导率张量我们考虑一个广义的各向异性介质，其导电性是由一个对称3×3张量描述。为了获得每个矩形单元的均电导率张量，单元首先在笛卡尔坐标系中被分成N x× N y×Nz个子单元，如图4所示。每个子单元假设有恒定的电导率张量，表示为 。不失一般性，我们可以推导出 < xx>和< yx>的平均电导率，其他的单元可以得到下面的一个类似的过程。为了得出< xx>，我们在单元的X方向施加电压V0，如图4所示。假设电场在每个子

14、单是均匀的元和电流在垂直于X方向子单元边界连续的，我们可以写出x方向每个子单元的电场如下单元X方向平均电流密度可以表示为具有相同上标 i 的子单元密度和，如果我们选择i1，平均电流为：平均电导率< xx>可以写为：方程（12）意味着平均系数< xx>是由一系列具有相同上标N x子单元组合获得，然后结合并行N x× N y子单元线。接下来，我们得到的非对角项的表达< yx>。对于一个广义各向异性介质，每个子单元的导电率不为零，这意味着任何施加在X方向的电压都可能导致电流在Y方向流动。电流密度J y在每个子单元由于E x由下式给出: (13)y方向平均

15、电流密度是所有N x ×N y×Z N个子单元的平均值，表示为 （14）将（15）代入（14），平均系数< yx>可得：(15）根据以上过程，我们可以得到其他的平均电导率< yy>，< zz>，< xz>和< yz>。为了保持平均电导率张量的对称性，我们设置平均电导率< >为计算的< >和< >的平均值，这就是：(16)图5。两个模型构建了内部一致性检查。（a）模型1。（b）模型23.数值结果基于上述理论，我们开发了一个代码来模拟在任意取向的线圈天线在复杂的各向异性介质的电磁场响

16、应，并进一步扩展到三轴感应工具响应的模拟。为了验证代码，我们首先提出一个内部一致性检查。考虑一个一维层状结构如图5所示，1层，3层和5层是分别具有 50ohm-m, 0.5ohm-m和1.0ohm-m的等方向性介质。2层是各向异性介质，具有11ohm-m水平电阻率和1.9ohm-m垂直电阻率。4层具有2.0ohm-m水平电阻率和 1.0ohm-m垂直电阻率。每一层的深度图所示。在第一个模型，我们假设介质的电阻率张量主轴与X-Y-Z坐标系统的发射器和接收器线圈倾斜角度为60°，如图5所示（一）。在第二个模型中，我们平铺介质的主轴和Z方向的对准源，而多层结构是平铺在X -0 Y -0 Z

17、 -0参考框架，如图5（b）。这两个模型将得到相同的结果，虽然代码的实现是不同的，因此提供了一个用于内部一致性检查的开发代码。在第一个模型中，由于网格图层的边界重合，不需要模拟平均电导率张量。但磁场的轴向分量由x方向源和z-向源的组合所产生的磁场。第二个模型，网格不一致，边界将形成更长时间，所以我们需要计算所哟单元的平均电导率张量。应该指出的是，在得出他们平均电导率张量之前，原有的工作坐标系X-Y-Z每个单元的电导率张量应转换到新的坐标系X'-Y'-Z'。转换可以由原来的电导率张量乘以一个旋转矩阵R实现：（17）旋转矩阵R可以表示为：（18）图6。从两个模型计算出的磁场

18、。(a)图(H xx )，（b）图 (H xx )，（c）图（(H zz ），（d）图 (H zz )。图7 7层非等方向性模型其中，两个欧拉角 和 分别对应叠片的倾斜角和方向角。旋转后，电导率张量将在X '-Y'- Z '坐标系中为全张量，从而代替在原坐标系X-Y-Z中对角张量。图6显示在不同的垂直深度一对发射机和接收机工作在20kHz间隔1.016m时，计算出的磁场强度Hxx（包括在x方向的发射线圈和接收线圈）和H zz（包括在Y方向发射器和接收器线圈。从两个不同的模型的结果可以看出，在验证有限差分方法的实现和平均电导率张量的计算之间显示出完美的一致性。在图6中，除

19、了各向异性的情况下的结果，我们还提出了H xx、H zz为各向同性的情况下，即2层和4层各向同性电阻率分别为1.0ohm-m和1.9ohm-m，从图6（c）和（d）我们可以看到，在各向同性和各向异性时z方向耦合H zz是相同的，这意味着H zz对垂直电阻率影响很小。接下来，我们考虑一个如图7所示的7层各向异性地层模型。介质层1，3，5和7是各向异性的，具有10ohm-m水平电阻率和 1ohm-m垂直电阻率。2层，4层和6层有。钻孔的直径和分别是21.59cm 和30.48cm。钻孔中的泥浆电阻率为0.4ohm-m。其他参数在图中给出。我们使用的工具包括三个并列的正交发射线圈和用来测量磁场分布的

20、三配置正交模型接收线圈。发射器和接收器线圈定位在X，图8。100kHz的2线圈轴工具的磁场响应虚部的比较。（a）H xx（b）H yy（C）H zz.Y和Z方向。倾角（工具轴线和层边界的法线之间的夹角）为30°。发射器和接收器之间的间距为1.8m。图8显示计算得到的H xx，Hyy和H zz在100kHz时响应的虚部。 15 的结果也在图中做了比较。观察到的良好的一致性，验证了目前的三维有限差分程序。在这些数据中，在忽略钻孔和入侵情况下我们还提出了一维计算结果的分析方法 22 。可以看出，在没有入侵，1D和3D结果非常接近彼此，因为在这种情况下，钻孔的影响不明显。然而，在层2，4和6

21、有入侵存在时，1D和3D效果有差异，这意味着3D仿真要得到准确的结果，需要更多对周围介质的信息。图9 有/无预处理时收敛性的比较。在图9中，我们比较了线性系统在有无ILU预调节器时的收敛速度。据观察，没有预处理，归一化残差范数（NRN）的1%即使经过200次迭代也不能实现。然而，在经过预处理之后，一个NRN 的0.1%只需要6次迭代就可以使Ifil= 400=或者只需要9次迭代就可以使Ifil=100。Ifil代表的矩阵L和U的第i行元素的数量，是为了保持建立预条件矩阵。大的Ifil代表更快的收敛速度，但对内存要求更高。因此，选择Ifil必须在存储和收敛性之间做良好的折衷。作为一个任意取向的线

22、圈天线的应用程序，我们考虑一个如图10所示例子。一个位于三层介质之间，间隔2.4384m的发射器-接收器对。中心层厚6.096m，电阻率为20ohm-m。上、下层的电阻率分别为2和5 ohm-m。图10。一对在3层介质之间倾斜的发射机-接收机图11。一个倾斜的发射器-接收器对一个三层各向同性各向异性介质的响应。（一）振幅衰减（B）相移。该工具（发射器-接收器对）的取向平行于边界，模拟90°水平井。发射线圈指出沿Z 0轴而接收器对准相对于Z 0轴和绕Z轴45°的测量方向。图11显示了振幅衰减和相移的响应（接收机指出45°-45°时接收信号的比率）作为一种在

23、100kHz和400kHz的垂直深度函数。我们可以从图11观察下面的事实，首先，两者的振幅比和相位移的峰值出现在边界，这就可以确定不同媒体的界限；第二，峰值随频率的增加而增加；此外，当工具到达导电层下面时，方向的相移和衰减是正值。相反，当工具在导电层上面时，方向的相移和衰减是负的。这是一个非常重要的地质应用。极性可用于在钻井中作为一个简单的指标来确定定向钻机应向上或者向下移动。在图11中，我们还提出了一个在中心层的水平电阻率和垂直电阻率为20ohm-m 和4ohm-m各向异性的情况下的比较结果。从图中，我们可以看到，在中心层的各向异性的存在会改变相移和衰减的极性，如图虚线曲线68m之间所示。4

24、结论在本文中，我们对任意取向的线圈天线在复杂地下介质产生的电磁场提出了一个三维有限差分模拟方法。该方法是用于预测和研究三轴感应工具在多层次倾斜各向异性地层的电磁响应。通过本方法得到的结果与从传输线理论得到的一维结果做了对比。它的结论是一维模拟结果与三维模拟结果在地层的和地层差异变化时差异明显，并且随着地层的变化差异将增大。当顶部介质是各向异性时这种差异也将增加了。因此，为了充分了解在一个复杂的多层各向异性地层钻孔的影响，全三维模拟是必要的，尽管它比一维分析模拟慢。该方法还可以被用来模拟和设计地质导向钻井工具。参考文献1. Moran, J. H. and S. Gianzero, “Effec

25、ts of formation anisotropy on resistivity-logging measurements,” Geophysics, Vol. 44, 12661286, 1984.2. Anderson, B. and S. K. Chang, “Synthetic induction logs by the finite element method,” The Log Analyst, Vol. 23, 1726, 1982.3. Anderson, B., “Simulation of induction logging by the finite- element

26、 method,” Geophysics, Vol. 49, 19431958, 1984.4. Everett, M. E., E. A. Badea, L. C. Shen, G. A. Merchant, and C. J. Weiss, “3-D finite element analysis of induction logging in a dipping formation,” IEEE Trans. Geoscience and Remote Sensing,Vol. 39, 22442252, 2001.5. Han, W. S., “3D finite element si

27、mulation method of induction andMWD tools,” Ph.D. Dissertation, University of Houston, 2004. Progress In Electromagnetics Research B, Vol. 44, 2012 2816. Zhdanov, M. S., S. K. Lee, and K. Yoshioka, “Integral equa-tion method for 3D modeling of electromagnetic fields in complex structures within homo

28、geneous background conductivity,” Geo-physics, Vol. 67, 333345, 2006.7. Fang, S., G. Z. Gao, and C. Torres-Verdin, “Efficient 3D elec-tromagnetic modeling in the presence of anisotropic conductive media, using integral equations,” Exploration Geophysics, Vol. 37,239244, 2006.8. Nie, X. C., N. Yuan,

29、and R. Liu, “Simulation of LWD tool response using a fast integral equation method,” IEEE Trans. On Geoscience and Remote Sensing, Vol. 48, No. 1, 7281, Jan. 2010.9. Weiss, C. J. and G. A. Newman, “Electromagnetic induction in a fully 3-D anisotropic earth,” Geophysics, Vol. 67, No. 4, 11041114, Jul

30、.Aug. 2002.10. Wang, T. and S. Fang, “3-D electromagnetic anisotropic modeling using finite difference,” Geophysics, Vol. 66, No. 5, 13861398,Sept.Oct. 2001.11. Davydycheva, S., V. Druskin, and T. Habashy, “An efficient finite- difference scheme for electromagnetic logging in 3D anisotropic inhomoge

31、neous media,” Geophysics, Vol. 68, No. 5, 15251536,Sept.Oct. 2003.12. Newman, G. A. and D. L. Alumbaugh, “Three-dimensional induction logging problems, Part 2: A finite-difference solution,” Geophysics, Vol. 67, No. 2, 484491, Mar.Apr. 2002.13. Streich, R., “3D finite-difference frequency-domain mod

32、eling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy,” Geophysics, Vol. 74, F95F105,2009.14. Davydycheva, S., D. Homan, and G. Minerbo, “Triaxial induction tool with electrode sleeve: FD modeling in 3D geometries,” Journal of Applied Geophysics, Vol. 67

33、, 98108, 2009.15. Hou, J. S., R. K. Mallan, and C. Torres-Verdin, “Finite-difference simulation of borehole EM measurements in 3D anisotropic media using coupled scalar-vector potentials,” Geophysics, Vol. 71, No. 5,G225G233, Sept.Oct. 2006.16. Novo, M. S., L. C. da Silva, and F. L. Teixeira, “Finit

34、e volume modeling of borehole electromagnetic logging in 3-D anisotropic formations using coupled scalar-vector potentials,”IEEE Antennas and Wireless Propagation Letters, 549552, 2007.17. Novo, M. S., L. C. da Silva, and F. L. Teixeira, “Comparisonof coupled-potentials and field-based finite-volume techniques for 282 Yuan, Liu, and Nie modeling of b