CN107092033A - It is a kind of to be used to regulate and control the method for measurement signal intensity - Google Patents
It is a kind of to be used to regulate and control the method for measurement signal intensity Download PDFInfo
- Publication number
- CN107092033A CN107092033A CN201710258059.4A CN201710258059A CN107092033A CN 107092033 A CN107092033 A CN 107092033A CN 201710258059 A CN201710258059 A CN 201710258059A CN 107092033 A CN107092033 A CN 107092033A
- Authority
- CN
- China
- Prior art keywords
- field
- component
- axles
- transmitting
- intensity
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/084—Sources
Abstract
The present invention provides a kind of method for regulating and controlling measurement signal intensity, the regularity of distribution of each electromagnetic field component after the different electric currents of analysis two, different angles, the transmitting field source synthesis of different length, by the difference for adjusting the starting phase angle between two field source electric currents, obtain initial phase angular difference between suitable two field source electric currents, strengthen region electromagnetic field intensity to be measured, or expand the virtual value of regional signal to be measured, including:Using two transmitting field sources, current signal is sent to position to be measured respectively;The position to be measured is fed back to the current signal of described two transmitting field sources and forms two different electromagnetic fields respectively, based on the field strength component of described two different electromagnetic fields, to build the first relation function;Based on first relation function, the optimal value of the superposition field strength formed after described two different electromagnetic field superpositions is obtained.Method of the invention, it is possible to effectively improve the laying flexibility ratio and efficiency of signal intensity and transmitting field source.
Description
Technical field
The present invention relates to geophysical exploration technology, it is used to regulate and control the side of measurement signal intensity more particularly, to a kind of
Method.
Background technology
Artificial source's electromagnetic method is developed into from natural source electromagnetic method, the introducing of artificial source enhances signal quality, improves letter
Make an uproar and compare, reduce the randomness of signal.Electric dipole and magnetic dipole are conventional transmitting field source.
In terms of electric dipole, generally from a horizontal electric dipole, using easy, economic scalar measurement.Due to water
Flat Electric Dipole length generally reaches 1-2km, arranges that orthogonal horizontal electric dipole difficulty is larger.It is mostly single galvanic couple for emission source
Pole component, the radiation pattern chart of theoretical calculation is bright, and each component of electromagnetic field all has necessarily weak in 360 degree of subtended angle scopes
Area.In order to ensure the quality of data, measured zone has strict limitation, general to require launching each 45 ° of dipole perpendicular bisector both sides
Measure, received and dispatched away from generally 10~15km with the subtended angle scope of each 30 ° of axial both sides.Rugged mountain area, is sometimes difficult to meet this
A little harsh conditions, so as to can not carry out normally measuring work, even if work reluctantly, it is also difficult to ensure the quality of data.
In terms of magnetic dipole, the form of magnetic dipole determines that it can not arrange two orthogonal field sources, it is difficult to realize
Tensor measuring.The signal attenuation of horizontal magnetic source transmitting is very fast, and less transmitting-receiving is away from can meet far zone condition, orthogonal horizontal magnetic
Even source arrangement is convenient, easily carries out quick shallow-layer and detects, but the electromagnetic field measurementses produced by orthogonal horizontal magnetic couple source are still
There is weak output signal region, it is impossible to meet whole district's measurement.
Either still there is signal attenuation in the measurement process of magnetic field intensity, soon in electric-field intensity, it is to be measured
The feedback signal strength that position is produced is weak, the problems such as measurement accuracy is restricted.Also, current generating apparatus or position to be measured
The laying of reception device needs to meet specific requirement, and adds application task difficulty.
The content of the invention
The present invention provides a kind of regulation and control that are used for for overcoming above mentioned problem or solving the above problems at least in part and measures letter
The method of number intensity, to solve, measurement signal decay is fast, can not meet the technical problem of whole district's measurement.
It is used to regulate and control the method for measurement signal intensity there is provided a kind of according to an aspect of the present invention, analyzes two differences
The regularity of distribution of each electromagnetic field component after electric current, different angles, the field source synthesis of different length, by adjusting two field source electricity
The difference of starting phase angle between stream, obtains initial phase angular difference between suitable two field source electric currents, strengthens region to be measured
Electromagnetic field intensity, or expand the virtual value of regional signal to be measured, including:
Step S1, using two transmitting field sources, respectively to position to be measured send current signal;
Step S2, the position to be measured are fed back to the current signal of described two transmitting field sources and form two differences respectively
Electromagnetic field, based on the field strength component of described two different electromagnetic fields, to build the first relation function;
Step S3, based on first relation function, obtain formed after described two different electromagnetic fields superpositions folded
The strong optimal value of extra show.
Further, the superposition field strength is electric-field intensity or magnetic field intensity.
Further, the angular range described in rapid S1 between two transmitting field sources is 0 ° -360 °.
Further, the field strength component of two different electromagnetic fields described in step S2, is specifically included:Described two transmittings
Field source includes the first transmitting field source and the second transmitting field source, and the intersection point using the described first transmitting field source and the second transmitting field source is seat
Mark origin,
The first rectangular coordinate system is set up by x1 axles of the described first transmitting field source, and by the feedback of the described first transmitting field source
Signal calculates in the radial field component and tangential field component of the position to be measured and obtains the first transmitting field source described
Position to be measured is along x1 axles and/or the first field strength component of y1 axles;
The second rectangular coordinate system is set up by x2 axles of the described second transmitting field source, and by the feedback of the described second transmitting field source
Signal calculates in the radial field component and tangential field component of the position to be measured and obtains the second transmitting field source described
Position to be measured is along x2 axles and/or the second field strength component of y2 axles.
Further, the method that the second field strength component is rotated by reference axis is converted into sits at first right angle
Along x1 axles and/or the 3rd field strength component of y1 axles in mark system.
Further, the particular content of the reference axis rotation is:By the x2 axles/y2 axles corresponding rotation to the x1
Coordinate in axle/y1 overlapping of axles, second rectangular coordinate system and the coordinate are corresponding in the first rectangular coordinate system new to be sat
Corresponding relation is set up between mark.
Further, the particular content of first relation function is built in step S2 to be included:
Step S21, using the described first transmitting field source as x-axis, set up overlapped with first rectangular coordinate system the 3rd straight
Angular coordinate system;
Step S22, by the first field strength component and the 3rd field strength component, obtain it is described first transmitting field source and second hair
Penetrate field source and the superimposed field strong component to be formed is superimposed in the position to be measured, the superimposed field strong component is in the 3rd rectangular coordinate system
The middle field strength component along x-axis and/or y-axis, to form the first relation function.
Further, first relation function is based in step S3, obtains and is superimposed by described two different electromagnetic fields
The optimal value of the superposition field strength formed afterwards, is specifically included:
Step S31, set up it is described first transmitting field source the first electric current with described second transmitting field source the second electric current it
Between the second relation function on phase difference;
Step S32, based on first relation function and the second relation function, take different phase difference values, calculating is obtained
Different superposition field intensity value, until obtaining the optimal value of the superposition field strength.
Further, the scope of phase difference described in step S31 is 0 ° -360 °.
Further, first relation function is expressed as:
Ex=Ex1+E′x1=Ex1+Ex2cosθ-Ey2sinθ
Ey=Ey1+E′y1=Ey1+Ex2sinθ+Ey2cosθ
Hx=Hx1+H′x1=Hx1+Hx2cosθ-Hy2sinθ
Hy=Hy1+H′y1=Hx1+Hx2sinθ+Hy2cosθ
In formula, Ex、EyRespectively component of the superposition of electric field intensity along x-axis/y-axis, Ex1、Ey1Respectively the first electric-field intensity edge
The component of x1 axles/y1 axles, E 'x1、E′y1Component of respectively the 3rd electric-field intensity along x1 axles/y1 axles, Ex2、Ey2Respectively second electricity
Component of the field intensity along x2 axles/y2 axles;Hx、HyRespectively it is superimposed component of the magnetic field intensity along x-axis/y-axis, Hx1、Hy1Respectively
Component of one magnetic field intensity along x1 axles/y1 axles, H 'x1、H′y1Component of respectively the 3rd magnetic field intensity along x1 axles/y1 axles, Hx2、Hy2
Component of respectively the second electromagnetic field intensity along x2 axles/y2 axles;θ is the angle between the first transmitting field source and the second transmitting field source.
A kind of method for regulating and controlling measurement signal intensity proposed by the present invention, its advantage is mainly as follows:
(1) transmitting field source uses two, and based on the first constructed relation function, to obtain the optimal of superposition field strength
Value, strengthens the feedback signal strength of position to be measured, improves measurement accuracy and reliability;
Angle between (2) two transmitting field sources is arbitrary, reduces the layout difficulty of transmitting field source;
(3) produced by the phase difference value between the current signal by introducing two transmitting field sources, two transmitting field sources of elimination
Weak output signal region, improve Measurement reliability, also, be adjusted merely by the phase difference value, just can adjust the strong of feedback signal
Degree, control methods are simple, effective.
Brief description of the drawings
Fig. 1 is a kind of schematic flow sheet of method for being used to regulate and control measurement signal intensity according to the embodiment of the present invention;
Fig. 2 is according to showing that a kind of reference axis of method for being used to regulate and control measurement signal intensity of the embodiment of the present invention rotates
It is intended to;
Fig. 3 is a kind of the 3rd rectangular coordinate system of method for being used to regulate and control measurement signal intensity according to the embodiment of the present invention
Schematic flow sheet;
Fig. 4 a do not carry out phase difference for a kind of method for being used to regulate and control measurement signal intensity according to the embodiment of the present invention
Electric-field intensity component, the radiation pattern figure of magnetic field strength component during value regulation;
Fig. 4 b are a kind of carry out phase difference value of method for being used to regulate and control measurement signal intensity according to the embodiment of the present invention
The radiation pattern figure of electric-field intensity component, magnetic field strength component after regulation.
Embodiment
With reference to the accompanying drawings and examples, the embodiment to the present invention is described in further detail.Implement below
Example is used to illustrate the present invention, but is not limited to the scope of the present invention.
It is shown in Figure 1, it is a kind of to be used to regulate and control the method for measurement signal intensity, analyze two different electric currents, different angles,
The regularity of distribution of each electromagnetic field component after the field source synthesis of different length, by adjusting the initial phase between two field source electric currents
The difference of parallactic angle, obtains initial phase angular difference between suitable two field source electric currents, strengthens region electromagnetic field intensity to be measured, or expand
The virtual value of big regional signal to be measured.It includes:
Step S1, using two transmitting field sources, respectively to position to be measured send current signal;
Step S2, the position to be measured are fed back to the current signal of described two transmitting field sources and form two differences respectively
Electromagnetic field, based on the field strength component of described two different electromagnetic fields, to build the first relation function;
Step S3, based on first relation function, obtain formed after described two different electromagnetic fields superpositions folded
The strong optimal value of extra show.
In the wild in actual measurement process, due to the influence of geographical position and geological environment, the electric and magnetic responses of some regions
Signal is weaker, using two transmitting field sources, and enhancing treats location and puts the current signal sent, to reach enhancing position response to be measured
The purpose of signal, the precision of analysis that location puts the mt regularity of distribution is treated so as to improve.
The current signal that two transmitting field sources are sent, the different electromagnetic fields that can be formed respectively in position to be measured, two are not
Same electromagnetic field is superimposed in position to be measured and forms superposition field strength.Based on two different field strength, structure obtains first
Relation function.First relation function reflects that the field strength that two transmitting field sources are individually formed is formed after being superimposed with two transmitting field sources
Superimposed field Qianghian relation.
Based on first relation function to obtain the superposition field strength of the position to be measured, by adjusting involved by the first relation function
And parameter or associated arguments are further associated by first relation function, with adjust superposition after superposition field strength intensity.
By the optimization for the superposition performance that two energy are launched with field strength, to regulate and control the signal intensity ultimately formed in position to be measured, with
Make up the problem of mono signal source signal intensity is difficult to be substantially improved.
In a specific embodiment, two transmitting field sources send current signal, position difference to be measured to position to be measured
Feedback signal is produced to two current signals, two different electromagnetic fields are formed by two produced feedback signals.According to
For the difference for the equipment for receiving the field intensity signal, the field strength can be electric-field intensity or magnetic field intensity.
Specifically, when the equipment for received field strength signal is potential difference measurement instrument, received feedback signal
The two different field strength formed are electric-field intensity, correspondingly, and what is formed after this two different field intensity signal superpositions is folded
Extra show is electric-field intensity by force;When the equipment for received field strength signal is Magnetic Sensor, received feedback signal institute shape
Into two different field strength be magnetic field intensity, correspondingly, the superimposed field formed after this two different field intensity signals superpositions
It is magnetic field intensity by force.
For same group of transmitting field source, using different signal receivers, you can for the earth target area to be measured
The research of magnetic distribution rule.It is that electric-field intensity and magnetic field are strong to the feedback signal received below with different embodiments
Degree is illustrated respectively.
Embodiment 1
In another specific embodiment, two transmitting field sources send current signal, two hairs to position to be measured respectively
The angle penetrated between field source is 0 ° -360 °.It is, in actual measurement process, the positions of two transmitting field sources can random cloth
If each other without special, relative status requirement, greatly reducing the layout difficulty of transmitting field source, improving work
Industry efficiency.
In another specific embodiment, two transmitting field sources include the first transmitting field source and the second transmitting field source, with
The intersection point of the intersection point or its extended line of first transmitting field source and the second transmitting field source is the origin of coordinates, and rectangular co-ordinate is set up respectively
System.Wherein, the line of position to be measured and the origin of coordinates is designated as the first line.
Specifically, the first rectangular coordinate system is set up using the first transmitting field source as x1 axles.Position to be measured is to the first transmitting field source
The first current signal produce feedback signal, the feedback signal position to be measured formed the first electric field intensity signal, by this first
Electric field intensity signal calculates in the radial electric field strength component and tangential electric-field intensity component of position to be measured and obtains first electric field
Strength signal is along x1 axles and/or the first electric-field intensity component of y1 axles.
On first electric field intensity signal position to be measured radial electric field strength component Er1With tangential electric-field intensity point
AmountCalculating be expressed as follows:
In formula, i1For the first electric current of the first transmitting field source, ρ is homogeneous half space resistivity, and r is position to be measured to coordinate
The distance of origin, k is wave number, and a is the length of the first transmitting field source Electric Dipole, and φ is between the first line and the first transmitting field source
Angle.
By the first electric field intensity signal position to be measured radial electric field strength component Er1With tangential electric-field intensity componentObtaining specific formula for calculation of first electric field intensity signal along x1 axles and/or the first electric-field intensity component of y1 axles is:
In formula, Ex1Component for the first electric-field intensity along x1 axles, Ey1Component for the first electric-field intensity along y1 axles, Er1、Respectively the first electric-field intensity is in the radial electric field strength component and tangential electric-field intensity component of position to be measured, and φ is first
Angle between line and the first transmitting field source.
Specifically, the second rectangular coordinate system is set up using the second transmitting field source as x2 axles.Position to be measured is to the second transmitting field source
The second current signal produce feedback signal, the feedback signal position to be measured formed the second electric field strength signal, by this second
Electric field strength signal calculates in the radial electric field strength component and tangential electric-field intensity component of position to be measured and obtains second electric field
Strength signal is along x2 axles and/or the second electric-field intensity component of y2 axles.
On second electric field intensity signal position to be measured radial electric field strength component Er2With tangential electric-field intensity point
AmountCalculating be expressed as follows:
In formula, i2For the second electric current of the second transmitting field source, ρ is homogeneous half space resistivity, and r is position to be measured to coordinate
The distance of origin, k is wave number, and b is the length of the second transmitting field source Electric Dipole, and θ is the second transmitting field source and the first transmitting field source
Between angle, φ be the first line and first transmitting field source between angle.
By the second electric field intensity signal position to be measured radial electric field strength component Er2With tangential electric-field intensity componentObtaining specific formula for calculation of second electric field intensity signal along x2 axles and/or the second electric-field intensity component of y2 axles is:
In formula, Ex2Component for the second electric-field intensity along x2 axles, Ey2Component for the second electric-field intensity along y2 axles, Er2、Respectively the second electric-field intensity is in the radial electric field strength component and tangential electric-field intensity component of position to be measured, and θ is the second hair
The angle penetrated between field source and the first transmitting field source, φ is the angle between the first line and the first transmitting field source.It is shown in Figure 2,
In another specific embodiment, for ease of calculate the first electric-field intensity be superimposed with the second electric-field intensity formed by superposition it is electric
Field intensity, by the second rectangular coordinate system by way of rotation, makes the x2 axles and the first rectangular co-ordinate of former second rectangular coordinate system
The x1 overlapping of axles of system, the y2 axles of former second rectangular coordinate system and the y1 overlapping of axles of the first rectangular coordinate system.Second electric-field intensity point
Amount is converted into the 3rd electric-field intensity component after the rotation of the second rectangular coordinate system.
Specifically, the second rectangular coordinate system is rotated to after being overlapped with the first rectangular coordinate system, in former second rectangular coordinate system
Position relationship do not change, it is, being merely the rotation of rectangular coordinate system.
In another specific embodiment, the second electric-field intensity component is converted into after the 3rd electric-field intensity component, the
Former coordinate of the two electric-field intensity components in the second rectangular coordinate system is converted into the new coordinate in the first rectangular coordinate system, and
Exist between the new coordinate and former coordinate and exist between corresponding relation, i.e. the 3rd electric-field intensity component and the second electric-field intensity component
Corresponding relation.
Specifically, former coordinate representation of the second electric-field intensity component in the second rectangular coordinate system is along x2 axles/y2 axles
Component, respectively Ex2、Ey2, the second electric-field intensity is converted into after the 3rd electric-field intensity, and the 3rd electric-field intensity component is straight first
Coordinate representation in angular coordinate system is the component E' along x1 axles/y1 axlesx1、E'y1。
Corresponding relation statement between 3rd electric-field intensity component and the second electric-field intensity component is as follows:
In formula, E 'x1、E′y1Component for the 3rd electric-field intensity component in the first rectangular coordinate system along x1 axles/y1 axles,
Ex2、Ey2Component for the second electric-field intensity component in the second rectangular coordinate system along x2 axles/y2 axles, θ be second transmitting field source with
Angle between first transmitting field source.
It is shown in Figure 3, in another specific embodiment, with the first electric-field intensity component and the second electricity in step S2
Based on field intensity component, to build the first relation function.The step of building first relation function be:
Step S21, using the described first transmitting field source as x-axis, set up overlapped with first rectangular coordinate system the 3rd straight
Angular coordinate system;
Step S22, by the first electric-field intensity component and the 3rd electric-field intensity component, obtain it is described first transmitting field source
The superposition of electric field strength component to be formed is superimposed in the position to be measured with the second transmitting field source, the superposition of electric field strength component is
Along the component of x-axis and/or y-axis in the 3rd rectangular coordinate system, to form the first relation function.
The 3rd rectangular coordinate system set up is overlapped with the first rectangular coordinate system, for statement superposition of electric field intensity.
First electric-field intensity and the second electric-field intensity are superimposed at position to be measured, the superposition of electric field formed after superposition, to increase
Signal intensity of the heavy current signal at position to be measured.Therefore, based on first electric-field intensity and the second electric-field intensity, it can obtain
To the superposition of electric field intensity on position to be measured.Correspondingly, can by the first electric-field intensity component and the second electric-field intensity component
The component of the superposition of electric field intensity on position to be measured is obtained, that is, forms the first pass for calculating the superposition of electric field strength component
It is function.
In another specific embodiment, have on the first relation function for calculating the superposition of electric field strength component
Body is expressed as follows:
In formula, Ex、EyRespectively component of the superposition of electric field intensity along x-axis/y-axis, Ex1、Ey1Respectively the first electric-field intensity edge
The component of x1 axles/y1 axles, E 'x1、E′y1Component of respectively the 3rd electric-field intensity along x1 axles/y1 axles, Ex2、Ey2Respectively second electricity
Component of the field intensity along x2 axles/y2 axles, θ is the angle between the first transmitting field source and the second transmitting field source.
In another specific embodiment, in step S3, the first pass obtained for stating superposition of electric field intensity is built
It is that after function, based on first relation function, the optimal value of superposition of electric field intensity can be acquired, its specific steps includes:
Step S31, set up it is described first transmitting field source the first electric current with described second transmitting field source the second electric current it
Between the second relation function on phase difference;
Step S32, based on first relation function and the second relation function, take different phase difference values, calculating is obtained
Different superposition of electric field strength component, until obtaining the optimal value of the superposition of electric field strength component.
In actual measurement process, when sending current signal to position to be measured, one to adjust feedback signal strength
Importance, is the attribute by adjusting the current signal that transmitting field source is provided, and the attribute includes the size and Orientation of electric current.
After the generation equipment for launching field source is determined, the regulation of size of current is more inconvenient, and adjusts current generation device conveying electric current
Phase method, then be easily achieved.
During using two transmitting field sources, by adjusting the phase difference value between its corresponding first electric current and the second electric current,
To adjust the superposition of electric field intensity after the first electric-field intensity and the second electric-field intensity are superimposed in position to be measured to optimal value, so as to reach
To the purpose for eliminating small-signal area, make the superposition of electric field strength signal intensity finally given high.
Only by adjusting the phase difference between the first electric current and the second electric current, the strong of position feed back signal to be measured can be optimized
Degree, substantially increases the laying flexibility of transmitting field source, reduces transmitting field source addressing and layout difficulty, can be according to practically
Matter environment, avoids unfavorable terrain, improves the flexibility of field work efficiency and construction, meanwhile, also improve the accuracy of measurement
And reliability.
In another specific embodiment, the second relation function between the first electric current and the second electric current on phase difference
Specifically it is expressed as follows:
In formula, i2For the second electric current of the second transmitting field source, i1For the first electric current of the first transmitting field source, λ is 0
Real number,For the first electric current and the phase difference value of the second electric current, i is imaginary unit.
It is understood that also formula (7) can be expressed as into i1On i2Relation function, i.e.,In formula, λ '
Or not 0 real number,For the first electric current and the phase difference value of the second electric current, i is imaginary unit.
Specifically, different phase differences are taken respectivelyValue, based on formula (6) and (7), can calculate and obtain different Ex、Ey
Value.Each phase differenceValue i.e. one group of E of correspondencex、EyValue.In resulting array Ex、EyIn value, optimal or numerical value is chosen maximum
One group of Ex、EyValue, it is now correspondingValue is optimum angle difference.
In another specific embodiment, the span of phase difference value is 0 ° -360 °.Due to the change of current phase
When showing numerically, there is certain repeatability, therefore, in actual mechanical process, the span of phase difference value is limited
At 0 ° -90 °.Take different phase differences successively in the range of 0 ° -90 °Value, is calculated by formula (6) and (7) and obtained
Corresponding ExValue, EyValue, until obtaining optimal or maximum superposition of electric field intensity level.
It is understood that in actual mechanical process, when it is determined that the first transmitting field source, the second transmitting field source and treating location
Postpone, when laying two transmitting field sources, the size of its first electric current, first launch field source Electric Dipole length, in midair
Between resistivity, the distance of position to be measured to the first rectangular coordinate system origin be known quantity;At the same time, the size of the second electric current,
The length of the Electric Dipole of second transmitting field source, half space resistivity, position to be measured are to the distance of the second rectangular coordinate system origin
Known quantity;In addition, the angle between the angle between two transmitting field sources, and position to be measured and the first transmitting field source is also known
Amount, i.e. only input different phase differencesValue, just can calculate and obtain corresponding ExValue, EyValue, so that it is poor to obtain optimum angle
Value.Obtain the optimum angle poorAfter value, the actual phase difference of the first transmitting field source and the second transmitting field source is adjusted to this most
Excellent phase differenceValue, i.e., can obtain the good feedback signal of signal condition in position to be measured in actual measurement process, so that
Improve measurement accuracy and reliability.
Further, it is to be appreciated that by formula (6) and (7), and take different phase differencesValue, obtains array pair
The E answeredxValue, EyValue, is carrying out the optimal E of selectionxValue, EyWhile value, the E with special characteristic can be also selectedxValue, Ey
Value, to be analyzed accordingly.
Embodiment 2
In another specific embodiment, two transmitting field sources send current signal, two hairs to position to be measured respectively
The angle penetrated between field source is 0 ° -360 °.It is, in actual measurement process, the positions of two transmitting field sources can random cloth
If each other without special, relative status requirement, greatly reducing the layout difficulty of transmitting field source, improving work
Industry efficiency.
In another specific embodiment, two transmitting field sources include the first transmitting field source and the second transmitting field source, with
The intersection point of the intersection point or its extended line of first transmitting field source and the second transmitting field source is the origin of coordinates, and rectangular co-ordinate is set up respectively
System.Wherein, the line of position to be measured and the origin of coordinates is designated as the first line.
Specifically, the first rectangular coordinate system is set up using the first transmitting field source as x1 axles.Position to be measured is to the first transmitting field source
The first current signal produce feedback signal, the feedback signal position to be measured formed the first magnetic field intensity signal, by this first
Magnetic field intensity signal calculates in the radial magnetic field strength component and tangential magnetic field strength component of position to be measured and obtains first electromagnetism
First magnetic field strength component of the field intensity signal along x1 axles and/or y1 axles.
On first magnetic field intensity signal position to be measured radial magnetic field strength component Hr1With tangential magnetic field intensity point
AmountCalculating be expressed as follows:
In formula, i1For first transmitting field source the first electric current, r be position to be measured to the distance of the origin of coordinates, k is wave number, a
For the length of the first transmitting field source Electric Dipole, φ is the angle between the first line and the first transmitting field source, and i is imaginary unit,WithRespectively withFor the first and second class Bessel function of imaginary argument of argument, I0、
I1、K0And K1Subscript " 0 " or " 1 " represent Bessel function exponent number.
By the first electromagnetic field intensity signal position to be measured radial magnetic field strength component Hr1With tangential magnetic field strength componentThe first electromagnetic field intensity signal is obtained along x1 axles and/or the specific formula for calculation of the first electromagnetic field intensity component of y1 axles
For:
In formula, Hx1Component for the first magnetic field intensity along x1 axles, Hy1Component for the first magnetic field intensity along y1 axles, Hr1、Respectively the first magnetic field intensity is in the radial magnetic field strength component and tangential magnetic field strength component of position to be measured, and φ is first
Angle between line and the first transmitting field source.
Specifically, the second rectangular coordinate system is set up using the second transmitting field source as x2 axles.Position to be measured is to the second transmitting field source
The second current signal produce feedback signal, the feedback signal position to be measured formed the second magnetic field intensity signal, by this second
Magnetic field intensity signal calculates in the radial magnetic field strength component and tangential magnetic field strength component of position to be measured and obtains second magnetic field
Second magnetic field strength component of the strength signal along x2 axles and/or y2 axles.
On second magnetic field intensity signal position to be measured radial magnetic field strength component Hr2With tangential magnetic field intensity point
AmountCalculating be expressed as follows:
In formula, i2For second transmitting field source the second electric current, r be position to be measured to the distance of the origin of coordinates, k is wave number, b
For the length of the second transmitting field source Electric Dipole, θ is the angle between the second transmitting field source and the first transmitting field source, and i is imaginary unit,WithRespectively withFor the first and second class Bessel function of imaginary argument of argument, I0、
I1、K0And K1Subscript " 0 " or " 1 " represent Bessel function exponent number.
By the second magnetic field intensity signal position to be measured radial magnetic field strength component Hr2With tangential magnetic field strength componentObtaining specific formula for calculation of second magnetic field intensity signal along x2 axles and/or the second magnetic field strength component of y2 axles is:
In formula, Hx2Component for the second magnetic field intensity along x2 axles, Hy2Component for the second magnetic field intensity along y2 axles, Hr2、Respectively the second magnetic field intensity is in the radial magnetic field strength component and tangential magnetic field strength component of position to be measured, and θ is second
Launch the angle between field source and the first transmitting field source, φ is the angle between the first line and the first transmitting field source.
It is shown in Figure 2, it is strong for ease of calculating the first magnetic field intensity and the second magnetic field in another specific embodiment
Magnetic field intensity is superimposed formed by degree superposition, by the second rectangular coordinate system by way of rotation, makes former second rectangular coordinate system
X2 axles and the first rectangular coordinate system x1 overlapping of axles.Second magnetic field strength component is turned after the rotation of the second rectangular coordinate system
It is changed to the 3rd magnetic field strength component.
Specifically, the second rectangular coordinate system is rotated to after being overlapped with the first rectangular coordinate system, in former second rectangular coordinate system
Position relationship do not change, it is, being merely the rotation of rectangular coordinate system.
In another specific embodiment, the second magnetic field strength component is converted into after the 3rd magnetic field strength component, the
Former coordinate of two magnetic field strength components in the second rectangular coordinate system is converted into the new coordinate in the first rectangular coordinate system, and
There is corresponding relation between the new coordinate and former coordinate, i.e., exist between the 3rd magnetic field strength component and the second magnetic field strength component
Corresponding relation.
Specifically, former coordinate representation of second magnetic field strength component in the second rectangular coordinate system is along x2 axles/y2 axles
Component, respectively Hx2、Hy2, the second magnetic field intensity is converted into after the 3rd magnetic field intensity, and the 3rd magnetic field strength component is straight first
Coordinate representation in angular coordinate system is the component H' along x1 axles/y1 axlesx1、H'y1。
Corresponding relation statement between 3rd magnetic field strength component and the second magnetic field strength component is as follows:
In formula, H'x1、H'y1Component for the 3rd magnetic field strength component in the first rectangular coordinate system along x1 axles/y1 axles,
Hx2、Hy2Component for the second magnetic field strength component in the second rectangular coordinate system along x2 axles/y2 axles, θ be second transmitting field source with
Angle between first transmitting field source.
It is shown in Figure 3, in another specific embodiment, with the first magnetic field strength component and the second magnetic in step S2
Based on field intensity component, to build the first relation function.The step of building first relation function be:
Step S21, using the described first transmitting field source as x-axis, set up overlapped with first rectangular coordinate system the 3rd straight
Angular coordinate system;
Step S22, by first magnetic field strength component and the 3rd magnetic field strength component, obtain it is described first transmitting field source
The component for the superposition magnetic field intensity to be formed, the superposition magnetic field strength component are superimposed in the position to be measured with the second transmitting field source
For the component in the 3rd rectangular coordinate system along x-axis and/or y-axis.
The 3rd rectangular coordinate system set up is overlapped with the first rectangular coordinate system, for statement superposition magnetic field intensity.
First magnetic field intensity and the second magnetic field intensity are superimposed at position to be measured, the magnetic field intensity formed after superposition, to increase
Signal intensity of the heavy current signal at position to be measured.Therefore, based on first magnetic field intensity and the second magnetic field intensity, it can obtain
To the superposition magnetic field intensity on position to be measured.Correspondingly, can by the first magnetic field strength component and the second magnetic field strength component
The component of the superposition magnetic field intensity on position to be measured is obtained, that is, forms the first pass for calculating the superposition magnetic field strength component
It is function.
In another specific embodiment, for being consistent property, make calculating superposition of electric field intensity and superposition magnetic field intensity
Relation function be expressed as the first relation function.Have on the first relation function for calculating the superposition magnetic field strength component
Body is expressed as follows:
In formula, Hx、HyRespectively it is superimposed component of the magnetic field intensity along x-axis/y-axis, Hx1、Hy1Respectively the first magnetic field intensity edge
The component of x1 axles/y1 axles, H'x1、H'y1Component of respectively the 3rd magnetic field intensity along x1 axles/y1 axles, Hx2、Hy2Respectively the second magnetic
Component of the field intensity along x2 axles/y2 axles, θ is the angle between the first transmitting field source and the second transmitting field source.
In another specific embodiment, in step S3, structure is obtained after the first relation function, based on first relation
Function, can acquire the optimal value of superposition magnetic field intensity, and its specific steps includes:
Step S31, set up it is described first transmitting field source the first electric current with described second transmitting field source the second electric current it
Between the second relation function on phase difference;
Step S32, by first relation function and the second relation function, take different phase difference values, calculating is obtained not
Same superposition magnetic field strength component, until obtaining the optimal value of the superposition magnetic field strength component.
In the actual measurement process in scene, when sending current signal to position to be measured, to adjust feedback signal strength
One importance, is the attribute by adjusting the current signal that transmitting field source is provided, the size of the attribute including electric current and
Direction.After the generation equipment for launching field source is determined, the regulation of size of current is more inconvenient, and adjusts current generation device conveying
The method of the phase of electric current, then be easily achieved.
During using two transmitting field sources, by adjusting the phase difference value between its corresponding first electric current and the second electric current,
To adjust the superposition magnetic field intensity after the first magnetic field intensity and the second magnetic field intensity are superimposed in position to be measured to optimal value, so as to reach
To the purpose for eliminating small-signal area, make the superposition magnetic field intensity signal intensity finally given high.
Only by adjusting the phase difference between the first electric current and the second electric current, the strong of position feed back signal to be measured can be optimized
Degree, substantially increases the laying flexibility of transmitting field source, reduces transmitting field source addressing and layout difficulty, can be according to practically
Matter environment, avoids unfavorable terrain, improves the flexibility of field work efficiency and construction, meanwhile, also improve the accuracy of measurement
And reliability.
In another specific embodiment, the second relation function between the first electric current and the second electric current on phase difference
Referring to formula (7).Specifically, different phase differences are taken respectivelyValue, can be calculated by formula (13) and (7) and obtain array edge
The component H of x-axis and/or y-axisx、Hy.Each phase differenceValue i.e. one group of H of correspondencexValue, HyValue.In resulting array HxValue,
HyIn value, one group of optimal or maximum numerical value H is chosenxValue, HyValue, it is now correspondingValue is optimum angle difference.
In another specific embodiment, the span of phase difference value is 0 ° -360 °.Due to the change of current phase
When showing numerically, there is certain repeatability, therefore, in actual mechanical process, the span of phase difference value is limited
At 0 ° -90 °.Take different phase differences successively in the range of 0 ° -90 °Value, is calculated by formula (13) and (7) and obtained
Corresponding Hx、HyValue, until obtaining optimal or maximum superposition field strength values.
It is understood that in actual mechanical process, when it is determined that the first transmitting field source, the second transmitting field source and treating location
Postpone, when laying two transmitting field sources, the size of its first electric current, first launch field source Electric Dipole length, in midair
Between resistivity, the distance of position to be measured to the first rectangular coordinate system origin be known quantity;At the same time, the size of the second electric current,
The length of the Electric Dipole of second transmitting field source, half space resistivity, position to be measured are to the distance of the second rectangular coordinate system origin
Known quantity;In addition, the angle between the angle between two transmitting field sources, and position to be measured and the first transmitting field source is also known
Amount, i.e. only input different phase differencesValue, just can calculate and obtain corresponding Hx、HyValue, so that it is poor to obtain optimum angleValue.
Obtain the optimum angle poorAfter value, the actual phase difference of the first transmitting field source and the second transmitting field source is adjusted to the optimal phase
Potential differenceValue, i.e., can obtain the good feedback signal of signal condition in position to be measured, so as to improve in actual measurement process
Measurement accuracy and reliability.
Further, it is to be appreciated that by formula (13) and (7), and take different phase differencesValue, obtains array pair
The H answeredx、HyValue, is carrying out the optimal H of selectionx、HyWhile value, the H with special characteristic can be also selectedx、HyValue, to enter
The corresponding analysis of row.
The a kind of of the present invention is used to regulate and control the method for measurement signal intensity, and field sources are launched using at an angle to each other two, to
Position to be measured sends current signal, the phase difference value between current signal by adjusting two transmitting field sources, eliminates faint letter
Number area, is superimposed field intensity signal, so as to improve measurement to strengthen formed by two feedback signals superposition that position to be measured is formed
Accuracy and reliability.Meanwhile, two transmitting field source angle and position to be measured between relative position relation do not have it is special will
Ask, greatly reduce layout difficulty during measurement, improve efficiency of construction.
It is shown in Figure 4, the comparison diagram respectively between adjustment current signal before and after phase difference value, Fig. 4 a are not enter line phase
Electric-field intensity component, the radiation pattern figure of magnetic field strength component when difference is adjusted, after Fig. 4 b is progress phase difference value regulations
The radiation pattern figure of electric-field intensity component, magnetic field strength component.
Finally, method of the invention is only preferably embodiment, is not intended to limit the scope of the present invention.It is all
Within the spirit and principles in the present invention, any modification, equivalent substitution and improvements made etc. should be included in the protection of the present invention
Within the scope of.
Claims (10)
1. a kind of be used to regulate and control the method for measurement signal intensity, for analyzing the nagneto-telluric field regularity of distribution, it is characterised in that bag
Include:
Step S1, using two transmitting field sources, respectively to position to be measured send current signal;
Step S2, the position to be measured are fed back to the current signal of described two transmitting field sources and form two different electricity respectively
Magnetic field, based on the field strength component of described two different electromagnetic fields, to build the first relation function;
Step S3, the superimposed field formed based on first relation function, acquisition after described two different electromagnetic fields superpositions
Strong optimal value.
2. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 1, it is characterised in that:The superposition field strength
For electric-field intensity or magnetic field intensity.
3. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 2, it is characterised in that:Described in step S1
Angular range between two transmitting field sources is 0 ° -360 °.
4. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 2, it is characterised in that described in step S2
The field strength component of two different electromagnetic fields, is specifically included:Described two transmitting field sources include the first transmitting field source and the second hair
Penetrate field source, using described first transmitting field source with second launch field source intersection point as the origin of coordinates,
The first rectangular coordinate system is set up by x1 axles of the described first transmitting field source, and by the feedback signal of the described first transmitting field source
Radial field component and tangential field component in the position to be measured, calculating obtain the first transmitting field source described to be measured
Position is along x1 axles and/or the first field strength component of y1 axles;
The second rectangular coordinate system is set up by x2 axles of the described second transmitting field source, and by the feedback signal of the described second transmitting field source
Radial field component and tangential field component in the position to be measured, calculating obtain the second transmitting field source described to be measured
Position is along x2 axles and/or the second field strength component of y2 axles.
5. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 4, it is characterised in that:Second field strength
The method that component is rotated by reference axis is converted into first rectangular coordinate system the 3rd along x1 axles and/or y1 axles
Strong component.
6. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 5, it is characterised in that the reference axis rotation
Turn particular content be:By the x2 axles/y2 axles corresponding rotation to the x1 axles/y1 overlapping of axles, second rectangular co-ordinate
Coordinate in system sets up corresponding relation in the first rectangular coordinate system with the coordinate between corresponding new coordinate.
7. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 5, it is characterised in that built in step S2
The particular content of first relation function includes:
Step S21, using described first transmitting field source as x-axis, set up overlapped with first rectangular coordinate system the 3rd right angle seat
Mark system;
Step S22, by the first field strength component and the 3rd field strength component, obtain the first transmitting field source and the second launching site
Source is superimposed the superimposed field strong component to be formed in the position to be measured, and the superimposed field strong component is the edge in the 3rd rectangular coordinate system
The field strength component of x-axis and/or y-axis, to form the first relation function.
8. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 4, it is characterised in that be based in step S3
First relation function, obtains the optimal value of the superposition field strength formed after described two different electromagnetic field superpositions, specifically
Including:
Closed between step S31, the first electric current for setting up the first transmitting field source and the second transmitting field source the second electric current
In the second relation function of phase difference;
Step S32, based on first relation function and the second relation function, take different phase difference values, calculating obtains difference
Superposition field intensity value, until obtain it is described superposition field strength optimal value.
9. a kind of method for regulating and controlling measurement signal intensity as claimed in claim 8, it is characterised in that:Institute in step S31
The scope for stating phase difference is 0 ° -360 °.
10. a kind of method for regulating and controlling measurement signal intensity as described in claim 6 or 8, it is characterised in that described first
Relation function is expressed as:
Ex=Ex1+E′x1=Ex1+Ex2cosθ-Ey2sinθ
Ey=Ey1+E′y1=Ey1+Ex2sinθ+Ey2cosθ
Hx=Hx1+H′x1=Hx1+Hx2cosθ-Hy2sinθ
Hy=Hy1+H′y1=Hx1+Hx2sinθ+Hy2cosθ
In formula, Ex、EyRespectively component of the superposition of electric field intensity along x-axis/y-axis, Ex1、Ey1Respectively the first electric-field intensity is along x1
The component of axle/y1 axles, Ex′1、Ey′1Component of respectively the 3rd electric-field intensity along x1 axles/y1 axles, Ex2、Ey2Respectively the second electric field
Component of the intensity along x2 axles/y2 axles;Hx、HyRespectively it is superimposed component of the magnetic field intensity along x-axis/y-axis, Hx1、Hy1Respectively first
Component of the magnetic field intensity along x1 axles/y1 axles, Hx′1、Hy′1Component of respectively the 3rd magnetic field intensity along x1 axles/y1 axles, Hx2、Hy2Point
Wei not component of second electromagnetic field intensity along x2 axles/y2 axles;θ is the angle between the first transmitting field source and the second transmitting field source.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710258059.4A CN107092033B (en) | 2017-04-19 | 2017-04-19 | A method of for regulating and controlling measuring signal intensity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710258059.4A CN107092033B (en) | 2017-04-19 | 2017-04-19 | A method of for regulating and controlling measuring signal intensity |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107092033A true CN107092033A (en) | 2017-08-25 |
CN107092033B CN107092033B (en) | 2019-06-11 |
Family
ID=59638252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710258059.4A Active CN107092033B (en) | 2017-04-19 | 2017-04-19 | A method of for regulating and controlling measuring signal intensity |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107092033B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781835A (en) * | 2019-01-28 | 2019-05-21 | 贵州电网有限责任公司 | A kind of detection method of electric pole tower grounded screen mesh shape |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5646524A (en) * | 1992-06-16 | 1997-07-08 | Elbit Ltd. | Three dimensional tracking system employing a rotating field |
CN102551711A (en) * | 2010-12-09 | 2012-07-11 | 同济大学 | Method for measuring magnetic field by aid of superconducting quantum interference device and positioning current source of magnetic field |
CN103389514A (en) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | MT (magnetotelluric) denoising device and method |
CN105204073A (en) * | 2015-09-18 | 2015-12-30 | 中南大学 | Tensor apparent conductivity measurement method |
CN105866844A (en) * | 2016-06-12 | 2016-08-17 | 何继善 | Device and method for measuring whole region measuring current source frequency domain electric field horizontal x component |
CN105891892A (en) * | 2016-06-12 | 2016-08-24 | 何继善 | Device and method for whole-region measurement of vertical components of current source frequency domain magnetic fields |
-
2017
- 2017-04-19 CN CN201710258059.4A patent/CN107092033B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5646524A (en) * | 1992-06-16 | 1997-07-08 | Elbit Ltd. | Three dimensional tracking system employing a rotating field |
CN102551711A (en) * | 2010-12-09 | 2012-07-11 | 同济大学 | Method for measuring magnetic field by aid of superconducting quantum interference device and positioning current source of magnetic field |
CN103389514A (en) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | MT (magnetotelluric) denoising device and method |
CN105204073A (en) * | 2015-09-18 | 2015-12-30 | 中南大学 | Tensor apparent conductivity measurement method |
CN105866844A (en) * | 2016-06-12 | 2016-08-17 | 何继善 | Device and method for measuring whole region measuring current source frequency domain electric field horizontal x component |
CN105891892A (en) * | 2016-06-12 | 2016-08-24 | 何继善 | Device and method for whole-region measurement of vertical components of current source frequency domain magnetic fields |
Non-Patent Citations (1)
Title |
---|
李帝铨: "E-Ex广域电磁法三维数值模拟", 《中国有色金属学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781835A (en) * | 2019-01-28 | 2019-05-21 | 贵州电网有限责任公司 | A kind of detection method of electric pole tower grounded screen mesh shape |
CN109781835B (en) * | 2019-01-28 | 2022-06-28 | 贵州电网有限责任公司 | Method for detecting grid shape of grounding grid of pole tower |
Also Published As
Publication number | Publication date |
---|---|
CN107092033B (en) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2829617C (en) | Multi-component electromagnetic prospecting apparatus and method of use thereof | |
CN101116010B (en) | A buried object locating and tracing method and system employing principal components analysis for blind signal detection | |
CN108873083B (en) | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method | |
CN105759316B (en) | A kind of method and apparatus of rectangular loop source transient electromagnetic detecting | |
CN109283499B (en) | Radar equation-based three-dimensional visualization method for detection range under active interference | |
WO2006135510A1 (en) | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys | |
CN104854480A (en) | Apparatus and methods to find a position in an underground formation | |
Falco et al. | Fracture characterisation using geoelectric null-arrays | |
CN106855635A (en) | A kind of method and apparatus of the bin attributes information for determining orthogonal observation system | |
CN105549100B (en) | Transient Electromagnetic Method in Mine apparatus and method based on U-shaped helical source | |
CN112433257A (en) | Vector electrical method scanning method based on frequency division electrical method | |
CN107092033A (en) | It is a kind of to be used to regulate and control the method for measurement signal intensity | |
CN103995286B (en) | A kind of acoustic receiver transducer and orientation, down-hole noise leakage detection apparatus and leak hunting method | |
Antti et al. | Separation of the geomagnetic variation field on the ground into external and internal parts using the spherical elementary current system method | |
WO2013063659A1 (en) | Improved metal detection method | |
CN108957589A (en) | More exploration geophysical field signal receiving sensors and sensor string, observation system | |
CN106873044B (en) | Array SQUID nuclear magnetic resonance underground water detecting device and imaging method | |
Smith | An airborne electromagnetic system with a three-component transmitter and three-component receiver capable of detecting extremely conductive bodies | |
Xue et al. | Comparison of frequency-domain electromagnetic responses at different offsets based on exact calculation of grounded-wire source | |
WO2021098667A1 (en) | Full-coverage anti-interference controlled-source electromagnetic detection method | |
AU2015249137B2 (en) | Multi-Component Electromagnetic Prospecting Apparatus and Method of Use Thereof | |
CN208654347U (en) | More exploration geophysical field signal receiving sensors and sensor string, observation system | |
Wu et al. | High-precision inversion of buried depth inurban underground iron pipelines based on AM-PSO algorithmfor magnetic anomaly | |
CN108519622B (en) | Underground electric target detection method and device based on natural field source excitation | |
Liu et al. | Research on ground-airborne frequency-domain electromagnetic rapid imaging method based on space magnetic gradient anomaly |
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 |