CN112083506A - Method for detecting interface orientation of medium layer based on fractional turn coil - Google Patents

Method for detecting interface orientation of medium layer based on fractional turn coil Download PDF

Info

Publication number
CN112083506A
CN112083506A CN202010964345.4A CN202010964345A CN112083506A CN 112083506 A CN112083506 A CN 112083506A CN 202010964345 A CN202010964345 A CN 202010964345A CN 112083506 A CN112083506 A CN 112083506A
Authority
CN
China
Prior art keywords
receiving coil
coil
electromotive force
electric field
medium
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.)
Withdrawn
Application number
CN202010964345.4A
Other languages
Chinese (zh)
Inventor
唐章宏
许月晨
王芬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Weizhi Jiachen Technology Development Co ltd
Original Assignee
Beijing Weizhi Jiachen Technology Development Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Weizhi Jiachen Technology Development Co ltd filed Critical Beijing Weizhi Jiachen Technology Development Co ltd
Priority to CN202010964345.4A priority Critical patent/CN112083506A/en
Publication of CN112083506A publication Critical patent/CN112083506A/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G01V3/28Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Mining & Mineral Resources (AREA)
  • Electromagnetism (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a method for detecting the interface orientation of a medium layer based on a fractional turn coil, which comprises the steps of adopting an axisymmetric winding round transmitting coil and an asymmetrically wound receiving coil which is coaxial and parallel with the transmitting coil, wherein the receiving coil can completely eliminate the direct coupling electromotive force generated by the transmitting coil to the receiving coil under an infinite uniform medium and only receives the signal of a transmitting signal reflected by a stratum interface; calculating a reflection field received by the receiving coil by adopting a mirror image method, obtaining a total electric field at the receiving coil according to a field superposition theory, and calculating induced electromotive force of the receiving coil by adopting Gaussian integration; the method has the advantages that high-precision measurement of the layer interface position is achieved according to the induced electromotive force amplitude and the phase angle change rule, the reflection field at the layer interface is calculated by adopting a mirror image method, the non-axisymmetric model is converted into the axisymmetric model, the calculated model can adopt a fast algorithm of the axisymmetric model, and the method has important significance for achieving high-precision measurement of the layer interface position and constantly inverting the information of the medium layer.

Description

Method for detecting interface orientation of medium layer based on fractional turn coil
Technical Field
The invention relates to the technical field of well logging, in particular to a method for detecting the interface orientation of a medium layer based on a fractional turn coil.
Background
With the rapid development of detection technology, remote detection and imaging systems have become a hot point problem, and especially in the fields of airplanes, radars, underwater submarines, geological exploration and the like, the development of imaging and medium detection of three-dimensional vector fields becomes very important research content, especially in complex environments, signals received by a measurement system are seriously interfered by background signals, and large errors exist when distance and direction information of an object/medium are remotely detected, so that a rapid response system cannot accurately identify the detected object. In the process of remote edge detection, the reflected signal of the medium layer interface is very weak, most of the signals received by the receiver come from background signals, serious misjudgment and missed judgment occur in the process of edge detection, the distance and the azimuth information of the layer interface cannot be effectively detected, and great challenges are met in the aspects of measuring the layer interface distance, electromagnetic information imaging and the like.
In the petroleum exploration and development engineering, the logging instrument can effectively detect the apparent conductivity of the layered stratum in a vertical well, the true conductivity of the stratum can be obtained through inversion, and the detection technology for the position of the layer interface is slowly developed. For the detection of the position of the layer interface, the traditional logging instrument can not completely eliminate the direct-coupled electromotive force in a medium due to the fact that the shielding coil can not completely eliminate the direct-coupled electromotive force, so that the proportion of direct-coupled signals in induced electromotive force generated by the receiving coil is large, the proportion of reflected signals is small, particularly when the layer interface is far, the layer interface reflected signals are seriously attenuated, the direction and the distance of the layer interface can not be effectively analyzed through the amplitude phase angle of the induced electromotive force, and the conductivity of the stratum can be solved through the induced electromotive force.
On the other hand, due to the existence of the layer interface, the layer interface reflects the electromagnetic wave emitted by the instrument transmitting coil, so that the induced electromotive force generated by the traditional circular symmetric receiving coil can generate a phenomenon of signal attenuation.
Aiming at the problems, the invention provides a method for detecting the orientation of a medium layer interface based on a fractional turn coil, which adopts an asymmetrically wound receiving coil, can effectively analyze the intensity of an electric field through rotation, realizes the measurement of the electric fields with different intensities at the receiving coil, and further realizes the orientation measurement of the layer interface through the amplitude of induced electromotive force.
Disclosure of Invention
Objects of the invention
In order to overcome at least one defect in the prior art, the invention provides a method for detecting the interface orientation of a medium layer based on a fractional turn coil. The receiving coil wound asymmetrically is adopted, the intensity of electric field can be effectively analyzed through rotation, electric fields with different intensities at the receiving coil are measured, and then the azimuth measurement of a layer interface is realized through the amplitude of induced electromotive force.
(II) technical scheme
As a first aspect of the invention, the invention discloses a method for detecting the interface orientation of a medium layer based on a fractional turn coil, which comprises the following steps:
step 1, completely eliminating direct coupling electromotive force generated by a transmitting coil to a receiving coil under an infinite uniform medium according to the receiving coil;
step 2, calculating a reflection field of the transmitting coil at the receiving coil by adopting a mirror image method according to the receiving coil and the transmitting coil;
step 3, calculating the total electric field generated by the transmitting coil at the receiving coil by adopting a field superposition theory according to the reflected field;
step 4, rotating the receiving coil according to the total electric field generated at the receiving coil and calculating induced electromotive force when the receiving coil rotates at different angles by adopting Gaussian integration;
and 5, measuring the layer interface orientation according to the induced electromotive force through the amplitude of the induced electromotive force and the change rule of the phase angle.
In a possible implementation manner, in the step 1, the transmitting coil is a circular transmitting coil wound in an axisymmetrical manner; the receiving coil is a fractional turn far detection receiving coil which is coaxial and parallel with the transmitting coil and is wound asymmetrically.
In a possible implementation manner, in the step 1, the asymmetrically wound fractional turn far detection receiving coil is specifically represented by:
a disc for winding a coil is equally divided into even fan-shaped areas passing through the center of a circle, and the even fan-shaped areas are clockwise arranged from 1 to N; n > = 4;
based on the sector areas, each turn of coil is divided into N-1 fractional sub-turns, each fractional sub-turn is wound on one sector area, the winding directions of the fractional sub-turns of two adjacent sector areas are opposite, and at the moment, the sector area N is empty;
based on the fractional sub-turns, the sub-turn coefficient ratio of the sector region from 1 to N-1 winding fractional sub-turns satisfies:
Figure 100002_DEST_PATH_IMAGE001
and is
Figure 479558DEST_PATH_IMAGE002
Wherein the content of the first and second substances,Nthe number of the sector areas;J iis the sector areaiThe corresponding sub-turn coefficient;
based on the fractional sub-turns, each fractional sub-turn comprises a circumferential arc part and a radial straight line part, the circumferential arc part is wound on the circumferential part of the sector area, and the radial straight line part is wound on the radial tangent plane part passing through the center of a circle of the sector area.
In a possible implementation manner, in step 1, the completely eliminating, by the receiving coil, the direct-coupled electromotive force generated by the transmitting coil to the receiving coil under an infinite homogeneous medium is specifically:
the transmitting coil and the receiving coil are coaxial and parallel, the transmitting coil generates an electric field in a concentric circle at the receiving coil, the integral result of the electric field along the receiving coil is zero, and the sum of magnetic fluxes passed by the receiving coil is zero, namely the induced electromotive force generated by the receiving coil is zero.
In a possible implementation manner, in step 2, the calculating the reflected field of the transmitting coil at the receiving coil by using a mirror image method specifically includes:
because the existence of the second layer medium causes the medium model not to be axisymmetric, the mirror image method eliminates the influence of the second layer medium, so that the calculated model adopts the algorithm of an axisymmetric model to realize the process of calculating the reflection field.
In a possible implementation manner, in the step 2, a dummy transmitting coil generated by the mirror image method is defined as a mirror image source, and a real transmitting coil is defined as a real source; the relational expression of the mirror image source and the real source is as follows:
Figure 100002_DEST_PATH_IMAGE003
wherein
Figure 508694DEST_PATH_IMAGE004
Is the current of the mirror image source,
Figure 100002_DEST_PATH_IMAGE005
is the current of the real source of the current,
Figure 952445DEST_PATH_IMAGE006
the dielectric layer interface reflection coefficient;
Figure 371925DEST_PATH_IMAGE006
the expression of (a) is:
Figure 100002_DEST_PATH_IMAGE007
Figure 300436DEST_PATH_IMAGE008
Figure 100002_DEST_PATH_IMAGE009
wherein
Figure 715236DEST_PATH_IMAGE010
And
Figure 100002_DEST_PATH_IMAGE011
is medium magnetic permeability;
Figure 329888DEST_PATH_IMAGE012
and
Figure DEST_PATH_IMAGE013
is a dielectric complex dielectric constant expressed as
Figure 813828DEST_PATH_IMAGE014
Figure DEST_PATH_IMAGE015
Figure 31183DEST_PATH_IMAGE016
Figure DEST_PATH_IMAGE017
Is a dielectric constant of a medium, and is,
Figure 51223DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE019
in order to be the conductivity of the medium,
Figure 335311DEST_PATH_IMAGE020
is the angular frequency of the electromagnetic wave,
Figure DEST_PATH_IMAGE021
as the incident angle of the electromagnetic wave,jin the form of imaginary unit, the imaginary part,
Figure 57280DEST_PATH_IMAGE022
=-1。
in a possible implementation, in the step 3, the field superposition theory is used to calculate the total electric field generated by the real source and the mirror source at the receiving coil after the mirror image method is used.
In a possible implementation manner, in step 4, the rotating the receiving coil and calculating the induced electromotive force when the receiving coil rotates by different angles by using gaussian integration specifically include:
calculating the induced electromotive force of the receiving coil using the Gaussian integration because the total electric field generated at the receiving coil is not an axisymmetric electric field due to the presence of a layer interface; and rotating the receiving coil, and calculating induced electromotive force when the receiving coil rotates at different angles according to the asymmetry of the receiving coil and the asymmetry of the electric field.
In a possible implementation manner, in the step 4, the induced electromotive force is calculated by using the gaussian integral, specifically represented as:
Figure 750429DEST_PATH_IMAGE024
Vis an induced electromotive force of the receiving coil,Efor the strength of the electric field,dlintegration units for the receiving coil, converted with respect to angleθFunction of (2)
Figure DEST_PATH_IMAGE025
Wherein
Figure 484030DEST_PATH_IMAGE026
For the receiving coil to be at an angleθThe vector mode of the tangential component of the electric field strength,Ris the radius of the coil, here,tdenoted as receiving coil atθThe tangential direction of (A);
using integral transformation
Figure DEST_PATH_IMAGE027
The result after transformation is
Figure 814386DEST_PATH_IMAGE028
Figure DEST_PATH_IMAGE029
For integral variables, using Gaussian integration
Figure 23650DEST_PATH_IMAGE030
Wherein
Figure DEST_PATH_IMAGE031
Is a weight factor of the gaussian integral,
Figure 864699DEST_PATH_IMAGE032
is the vector mode of the tangential component of the electric field at the gaussian point that needs to be transformed to the electric field at the actual coordinates, i.e. the
Figure DEST_PATH_IMAGE033
Wherein
Figure 108598DEST_PATH_IMAGE034
Is the Gaussian point
Figure DEST_PATH_IMAGE035
And the coordinates of the sampling point at the position corresponding to the receiving coil.
In a possible embodiment, the step 5 is specifically:
in the rotation process of the receiving coil, the induced electromotive force which changes along with the rotation angle is generated in a non-uniform electric field by the receiving coil due to an asymmetric structure, a plurality of peak values are formed, and the judgment of the layer interface direction is realized by analyzing the amplitude and the phase angle change of the induced electromotive force.
(III) advantageous effects
The invention provides a method for detecting the interface orientation of a medium layer based on a fractional turn coil, which is characterized in that a novel asymmetrically wound far detection receiving coil is designed, a coaxial and parallel arrangement mode of a transmitting coil and the receiving coil is adopted, so that the transmitting coil generates a concentric electric field at the receiving coil, the integral result of the electric field along the receiving coil is zero, namely the induced electromotive force generated by the receiving coil is zero, the direct coupling electromotive force generated by the transmitting coil under an infinite uniform medium is completely eliminated in the medium, and the direct coupling signal has zero interference on a reflected signal; because the existence of the second layer of medium causes the medium model not to be axisymmetric, the mirror image method can eliminate the influence of the second layer of medium, so that the calculated model can adopt the fast algorithm of the axisymmetric model to realize fast and high-precision calculation of the reflection field; calculating the total electric field generated by the real source and the mirror image source at the receiving coil after using the mirror image method by adopting a field superposition theory, and analyzing the influence of the layer interface on the electric field at the receiving coil through the total electric field at the receiving coil; the total electric field at the receiving coil is not an axisymmetric electric field due to the existence of the layer interface, induced electromotive force generated when the receiving coil rotates at different angles is calculated by rotating the receiving coil through Gaussian integration, and the layer interface azimuth information is judged through the variation rule of the amplitude and the phase angle of the induced electromotive force, so that the electromagnetic information of the layer interface is rapidly and highly accurately calculated.
Drawings
The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining and illustrating the present invention and should not be construed as limiting the scope of the present invention.
Fig. 1 is a flowchart of a method for detecting the orientation of a medium layer interface based on a fractional turn coil according to the present invention.
Fig. 2 is a model schematic diagram of a novel asymmetric-winding far-detection circular receiving coil (N = 4) based on a fractional-turn coil method for detecting the interface orientation of a medium layer provided by the invention.
FIG. 3 is a schematic diagram of a constructed medium model of a method for detecting the interface orientation of a medium layer based on a fractional turn coil provided by the invention.
Fig. 4 is a schematic diagram of a medium model without considering the medium 2 of the method for detecting the interface orientation of the medium layer based on the fractional turn coil provided by the invention.
Fig. 5 is a schematic model diagram of a transmission coil after a medium 2 is considered in the method for detecting the interface orientation of the medium layer based on the fractional turn coil provided by the invention.
Fig. 6 is a distribution diagram of the total electric field at the receiving coil obtained by the field superposition theory of the method for detecting the interface orientation of the medium layer based on the fractional turn coil provided by the invention.
Fig. 7 is a schematic diagram of changes of induced electromotive force amplitude of the receiving coil and a rotation angle of the receiving coil in the method for detecting the interface orientation of the dielectric layer based on the fractional turn coil according to the present invention.
Fig. 8 is a schematic diagram illustrating changes of a phase angle of induced electromotive force of a receiving coil and a rotation angle of the receiving coil in the method for detecting the interface orientation of the dielectric layer based on the fractional turn coil according to the present invention.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in more detail below with reference to the accompanying drawings in the embodiments of the present invention.
It should be noted that: in the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described are some embodiments of the present invention, not all embodiments, and features in embodiments and embodiments in the present application may be combined with each other without conflict. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it is to be understood that the terms "central," "longitudinal," "lateral," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like are used in the orientation or positional relationship indicated in the drawings, which are used for convenience in describing the invention and for simplicity in description, and are not intended to indicate or imply that the referenced device or element must have a particular orientation, be constructed and operated in a particular orientation, and are not to be considered limiting of the scope of the invention.
A first embodiment of a method for detecting the orientation of a medium layer interface based on a fractional turn coil according to the present invention is described in detail below with reference to fig. 1-8. As shown in fig. 1 to 8, the method for detecting the interface orientation of the medium layer provided in this embodiment mainly includes: step 1, step 2, step 3, step 4 and step 5.
With the wide use of inclined wells and horizontal wells, the existence of layer interfaces causes that a computed logging model is not axisymmetric any more, and for computing the azimuth distance of the layer interfaces, a three-dimensional numerical computation method is needed, so that the computation precision is not high and the computation amount is large, thereby increasing the logging response time. The invention adopts a mirror image method to calculate the reflection field of the transmitting coil at the layer interface aiming at the inclined well and the horizontal well, converts the three-dimensional non-axisymmetric model into the axisymmetric model, greatly reduces the calculated amount, and realizes the measurement of the layer interface and the stratum conductivity through the electromagnetic wave phase angle and the attenuation rule. The invention enables the logging instrument to detect the position and the orientation of the boundary of the medium layer in the reservoir, and the logging instrument advances along the parallel direction of the reservoir, thereby accurately evaluating the size of the oil and gas reservoir.
Step 1, completely eliminating direct coupling electromotive force generated by a transmitting coil to a receiving coil under an infinite uniform medium according to the receiving coil;
in the step 1, the transmitting coil is a circular transmitting coil wound in an axisymmetrical manner; the receiving coil is a novel fractional-turn far detection receiving coil which is wound asymmetrically and is coaxial and parallel to the transmitting coil.
Step 2, calculating a reflection field of the transmitting coil at the receiving coil by adopting a mirror image method according to the receiving coil and the transmitting coil;
step 3, calculating the total electric field generated by the transmitting coil at the receiving coil by adopting a field superposition theory according to the reflected field;
step 4, rotating the receiving coil according to the total electric field generated at the receiving coil and calculating induced electromotive force when the receiving coil rotates at different angles by adopting Gaussian integration;
and 5, realizing high-precision measurement of the layer interface orientation by the amplitude of the induced electromotive force and the change rule of the phase angle according to the induced electromotive force.
In step 1, the asymmetrically wound fractional-turn far detection receiving coil is specifically represented by a model:
a disc for winding a coil is equally divided into even fan-shaped areas passing through the center of a circle, and the even fan-shaped areas are clockwise arranged from 1 to N; n > = 4;
based on the sector areas, each turn of coil is divided into N-1 fractional sub-turns, each fractional sub-turn is wound on one sector area, the winding directions of the fractional sub-turns of two adjacent sector areas are opposite, and at the moment, the sector area N is empty;
based on the fractional sub-turns, the sub-turn coefficient ratio of the fractional sub-turns wound in the sector areas from 1 to N-1 is
Figure DEST_PATH_IMAGE037
And the ratio satisfies the equation:
Figure 937752DEST_PATH_IMAGE038
and is
Figure DEST_PATH_IMAGE039
Wherein the content of the first and second substances,Nthe number of the sector areas;J iis the sector areaiThe coefficient of the corresponding sub-turn is,J iis a fraction;
on the basis of the fractional sub-turns, a sub-turn lead divides the fractional sub-turns into a circumferential arc part and a radial straight line part, each fractional sub-turn comprises a circumferential arc part and a radial straight line part, the circumferential arc part is wound on the circumferential part of the sector area, and the radial straight line part is wound on the radial tangent plane part of the sector area passing through the center of a circle; for the purpose of intuitively describing a model of coil winding completion, fig. 2 shows a model schematic diagram of a receiving coil when N =4, and the sub-turn coefficients of the sector areas 1 to 3 are 0.25, 0.5, and 0.25, respectively.
In step 1, the completely eliminating, according to the receiving coil, the direct-coupled electromotive force generated by the transmitting coil to the receiving coil under the infinite homogeneous medium includes: the transmitting coil and the receiving coil are coaxial and parallel, the transmitting coil generates an electric field in a concentric circle at the receiving coil, the integral result of the electric field along the receiving coil is zero, and the sum of magnetic fluxes passed by the receiving coil is zero, namely the induced electromotive force generated by the receiving coil is zero.
According to the asymmetric winding receiving coil, the direct coupling electromotive force of the transmitting coil is eliminated, and the method is specifically expressed as follows:
in an electric field generated by a uniform electric field or a coaxial parallel circular transmitting coil in a uniform medium, magnetic fluxes passing through N fan-shaped areas are all equal, and according to the winding direction of a receiving coil, the induced electromotive force received by the receiving coil is 0 under the condition;
the electric field generated by the coaxial parallel circular transmitting coil in the uniform medium is concentric on the plane of the receiving coil, so that all radial straight line parts of the receiving coil passing through the center of a circle are always vertical to the direction of the electric field, all circumferential arc parts on the circumference are always parallel to the direction of the electric field, and according to the winding direction of the receiving coil, the induced electromotive force received by the receiving coil is 0 under the condition;
under the condition that a stratum interface exists near the receiving coil, the interface reflects electromagnetic waves generated by the transmitting coil, the electromagnetic waves are relative to an electric field generated by a non-uniform or non-concentric circle of the receiving coil, and therefore induced electromotive force is generated on the receiving coil.
In the step 2, in the horizontal well, the medium model is a non-axisymmetric model, and for convenience of calculation, the medium model is configured as shown in fig. 3; the reflection field of the transmitting coil at the receiving coil is calculated by using a mirror image method (i.e. the reflection field of the transmitting coil at the layer interface is calculated by using the mirror image method), and an equivalent model of the mirror image method is shown in fig. 4 and 5, and specifically includes: because the existence of the second layer of medium causes the medium model not to be axisymmetric, the mirror image method eliminates the influence of the second layer of medium, so that the calculated model can adopt a fast algorithm of an axisymmetric model to realize fast and high-precision calculation of the reflection field.
In the step 2, a virtual transmitting coil generated by adopting the mirror image method is defined as a mirror image source, and a real transmitting coil is defined as a real source;
because the existence of the second layer of medium causes the (horizontal well) medium model not to be axisymmetric, the mirror image method can eliminate the influence of the second layer of medium, and the specific method is as follows:
to visualize the application of the mirror method to logging while drilling, a dummy transmitter coil (i.e., the mirror transmitter coil in fig. 5) is introduced into the medium 2, which is mirror-symmetric to the transmitter coil in the medium 1, and is defined as a mirror source, while a real transmitter coil is defined as a real source. The regions in medium 2 were completely replaced by the medium of medium 1.
The relational expression of the mirror image source and the real source is as follows:
Figure 385045DEST_PATH_IMAGE040
wherein
Figure DEST_PATH_IMAGE041
Is the current of the mirror image source,
Figure 13472DEST_PATH_IMAGE042
is the current of the real source of the current,
Figure DEST_PATH_IMAGE043
is the interface reflection coefficient of medium 1 and medium 2;
Figure 361146DEST_PATH_IMAGE043
the expression of (a) is:
Figure 49616DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE045
Figure 249784DEST_PATH_IMAGE046
wherein
Figure DEST_PATH_IMAGE047
And
Figure 681903DEST_PATH_IMAGE048
the magnetic permeability of medium 1 and medium 2;
Figure DEST_PATH_IMAGE049
and
Figure 618504DEST_PATH_IMAGE050
is the complex dielectric constant of medium 1 and medium 2, and has the expression
Figure DEST_PATH_IMAGE051
Figure 87662DEST_PATH_IMAGE052
Figure DEST_PATH_IMAGE053
Figure 493236DEST_PATH_IMAGE054
Is the dielectric constant of medium 1 and medium 2,
Figure DEST_PATH_IMAGE055
Figure 712733DEST_PATH_IMAGE056
which is the electrical conductivity of the medium 1 and the medium 2,
Figure DEST_PATH_IMAGE057
being angular frequency of electromagnetic waves,
Figure 520152DEST_PATH_IMAGE058
As the incident angle of the electromagnetic wave,jin the form of imaginary unit, the imaginary part,
Figure DEST_PATH_IMAGE059
=-1。
to further illustrate the calculation process, the calculation model shown in fig. 3 is used for demonstration: setting the relative dielectric constants of the medium 1 and the medium 2 to be 5, and the relative magnetic conductivities of the medium 1 and the medium 2 to be 1; the conductivity of the medium 1 is 0.01S/m, and the conductivity of the medium 2 is 0.5S/m; the circle center coordinates of the transmitting coil (0, 0, 0) and the receiving coil (0, 0, 0.5), the radiuses of the transmitting coil and the receiving coil are both 0.03m, the number of turns of the transmitting coil is 100, and the number of turns of the receiving coil is 50; the layer interface distance is 1m from the z axis; the current of the transmitting coil is 1A sine alternating current with the frequency of 100 KHz.
The non-axisymmetric model adopts a mirror image method, so that the calculated model can adopt a quick algorithm of the axisymmetric model to realize quick and high-precision calculation of the reflection field.
In the step 3, the total electric field generated by the real source and the mirror source at the receiving coil after the mirror image method is used is calculated by adopting the field superposition theory. The effect of the layer interface on the electric field at the receiver coil can be analyzed by the total electric field at the receiver coil, which is shown in fig. 6.
In step 4, the rotating the receiving coil and calculating the induced electromotive force when the receiving coil rotates by different angles by using gaussian integration specifically include: the total electric field generated at the receiving coil is not an axisymmetric electric field due to the existence of a layer interface, the induced electromotive force generated by the asymmetric receiving coil is nonzero, and the induced electromotive force of the receiving coil is calculated by adopting the Gaussian integral; and rotating the receiving coil, and calculating the induced electromotive force of the receiving coil when the receiving coil rotates at different angles according to the asymmetry of the receiving coil and the asymmetry of the electric field, so as to realize the rapid and high-precision calculation of the induced electromotive force of the receiving coil in the non-uniform electric field.
In step 4, the specific gaussian integration method is as follows (only the algorithm of gaussian integrals of circumferential arc sections 0 to pi/2 of the model of the N =4 receiving coil is given here, which is the same as the algorithm of gaussian integrals of other circumferential arc sections and radial straight line sections, and is not stated again):
calculating the induced electromotive force by using the Gaussian integral is specifically expressed as:
Figure 301158DEST_PATH_IMAGE061
Vis an induced electromotive force of the receiving coil,Efor the strength of the electric field,dlintegration units for the receiving coil, converted with respect to angleθFunction of (2)
Figure DEST_PATH_IMAGE062
Wherein
Figure 567929DEST_PATH_IMAGE063
For the receiving coil to be at an angleθThe vector mode of the tangential component of the electric field strength,Ris the radius of the coil, here,tdenoted as receiving coil atθThe tangential direction of (A);
using integral transformation
Figure DEST_PATH_IMAGE064
The result after transformation is
Figure 341850DEST_PATH_IMAGE065
Figure DEST_PATH_IMAGE066
For integral variables, using Gaussian integration
Figure 754507DEST_PATH_IMAGE067
Wherein
Figure 221261DEST_PATH_IMAGE031
Is a weight factor of the gaussian integral,
Figure DEST_PATH_IMAGE068
is the vector mode of the tangential component of the electric field at the gaussian point that needs to be transformed to the electric field at the actual coordinates, i.e. the
Figure 116273DEST_PATH_IMAGE069
Wherein
Figure DEST_PATH_IMAGE070
Is the Gaussian point
Figure 428306DEST_PATH_IMAGE071
Corresponding to the coordinates of the sampling point at the actual receiving coil.
Wherein, the step 5 specifically comprises the following steps: the method is characterized in that a receiving coil is wound asymmetrically, the receiving coil generates induced electromotive force which changes along with a rotation angle in a non-uniform electric field due to an asymmetric structure in the rotation process of the receiving coil, a plurality of peak values are formed, and the judgment of the layer interface direction is realized by analyzing the amplitude and phase angle change of the induced electromotive force.
The induced electromotive force generated by the sending coil under the infinite uniform medium can be completely eliminated by adopting the receiving coil wound asymmetrically, so that signals received by the receiving coil all come from reflected waves, and the direction of the layer interface is calculated through the amplitude and the phase angle change rule of the induced electromotive force. FIG. 7 is an image of the induced electromotive force amplitude change of 360 ° rotation of the receiving coil when the layer interface is 1m from the z-axis; FIG. 8 is an image of the phase angle change of the induced electromotive force induced by a 360 ° rotation of the receiving coil when the layer interface is at a distance of 1m from the z-axis.
As shown in fig. 7, for different rotation angles, 4 peaks of induced electromotive force exist, wherein the first and third peaks are two parts with less reverse winding, the maximum peak is the part with the most forward winding turns, and the minimum peak is the part without winding the receiving coil, and by rotating the receiving coil, the position of the layer interface can be effectively judged according to the peak change of the induced electromotive force; when N =4, the peak value is maximum when the receiving coil rotates 135 ° as judged by the amplitude of the coil, and the angular bisector of the sector region 2 is the x axis at this time, that is, the direction in which the electric field is maximum is the positive direction of the x axis, so the layer interface is in the positive direction of the x axis.
When the receiving coil has N sectors, it is assumed that sector 1 is located
Figure DEST_PATH_IMAGE072
Rotating the receiving coil within 0 deg. to obtain different peak values of induced electromotive force, and obtaining the rotation angle when the peak value is maximum by image, i.e. the maximum position of electric field is the rotation of angular bisector of sector region N/2
Figure 695471DEST_PATH_IMAGE073
The position of the layer interface can be obtained, and the obtained direction of the layer interface is
Figure DEST_PATH_IMAGE074
From fig. 8, it can be obtained that the phase angle of the induced electromotive force in the rotation process of the instrument jumps between two constant positive and negative values, the phase angle corresponding to the part with the most winding turns keeps the angle the largest, and the phase angles corresponding to the other three peaks keep the angle smaller; the azimuth angle of the layer interface can be judged according to the rotating angle corresponding to the phase angle, so that the high-precision measurement of the azimuth information of the layer interface can be realized by analyzing the amplitude of the induced electromotive force of the receiving coil and the change rule of the phase angle.
The invention relates to a method for detecting the interface orientation of a medium layer based on a fractional turn coil, which is characterized in that a novel asymmetrically wound far detection receiving coil is designed, a coaxial and parallel arrangement mode of a transmitting coil and the receiving coil is adopted, so that the transmitting coil generates a concentric electric field at the receiving coil, the integral result of the electric field along the receiving coil is zero, namely the induced electromotive force generated by the receiving coil is zero, the direct coupling electromotive force generated by the transmitting coil under an infinite uniform medium is completely eliminated in the medium, and the direct coupling signal has zero interference on a reflected signal; because the existence of the second layer of medium causes the medium model not to be axisymmetric, the mirror image method can eliminate the influence of the second layer of medium, so that the calculated model can adopt the fast algorithm of the axisymmetric model to realize fast and high-precision calculation of the reflection field; calculating the total electric field generated by the real source and the mirror image source at the receiving coil after using the mirror image method by adopting a field superposition theory, and analyzing the influence of the layer interface on the electric field at the receiving coil through the total electric field at the receiving coil; the total electric field at the receiving coil is not an axisymmetric electric field due to the existence of the layer interface, induced electromotive force generated when the receiving coil rotates at different angles is calculated by rotating the receiving coil through Gaussian integration, and the layer interface azimuth information is judged through the variation rule of the amplitude and the phase angle of the induced electromotive force, so that the electromagnetic information of the layer interface is rapidly and highly accurately calculated.
The above description is only for the specific embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (10)

1. A method for detecting the interface orientation of a medium layer based on a fractional turn coil is characterized by comprising the following steps:
step 1, completely eliminating direct coupling electromotive force generated by a transmitting coil to a receiving coil under an infinite uniform medium according to the receiving coil;
step 2, calculating a reflection field of the transmitting coil at the receiving coil by adopting a mirror image method according to the receiving coil and the transmitting coil;
step 3, calculating the total electric field generated by the transmitting coil at the receiving coil by adopting a field superposition theory according to the reflected field;
step 4, rotating the receiving coil according to the total electric field generated at the receiving coil and calculating induced electromotive force when the receiving coil rotates at different angles by adopting Gaussian integration;
and 5, measuring the layer interface orientation according to the induced electromotive force through the amplitude of the induced electromotive force and the change rule of the phase angle.
2. The method for detecting the interface orientation of the medium layer according to claim 1, wherein in the step 1, the transmitting coil is a circular transmitting coil wound with axial symmetry; the receiving coil is a fractional turn far detection receiving coil which is coaxial and parallel with the transmitting coil and is wound asymmetrically.
3. The method for detecting the interface orientation of a medium layer as claimed in claim 2, wherein in the step 1, the asymmetrically wound fractional turn far detection receiving coil is specifically represented by:
a disc for winding a coil is equally divided into even fan-shaped areas passing through the center of a circle, and the even fan-shaped areas are clockwise arranged from 1 to N; n > = 4;
based on the sector areas, each turn of coil is divided into N-1 fractional sub-turns, each fractional sub-turn is wound on one sector area, the winding directions of the fractional sub-turns of two adjacent sector areas are opposite, and at the moment, the sector area N is empty;
based on the fractional sub-turns, the sub-turn coefficient ratio of the sector region from 1 to N-1 winding fractional sub-turns satisfies:
Figure DEST_PATH_IMAGE001
and is
Figure 709362DEST_PATH_IMAGE002
Wherein the content of the first and second substances,Nthe number of the sector areas;J iis the sector areaiThe corresponding sub-turn coefficient;
based on the fractional sub-turns, each fractional sub-turn comprises a circumferential arc part and a radial straight line part, the circumferential arc part is wound on the circumferential part of the sector area, and the radial straight line part is wound on the radial tangent plane part passing through the center of a circle of the sector area.
4. The method for detecting the interface orientation of the dielectric layer according to claim 1, wherein in the step 1, the direct-coupled electromotive force generated by the transmitting coil to the receiving coil under the infinite homogeneous medium is completely eliminated according to the receiving coil, and specifically:
the transmitting coil and the receiving coil are coaxial and parallel, the transmitting coil generates an electric field in a concentric circle at the receiving coil, the integral result of the electric field along the receiving coil is zero, and the sum of magnetic fluxes passed by the receiving coil is zero, namely the induced electromotive force generated by the receiving coil is zero.
5. The method for detecting the interface orientation of the dielectric layer according to claim 1, wherein in the step 2, the reflected field of the transmitting coil at the receiving coil is calculated by a mirror image method, specifically:
because the existence of the second layer medium causes the medium model not to be axisymmetric, the mirror image method eliminates the influence of the second layer medium, so that the calculated model adopts the algorithm of an axisymmetric model to realize the process of calculating the reflection field.
6. The method for detecting the interface orientation of the dielectric layer according to claim 5, wherein in the step 2, a dummy transmitting coil generated by the mirror image method is defined as a mirror image source, and a real transmitting coil is defined as a real source; the relational expression of the mirror image source and the real source is as follows:
Figure DEST_PATH_IMAGE003
wherein
Figure 157661DEST_PATH_IMAGE004
Is the current of the mirror image source,
Figure DEST_PATH_IMAGE005
is the current of the real source of the current,
Figure 514562DEST_PATH_IMAGE006
the dielectric layer interface reflection coefficient;
Figure DEST_PATH_IMAGE007
the expression of (a) is:
Figure 50717DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE009
Figure 927406DEST_PATH_IMAGE010
wherein
Figure DEST_PATH_IMAGE011
And
Figure 327032DEST_PATH_IMAGE012
is medium magnetic permeability;
Figure 62907DEST_PATH_IMAGE013
and
Figure DEST_PATH_IMAGE014
is a dielectric complex dielectric constant expressed as
Figure 402752DEST_PATH_IMAGE015
Figure DEST_PATH_IMAGE016
Figure 399527DEST_PATH_IMAGE017
Figure DEST_PATH_IMAGE018
Is a dielectric constant of a medium, and is,
Figure 907738DEST_PATH_IMAGE019
Figure DEST_PATH_IMAGE020
in order to be the conductivity of the medium,
Figure 989963DEST_PATH_IMAGE021
is the angular frequency of the electromagnetic wave,
Figure DEST_PATH_IMAGE022
as the incident angle of the electromagnetic wave,jin the form of imaginary unit, the imaginary part,
Figure 540025DEST_PATH_IMAGE023
=-1。
7. the method for detecting the interface orientation of a medium layer as claimed in claim 6, wherein in the step 3, the total electric field generated at the receiving coil by the real source and the mirror source after the mirror image method is used is calculated by adopting the field superposition theory.
8. The method for detecting the interface orientation of the medium layer as claimed in claim 1, wherein in the step 4, the step of rotating the receiving coil and calculating the induced electromotive forces of the receiving coil when the receiving coil rotates at different angles by using gaussian integration specifically comprises:
calculating the induced electromotive force of the receiving coil using the Gaussian integration because the total electric field generated at the receiving coil is not an axisymmetric electric field due to the presence of a layer interface; and rotating the receiving coil, and calculating induced electromotive force when the receiving coil rotates at different angles according to the asymmetry of the receiving coil and the asymmetry of the electric field.
9. The method for detecting an interface orientation of a dielectric layer as claimed in claim 8, wherein in the step 4, the induced electromotive force is calculated by the gaussian integration as:
Figure 656885DEST_PATH_IMAGE025
Vis an induced electromotive force of the receiving coil,Efor the strength of the electric field,dlintegration units for the receiving coil, converted with respect to angleθFunction of (2)
Figure DEST_PATH_IMAGE026
Wherein
Figure 601576DEST_PATH_IMAGE027
For the receiving coil to be at an angleθThe vector mode of the tangential component of the electric field strength,Ris the radius of the coil, here,tdenoted as receiving coil atθThe tangential direction of (A);
using integral transformation
Figure DEST_PATH_IMAGE028
The result after transformation is
Figure 905519DEST_PATH_IMAGE029
Figure DEST_PATH_IMAGE030
For integral variables, using Gaussian integration
Figure 259271DEST_PATH_IMAGE031
Wherein
Figure DEST_PATH_IMAGE032
Is a weight factor of the gaussian integral,
Figure 542222DEST_PATH_IMAGE033
is a Gaussian point
Figure DEST_PATH_IMAGE034
The vector mode of the tangential component of the electric field at which the electric field at said gaussian point needs to be transformed to the electric field at the actual coordinates, i.e. the electric field at the actual coordinates
Figure 939705DEST_PATH_IMAGE035
Wherein
Figure DEST_PATH_IMAGE036
Is the Gaussian point
Figure 747256DEST_PATH_IMAGE034
And the coordinates of the sampling point at the position corresponding to the receiving coil.
10. The method for detecting the interface orientation of the medium layer according to claim 1, wherein the step 5 is specifically:
in the rotation process of the receiving coil, the induced electromotive force which changes along with the rotation angle is generated in a non-uniform electric field by the receiving coil due to an asymmetric structure, a plurality of peak values are formed, and the judgment of the layer interface direction is realized by analyzing the amplitude and the phase angle change of the induced electromotive force.
CN202010964345.4A 2020-09-15 2020-09-15 Method for detecting interface orientation of medium layer based on fractional turn coil Withdrawn CN112083506A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010964345.4A CN112083506A (en) 2020-09-15 2020-09-15 Method for detecting interface orientation of medium layer based on fractional turn coil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010964345.4A CN112083506A (en) 2020-09-15 2020-09-15 Method for detecting interface orientation of medium layer based on fractional turn coil

Publications (1)

Publication Number Publication Date
CN112083506A true CN112083506A (en) 2020-12-15

Family

ID=73737993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010964345.4A Withdrawn CN112083506A (en) 2020-09-15 2020-09-15 Method for detecting interface orientation of medium layer based on fractional turn coil

Country Status (1)

Country Link
CN (1) CN112083506A (en)

Similar Documents

Publication Publication Date Title
US8055193B2 (en) Underwater remote sensing
US7990153B2 (en) Compensated directional resistivity measurements
US8159227B2 (en) Methods for making directional resistivity measurements
EA007644B1 (en) Method of imaging subsurface formations using a virtual source array
JP2008538042A (en) Pointing vector based method for determining the direction and position of electromagnetic sources
JP4527484B2 (en) Condition detection device
CN111025231B (en) Magnetic induction through-the-earth positioning method based on signal direction
WO2006134329A2 (en) Underwater remote sensing
WO2019245487A1 (en) Operating method of a metal detector capable of measuring target depth
CN109188534A (en) A kind of underwater form of metal detection method and device based on active electric-field principle
CN112083506A (en) Method for detecting interface orientation of medium layer based on fractional turn coil
Liu et al. Synthesized magnetic field focusing for the non-destructive testing of oil and gas well casing pipes using pulsed eddy-current array
CN111929737B (en) Method and device for remotely detecting layer interface position and electromagnetic information
WO2016123803A1 (en) Antenna circuit for logging, antenna system and corresponding logging method
Li et al. Calibration method of azimuth electromagnetic boundary detection logging while drilling based on air–seawater double-layer dielectric as reflection interface
CN111220947A (en) Magnetic induction through-the-earth positioning method based on path loss
CN111929736B (en) Asymmetric fractional turn coil winding method for measuring alternating magnetic field
NO166897B (en) PROCEDURE AND APPARATUS FOR MEASURING THE INTERNAL DIAMETER OF A METAL TUBE.
CN109814157B (en) Method and device for identifying relative gesture and computer storage medium
EP3655769A1 (en) Multi-coil tool for attenuation of motion-induced noise
Liu et al. A study on directional resistivity logging-while-drilling based on self-adaptive hp-FEM
CN113341471A (en) Wide-width electromagnetic induction logging method and device, storage medium and electronic equipment
GB2459218A (en) Underwater signal direction determination
CN112081586B (en) Multi-spectrum stratum boundary far detection method and device
CN109952519B (en) Application of E-field antenna to resistivity logging tool

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20201215