CN106814403A - A kind of method for compensating transient electromagnetic signal negative value - Google Patents

A kind of method for compensating transient electromagnetic signal negative value Download PDF

Info

Publication number
CN106814403A
CN106814403A CN201710031598.4A CN201710031598A CN106814403A CN 106814403 A CN106814403 A CN 106814403A CN 201710031598 A CN201710031598 A CN 201710031598A CN 106814403 A CN106814403 A CN 106814403A
Authority
CN
China
Prior art keywords
negative value
fitting
electromagnetic signal
transient electromagnetic
section
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.)
Granted
Application number
CN201710031598.4A
Other languages
Chinese (zh)
Other versions
CN106814403B (en
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.)
Shanghai Institute of Microsystem and Information Technology of CAS
Original Assignee
Shanghai Institute of Microsystem and Information Technology of CAS
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 Shanghai Institute of Microsystem and Information Technology of CAS filed Critical Shanghai Institute of Microsystem and Information Technology of CAS
Priority to CN201710031598.4A priority Critical patent/CN106814403B/en
Publication of CN106814403A publication Critical patent/CN106814403A/en
Application granted granted Critical
Publication of CN106814403B publication Critical patent/CN106814403B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/38Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention provides a kind of method for compensating transient electromagnetic signal negative value, wherein, the method for the compensation transient electromagnetic signal negative value at least comprises the following steps:Measurement data curve according to transient electromagnetic signal, chooses negative value section as fitting section;E index fitting is carried out to the fitting section, to obtain fitting data;The measurement data is poor with the fitting data, to compensate the transient electromagnetic signal negative value.The method of compensation transient electromagnetic signal negative value of the invention, has the advantages that:Using the negative value method of the invention, it is possible to effectively treatment transient electromagnetic signal, so as to obtain the valid data of long period, geology detecting depth is effectively improved.Through the transient electromagnetic signal after overcompensation, it is possible to process means of interpretation by traditional ripe TEM data and it is explained, improve the accuracy of resistivity explanation.

Description

