CN117270063A - Signal frequency acquisition method and device - Google Patents
Signal frequency acquisition method and device Download PDFInfo
- Publication number
- CN117270063A CN117270063A CN202311180437.3A CN202311180437A CN117270063A CN 117270063 A CN117270063 A CN 117270063A CN 202311180437 A CN202311180437 A CN 202311180437A CN 117270063 A CN117270063 A CN 117270063A
- Authority
- CN
- China
- Prior art keywords
- autocorrelation
- frequency
- sampling
- signal
- data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005070 sampling Methods 0.000 claims abstract description 98
- 238000004364 calculation method Methods 0.000 claims abstract description 48
- 238000005553 drilling Methods 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 13
- 238000004891 communication Methods 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 3
- 230000001360 synchronised effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
-
- 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/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/30—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Electromagnetism (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The application discloses a signal frequency acquisition method and device. The signal frequency acquisition method comprises the following steps: acquiring sampling point data; grouping the acquired sampling point data to form at least four groups of sampling groups; performing autocorrelation calculation on each group of sample group data so as to obtain autocorrelation parameters; and calculating the frequency and phase information of the received signal according to the autocorrelation parameters. The signal frequency acquisition method can solve the problem that the frequency of the signal received by the receiving coil is not synchronous with the frequency of the signal transmitted by the transmitting coil in the logging process of the azimuth electromagnetic wave resistivity instrument while drilling. Through frequency synchronization algorithm research, accurate frequency of the received signal can be obtained so as to acquire amplitude and phase information of the received signal subsequently.
Description
Technical Field
The application relates to the technical field of azimuth electromagnetic wave resistivity instruments while drilling, in particular to a signal frequency acquisition method and a signal frequency acquisition device.
Background
The logging instrument comprises a transmitting coil and a receiving coil, wherein a transmitting antenna transmits signals with fixed frequency (such as 400kHz and 2 MHz) in the application process, and the receiving antenna receives corresponding transmitting signals. However, when the receiving antenna receives a signal, the frequency of the signal may be shifted, so that the frequency of the received signal and the frequency of the transmitted signal at this time are not uniform. The frequency of the signal needs to be determined by calculation with a corresponding algorithm.
It is therefore desirable to have a solution that overcomes or at least alleviates at least one of the above-mentioned drawbacks of the prior art.
Disclosure of Invention
It is an object of the present invention to provide a signal frequency acquisition method that overcomes or at least alleviates at least one of the above-mentioned drawbacks of the prior art.
In one aspect of the present invention, a signal frequency acquisition method is provided for an azimuth while drilling electromagnetic wave resistivity instrument, the signal frequency acquisition method comprising:
acquiring sampling point data;
grouping the acquired sampling point data to form at least four groups of sampling groups;
performing autocorrelation calculation on each group of sample group data so as to obtain autocorrelation parameters;
and calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
Optionally, the calculating the frequency and the phase information of the received signal according to the autocorrelation parameters is obtained by adopting the following formula:
wherein,
f 0 ' is the frequency of the received signal, nf 0 For the sampling frequency, angle represents a function of complex phase angle, a represents a first parameter directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents the direct phase with the phase of the received signalA fourth parameter of,For receiving signal phase angle information, the a, b, c, d constitutes the autocorrelation parameters.
Optionally, the autocorrelation parameters include a first parameter, a second parameter, a third parameter, and a fourth parameter;
the performing autocorrelation calculation on each set of samples, respectively, so as to obtain autocorrelation parameters includes:
the first parameter, the second parameter, the third parameter and the fourth parameter are obtained by adopting the following formula:
wherein,
z 1 ＝2N-4y 1 ,z 2 ＝2N-4y 2 ,z 3 ＝2N-4y 3 ,z 4 ＝2N-4y 4 ；
n represents the data in each set of sample groups, y1 is the result of the autocorrelation calculation of the first set of sample group data, reflects the autocorrelation of the first set of sample group data, y2 is the result of the autocorrelation calculation of the second set of sample group data, reflects the autocorrelation of the second set of sample group data, y3 is the result of the autocorrelation calculation of the third set of sample group data, reflects the autocorrelation of the third set of sample group data, y4 is the result of the autocorrelation calculation of the fourth set of sample group data, reflects the autocorrelation of the fourth set of sample group data, and a represents the first parameter that is directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents a fourth parameter directly related to the phase of the received signal, said a, b, c, d comprising said autocorrelation parameters, z 1 Representing the result (y) of the autocorrelation calculation with the first set of samples 1 ) The parameter of the correlation is set to be,z 2 representing the result (y) of the autocorrelation calculation with the first set of samples 2 ) Related parameters, z 3 Representing the result (y) of the autocorrelation calculation with the first set of samples 3 ) Related parameters, z 4 Representing the result (y) of the autocorrelation calculation with the first set of samples 4 ) The related parameters, a, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 ) The parameters of the correlation, B, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 、y 4 ) The parameter of correlation, C, represents the result of the autocorrelation calculation (y 2 、y 3 、y 4 ) Related parameters, z 1 、z 2 、z 3 、z 4 Both A, B, C are intermediate parameters introduced to simplify the expression of a, b, c, d.
Optionally, the y1, y2, y3, y4 are respectively obtained by adopting the following formulas:
wherein,
n represents the number of sampling points of each set of sampling group data,
Nf 0 For sampling frequency f 0 ' is the frequency of the received signal, n represents the position mark of the sampling point in the sampled data, pi represents the circumference ratio, ">For receiving signal phase angle information.
Optionally, the acquiring the sampling point data includes:
using Nf 0 Sampling the receiving end for sampling frequency so as to obtain sampling data, wherein the sampling point data are as follows:
(n 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 s (n+2)); wherein,
n=0, n=1, … …, n=n-1, n=n, n=n+1, n=n+2 denote position marks of sampling points in the sampling data.
Optionally, grouping the acquired sample point data to form a plurality of sample groups includes:
the (n) 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 S (N+2)) are grouped, wherein the first group is from the 0 th sampling point to the N-1 th sampling point; the second group is from the 1 st sampling point to the N th sampling point; the third group is from the 2 nd sampling point to the (n+1) th sampling point; the fourth group is the 3 rd sample point to the n+2 th sample point.
Optionally, the sampling point data is obtained using the following formula:
wherein,
Nf 0 for sampling frequency f 0 ' is the frequency of the received signal, at which point f 0 ′≠f 0 Pi represents the circumference ratio,For receiving signal phase angle information.
The application also provides a signal frequency acquisition device, the signal frequency acquisition device includes:
the sampling point data acquisition module is used for acquiring sampling point data;
the grouping module is used for grouping the acquired sampling point data so as to form at least four groups of sampling groups;
the autocorrelation technology module is used for carrying out autocorrelation calculation on each group of sampling group data so as to obtain autocorrelation parameters;
and the signal frequency and phase information acquisition module is used for calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
The application also provides an electronic device, which comprises: the device comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory are communicated with each other through the communication bus; the memory has stored therein a computer program which, when executed by the processor, causes the processor to perform the steps of the signal frequency acquisition method as described above.
The present application also provides a computer readable storage medium storing a computer program executable by an electronic device, which when run on the electronic device is capable of implementing the steps of the signal frequency acquisition method as described above.
The beneficial effects are that:
the signal frequency acquisition method can solve the problem that the frequency of the signal received by the receiving coil is not synchronous with the frequency of the signal transmitted by the transmitting coil in the logging process of the azimuth electromagnetic wave resistivity instrument while drilling. Through frequency synchronization algorithm research, accurate frequency of the received signal can be obtained so as to acquire amplitude and phase information of the received signal subsequently.
Drawings
Fig. 1 is a flowchart of a signal frequency acquisition method according to an embodiment of the present application.
Fig. 2 is an electronic device for implementing the signal frequency acquisition method shown in fig. 1.
Fig. 3 is a schematic structural diagram of an antenna of an electromagnetic wave resistivity instrument in azimuth while drilling according to an embodiment of the present application.
Fig. 4 is a schematic diagram of a signal collected by a receiving coil according to a method of an embodiment of the present application.
Detailed Description
In order to make the purposes, technical solutions and advantages of the implementation of the present application more clear, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the accompanying drawings in the embodiments of the present application. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are some, but not all, of the embodiments of the present application. The embodiments described below by referring to the drawings are exemplary and intended for the purpose of explaining the present application and are not to be construed as limiting the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure. Embodiments of the present application are described in detail below with reference to the accompanying drawings.
It should be noted that in the description of the present invention, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
Fig. 1 is a flowchart of a signal frequency acquisition method according to an embodiment of the present application.
The signal frequency acquisition method as shown in fig. 1 is used for an azimuth while drilling electromagnetic wave resistivity instrument, and is characterized in that the signal frequency acquisition method comprises the following steps:
acquiring sampling point data;
grouping the acquired sampling point data to form at least four groups of sampling groups;
performing autocorrelation calculation on each sampling group so as to obtain autocorrelation parameters;
and calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
The signal frequency acquisition method can solve the problem that the frequency of the signal received by the receiving coil is not synchronous with the frequency of the signal transmitted by the transmitting coil in the logging process of the azimuth electromagnetic wave resistivity instrument while drilling. Through frequency synchronization algorithm research, accurate frequency of the received signal can be obtained so as to acquire amplitude and phase information of the received signal subsequently.
In this embodiment, the frequency and phase information of the received signal calculated according to the autocorrelation parameters are obtained by using the following formula:
wherein,
f 0 ' is the frequency of the received signal, nf 0 For the sampling frequency, angle represents a function of complex phase angle, a represents a first parameter directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents a fourth parameter directly related to the phase of the received signal,For receiving signal phase angle information, the a, b, c, d constitutes the autocorrelation parameters.
In this embodiment, the autocorrelation parameters include a first parameter, a second parameter, a third parameter, and a fourth parameter;
the performing autocorrelation calculation on each set of samples, respectively, so as to obtain autocorrelation parameters includes:
the first parameter, the second parameter, the third parameter and the fourth parameter are obtained by adopting the following formula:
wherein,
z 1 ＝2N-4y 1 ,z 2 ＝2N-4y 2 ,z 3 ＝2N-4y 3 ,z 4 ＝2N-4y 4 ；
n represents the data in each set of sample sets, y1 is the autocorrelation calculation of the first set of sample set data, reflects the autocorrelation of the first set of sample set data, y2 is the autocorrelation calculation of the second set of sample set data, reflects the autocorrelation of the second set of sample set data, y3 is the autocorrelation calculation of the third set of sample set data, and reflects the third set of sample set dataThe autocorrelation, y4 is the result of the calculation of the autocorrelation of the fourth set of sample set data, reflecting the autocorrelation of the fourth set of sample set data, a representing the first parameter directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents a fourth parameter directly related to the phase of the received signal, said a, b, c, d comprising said autocorrelation parameters, z 1 Representing the result (y) of the autocorrelation calculation with the first set of samples 1 ) Related parameters, z 2 Representing the result (y) of the autocorrelation calculation with the first set of samples 2 ) Related parameters, z 3 Representing the result (y) of the autocorrelation calculation with the first set of samples 3 ) Related parameters, z 4 Representing the result (y) of the autocorrelation calculation with the first set of samples 4 ) The related parameters, a, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 ) The parameters of the correlation, B, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 、y 4 ) The parameter of correlation, C, represents the result of the autocorrelation calculation (y 2 、y 3 、y 4 ) Related parameters, z 1 、z 2 、z 3 、z 4 Both A, B, C are intermediate parameters introduced to simplify the expression of a, b, c, d.
In this embodiment, the y1, y2, y3, y4 are obtained by the following formulas:
wherein,
n represents the number of sampling points of each set of sampling group data,Nf 0 For sampling frequency f 0 ' position mark for the frequency of the received signal, n representing the sampling point in the sampled data。
In this embodiment, the acquiring the sampling point data includes:
sampling is carried out on a receiving end so as to obtain sampling point data, wherein the sampling point data are as follows:
(n 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 ,S(N+2))。
in this embodiment, grouping the acquired sampling point data, thereby forming a plurality of sampling groups includes:
the (n) 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 S (N+2)) are grouped, wherein the first group is from the 0 th sampling point to the N-1 th sampling point; the second group is from the 1 st sampling point to the N th sampling point; the third group is from the 2 nd sampling point to the (n+1) th sampling point; the fourth group is the 3 rd sample point to the n+2 th sample point.
In this embodiment, the following formula is adopted to obtain sampling point data:
wherein,
Nf 0 for sampling frequency f 0 ' is the frequency of the received signal, at which point f 0 ′≠f 0 。
The present application is described in further detail below by way of examples, which are not to be construed as limiting the present application in any way.
Referring to fig. 3, in the present embodiment, the antenna structure of the azimuth electromagnetic wave resistivity instrument while drilling is shown in fig. 3, T1-T4 are transmitting antennas, R1-R4 are receiving antennas, wherein R1 and R2 are axial receiving antennas, and R3 and R4 are azimuth receiving antennas. In the figure, R1/R2 receives signals from T1 and T4 respectively, and R1/R2/R3/R4 receives signals from T1/T2/T3/T4 respectively in the practical application process. It can be appreciated that the antenna structure of the electromagnetic wave resistivity instrument while drilling azimuth of the present application is in the prior art, for example, the patent No. ZL202210621327.5, and the patent name is a measurement while drilling device described in the calibration method of the measurement while drilling device.
In the present embodiment, it is assumed that the transmission signal is:
wherein f 0 For the frequency of the transmitted signal,for the initial phase of the transmit/receive signal, the receiving end can obtain the discrete signal by sampling as follows:
wherein Nf is 0 For sampling frequency f 0 ' is the frequency of the received signal, at which point f 0 ′≠f 0 . Discrete signals are obtained by sampling, the sampled signals being expressed as:
(n 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 s (N+2)), dividing the sampling points into four groups, wherein the first group is 0-N-1 sampling points; the second group is 1-N sampling points; the third group is 2-N+1 sampling points; the fourth group is 3-N+2 sampling points. Expressed using an e-exponential function as:
using the calculation of the auto-correlation,
wherein:
and (3) summing and calculating to obtain:
and let z 1 ＝2N-4y 1 ,z 2 ＝2N-4y 2 ,z 3 ＝2N-4y 3 ,z 4 ＝2N-4y 4 Then it is simplified as:
the formula can be obtained by the primordial elimination method:
wherein,the values of a, b, c, d in the middle can be obtained by sampling the original signal and integrating the sampled values. It can be seen that:
in the above formula, angle represents a phase angle function for obtaining a complex number.
In practical use, the received signal is set to be sinusoidal, the amplitude is set to be 1V, and the frequency of the received signal is f 0 ' the received signal frequency taken here is the verification frequency set according to the requirement in the present embodiment (specifically, since simulation verification is being performed at this time, i amThe frequency of the receiver signal is set to 15Hz, here only a signal with a frequency of 15Hz is generated. At this time, the frequency of the signal at the transmitting end is 14Hz, and the frequency at the receiving end is only known before the receiving end receives the signal, and the frequency at the receiving end is not clear to be 15Hz. So 15Hz is set here only for verification use. ) Initial phase isThe above information is the signal to be received at the receiving end, and when the signal is received, we do not know the frequency of the received signal, and the frequency according to the transmitted signal is f 0 =14 Hz, so the sampling signal frequency chosen according to the Nyquist sampling theorem is f s Specifically, 16 sampling points in one period are generally selected, i.e., 16×14=224 Hz. Using the sampling frequency, 19 points were acquired (n=16), and the signal was acquired as shown in fig. 4.
By integrating the autocorrelation calculation and combining (9), it is possible to obtain
From the calculation of (10), it can be seen that:
through the calculation process, the frequency and phase information of the received signal can be accurately obtained.
The application also provides a signal frequency acquisition device, which comprises a sampling point data acquisition module, a grouping module, an autocorrelation technology module and a signal frequency and phase information acquisition module, wherein,
the sampling point data acquisition module is used for acquiring sampling point data;
the grouping module is used for grouping the acquired sampling point data so as to form at least four groups of sampling groups;
the autocorrelation technology module is used for carrying out autocorrelation calculation on each group of sampling group data so as to obtain autocorrelation parameters;
the signal frequency and phase information acquisition module is used for calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
The application also provides an electronic device, which comprises: the device comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory are communicated with each other through the communication bus; the memory has stored therein a computer program which, when executed by the processor, causes the processor to perform the steps of the signal frequency acquisition method as described above.
The present application also provides a computer readable storage medium storing a computer program executable by an electronic device, which when run on the electronic device is capable of implementing the steps of the signal frequency acquisition method as described above.
It will be appreciated that the description of the method above applies equally to the description of the device.
Fig. 2 is an exemplary block diagram of an electronic device capable of implementing the signal frequency acquisition method provided according to one embodiment of the present application.
As shown in fig. 2, the electronic device includes an input device 501, an input interface 502, a central processor 503, a memory 504, an output interface 505, and an output device 506. The input interface 502, the central processing unit 503, the memory 504, and the output interface 505 are connected to each other through a bus 507, and the input device 501 and the output device 506 are connected to the bus 507 through the input interface 502 and the output interface 505, respectively, and further connected to other components of the electronic device. Specifically, the input device 501 receives input information from the outside, and transmits the input information to the central processor 503 through the input interface 502; the central processor 503 processes the input information based on computer executable instructions stored in the memory 504 to generate output information, temporarily or permanently stores the output information in the memory 504, and then transmits the output information to the output device 506 through the output interface 505; the output device 506 outputs the output information to the outside of the electronic device for use by the user.
That is, the electronic device shown in fig. 2 may also be implemented to include: a memory storing computer-executable instructions; and one or more processors that, when executing the computer-executable instructions, implement the signal frequency acquisition method described in connection with fig. 1.
In one embodiment, the electronic device shown in FIG. 2 may be implemented to include: a memory 504 configured to store executable program code; the one or more processors 503 are configured to execute the executable program code stored in the memory 504 to perform the signal frequency acquisition method in the above-described embodiments.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.
Computer-readable media include both permanent and non-permanent, removable and non-removable media, and the media may be implemented in any method or technology for storage of information. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
Furthermore, it is evident that the word "comprising" does not exclude other elements or steps. A plurality of units, modules or means recited in the apparatus claims can also be implemented by means of software or hardware by means of one unit or total means. The terms first, second, etc. are used to identify names, and not any particular order.
The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The processor referred to in this embodiment may be a central processing unit (Central Processing Unit, CPU), or other general purpose processor, digital signal processor (Digital Signal Processor, DSP), application specific integrated circuit (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may be used to store computer programs and/or modules, and the processor may perform various functions of the apparatus/terminal device by executing or executing the computer programs and/or modules stored in the memory, and invoking data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program (such as a sound playing function, an image playing function, etc.) required for at least one function, and the like; the storage data area may store data (such as audio data, phonebook, etc.) created according to the use of the handset, etc. In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as a hard disk, memory, plug-in hard disk, smart Media Card (SMC), secure Digital (SD) Card, flash Card (Flash Card), at least one disk storage device, flash memory device, or other volatile solid-state storage device.
In this embodiment, the modules/units of the apparatus/terminal device integration may be stored in a computer readable storage medium if implemented in the form of software functional units and sold or used as a separate product. Based on such understanding, the present invention may implement all or part of the flow of the method of the above embodiment, or may be implemented by hardware related to the instructions of a computer program, where the computer program may be stored in a computer readable storage medium, and when executed by a processor, the computer program may implement the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, executable files or in some intermediate form, etc. The computer readable medium may include: any entity or device capable of carrying computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth.
It should be noted that the content of the computer readable medium can be appropriately increased or decreased according to the requirements of the legislation and the practice of the patent in the jurisdiction. While the preferred embodiments have been described, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention, and it is intended that the scope of the invention shall be limited only by the claims appended hereto.
While the invention has been described in detail in the foregoing general description and with reference to specific embodiments thereof, it will be apparent to one skilled in the art that modifications and improvements can be made thereto. Accordingly, such modifications or improvements may be made without departing from the spirit of the invention and are intended to be within the scope of the invention as claimed.
Claims (10)
1. A signal frequency acquisition method for an azimuth while drilling electromagnetic wave resistivity instrument, the signal frequency acquisition method comprising:
acquiring sampling point data;
grouping the acquired sampling point data to form at least four groups of sampling groups;
performing autocorrelation calculation on each group of sample group data so as to obtain autocorrelation parameters;
and calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
2. The signal frequency acquisition method of claim 1 wherein the calculating the received signal frequency and phase information from the autocorrelation parameters is obtained using the formula:
wherein,
f 0 ' is the frequency of the received signal, nf 0 For samplingFrequency, angle represents a function for determining a complex phase angle, and a represents a first parameter directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents a fourth parameter directly related to the phase of the received signal,For receiving signal phase angle information, the a, b, c, d constitutes the autocorrelation parameters.
3. The signal frequency acquisition method of claim 2, wherein the autocorrelation parameters include a first parameter, a second parameter, a third parameter, and a fourth parameter;
the performing autocorrelation calculation on each set of samples, respectively, so as to obtain autocorrelation parameters includes:
the first parameter, the second parameter, the third parameter and the fourth parameter are obtained by adopting the following formula:
wherein,
B＝z 1 z 4 -z 2 z 3 ,/>
z 1 ＝2N-4y 1 ,z 2 ＝2N-4y 2 ,z 3 ＝2N-4y 3 ,z 4 ＝2N-4y 4 ；
n represents the data in each set of sample groups, y1 is the autocorrelation calculation of the first set of sample group data, reflects the autocorrelation of the first set of sample group data, y2 is the autocorrelation calculation of the second set of sample group data, and reflects the second set of sample group dataY3 is the result of the autocorrelation calculation of the third set of sample set data reflecting the autocorrelation of the third set of sample set data, y4 is the result of the autocorrelation calculation of the fourth set of sample set data reflecting the autocorrelation of the fourth set of sample set data, a represents a first parameter directly related to the frequency and phase of the received signal; b represents a second parameter directly related to the frequency and phase of the received signal; c represents a third parameter directly related to the phase of the received signal; d represents a fourth parameter directly related to the phase of the received signal, said a, b, c, d comprising said autocorrelation parameters, z 1 Representing the result (y) of the autocorrelation calculation with the first set of samples 1 ) Related parameters, z 2 Representing the result (y) of the autocorrelation calculation with the first set of samples 2 ) Related parameters, z 3 Representing the result (y) of the autocorrelation calculation with the first set of samples 3 ) Related parameters, z 4 Representing the result (y) of the autocorrelation calculation with the first set of samples 4 ) The related parameters, a, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 ) The parameters of the correlation, B, represent the results of the autocorrelation calculations (y 1 、y 2 、y 3 、y 4 ) The parameter of correlation, C, represents the result of the autocorrelation calculation (y 2 、y 3 、y 4 ) Related parameters, z 1 、z 2 、z 3 、z 4 Both A, B, C are intermediate parameters introduced to simplify the expression of a, b, c, d.
4. A signal frequency acquisition method according to claim 3, wherein y1, y2, y3, y4 are each obtained using the following formula:
n represents the number of sampling points of each set of sampling group data,
Nf 0 For sampling frequency f 0 ' is the frequency of the received signal, n represents the position mark of the sampling point in the sampled data, pi represents the circumference ratio, ">For receiving signal phase angle information.
5. The signal frequency acquisition method of claim 4, wherein the acquiring sample point data comprises:
using Nf 0 Sampling the receiving end for sampling frequency so as to obtain sampling data, wherein the sampling point data are as follows:
(n 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 s (n+2)); wherein,
n=0, n=1, … …, n=n-1, n=n, n=n+1, n=n+2 denote position marks of sampling points in the sampling data.
6. The signal frequency acquisition method of claim 5 wherein grouping the acquired sample point data to form a plurality of groups of samples comprises:
the said
(n 0 ,S(0)),(n 1 ,S(1)),…(n N-1 ,S(N-1)),(n N ,S(N)),(n N+1 ,S(N+1)),(n N+2 S (N+2)) are grouped, wherein the first group is from the 0 th sampling point to the N-1 th sampling point; the second group is from the 1 st sampling point to the N th sampling point; the third group is from the 2 nd sampling point to the (n+1) th sampling point; the fourth group is the 3 rd sample point to the n+2 th sample point.
7. The signal frequency acquisition method of claim 6 wherein the sample point data is acquired using the formula:
wherein,
Nf 0 for sampling frequency f 0 ' is the frequency of the received signal, at which point f 0 ′≠f 0 Pi represents the circumference ratio,For receiving signal phase angle information.
8. A signal frequency acquisition device, characterized in that the signal frequency acquisition device comprises:
the sampling point data acquisition module is used for acquiring sampling point data;
the grouping module is used for grouping the acquired sampling point data so as to form at least four groups of sampling groups;
the autocorrelation technology module is used for carrying out autocorrelation calculation on each group of sampling group data so as to obtain autocorrelation parameters;
and the signal frequency and phase information acquisition module is used for calculating the frequency and phase information of the received signal according to the autocorrelation parameters.
9. An electronic device, the electronic device comprising: the device comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory are communicated with each other through the communication bus; a computer program is stored in a memory, which when executed by a processor causes the processor to perform the steps of the signal frequency acquisition method according to any one of claims 1 to 7.
10. A computer readable storage medium, characterized in that it stores a computer program executable by an electronic device, which, when run on the electronic device, is capable of implementing the steps of the signal frequency acquisition method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311180437.3A CN117270063B (en) | 2023-09-13 | 2023-09-13 | Signal frequency acquisition method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311180437.3A CN117270063B (en) | 2023-09-13 | 2023-09-13 | Signal frequency acquisition method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117270063A true CN117270063A (en) | 2023-12-22 |
| CN117270063B CN117270063B (en) | 2024-08-27 |
Family
ID=89199991
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311180437.3A Active CN117270063B (en) | 2023-09-13 | 2023-09-13 | Signal frequency acquisition method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117270063B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101223747A (en) * | 2005-07-13 | 2008-07-16 | 文卡达·古鲁普拉赛德 | Distance-dependent spectra with uniform sampling spectrometry |
| CN111947957A (en) * | 2020-08-17 | 2020-11-17 | 浙江中控技术股份有限公司 | Vibration signal sampling method and device, storage medium and electronic equipment |
| CN115685355A (en) * | 2022-11-01 | 2023-02-03 | 中国科学院地质与地球物理研究所 | Amplitude-phase extraction method and device for while-drilling azimuth electromagnetic wave resistivity instrument |
| US20230144980A1 (en) * | 2020-04-07 | 2023-05-11 | Zte Corporation | Method and apparatus for estimating frequency offset, electronic device and computer-readable medium |
| CN116299734A (en) * | 2023-04-13 | 2023-06-23 | 中南大学 | Method, device, medium and equipment for acquiring excitation parameters by using spread spectrum signals |
-
2023
- 2023-09-13 CN CN202311180437.3A patent/CN117270063B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101223747A (en) * | 2005-07-13 | 2008-07-16 | 文卡达·古鲁普拉赛德 | Distance-dependent spectra with uniform sampling spectrometry |
| US20230144980A1 (en) * | 2020-04-07 | 2023-05-11 | Zte Corporation | Method and apparatus for estimating frequency offset, electronic device and computer-readable medium |
| CN111947957A (en) * | 2020-08-17 | 2020-11-17 | 浙江中控技术股份有限公司 | Vibration signal sampling method and device, storage medium and electronic equipment |
| CN115685355A (en) * | 2022-11-01 | 2023-02-03 | 中国科学院地质与地球物理研究所 | Amplitude-phase extraction method and device for while-drilling azimuth electromagnetic wave resistivity instrument |
| CN116299734A (en) * | 2023-04-13 | 2023-06-23 | 中南大学 | Method, device, medium and equipment for acquiring excitation parameters by using spread spectrum signals |
Non-Patent Citations (1)
| Title |
|---|
| 和丽真;杨国华;王珺;郭庆明;: "基于随钻电磁波仪器的欠采样DPSD设计与实现", 电子测量技术, no. 12, 23 June 2020 (2020-06-23) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117270063B (en) | 2024-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110727004A (en) | Weil code generation method and device and storage medium | |
| CN106546817A (en) | A kind of Frequency Estimation with feedback function and energy state postulate | |
| US8615031B2 (en) | Signal processing method, correlator, software signal receiver by using code cycle | |
| CN116015492A (en) | Radio station fault diagnosis method, system and equipment | |
| CN106199607B (en) | A kind of the Sounnd source direction localization method and device of microphone array | |
| CN117270063B (en) | Signal frequency acquisition method and device | |
| CN107708145B (en) | Synchronous detection method and synchronous detection equipment | |
| CN113031026A (en) | Ranging code generation method, device, equipment and storage medium | |
| CN116125502B (en) | Navigation signal generation method, device, equipment and storage medium | |
| CN114690214A (en) | GNSS signal capturing method and device and computer storage medium | |
| CN109736793B (en) | Method and equipment for evaluating pore structure of reservoir | |
| US11740270B2 (en) | Pattern generator and built-in-self test device including the same | |
| CN116319219A (en) | Synchronization method based on duplicate correlation in unmanned aerial vehicle and satellite communication and related equipment | |
| CN102307164B (en) | digital frequency estimation method and system | |
| CN112816940A (en) | Target distance estimation method and device based on sound pressure and particle vibration velocity | |
| CN117270075A (en) | Method and device for correcting azimuth signal background value of azimuth electromagnetic wave instrument while drilling | |
| CN117849509B (en) | Method for determining frequency consistency level between channels of frequency converter and related equipment | |
| CN117031913B (en) | Time service method and device based on sine wave signals | |
| CN114018354B (en) | Method, system, equipment and medium for converting dial plate image of pointer type metering device | |
| US10425122B1 (en) | System-on-a-Chip for reception of telemetry messages over a radio frequency channel | |
| CN116988782B (en) | Deep well power supply and data transmission method and system based on single-core cable | |
| RU2303786C2 (en) | Method and device for evaluating bearing frequency of signal | |
| CN113094176A (en) | Hybrid cloud management method and device | |
| CN107290761B (en) | Parallel code phase searching device and method for realizing parallel code phase searching | |
| CN116430114A (en) | High-cost-performance high-precision frequency detection method and system |
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 |