CN102265188B - 电磁感应和诱导极化效应的量化分解方法 - Google Patents
电磁感应和诱导极化效应的量化分解方法 Download PDFInfo
- Publication number
- CN102265188B CN102265188B CN200880132434.XA CN200880132434A CN102265188B CN 102265188 B CN102265188 B CN 102265188B CN 200880132434 A CN200880132434 A CN 200880132434A CN 102265188 B CN102265188 B CN 102265188B
- Authority
- CN
- China
- Prior art keywords
- field
- time
- input function
- polarization
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000005674 electromagnetic induction Effects 0.000 title claims abstract description 40
- 230000000694 effects Effects 0.000 title claims description 11
- 230000001052 transient effect Effects 0.000 claims abstract description 38
- 238000012905 input function Methods 0.000 claims abstract description 36
- 230000006698 induction Effects 0.000 claims description 13
- 238000011002 quantification Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 abstract description 27
- 238000011160 research Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000005672 electromagnetic field Effects 0.000 description 6
- 230000005284 excitation Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000010606 normalization Methods 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000013016 damping Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 125000001183 hydrocarbyl group Chemical group 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 230000005520 electrodynamics Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001502 supplementing effect Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/086—Processing
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于电法勘探研究领域。其要点是：利用偶极子轴向装置测量极化介质在通过电流脉冲时的瞬变过程。建立同电磁感应场和诱导极化场有不同依赖关系的若干个函数。其中一个函数同DU(t)＝ΔU(t)/ΔU0相比，可增大电磁感应同诱导极化的比值，式中ΔU0——在电流通过时测量的电势差ΔU。另外一个函数同DU(t)相比，可缩小电磁感应同诱导极化的比值。第三个函数的构成则是瞬变场时间导数和空间导数的组合。对同一记录点录得的所有输入函数进行同时反演，包括DU(t)，从而建立介质剖面的地电模型。在建立的模型中使每一层的极化率等于零，计算正解问题求出电磁感应场。在同一模型中使波数等于零，正解问题求出诱导极化的电流分量场。评估诱导极化电流分量的局域变化，对其进行地质解释。技术成果：可量化确定电磁感应和诱导极化场强，误差不超过0.5％，并且可在0.1μs至数秒的较大时间间隔内分解瞬变场，其中包括针对高导纳剖面。
Description
技术领域
本发明属于电法勘探研究领域，是一种地下岩层性质的测定方法，同时还对地层的电学性质参数进行测量，并分解测得参数。
本发明用于陆上及海洋中的可控电磁场源电法勘探，适用于寻找并圈定油气藏的综合性油气勘探工作。它对次生蚀变岩在油气迁移作用下的回应进行分离，并以其为基础将测得信号分解成某些分量。本发明可保证一系列用于分层确定介质的地电参数值、诱导极化和电磁感应过程特征的技术方法。
背景技术
已知发明《地电勘探方法(方案)》，专利号2235347RU，发表日期2004年8月27日，国际专利分类号G01V3/06。根据该发明，在被研究介质的深层周期性依次通过矩形电流脉冲，并在每两次脉冲之间保持间歇，从而激发电磁场。在整个过程中，于每次电流脉冲的结尾和两次脉冲之间的间隙测量轴向和正交电势差的瞬时值，确定电势差。按公式计算n个归一化参数，根据极化导电介质中的偶极子源场强的数学物理差分波动方程运算反解问题，建立其几何结构和电学参数最接近被研究介质的介质模型根据该方程中包含的电物理参数，建立该模型的时间剖面图。利用该发明可发现并圈定油气藏、评估其饱和质量，并可分解电导率参数和诱导极化参数。然而，该发明不能解决分别量化研究电磁感应过程和诱导极化过程的重要问题，未考虑所建立模型的电物理参数的等值性，不能确定这些参数的信息容量。
已知发明《使用聚焦电流的海洋地电勘探方法(方案)》，专利号RU2284555，发表日期2006年9月27日，国际专利分类号G01V3/06。根据该发明，在被研究介质的深层激发电磁场，测量第一和第二电势差的瞬时值，使沿剖面的轴向合成电势差等于零并从等式中求得因数，计算三个独立的归一化电学参数集合，进行反解数学运算，求得三个电物理参数：电导率、诱导极化度和诱导极化电势差的衰减时间常数，并按这些参数建立三个时间剖面图。该发明的本质在于分解电导率参数和诱导极化参数，作为补充，还可确定与前两个参数同样重要的第三个参数——诱导极化电势差的衰减时间常数。然而，该发明仅能获得有关被研究介质所有结构元素的总和信息，这是因为电磁场源的被测电流I在介质中的空间分布不受任何监测，没有在真实存在的三维非均质介质中关于该分布的信息。该发明不能解决分别量化研究电磁感应过程和诱导极化过程的重要问题，未考虑所建立模型电物理参数的等值性，不能确定这些参数的信息容量。
已知发明《地电勘探方法(方案)》，专利号2231089，发表日期2003年7月8日，国际专利分类号G01V3/06。根据该发明，在被研究介质的深层周期性依次通过矩形电流脉冲，并在每次脉冲之后保持间歇，从而激发电磁场。在观测点测量电势差，使用归一化参数值和电化学极化导电介质中的偶极子源电场场强的数学物理差分波动方程，运算反解问题，得出每个介质单元所固有的电物理参数：电导率、诱导极化和诱导极化电势差的衰减时间常数，并按这些参数建立时间剖面。利用该发明可发现并圈定油气藏、评估其饱和质量，并可分解电导率参数和诱导极化参数。然而，该发明不能将地质介质单元所固有的诱导极化与瞬变电动力学过程彻底分解开来，该电动力学过程同所示地质岩石的剖分单元的电导率相关。该发明亦不能解决分别量化研究电磁感应过程和诱导极化过程的重要问题，未考虑所建立模型电物理参数的等值性，不能确定这些参数的信息容量。
已知专利《地电勘探信号的分析方法》，专利号2253137RU，发表日期2005年1月10日，国际专利分类号G01V3/38。根据该发明，将每个测量信号同时用至少两个调制器转换，然后找出最小的步距间隔。该发明通过对信号(固定分量)和干扰(随机分量)的同时分析，能够提高准确性和信息容量。然而，该发明的目的是利用该方法在计算机中对比增量调制器集成电路的输出信号和相应通道的输入信号。而这种分析方法的使用，却必须具备至少两个量子化步距差别很大的增量调节器、一个接收电极系统、一套磁场传感器或者具有足够宽频带性和空间角度选择性的任何其他传感器。该方法不能用于分析在使用诱导极化过程寻找和勘探油气藏时获得的参数。
已知发明《充满井内流体的井筒周围地质构造的电阻率测定方法》，专利号2209451RU，发表日期2003年7月27日，国际专利分类号G01V3/38。该发明包括获取某些图表、选择模拟剖面、引入附加模拟的剖面、重复修正模拟剖面直至图表和相应模拟图之间的差值小于选定的阈值、将每个模拟剖面作为值的关系函数进行修正等。该发明考虑了井内流体对井筒周围构造的入侵，完善了地质构造电阻率的测定方法。然而，该发明不能发现并圈定油气藏、评估其饱和质量，并分解电导率参数和诱导极化参数。此外，该发明亦不能解决分别量化研究电磁感应过程和诱导极化过程的重要问题，未考虑所建立模型电物理参数的等值性，不能确定这些参数的信息容量。
同本发明在地电性征的处理方法上最接近的是《地下构造的性质测定方法》发明，专利号2294547RU，发表日期2007年2月27日，国际专利分类号G01V3/38。该发明包括构造参数模型的拟定、模拟电学性质的测量、将模拟测量同已测得电学性质相比较、以各种测量为基础确定构造的物理/石油物理性质以保证联合反演、利用电学性质模型和各种测量评估构造的物理/石油物理性质，还包括使用综合模型。该发明可保证在应用电力对钻井进行剖分时能够实施更完善的数据处理方法，并且更准确地确定Rt值以评价构造特征，在此过程中使用将初始信息引入反演图的构造层系测定。初始信息即稳定未精确设置的“可能性扩展”任务的几何信息。然而，该发明不能发现并圈定油气藏、评估其饱和质量，并分解电导率参数和诱导极化参数。此外，该发明亦不能解决分别量化研究电磁感应过程和诱导极化过程的重要问题，未考虑所建立模型的电物理参数的等值性，不能确定这些参数的信息容量。该发明有另一个目的，即当介质模型可通过迭代反演方法进行简单明确时实现该方法，不将测得信号分解成某些分量。而在我们的发明中，通过反解运算得出的介质模型用于将测得信号ΔU(t)分解成分量EM(t)和IPg(t)。
发明内容
当介质中通过电流时，会产生一些由不同物理现象所决定的效应。此时会发出电磁感应，其特点是符合波动理论的感应场，并且发生由介质中的电化学过程所决定的诱导极化。
到目前为止，对这两种效应的分解办法是增加记录介质中激发脉冲瞬变的时间，而电磁感应和极化正是由该脉冲所引起的。然而，随着记录时间的增加，电磁感应场在测量之前几经衰减到了几乎可以忽略不计的微小数值。相反，如果缩短测量时间，在激发脉冲之后的初始时间内，极化效应几乎没有。在测量沉积盖层较厚的地台时，由于沉积盖层的导纳较大(数百或数千西门子)，电磁波衰减会比较慢。该问题的存在及其解决方法的缺失在J.维特的书中有提到(《大地电磁学》，James R.Wait，莫斯科，地址出版社，1987年，235页)。
本方法的技术成果是它可以解决分别量化研究电磁感应场和诱导极化的重要问题，考虑了所建立模型的电物理参数等值性，能够确定其信息容量。
本方法能够：
-量化测定直流电的电磁感应场强和诱导极化场强，误差不超过0.5％；
-保证电磁场分解从1微秒至数秒的较大时间间隔，包括针对高导纳的地质剖分。
由于函数ΔU总是解析性的，即总是连续的，并且其二阶导数亦为连续，因而本方法是可用的。此外，由于部分诱导极化场由感应电流引起，而部分感应场反过来又由极化场引起，所以电磁感应场和诱导极化场是非相加性的。我们将只讨论直流电在通过介质的时间段内所引起的诱导极化。
对所声明技术成果的保证，是因为电磁感应和诱导极化效应的量化分解方法包括：
-测量极化介质的瞬变过程；
-水平层状极化介质范围内的电磁测量模拟；
-从瞬变过程中分离出电动力学分量；
-从瞬变过程中分离出诱导极化的电流分量；
-评估诱导极化电流分量的局域变化，对其进行地质解释。
本方法的特点是记录几个取决于介质参数的、在介质表面测量的输入函数，这些函数对电磁感应场和诱导极化场有不同的依赖关系，同时其中一个函数(例如第一输入函数Ps)同电势差输入函数DU(t)相比，可增大电磁感应和诱导极化的比值，而另一个函数(例如第二输入函数P1)同DU(t)相比，则可缩小电磁感应和诱导极化的比值，还有一个函数(例如第三输入函数Dφ)则是瞬变场的时间导数和空间导数的组合。然后，同时反演所有输入函数，得出剖面的地电模型，并且此时表示诱导极化的介质参数数量明显过剩，因此得出的介质模型不是唯一的，而是位于某个等值区域内。在得出的模型中使每一层的极化率都等于零，通过正解问题计算出电磁感应场EM；使同一模型中的波数等于零，正解问题计算出诱导极化的电流分量场IP。此外，举例而言，建立参数P1，将其作为同ΔU(t)相比提高诱导极化/电动力比值的筛选参数，以瞬变场第二轴向电势差Δ2U(t)同瞬变场第一轴向电势差ΔU(t)在断电制式下的比值计算P1，计算公式为：P1(t)＝Δ2U(t)/ΔU(t)Itime off；或者，举例而言，建立参数Ps，将其作为同ΔU(t)相比降低诱导极化/电动力比值的筛选参数，以瞬变场第二轴向电势差Δ2U(t)同瞬变场第一轴向电势差ΔU(t)在通电制式下的比值计算Ps，计算公式为：Ps(t)＝Δ2U(t)/ΔU(t)Itime on；再或者建立参数作为按公式Dφ(t)＝It(t)-P1(t)计算的参数It(t)同参数P1(t)的差值，并且It(t)作为时间导数Δ2U(t)同ΔU(t)的比值，其计算公式为：It(t)＝Δt(Δ2U(t))/Δt(ΔU(t))。
在建立同一记录点测得的输入函数之后，将它们最小化，为此在输入时给出参数的最小化程序：
-DU(t)＝ΔU(t)/ΔU0，式中：ΔU0为电流通过时测量的电势差ΔU(t)；
-P1(t)＝Δ2U(t)/ΔU(t)，在断电制式下进行计算；
-Ps(t)＝Δ2U(t)/ΔU(t)，在通电制式下进行计算；
-Dφ(t)＝It(t)-P1(t)，式中：It(t)＝Δt(Δ2U(t))/Δt(ΔU(t))，
并对它们在水平分层极化介质范围内进行反演，并且此时表示模型极化性质的地电参数数量明显过剩。同时作为一种特殊情况，正解波数等于零求出IP场，正解极化率等于零求出EM场。
完成所有运算之后，检查解的可重复性。
为解决所设置的效应分解任务，需要完成以下几个步骤。
第I阶段：
记录几个可在介质表面测量并建立的输入函数。
函数的设计应使它们对电磁感应和诱导极化有不同的依赖关系。其中一个同电压降曲线DU相比，可增加电磁感应同诱导极化的比值，而另一个则缩小该比值。这些函数是瞬变场同其空间导数的组合。此外，还有一个函数是时间导数和空间导数的组合。
第II阶段：
1.同时反演(反解问题)所有输入函数，建立一个剖面模型。同时逐层求解介质参数：电阻率ρ，极化率η，弛豫时间τ，幂c。极化率、弛豫时间和幂根据Cole-Cole模型确定。利用这些参数对模型进行参数化时应使之明显过剩，意即获得等值性的某些解。
2.在建立的模型中使每层的极化率等于零，正解问题求得电磁感应场EM。
3.在同一模型中使波数等于零，正解问题求得直流电的诱导极化场IP。
重复运算(第II阶段)显示计算出的EM场和IP场重合，并精确至反解运算的解。所得出的模型是不相同的，但场是重合的。
其中，IP——使剖面任意点的波数等于零，对该模型进行正解运算得出的诱导极化场；EM——使剖面任意点的极化率等于零，对该模型进行正解运输得出的电磁感应场。
附图说明
本方法用图表进行说明和证明。
图1是ΔU1和ΔU2的曲线，即偶极子轴向装置在接收线OM和ON上的电势差。
这个图形显示存在着IP场和EM场不同组合的无穷集。对合成模型计算了衰减曲线DU并进行了两次不同零点逼近的反演，精确至0.1％。得到两个不同的等值(按DU参数)解，其中每一个符合显著不同的IP场和EM场的值。
图2显示的是在瞬时过程不同时期(不同相)的EM和IP电流分布间隔：初期(A)、过渡期(B)和晚期(C)。
此时EM场服从集肤效应，于瞬态过程中在下半空间内扩散，而由其引起的涡流则趋向于在下半空间内呈均匀分布。IP场则相反，总是呈空间不均匀分布，并且其强度随着同电流源(发射源)的距离增大而减弱。
根据所设置问题的涵义，需要为瞬变场引入可在同一点记录的转换参数，并且这些转换参数对诱导极化场和感应场具有不同的函数依赖。最自然的解次方法——引入瞬变场的归一化空间导数(更准确来讲是有限差分)，做为对瞬变场本身的补充。这是因为感应场和诱导极化场的空时结构是不同的。
在图3中显示的是：a).测量装置在轮船移动时记录DU(t)、P1(t)、Ps(t)、参数的示意图；b).瞬变激发电流在振荡线AB上的形状；c).在接收线上测量的信号Δ2U(t)иΔU(t)。
图4显示的是DU(t)曲线在地球表面上随时间的衰减，包括存在(B)和不存在(A)极化、在接收线MO(距离为1000m(1))和接收线ON(距离为2400米(2))上的情况。
图5显示的是参数P1(t)在地球表面上的曲线，包括存在(2)和不存在(1)极化的情况。条件是η＝1％，ρ＝50Ω·m和η＝0％。
图6显示的是参数Ps在地球表面上的曲线，相对于条件η＝4％和η＝0，其中η-在Cole-Cole模型术语中的极化系数。
图7显示的是在η＝0和η＝1％条件下计算的曲线，其中η-在Cole-Cole模型术语中的极化系数。曲线相对于的直线建立。除了空间性质外还使用时间导数，从而获得更详细的瞬时过程特征。
图8显示的是三维参考模型。利用它计算了符合实地测量所得到的参数综合曲线(DU(t)、P1(t)、Ps(t)、)，以及EM场和IP场的值。并且给出了与其等效的一维模型。
图9显示的是EM场和IP场，它们同图8中所给出的模型相吻合，精确度为0.15％。
图10显示的是在北古利亚耶夫斯卡娅区(巴伦支海大陆架)进行实地测量时的观测曲线及与其相符的模型。给出了等效地电模型之一。
图11显示的是在北古利亚耶夫斯卡娅区(巴伦支海大陆架)进行实地测量时的观测曲线及与其相符的模型。给出了与图10不同的等效地电模型之一。
图12显示的是从等效模型中得出的EM场和IP场的衰减曲线。它们彼此间的差值要小于测量误差，即小于0.5％。此时IP场从在大范围时间间隔内比其高2个数量级的EM场中分离出来。
图13显示的是已知油气田的IP场剖面等值线和油气藏沿IP场的划定。
具体实施方式
已知对地球表面进行测量获得的一系列输入数据。
数据的组成应使它们对诱导极化和电磁感应有不同的函数依赖关系。以这些数据为基础，在表示诱导极化过程的剖面参数数量过剩的条件下反演得出地电模型。
如果在这个模型中，令波数为零正解得出IP场，令极化率为零正解得出EM场，则此时的IP场和EM场符合可重复性条件。这就意味着对不同等效模型进行IP场和EM场计算所得出的偏差，不会超过实际观测误差(约0.5％)。
为实现本方法，
1.在电流脉冲切断制式下测量时建立参数P1(见图1)，将其作为同ΔU(t)相比增大诱导极化/电动力比值的筛选参数。
-为此以瞬变场第二轴向电势差Δ2U(t)同瞬变场第一轴向电势差ΔU(t)在断电制式下的比值计算P1，
即：P1(t)＝Δ2U(t)/ΔU(t)。
在图1中显示了ΔU1和ΔU2——偶极子轴向装置在接收线OM和ON上的电势差。参数P1对瞬变场的空间不均匀性作出反应。电磁场以扩散方程表示，因此在增大衰减时间时会趋向于在介质中均匀分布。诱导极化场正比于极化电流密度，因此它总是呈空间不均匀性(图2)[列格依多P.Yu.，《油气地球物理中研究极化剖面的差分归一化地电勘探理论和技术》。地质与矿物学正博士学位答辩论文。1998年]。
2.在电流脉冲接通制式下测量时建立参数Ps，将其作为同ΔU(t)相比降低诱导极化/电动力比值的筛选参数。
-为此计算瞬变场第二轴向电势差Δ2U同第一轴向电势差ΔU在通电制式下的比值。
在瞬变初期的诱导极化/电动力比值很小，而在晚期阶段Ps参数基本上由阻抗的分布来决定，极化率对它的影响很小。因此Ps参数可被视为同瞬变场ΔU相比较能够缩小诱导极化/电动力比值的筛选参数(图6)。
3.建立参数以便更详细地表示瞬变过程，因为这个参数以时间导数Δ2U(t)和ΔU(t)为基础计算。
-它是相对于作为瞬变时间对数的假定单位的参数差，即参数It同参数P1的差值：
-其中，It＝Δt(Δ2U(t))/Δt(U)，即时间导数Δ2U(t)同ΔU(t)的比值。诱导极化场和电磁感应场的衰退时间一般是不同的，因此有意识地使用时间参数，以更详细地描述瞬变过程(图7)。
4.对在同一点记录的参数进行最小化。
-为此在输入时给出参数的最小化程序(图3)；
-DU(t)＝ΔU(t)/ΔU0，式中：ΔU0为电流通过时测量的电势差ΔU(t)；
-P1(t)＝Δ2U(t)/ΔU(t)，在断电制式下进行计算；
-Ps(t)＝Δ2U(t)/ΔU(t)，在通电制式下进行计算；
-Dφ(t)＝It-P1，式中：It＝Δt(Δ2U(t))/Δt(ΔU)。
在模型的地电参数数量明显过剩的条件下，对它们在水平分层极化介质范围内进行反演，以尽可能准确地接近对非均质各向异性介质和梯度介质记录的瞬变场曲线。
5.按建立的介质模型计算IP场和EM场。
-使模型在介质所有区域的极化率等于零，正解问题求出电磁感应场EM，
-使模型在介质所有区域的波数等于零，正解问题求出诱导极化场IP。
6.检查解的可重复性。
-对EM场和IP场的重复性检查显示，它们实际上对所有种类的剖面都是可重复的，并且误差不超过测量误差(约0.5％)(见图10、13)。
感应场和极化场的分解方法，可以用以下示例进行展示。
在三维模型中获得符合实地测量的参数综合曲线和EM场与IP场的值(图8)。
对于综合曲线，已在水平分层极化介质范围内建立了等值模型，以其为基础计算了EM场和IP场(图8)。
对于从三维模型和水平分层模型中所得的场值进行了比较。它们之间的误差(0.15％)不超过三维模型综合曲线与水平分层极化介质模型曲线之间的偏离误差。由于偏离很小，曲线之间的差别几乎看不出来。
实践当中，对电磁感应效应和诱导极化效应的量化分解，在各种不同的地质和地球物理条件下进行了实验。
对已知的北古利亚耶夫斯卡娅区(巴伦支海大陆架)的油气田进行了电法勘探，测量了参数DU(t)、P1(t)、Ps(t)、为场曲线建立了水平分层极化模型，逼近误差小于0.5％(图11、12)。此时，为同一个测量标点建立了几个等效模型，并且计算出的EM场和IP场的图形偏离小于0.5％(图13)。在此条件下，成功分离出了比EM场小两个数量级的IP场。通过建立IP场的等值线，清除地显示出了油气藏(图14)。其中：
Ro-阻抗，Ω·m；
c-在Cole-Cole模型中的幂；
H-层的厚度，m；
η-Cole-Cole模型的极化率系数，％
t-Cole-Cole模型的弛豫时间，s。
例如以下参数的模型：
Ro＝10Ω·m，η＝5％，t＝0.5s，c＝0.5，H＝100m；
Ro＝50Ω·m，η＝2％，t＝0.5s，c＝0.5，H＝100m；
Ro＝20Ω·m，η＝2％，t＝1s，c＝0.5，H＝200m；
Ro＝5Ω·m，η＝5％，t＝1s，c＝0.5，H＝300m；
Ro＝200Ω·m，η＝7％，t＝1s，c＝0.5；
Ro＝500Ω·m，η＝0，H＝450m；
Ro＝5Ω·m，η＝0，H＝800m；
Ro＝1000Ω·m，η＝0，H＝∞。
对该模型计算了符合实地测量所得的输入数据值，并计算了IP场和EM场，为巴伦支海大陆架北古利亚耶夫斯卡娅油气田建立了曲线图(见图13)。换句话说，根据所设置问题的涵义，为瞬变场引入了转换参数，这些转换参数理论上可在同一点记录得到，并且对诱导极化场和感应场具有不同的函数依赖。同时引入瞬变场的归一化空间导数(更准确来讲是有限差分)做为对瞬变场本身的补充，因为感应场和诱导极化场的空时结构是不同的。从而得出了所声明的技术成果。
Claims (9)
1.一种电磁感应和诱导极化效应的量化分解方法，包括对极化介质进行瞬变场测量，在水平分层极化介质范围内模拟电磁测量，从瞬变场中分离出电磁分量，从瞬变场中分离出诱导极化的电流分量，评估诱导极化电流分量的局域变化，对其进行地质解释；
其特点是记录复数个依赖于介质参数的、在介质表面测得的输入函数，所述输入函数包括电势差输入函数DU(t)、第一输入函数、第二输入函数与第三输入函数，并且这些函数同电磁感应场和诱导极化场有着不同的依赖关系；其中第一输入函数同电势差输入函数DU(t)相比，可提高电磁感应同诱导极化的比值；第二输入函数同电势差输入函数DU(t)相比，可降低电磁感应同诱导极化的比值；第三输入函数则是瞬变场时间导数和空间导数的组合；其中，DU(t)＝ΔU(t)/ΔU0，ΔU0为电流通过时测量的电势差ΔU(t)，ΔU(t)为时间t时电流通过时测量的电势差；然后对同一记录点测得的所有输入函数进行同时反演，从而建立介质剖面的地电模型；此时表示诱导极化的介质参数数量明显过剩，因此建立的介质模型不是唯一的，而是位于某个等值区域内；在建立的模型中使每一层的极化率等于零，计算正解问题得出电磁感应场EM；同时在同一模型中使波数等于零，正解问题求得诱导极化场IP的电流分量。
2.根据权利要求1所述的量化分解的方法，其特点是建立一组输入函数，并且其中，所述第一输入函数为Ps(t)，通过与电势差输入函数DU(t)相比增大电磁感应同诱导极化的比值进行计算，其中，Ps(t)＝Δ2U(t)/ΔU(t)，ΔU(t)为在通电制式下，时间t时电流通过时测量的电势差ΔU(t)，Δ2U(t)为在通电制式下，时间t时瞬变场第二轴向电势差；所述第二输入函数是P1(t)，通过与电势差输入函数DU(t)相比缩小电磁感应同诱导极化的比值进行计算，其中，P1(t)＝Δ2U(t)/ΔU(t)，ΔU(t)为在断电制式下，时间t时电流通过时测量的电势差，Δ2U(t)为在断电制式下，时间t时瞬变场第二轴向电势差；所述第三输入函数是Dφ(t)，通过组合瞬变场时间导数和空间导数进行计算，其中，所述第三输入函数是Dφ(t)作为按公式Dφ(t)＝It(t)-P1(t)计算的参数It(t)同参数P1(t)的差值，并且It(t)作为时间导数Δ2U(t)同ΔU(t)的比值，其计算公式为It(t)＝Δt(Δ2U(t))/Δt(ΔU(t))，Δt(Δ2U(t))为Δ2U(t)的时间导数，Δt(ΔU(t))为ΔU(t)的时间导数。
3.根据权利要求1所述的量化分解的方法，其特点是建立所述第二输入函数P1(t)，将其作为同ΔU(t)相比增大诱导极化/电动力比值的筛选参数，以瞬变场第二轴向电势差Δ2U(t)同第一轴向电势差ΔU(t)在断电制式下的比值进行计算，计算公式为：P1(t)＝Δ2U(t)/ΔU(t)|time off。
4.根据权利要求1所述的量化分解的方法，其特点是建立所述第一输入函数Ps(t)，将其作为同ΔU(t)相比缩小诱导极化/电动力比值的筛选参数，以瞬变场第二轴向电势差Δ2U(t)同第一轴向电势差ΔU(t)在通电制式下的比值进行计算，计算公式为：Ps(t)＝Δ2U(t)/ΔU(t)|time on。
5.根据权利要求1所述的量化分解的方法，其特点是建立所述第三输入函数Dφ(t)，作为按公式Dφ(t)＝It(t)-P1(t)计算的参数It(t)同参数P1(t)的差值，并且It(t)作为时间导数Δ2U(t)同ΔU(t)的比值，其中，Δ2U(t)为瞬变场第二轴向电势差、ΔU(t)为第一轴向电势差，其计算公式为：It(t)＝Δt(Δ2U(t))/Δt(ΔU(t))，其中，Δt(Δ2U(t))为Δ2U(t)的时间导数，Δt(ΔU(t))为ΔU(t)的时间导数。
6.根据权利要求1所述的量化分解的方法，其特点是对在同一记录点测得的输入函数同时进行最小化，为此在输入时给出参数的最小化程序，其中，所述的第一输入函数Ps(t)、第二输入函数P1(t)与第三输入函数Dφ(t)为：
-P1(t)＝Δ2U(t)/ΔU(t)，式中：Δ2U(t)为时间t时瞬变场第二轴向电势差，ΔU(t)为第一轴向电势差，在断电制式下进行计算；
-Ps(t)＝Δ2U(t)/ΔU(t)，在通电制式下进行计算；
-Dφ(t)＝It(t)-P1(t)，式中：It(t)＝Δt(Δ2U(t))/Δt(ΔU(t))，其中，Δt(Δ2U(t))为Δ2U(t)的时间导数，Δt(ΔU(t))为ΔU(t)的时间导数；
并对它们在水平分层极化介质范围内进行反演，而且此时表示模型极化性质的地电参数数量明显过剩。
7.根据权利要求1所述的量化分解的方法，其特点是使波数等于零，计算对模型进行最小化得到的正解问题，从而算出诱导极化场IP。
8.根据权利要求1所述的量化分解的方法，其特点是使极化率等于零，计算对模型进行最小化得到的正解问题，从而算出电磁感应场EM。
9.根据权利要求1～8任一项所述的量化分解的方法，其特点是检查解的可重复性。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2008141889/28A RU2399931C2 (ru) | 2008-10-23 | 2008-10-23 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
| RU2008141889 | 2008-10-23 | ||
| PCT/RU2008/000711 WO2010047613A1 (ru) | 2008-10-23 | 2008-11-17 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102265188A CN102265188A (zh) | 2011-11-30 |
| CN102265188B true CN102265188B (zh) | 2014-07-30 |
Family
ID=42119492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200880132434.XA Expired - Fee Related CN102265188B (zh) | 2008-10-23 | 2008-11-17 | 电磁感应和诱导极化效应的量化分解方法 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8138763B2 (zh) |
| EP (1) | EP2352041A4 (zh) |
| CN (1) | CN102265188B (zh) |
| MX (1) | MX2011002163A (zh) |
| RU (1) | RU2399931C2 (zh) |
| WO (1) | WO2010047613A1 (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120099396A1 (en) * | 2010-10-22 | 2012-04-26 | Chevron U.S.A. Inc. | System and method for characterization with non-unique solutions of anisotropic velocities |
| RU2491580C1 (ru) * | 2012-02-07 | 2013-08-27 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Способ измерения геофизических характеристик с применением последующей инверсии геоэлектрических данных с дополнительным временным фильтром |
| RU2494419C1 (ru) * | 2012-03-11 | 2013-09-27 | Открытое акционерное общество Научно-производственное предприятие "Научно-исследовательский и проектно-конструкторский институт геофизических исследований геологоразведочных скважин" (ОАО НПП "ВНИИГИС") | Способ геоэлектроразведки |
| EA025583B1 (ru) * | 2012-09-05 | 2017-01-30 | Сергей Иванович Говорков | Способ геоэлектроразведки и устройство для его осуществления |
| RU2595327C1 (ru) * | 2015-03-03 | 2016-08-27 | Общество с ограниченной ответственностью "ЛУКОЙЛ-Приморьенефтегаз" | Устройство сейсморазведки 2d или 3d, электроразведки и гис для повышения надежности картирования кровли соли и для прогноза нефтегазоносности подсолевых отложений в районах с развитой солянокупольной тектоникой |
| WO2016176541A1 (en) * | 2015-04-30 | 2016-11-03 | Schlumberger Technology Corporation | Method to estimate water saturation in electromagnetic measurements |
| CN105068130B (zh) * | 2015-07-22 | 2018-02-16 | 中石化石油工程地球物理有限公司江汉分公司 | 全频段电性参数的统一数学表征形式及参数的计算方法 |
| CN105549101B (zh) * | 2016-01-28 | 2018-01-05 | 中国矿业大学 | 一种瞬变电磁数据微分电导解释方法 |
| CN108241781B (zh) * | 2018-01-04 | 2021-07-30 | 国家电网公司 | 一种实验室标准电场发生器的参数扫描优化设计方法 |
| CN110058316B (zh) * | 2019-05-10 | 2020-09-11 | 成都理工大学 | 一种基于电阻率等值性原理的电磁测深约束反演方法 |
| CN110673218A (zh) * | 2019-09-05 | 2020-01-10 | 中国科学院地质与地球物理研究所 | 一种接地导线源瞬变电磁响应中ip信息的提取方法 |
| CN114791633B (zh) * | 2022-06-21 | 2022-09-30 | 西安石油大学 | 监测页岩气压裂的方法、系统及介质 |
| CN115829001B (zh) * | 2022-11-08 | 2023-06-20 | 中国科学院地质与地球物理研究所 | 一种瞬变电磁-激电场分离及多参数信息提取方法和系统 |
| CN116699724B (zh) * | 2023-08-01 | 2023-10-10 | 北京建工环境修复股份有限公司 | 一种时间域激发极化数据质量评价方法、体系及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1122998A1 (ru) * | 1983-06-03 | 1984-11-07 | Всесоюзное Морское Научно-Производственное Объединение Инженерной Геологии "Союзморинжгеология" | Способ морской геоэлектроразведки |
| CA2299743A1 (en) * | 2000-02-22 | 2001-08-22 | Kenneth Duckworth | A method and apparatus for the detection of induced polarization and electrical resistivity in rocks and soils free of the effect of the electromagnetic coupling phenomenon |
| RU2235347C1 (ru) * | 2003-08-11 | 2004-08-27 | РЫХЛИНСКИЙ Николай Иванович | Способ геоэлектроразведки (варианты) |
| CN101166999A (zh) * | 2005-03-24 | 2008-04-23 | 电法勘探有限公司俄罗斯联邦 | 一种使用空间微分法将瞬变场分成若干区段的电法勘探方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2304310C (en) | 1997-10-08 | 2004-11-30 | Shell Canada Limited | Resistivity log correction method |
| DE10006083B4 (de) * | 2000-02-11 | 2004-01-22 | INSTITUT FüR MIKROTECHNIK MAINZ GMBH | Verfahren zur quantitativen und/oder qualitativen Bestimmung von Schichtdicken sowie ein Mikroreaktionsgefäß und eine Titerplatte |
| FR2815124A1 (fr) | 2000-09-30 | 2002-04-12 | Schlumberger Services Petrol | Procede de determination de la saturation en hydrocarbure d'une formation |
| RU2253137C2 (ru) | 2003-07-02 | 2005-05-27 | Федеральное государственное унитарное предприятие Забайкальский комплексный научно-исследовательский институт (ФГУП ЗабНИИ) | Способ анализа геоэлектроразведочного сигнала |
| RU2231089C1 (ru) | 2003-07-08 | 2004-06-20 | РЫХЛИНСКИЙ Николай Иванович | Способ геоэлектроразведки (варианты) |
| RU2284555C1 (ru) | 2005-06-01 | 2006-09-27 | Николай Иванович РЫХЛИНСКИЙ | Способ морской геоэлектроразведки с фокусировкой электрического тока (варианты) |
| US8447521B2 (en) * | 2009-03-13 | 2013-05-21 | Schlumberger Technology Corporation | Subsurface electrical conductivity contrast reconstruction using narrow-band electromagnetic data and fullwave imaging |
-
2008
- 2008-10-23 RU RU2008141889/28A patent/RU2399931C2/ru active
- 2008-11-17 EP EP08877586.1A patent/EP2352041A4/en not_active Withdrawn
- 2008-11-17 CN CN200880132434.XA patent/CN102265188B/zh not_active Expired - Fee Related
- 2008-11-17 WO PCT/RU2008/000711 patent/WO2010047613A1/ru active Application Filing
- 2008-11-17 MX MX2011002163A patent/MX2011002163A/es active IP Right Grant
-
2009
- 2009-10-19 US US12/588,507 patent/US8138763B2/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1122998A1 (ru) * | 1983-06-03 | 1984-11-07 | Всесоюзное Морское Научно-Производственное Объединение Инженерной Геологии "Союзморинжгеология" | Способ морской геоэлектроразведки |
| CA2299743A1 (en) * | 2000-02-22 | 2001-08-22 | Kenneth Duckworth | A method and apparatus for the detection of induced polarization and electrical resistivity in rocks and soils free of the effect of the electromagnetic coupling phenomenon |
| RU2235347C1 (ru) * | 2003-08-11 | 2004-08-27 | РЫХЛИНСКИЙ Николай Иванович | Способ геоэлектроразведки (варианты) |
| CN101166999A (zh) * | 2005-03-24 | 2008-04-23 | 电法勘探有限公司俄罗斯联邦 | 一种使用空间微分法将瞬变场分成若干区段的电法勘探方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102265188A (zh) | 2011-11-30 |
| US20100164500A1 (en) | 2010-07-01 |
| MX2011002163A (es) | 2011-10-06 |
| US8138763B2 (en) | 2012-03-20 |
| WO2010047613A1 (ru) | 2010-04-29 |
| EP2352041A1 (en) | 2011-08-03 |
| EP2352041A4 (en) | 2017-08-02 |
| RU2008141889A (ru) | 2010-04-27 |
| RU2399931C2 (ru) | 2010-09-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102265188B (zh) | 电磁感应和诱导极化效应的量化分解方法 | |
| Fiandaca et al. | Time-domain-induced polarization: Full-decay forward modeling and 1D laterally constrained inversion of Cole-Cole parameters | |
| Holtham et al. | Three-dimensional inversion of ZTEM data | |
| CN101166999B (zh) | 一种使用空间微分法将瞬变场分成若干区段的电法勘探方法 | |
| CN110673218A (zh) | 一种接地导线源瞬变电磁响应中ip信息的提取方法 | |
| CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| Behroozmand et al. | A comprehensive study of parameter determination in a joint MRS and TEM data analysis scheme | |
| Mozaffari et al. | 2.5 D crosshole GPR full-waveform inversion with synthetic and measured data | |
| Grombacher et al. | Steady‐State Surface NMR for Mapping of Groundwater | |
| Dalgaard et al. | Enhancing SNMR model resolution by selecting an optimum combination of pulse moments, stacking, and gating | |
| CN108169802A (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
| Yang et al. | 3-D DC resistivity modelling with arbitrary long electrode sources using finite element method on unstructured grids | |
| Kalscheuer et al. | A non-linear truncated SVD variance and resolution analysis of two-dimensional magnetotelluric models | |
| Mukanova et al. | The boundary element method in geophysical survey | |
| Wang et al. | 3-D parallel inversion of multichannel transient electromagnetic data using a moving footprint | |
| Danilovskiy et al. | Deep-learning-based noniterative 2D-inversion of unfocused lateral logs | |
| Cozzolino et al. | Synthetic focusing and simulation of dual laterolog tool in axisymmetric subsurface models | |
| Kang et al. | Using logging while drilling resistivity imaging data to quantitatively evaluate fracture aperture based on numerical simulation | |
| CN115929292A (zh) | 基于成像测井和常规测井融合的测井方法、装置及设备 | |
| Kang et al. | Using the multi-component induction logging data to evaluate the geometric parameters of the hydraulic fracture | |
| Stephan et al. | Adding realistic noise models to synthetic ground‐penetrating radar data | |
| Meng et al. | Time-domain induced polarization forward modeling and an apparent spectral parameter solution based on the Weibull growth model | |
| Li et al. | 2D cross-hole electromagnetic inversion algorithms based on regularization algorithms | |
| Hernández Contreras et al. | A computational package to compute the electrical resistivity tomography response for regular bodies immersed in a homogeneous half-space | |
| Zou et al. | A rapid NMR T2 inversion method based on norm smoothing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140730 Termination date: 20181117 |