A kind of method for compensating transient electromagnetic signal negative value
Technical field
The present invention relates to technical field of geological exploration, more particularly to a kind of method for compensating transient electromagnetic signal negative value.
Background technology
Magnetic source transient electromagnetic method (Transient electromagnetic method, abbreviation TEM method) is that one kind has The geophysical prospecting method of effect, TEM is owned by France in time-domain artificial source's electromagnetic method, is with the electric conductivity and magnetic conduction of rock in the earth (ore) Property be physical property premise, according to electromagnetic induction principle observation, research the electromagnetic field room and time regularity of distribution, to find, underground is good to lead A kind of exploration method of ore body or the related geological problem of solution.
The general principle of TEM methods is as follows:It is provided with bi-directional pulse current generation with earth-free loop line or ground connection wire long to swash Generating magnetic field, under the excitation of the electromagnetic field, conductive geologic body is sensed and is produced vortex current.Due to conductive geologic body right and wrong Linear, so pulse current jumps to zero from peak value, one time electromagnetic field disappears immediately, and vortex does not disappear immediately, there is one Individual transient process, the speed of this process is relevant with the electrical parameter of conductor.The electric conductivity of geologic body is better, the thermosteresis of vortex Smaller, transient process is then more long.This vortex transient process, corresponding transient magnetic field, pulse current shut-off phase are formed in space Between in ground observation transient magnetic field, that is, observe secondary electromagnetic field, so that it may find the presence of subsurface anomaly geologic body, so that it is determined that ground The electrical structure and spatial distribution state of lower conductor.
Due to the measurement of TEM methods is the transient process that is vortexed in conductor, and observation carries out during pulse interval, do not deposit In an interference for Electromagnetic Sources, observed parameter is pure secondary electromagnetic field, and being uniquely can be using the side with point device in electromagnetic method Method, detection target coupling is most tight, and the response of acquisition is most strong.Also, magnetic source excitation, is not limited by grounding requirement.
The transient electromagnetic signal that TEM method measurements are obtained be generally the decay of for changing over time approximate e index on the occasion of Curve, but signal reversal development is found in increasing measurement data, i.e. partial transient electromagnetic signal shows as negative value. Existing processing method can be attributed to two kinds:
1st, negative value is directly neglected, with remaining on the occasion of resistivity explanation is carried out, causes that the valid data time is short, geology is visited The higher problem of depth as shallow, resistivity is surveyed, is a kind of means of interpretation without scientific basis.
2nd, it is believed that effect of induced polarization is the main cause for causing transient electromagnetic signal to invert, and will Cole-cole models are incorporated into the middle of TEM method forward modellings, by the equivalent system as resistance composition in parallel with electric capacity of the earth. The theory is set up after VOF is launched, and secondary electromagnetic field charges to polarizable the earth, late period electric discharge, occurs in theory Negative value all very littles, it is difficult to observe in practice, are differed with band Cole-Cole model explanations resistivity out with actual conditions It is very big, and because cole-cole model parameters are more, inverting difficulty is very big, also fails to be applied to measured data explanation.
Based on problem above, a kind of method that can effectively process transient electromagnetic signal negative value is needed badly now, so as to To the valid data of long period, geology detecting depth is improved, and the accuracy of resistivity explanation can be improved.
The content of the invention
The shortcoming of prior art, bears it is an object of the invention to provide one kind compensation transient electromagnetic signal in view of the above The method of value, for solving the problems, such as processing method existing defects in the prior art to transient electromagnetic signal negative value.
In order to achieve the above objects and other related objects, the present invention provides a kind of side for compensating transient electromagnetic signal negative value Method, wherein, the method for the compensation transient electromagnetic signal negative value at least comprises the following steps:
Measurement data curve according to transient electromagnetic signal, chooses negative value section as fitting section;
E index fitting is carried out to the fitting section, to obtain fitting data;
The measurement data is poor with the fitting data, to compensate the transient electromagnetic signal negative value.
Preferably, the starting position as the fitting section at the negative peak of the negative value section is chosen.
Preferably, choose in the negative value section and be located at after the negative peak and the random time point before noise level As the scram position of the fitting section at negative value.
Preferably, the measurement data curve according to transient electromagnetic signal, chooses negative value section as fitting section, specific method For:
Electromagnetic field is launched to the earth by coil, the secondary electromagnetic field signal of feedback is then received, it is described to obtain The measurement data curve of transient electromagnetic signal;Wherein, the earth is equivalent to the circuit system of resistance and capacitances in series, so that The overall magnetic field change curve of the secondary electromagnetic field has identical trend with the charge and discharge conditional curve of the electric capacity;
The induced-current for choosing the secondary electromagnetic field is less than measurement data corresponding during the discharge current of the electric capacity Curved section as fitting section;Wherein, when the induced-current of the secondary electromagnetic field is less than the discharge current of the electric capacity, institute The curved section of corresponding measurement data is the negative value section of the measurement data curve of the transient electromagnetic signal.
Preferably, choose at measurement data corresponding during the discharge current maximum of the electric capacity as the fitting section Starting position;Wherein, when the discharge current of the electric capacity is maximum, corresponding measurement data is that the maximum of the negative value section is born Value.
Preferably, choose from the electric capacity there is maximum discharge current after and random time point institute before noise level As the scram position of the fitting section at corresponding measurement data.
Preferably, the electric discharge magnetic field change curve that the discharge current of the electric capacity is produced is suitable to enter by the fitting section The row e index fitting data that obtains of fitting is obtained, wherein, the fitting data is the electromotive force decay when electric capacity discharges Coefficient.
Preferably, the transient electromagnetic signal is suitable to by the way that the measurement data is poor with the fitting data, with reality The electric discharge magnetic field change curve that the overall magnetic field change curve of the existing secondary electromagnetic field is produced with the discharge current of the electric capacity Differ from, so as to compensate negative value.
Preferably, e index fitting is carried out to the fitting section, to obtain fitting data, specific method is:
The fitting segment data is taken the logarithm the bottom of by of e, formula Y=lny is obtained;Wherein, Y is logarithmic data, and y is described Fitting segment data;
Fitting a straight line is carried out to the logarithmic data by formula Y=a+bx using least square method, a and b is obtained;Its In, a, b are constant, and x is independent variable;
Fitting data is obtained by formula y=eaebx.
As described above, the method for compensation transient electromagnetic signal negative value of the invention, has the advantages that:Using this hair Bright method, can effectively process the negative value of transient electromagnetic signal, so as to obtain the valid data of long period, effectively improve ground Matter investigation depth.Through the transient electromagnetic signal after overcompensation, it is possible to process means of interpretation pair by traditional ripe TEM data It is explained, and improves the accuracy of resistivity explanation.
Brief description of the drawings
Fig. 1 is shown as the schematic flow sheet of the method for the compensation transient electromagnetic signal negative value of the embodiment of the present invention.
Fig. 2 is shown as the equivalent circuit diagram of the earth in the method for the compensation transient electromagnetic signal negative value of the embodiment of the present invention.
Fig. 3 is shown as the charge and discharge conditional curve of electric capacity in Fig. 2.
Fig. 4 is shown as the electromagnetic field launched in the method for the compensation transient electromagnetic signal negative value of the embodiment of the present invention Current waveform figure.
The method coil central data that Fig. 5 is shown as the compensation transient electromagnetic signal negative value of the embodiment of the present invention arrive two The overall magnetic field change curve of secondary electromagnetic field, wherein, 1. interior illustration is the Linear Amplifer figure of (I) section, and 2. interior illustration is (III) The Linear Amplifer figure of section.
Fig. 6 is shown as comparison diagram after the before processing of the method for the compensation transient electromagnetic signal negative value of the embodiment of the present invention.
Component label instructions
S1~S3 steps
Specific embodiment
Embodiments of the present invention are illustrated below by way of specific instantiation, those skilled in the art can be by this specification Disclosed content understands other advantages of the invention and effect easily.The present invention can also be by specific realities different in addition The mode of applying is embodied or practiced, the various details in this specification can also based on different viewpoints with application, without departing from Various modifications or alterations are carried out under spirit of the invention.
In view of for many disadvantages of transient electromagnetic signal negative value treatment explanation in background of invention, inventor is deep The principle of the regional the earth equivalent capacity discharge and recharge of polarization of having analyzed and researched, has newly to the negative value phenomenon that transient electromagnetic signal occurs Discovery:
The earth equivalent capacity actually would have been completed charging during an electromagnetism VOF, and start electric discharge;By The electromotive force caused in electromagnetism VOF is very big, and leading role is played in the charging to electric capacity, and electric capacity is in the secondary electromagnetic field phase Between discharge, and discharge current accounts for leading role in transient electromagnetic signal late period, and attenuation curve is presented reversion.Also, electric capacity Charge and discharge behavior occurs during an electromagnetic field, but unrelated with an electromagnetism VOF speed.
Based on above-mentioned discovery, inventor has createed the present invention, by being fitted the negative response of attenuation curve late period, from measurement The Negative Acknowledgment is subtracted in data, real response greatly is obtained, so as to realize the compensation to transient electromagnetic signal negative value.
Fig. 1~Fig. 6 is referred to, embodiments of the invention are related to a kind of method for compensating transient electromagnetic signal negative value.Need Illustrate, the diagram provided in the present embodiment only illustrates basic conception of the invention in a schematic way, only show in schema then Show the component relevant with the present invention rather than component count during according to actual implementation, shape and size are drawn, its actual implementation When the kenel of each component, quantity and ratio can be a kind of random change, and its assembly layout kenel be likely to it is increasingly complex.
As shown in figure 1, the method for the compensation transient electromagnetic signal negative value of the present embodiment at least comprises the following steps:
Step S1, the measurement data curve according to transient electromagnetic signal chooses negative value section as fitting section.
In step sl, the starting position as fitting section at the negative peak of negative value section is chosen, position in negative value section is chosen In after negative peak and before noise level at the negative value of random time point as the scram position of fitting section.Wherein, making an uproar Refer to that system does not start acquisition noise before sound level.
Step S2, carries out e index fitting, to obtain fitting data to fitting section.
Wherein, the specific method of step S2 is:
Step S201, takes the logarithm to fitting segment data the bottom of by of e, obtains formula Y=lny;Wherein, Y is logarithmic data, y It is fitting segment data.
Step S202, fitting a straight line is carried out using least square method by formula Y=a+bx to logarithmic data, obtain a and b;Wherein, a, b are constant, and x is independent variable.
Step S203, fitting data is obtained by formula y=eaebx.
Step S3, measurement data is poor with fitting data, to compensate transient electromagnetic signal negative value.
Additionally, in the present embodiment, the specific method of step S1 is:
Step S101, an electromagnetic field is launched by coil to the earth, then receives the secondary electromagnetic field signal of feedback, with Obtain the measurement data curve of transient electromagnetic signal;Wherein, the earth is equivalent to the circuit system of resistance and capacitances in series, so that The overall magnetic field change curve of secondary electromagnetic field has identical trend with the charge and discharge conditional curve of electric capacity.
Step S102, the induced-current for choosing secondary electromagnetic field is less than measurement data corresponding during the discharge current of electric capacity Curved section as fitting section;Wherein, when the induced-current of secondary electromagnetic field is less than the discharge current of electric capacity, corresponding survey It is the negative value section of the measurement data curve of transient electromagnetic signal to measure the curved section of data.
As a preferred scheme, choose at measurement data corresponding during the discharge current maximum of electric capacity as fitting section Starting position;Wherein, when the discharge current of electric capacity is maximum, corresponding measurement data is the negative peak of negative value section.And And, choose from electric capacity there is maximum discharge current after and measurement data before noise level corresponding to random time point at make To be fitted the scram position of section, wherein, referred to electric capacity still in discharge regime before noise level and system does not start collection Noise.The starting position and scram position that are fitted section are set according to the preferred scheme, more accurate transition can be realized Electromagnetic signal negative offset, obtains the valid data of longer time, so as to improve geology detecting depth.
Wherein, in step S101, the earth is equivalent to the circuit system that resistance R connects with electric capacity C, as shown in Fig. 2 so that The overall magnetic field change curve of secondary electromagnetic field has identical trend with the charge and discharge conditional curve of electric capacity.Its principle is as follows:
First, the electric current behavior of electric capacity in circuit analysis, the Main Analysis equivalent circuits is simulated to the equivalent circuit.
As shown in figure 3, a curves are the capacitance current under constant electromotive force, b curves are for electromotive forceElectric capacity electricity Stream.
When electromotive force is constant, its current expression is:
Wherein τc=RC, A are constant coefficient.
As shown in a curves of Fig. 3, the electric current by electric capacity is an e index attenuation curve, is terminated until charging, electric current It is zero.
When electromotive force isWhen, the electric current by electric capacity is an e index curve decayed with the time, and its ammeter reaches Formula is:
Wherein τ is the attenuation coefficient of electromotive force, and A, B are constant coefficient.
As shown in the b curves of Fig. 3, occurs reversal development (the i.e. dotted line of b curves in attenuation process by the electric current of electric capacity Section, the data of the phantom line segments are actual to be negative value, for the ease of expressing in the example shown, the data of the phantom line segments is taken absolute value table Reach, in fact, the phantom line segments curve should be the curve symmetrical relative to X-axis).Wherein, electric current is on the occasion of being expressed as filling for electric capacity Electric process, electromotive force of source ε (0) is more than capacitance voltage at the beginning, now starts to be that electric capacity charges, with subtracting for electromotive force of source Small accumulative with electric capacity electricity, charging current is gradually reduced.Until electromotive force of source is equal with capacitive e.m.f., completion was charged Journey, electric current is zero.With the further reduction of electromotive force of source, electric capacity starts electric discharge, because in the opposite direction, electric current shows as bearing Value (phantom line segments of the b curves of such as Fig. 3), when electromotive force of source decay levels off to zero, the discharge current of electric capacity reaches maximum, with Electric capacity continues to discharge afterwards, and with the reduction of electric capacity electricity, the discharge current of electric capacity is intended to zero.
Secondly, situation is surveyed using the Analysis of Equivalent Circuit.By taking the polarization the earth in Inner Mongolia Autonomous Region mining area as an example, When carrying out transient electromagnetic detecting to it, it is found that the data of all sensing points all have negative value phenomenon, from measurement data analysis The charge and discharge process of electric capacity.
It is as shown in Figure 4 a current waveform figure for electromagnetic field of actual transmission, the electromagnetic field tool of the actual transmission There are the turn-off time for opening electric time and 920ms of 80ms, and the actually trailing edge of 2ms about wide.In figure, (I) section is linear pass Disconnected, (II) section is slow shut-off, and (III) section is to complete switch off.And the variation tendency of the charge and discharge conditional curve of electric capacity is still such as Shown in Fig. 2.The overall magnetic field change curve of the secondary electromagnetic field that hub of a spool is collected is illustrated in figure 5, (I) in figure (II) (I) (II) (III) section in (III) section difference corresponding diagram 4.
Please continue to refer to Fig. 5, (I) section is a behavior for electromagnetic field in 1ms, and 1. its Linear Amplifer figure is shown in interior illustration, The partial magnetic field homogenous linear changes, and produces constant electromotive force, and the part is identical with a curve procedures in Fig. 3, while and b Curve front half section is also identical, is charging process.(II) section is 1ms to 2ms, and this is because the aftercurrent in transmitting coil is produced Magnetic field, its attenuation process is approximate e index, and the part is identical with the b curve procedures in Fig. 3, is discharge process after first charging. (III) after section is 2ms, i.e. overstriking curved section in Fig. 5, the data of the overstriking curved section are actual to be negative value, for the ease of Express in the example shown, the data of the phantom line segments are taken absolute value expression, in fact, the overstriking curved section should be relative to X-axis pair 2. the curve of title, its Linear Amplifer figure is shown in interior illustration;The part is electric capacity discharge process, the induced-current side with secondary electromagnetic field To conversely, when the induced-current of secondary electromagnetic field is more than the discharge current of electric capacity, magnetic field value B (i.e. magnetic induction intensity) shows as On the occasion of, the induced-current with secondary electromagnetic field gradually weakens, and less than the discharge current of electric capacity, now magnetic field value B appearance Negative value, next as the induced-current of secondary electromagnetic field decays to very little level so that this section of magnetic field of secondary electromagnetic field becomes Change curve to show and electric capacity discharge process curve identical trend.
By above-mentioned analysis and Fig. 3 and Fig. 5, the earth is equivalent to the circuit system that resistance R connects with electric capacity C, can To obtain the overall magnetic field change curve of the secondary electromagnetic field that there is same trend with the charge and discharge conditional curve of electric capacity.
In addition, the preferred scheme in step S102, chooses measurement number corresponding during the discharge current maximum of electric capacity According to place as the starting position for being fitted section;Wherein, when the discharge current of electric capacity is maximum, corresponding measurement data is negative value section Negative peak.Meanwhile, choose from electric capacity there is maximum discharge current after and before noise level corresponding to random time point Measurement data at as fitting section scram position.Continue by taking the polarization the earth in Inner Mongolia Autonomous Region mining area as an example, by scheming 3rd, Fig. 4 and Fig. 5 are visible, and the discharge current of electric capacity measurement data corresponding when maximum (negative peak of negative value section) place is Magnetic field value B during 5ms, therefore, the magnetic field value B when starting position for being fitted section is 5ms;Meanwhile, from after 5ms to noise level Still there is the measurement data of more than 100ms before, in order to ensure that compensation negative value more prepares, choose in the 20ms~100ms after 5ms As the scram position of fitting section at the magnetic field value B of random time point.
In addition, in the present embodiment, the electric discharge magnetic field change curve that the discharge current of electric capacity is produced is suitable to by fitting The Duan Jinhang e indexes fitting data that obtains of fitting is obtained, wherein, fitting data is that electromotive force decay when electric capacity discharge is Number.Also, transient electromagnetic signal is suitable to by the way that measurement data is poor with fitting data, to realize the overall magnetic of secondary electromagnetic field Field change curve is poor with the electric discharge magnetic field change curve that the discharge current of electric capacity is produced, so as to compensate negative value.Need explanation It is that, please continue to refer to Fig. 5, after (III) section 5ms, the discharge current of electric capacity is occupied an leading position, in can be to this period Curve carries out e index and fits τ values when electric capacity discharges, so as to be inferred to entirely (III) section according to above-mentioned formula (1) and (2) During electric capacity discharge current produce electric discharge magnetic field.The electric discharge magnetic field of the fitting is so subtracted by total magnetic field, so that it may With the magnetic responsiveness being inferred under pure resistor element circuit system greatly, then can just be processed by traditional ripe TEM data Means of interpretation is explained to it.
Still by taking the polarization the earth in Inner Mongolia Autonomous Region mining area as an example, using the compensation transient electromagnetic signal of the present embodiment The treatment effect that the method for negative value is obtained is as shown in fig. 6, in figure, solid line is represented the measurement of actually measured transient electromagnetic signal Data and curves blocked using the method that tradition blocks negative value after data and curves, dotted line is represented believes actually measured transient electromagnetic Number measurement data curve by the data and curves after the method compensation correction of the present embodiment.It can be seen that, block negative using traditional The method of value, the data that block curve only has 3ms can be utilized, and correspondence subterranean depth is only about 400 meters;And use the present embodiment Method compensation correction after, the utilizable data of dashed curve reach 10ms or so, and correspondence subterranean depth is about 800 meters.Cause This, the method for the compensation transient electromagnetic signal negative value of the present embodiment can effectively improve investigation depth.
Above the step of various methods divide, be intended merely to description it is clear, can be merged into when realizing a step or Some steps are split, multiple steps are decomposed into, as long as comprising identical logical relation, all in the protection domain of this patent It is interior;To adding inessential modification in algorithm or in flow or introducing inessential design, but its algorithm is not changed With the core design of flow all in the protection domain of the patent.
In sum, the method for compensation transient electromagnetic signal negative value of the invention, has the advantages that:Using this hair Bright method, can effectively process the negative value of transient electromagnetic signal, so as to obtain the valid data of long period, effectively improve ground Matter investigation depth.Through the transient electromagnetic signal after overcompensation, it is possible to process means of interpretation pair by traditional ripe TEM data It is explained, and improves the accuracy of resistivity explanation.So, the present invention effectively overcomes various shortcoming of the prior art And have high industrial utilization.
The above-described embodiments merely illustrate the principles and effects of the present invention, not for the limitation present invention.It is any ripe The personage for knowing this technology all can carry out modifications and changes under without prejudice to spirit and scope of the invention to above-described embodiment.Cause This, those of ordinary skill in the art is complete with institute under technological thought without departing from disclosed spirit such as Into all equivalent modifications or change, should be covered by claim of the invention.

Claims (9)

1. it is a kind of compensate transient electromagnetic signal negative value method, it is characterised in that the side of the compensation transient electromagnetic signal negative value Method at least comprises the following steps:
Measurement data curve according to transient electromagnetic signal, chooses negative value section as fitting section;
E index fitting is carried out to the fitting section, to obtain fitting data;
The measurement data is poor with the fitting data, to compensate the transient electromagnetic signal negative value.
2. the method for compensation transient electromagnetic signal negative value according to claim 1, it is characterised in that choose negative value section Negative peak at as it is described fitting section starting position.
3. the method for compensation transient electromagnetic signal negative value according to claim 2, it is characterised in that choose negative value section As the shut-off position of the fitting section at the negative value of the interior random time point after the negative peak and before noise level Put.
4. it is according to claim 1 compensation transient electromagnetic signal negative value method, it is characterised in that according to transient electromagnetic believe Number measurement data curve, choose negative value section as fitting section, specific method is:
Electromagnetic field is launched to the earth by coil, the secondary electromagnetic field signal of feedback is then received, to obtain the transition The measurement data curve of electromagnetic signal;Wherein, the earth is equivalent to the circuit system of resistance and capacitances in series, so that described The overall magnetic field change curve of secondary electromagnetic field has identical trend with the charge and discharge conditional curve of the electric capacity;
The induced-current for choosing the secondary electromagnetic field is less than the song of measurement data corresponding during the discharge current of the electric capacity Line segment is used as fitting section;Wherein, it is corresponding when the induced-current of the secondary electromagnetic field is less than the discharge current of the electric capacity Measurement data curved section be the transient electromagnetic signal measurement data curve negative value section.
5. the method for compensation transient electromagnetic signal negative value according to claim 4, it is characterised in that choose the electric capacity As the starting position of the fitting section at discharge current measurement data corresponding when maximum;Wherein, putting in the electric capacity When electric current is maximum, corresponding measurement data is the negative peak of the negative value section.
6. the method for compensation transient electromagnetic signal negative value according to claim 5, it is characterised in that choose from the electric capacity As the fitting at measurement data after appearance maximum discharge current and corresponding to the random time point before noise level The scram position of section.
7. it is according to claim 4 compensation transient electromagnetic signal negative value method, it is characterised in that the electric discharge of the electric capacity The electric discharge magnetic field change curve that electric current is produced be suitable to by the fitting section carried out the e index fitting data that obtains of fitting come Arrive, wherein, the fitting data is the electromotive force attenuation coefficient when electric capacity discharges.
8. the method for compensation transient electromagnetic signal negative value according to claim 7, it is characterised in that transient electromagnetic letter Number it is suitable to, by the way that the measurement data is poor with the fitting data, to change with the overall magnetic field for realizing the secondary electromagnetic field Curve is poor with the electric discharge magnetic field change curve that the discharge current of the electric capacity is produced, so as to compensate negative value.
9. it is according to claim 1 compensation transient electromagnetic signal negative value method, it is characterised in that to it is described fitting section enter Row e index is fitted, and to obtain fitting data, specific method is:
The fitting segment data is taken the logarithm the bottom of by of e, formula Y=lny is obtained;Wherein, Y is logarithmic data, and y is the fitting Segment data;
Fitting a straight line is carried out to the logarithmic data by formula Y=a+bx using least square method, a and b is obtained;Wherein, a, b Constant is, x is independent variable;
By formula y=ea·ebxObtain fitting data.
CN201710031598.4A 2017-01-17 2017-01-17 A method of compensation transient electromagnetic signal negative value Active CN106814403B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710031598.4A CN106814403B (en) 2017-01-17 2017-01-17 A method of compensation transient electromagnetic signal negative value

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710031598.4A CN106814403B (en) 2017-01-17 2017-01-17 A method of compensation transient electromagnetic signal negative value

Publications (2)

Publication Number Publication Date
CN106814403A true CN106814403A (en) 2017-06-09
CN106814403B CN106814403B (en) 2019-01-04

Family

ID=59111794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710031598.4A Active CN106814403B (en) 2017-01-17 2017-01-17 A method of compensation transient electromagnetic signal negative value

Country Status (1)

Country Link
CN (1) CN106814403B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285790A (en) * 2020-03-30 2021-01-29 中国科学院地质与地球物理研究所 Method and device for determining attenuation coefficient of electromagnetic wave field and medium
CN116449442A (en) * 2023-06-16 2023-07-18 安徽惠洲地质安全研究院股份有限公司 Three-dimensional unfolding processing method and system for drilling transient electromagnetic method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101882964A (en) * 2010-06-12 2010-11-10 桂林电子科技大学 De-noising method of transient electromagnetic detecting echo signal
US9310511B2 (en) * 2012-11-01 2016-04-12 Baker Hughes Incorporated Apparatus and method for deep transient resistivity measurement
CN105549097A (en) * 2015-12-22 2016-05-04 吉林大学 Transient electromagnetic signal power frequency and harmonic interference elimination method and apparatus thereof
CN105652325A (en) * 2016-01-05 2016-06-08 吉林大学 Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101882964A (en) * 2010-06-12 2010-11-10 桂林电子科技大学 De-noising method of transient electromagnetic detecting echo signal
US9310511B2 (en) * 2012-11-01 2016-04-12 Baker Hughes Incorporated Apparatus and method for deep transient resistivity measurement
CN105549097A (en) * 2015-12-22 2016-05-04 吉林大学 Transient electromagnetic signal power frequency and harmonic interference elimination method and apparatus thereof
CN105652325A (en) * 2016-01-05 2016-06-08 吉林大学 Exponential fit-adaptive Kalman-based ground-air electromagnetic data de-noising method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周磊 等: "激发极化介质的瞬变电磁一维正演", 《物探与化探》 *
王伟明: "考虑激电效应的瞬变电磁中心回线一维正反演", 《新疆有色金属》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285790A (en) * 2020-03-30 2021-01-29 中国科学院地质与地球物理研究所 Method and device for determining attenuation coefficient of electromagnetic wave field and medium
CN116449442A (en) * 2023-06-16 2023-07-18 安徽惠洲地质安全研究院股份有限公司 Three-dimensional unfolding processing method and system for drilling transient electromagnetic method
CN116449442B (en) * 2023-06-16 2023-08-29 安徽惠洲地质安全研究院股份有限公司 Three-dimensional unfolding processing method and system for drilling transient electromagnetic method

Also Published As

Publication number Publication date
CN106814403B (en) 2019-01-04

Similar Documents

Publication Publication Date Title
CN101650443B (en) Back-propagation network calculating method of apparent resistivity
CN112415615B (en) Time-frequency electromagnetic fracturing monitoring system and monitoring method based on distributed optical fiber sensing
CN103995301A (en) Method and device for evaluating total organic carbon content in shale gas reservoir
CN106199732A (en) A kind of transient electromagnetic multiple tracks covers observation device and method
CN107044866A (en) The soil moisture, moisture and conductivity measuring apparatus
CN106814403B (en) A method of compensation transient electromagnetic signal negative value
CN110488365A (en) A kind of multipolarization telluric electromagnetic sounding method
YAN et al. The probing depth of transient electromagnetic field method
Flores Orozco et al. Investigation of cable effects in spectral induced polarization imaging at the field scale using multicore and coaxial cables
CN105549100B (en) Transient Electromagnetic Method in Mine apparatus and method based on U-shaped helical source
JI et al. A study on solution of transient electromagnetic response during transmitting current turn‐off in the ATTEM system
Zhou et al. Induced polarization effect on grounded-wire transient electromagnetic data from transverse electric and magnetic fields
CN106869913B (en) A method of water injection well in oil fields waterflood front is detected using well-in-situ potential technology
CN108169802A (en) A kind of time domain electromagnetic data slow diffusion imaging method of harsh media model
CN108776357B (en) The bearing calibration of sedimentary formation transient electromagnetic method electromagnetic interference and device
CN114217354A (en) Electromagnetic data acquisition system and method based on optical fiber electromagnetic sensor
CN206573161U (en) A kind of soil moisture, moisture and conductivity measuring apparatus
CN105137495B (en) A kind of method and system of oil and gas detection
CN105700026B (en) A kind of buried target body conduction index extracting method and device
Zhang et al. Analysis of underground propagation effects of lightning electromagnetic fields in different geological environments
Li et al. Airborne transient electromagnetic simulation: detecting geoelectric structures for HVdc monopole operation
CN106248740B (en) A kind of soil resistivity humidity correcting method
CN104991281B (en) A kind of detection method and device of coal seam buried depth
CN104502986B (en) Physical prospecting induced electrical sounding data chromatography processing method
WANG et al. Research on the effect of 3D body between transmitter and receivers on CSAMT response using Integral Equation method

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
GR01 Patent grant
GR01 Patent grant