CN117675091A - Logging signal demodulation method and device, computer equipment and storage medium - Google Patents
Logging signal demodulation method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN117675091A CN117675091A CN202211023979.5A CN202211023979A CN117675091A CN 117675091 A CN117675091 A CN 117675091A CN 202211023979 A CN202211023979 A CN 202211023979A CN 117675091 A CN117675091 A CN 117675091A
- Authority
- CN
- China
- Prior art keywords
- function
- functions
- signal
- time domain
- determining
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000003860 storage Methods 0.000 title claims abstract description 20
- 230000006870 function Effects 0.000 claims abstract description 195
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 238000004590 computer program Methods 0.000 claims description 15
- 238000007781 pre-processing Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 13
- 238000004891 communication Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 238000009795 derivation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000003190 augmentative effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000003208 petroleum Substances 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
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The disclosure relates to the technical field of logging, and in particular relates to a logging signal demodulation method, a device, computer equipment and a storage medium. The method comprises the steps of responding to a received signal to be demodulated, and determining a first time domain function and a power spectrum image corresponding to the signal to be demodulated; determining a target harmonic for demodulating the signal from the power spectrum image; processing a plurality of preset second time domain functions and first time domain functions according to the target harmonic wave to obtain a plurality of filter functions and functions to be processed; and demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements. By using the embodiment of the disclosure, the frequency band of the signal to be demodulated, which is demodulated by the determined target harmonic wave, is least interfered by the power spectrum image, so that a code element obtained by demodulation is more accurate.
Description
Technical Field
The disclosure relates to the technical field of logging, and in particular relates to a logging signal demodulation method, a device, computer equipment and a storage medium.
Background
At present, in the use process of petroleum logging instruments, manchester codes are adopted for modulating and demodulating collected signals. In signal demodulation, only a filter is currently used to filter out high-frequency signals, and fundamental waves are used for detecting and demodulating fundamental frequency signals. However, because of the very complex environment downhole in oil, the fundamental frequency signals are often more disturbed than the high frequency signals. In this case, if the low-frequency signal of the signal is still demodulated using the fundamental wave, an accurate symbol cannot be demodulated.
How to determine the appropriate wave for demodulation to improve the signal demodulation accuracy is a problem in the prior art.
Disclosure of Invention
In order to solve the problems in the prior art, the disclosed embodiments provide a method, a device, a computer device and a storage medium for demodulating a logging signal, and the frequency band of the signal to be demodulated, which is demodulated by a determined target harmonic, is minimally interfered by a power spectrum image, so that a symbol obtained by demodulation is more accurate.
In order to solve the technical problems, the specific technical scheme of the disclosure is as follows:
in one aspect, embodiments of the present disclosure provide a method of demodulating a logging signal, comprising,
in response to receiving a signal to be demodulated, determining a first time domain function and a power spectrum image corresponding to the signal to be demodulated;
determining a target harmonic for demodulating a signal according to the power spectrum image;
processing a plurality of preset second time domain functions and the first time domain functions according to the target harmonic wave to obtain a plurality of filter functions and functions to be processed; and
demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements.
Further, determining a target harmonic for demodulation from the power spectral image further comprises:
according to the power spectrum image, determining energy values corresponding to a plurality of frequency segments to obtain a plurality of energy values;
for each frequency segment, determining at least one energy difference between the frequency segment and at least one target frequency segment, respectively; and
and determining the target harmonic according to the frequency segment corresponding to the at least one energy difference value meeting a preset condition.
Further, processing the preset second time domain functions according to the target harmonic wave to obtain a plurality of filtering functions further includes:
performing Fourier transform on each second time domain function to obtain a plurality of first Fourier series expansions;
determining a plurality of first zeroing terms in each first fourier series expansion according to the target harmonic; and
and determining that coefficients of the first zeroing terms in each first Fourier series expansion are zero, and obtaining a plurality of filtering functions.
Further, in the case where the plurality of second time domain functions is two, the plurality of second time domain functions further includes:
wherein the g 1 (t) and said g 2 (T) is a second time domain function, T is time of day, and T is period, respectively.
Further, according to the target harmonic, processing the first time domain function to obtain a function to be processed includes:
performing Fourier transform on the first time domain function to obtain a second Fourier series expansion;
determining a plurality of second zeroing terms in the second fourier series expansion according to the target harmonic;
determining that coefficients of the plurality of second zeroing terms in the second Fourier series expansion are zero to obtain a function to be processed; and
and carrying out inverse Fourier transform on the pre-processing function to obtain the pre-processing function.
Further, demodulating the signal to be demodulated according to the plurality of filter functions and the function to be processed, obtaining a plurality of symbols further includes:
performing convolution operation on the function to be processed and each filtering function to obtain a plurality of objective functions;
determining a plurality of target values for each moment according to the objective function; and
making decisions for a plurality of target values at each time to obtain symbols corresponding to each time, and
and obtaining a plurality of code elements corresponding to the signal to be demodulated according to the code elements corresponding to each moment.
Further, in the case where the plurality of target values is two, the deciding further includes:
d is then n ＝-1,d n+1 ＝-1
D is then n ＝1,d n+1 ＝1
D is then n ＝-1,d n+1 ＝-1
D is then n ＝-1,d n+1 ＝1
Wherein n is a positive integer, T is a period, and y is 1 () And said y 2 () Respectively as objective functions, d n For the symbol of the corresponding time instant, and d n+1 Is the symbol at the next time of the corresponding time.
In another aspect, embodiments of the present disclosure also provide a logging signal demodulation apparatus, comprising,
a first determining unit, configured to determine a first time domain function and a power spectrum image corresponding to a signal to be demodulated in response to receiving the signal to be demodulated;
a second determining unit for determining a target harmonic for demodulating a signal based on the power spectrum image;
the processing module is used for processing a plurality of preset second time domain functions and the first time domain functions according to the target harmonic wave to obtain a plurality of filter functions and functions to be processed; and
and the third determining unit is used for demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements.
In another aspect, an embodiment of the disclosure further provides a computer device, including a memory, a processor, and a computer program stored on the memory and capable of running on the processor, where the processor implements the method described above when executing the computer program.
In another aspect, embodiments of the present disclosure also provide a computer-readable storage medium having stored thereon computer instructions that, when executed by a processor, implement the above-described method.
By utilizing the embodiment of the disclosure, firstly, a target frequency segment with the minimum interference degree is determined according to the power spectrum image of the signal to be demodulated. And constructing a plurality of filter functions based on the target harmonics corresponding to the target frequency band. And demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed corresponding to the signal to be demodulated, thereby improving the accuracy of the code element obtained by demodulation.
Drawings
In order to more clearly illustrate the embodiments of the present disclosure or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present disclosure, and other drawings may be obtained according to these drawings without inventive effort to a person of ordinary skill in the art.
FIG. 1 is a schematic diagram of an embodiment of a method for demodulating a logging signal according to an embodiment of the disclosure;
FIG. 2 is a flow chart illustrating a method of demodulating a log signal according to an embodiment of the present disclosure;
FIG. 3A is a flow chart illustrating a method of demodulating a log signal according to another embodiment of the present disclosure;
FIG. 3B is a flow chart illustrating a method of demodulating a log signal according to another embodiment of the present disclosure;
FIG. 4 is a schematic diagram of an ideal power spectrum according to an embodiment of the present disclosure;
FIG. 5 is a schematic diagram illustrating demodulation of a logging signal according to another embodiment of the disclosure;
FIG. 6 is a schematic diagram of a well logging signal demodulation apparatus according to an embodiment of the disclosure;
fig. 7 is a schematic structural diagram of a computer device according to an embodiment of the disclosure.
[ reference numerals description ]
101. A collecting terminal;
102. a server;
301. a signal to be demodulated;
302. a target harmonic;
303. a first time domain function;
310. a plurality of second time domain functions;
311. a second time domain function;
312. a second time domain function;
320. a plurality of filter functions;
321. a filtering function;
322. a filtering function;
330. a function to be processed;
340. a plurality of objective functions;
341. an objective function;
342. an objective function;
350. a plurality of symbols;
351. a code element;
352. a code element;
501. a signal diagram to be demodulated;
502. a coded function map;
503. a first time domain function diagram;
504. a harmonic component map;
505. a symbol map;
610. a first determination unit;
620. a second determination unit;
630. a processing unit;
640. a third determination unit;
702. a computer device;
704. a processing device;
706. storing the resource;
708. a driving mechanism;
710. an input/output module;
712. an input device;
714. an output device;
716. a presentation device;
718. a graphical user interface;
720. a network interface;
722. a communication link;
724. a communication bus.
Detailed Description
The following description of the technical solutions in the embodiments of the present disclosure will be made clearly and completely with reference to the accompanying drawings in the embodiments of the present disclosure, and it is apparent that the described embodiments are only some embodiments of the present disclosure, not all embodiments. Based on the embodiments in this disclosure, all other embodiments that a person of ordinary skill in the art would obtain without making any inventive effort are within the scope of protection of this disclosure.
It should be noted that the terms "first," "second," and the like in the description and claims of the present disclosure and in the foregoing figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the disclosure described herein may be capable of operation in sequences other than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, apparatus, article, or device that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed or inherent to such process, method, article, or device.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is illustrated in the flowcharts, in some cases the steps illustrated or described may be performed in an order other than that illustrated herein.
Fig. 1 is a schematic diagram of an implementation system of a method for demodulating a logging signal according to an embodiment of the disclosure, which may include: the acquisition terminal 101 and the server 102 communicate with each other through a network, which may include a local area network (Local Area Network, abbreviated as LAN), a wide area network (Wide Area Network, abbreviated as WAN), the internet, or a combination thereof, and are connected to a website, user equipment (e.g., a computing device), and a backend system. The server 102 may demodulate the signal to be demodulated according to the signal to be demodulated obtained by the acquisition terminal 101, determine a plurality of code elements or target information corresponding to the signal to be demodulated, and then send the plurality of code elements or target information to the acquisition terminal 101 through a network, where the acquisition terminal 101 displays the plurality of code elements or target information provided by the server 102, so that a user corresponding to the acquisition terminal 101 can check on the acquisition terminal 101. Alternatively, the servers 102 may be nodes of a cloud computing system (not shown), or each server 102 may be a separate cloud computing system, including multiple computers interconnected by a network and operating as a distributed processing system.
In an alternative embodiment, the acquisition terminal 101 may include downhole instrumentation and uphole electronics. The downhole instrument may be, for example, any sensor that can collect a signal. The above-well electronic device and the oil well instrument communicate through a network to receive signals, and the above-well electronic device is not limited to a smart phone, an acquisition device, a desktop computer, a tablet computer, a notebook computer, a smart speaker, a digital assistant, an augmented Reality (AR, augmented Reality)/Virtual Reality (VR) device, an intelligent wearable device, and the like. Alternatively, the operating system running on the electronic device may include, but is not limited to, an android system, an IOS system, linux, windows, and the like.
In addition, it should be noted that, fig. 1 is only an application environment provided by the present disclosure, and in practical application, a plurality of acquisition terminals 101 may also be included, which is not limited in this specification.
Fig. 2 is a flowchart illustrating a method for demodulating a logging signal according to an embodiment of the present disclosure, and in order to solve the problems in the prior art, the embodiment of the present disclosure provides a method for demodulating a logging signal, which implements determining a suitable wave for demodulation, and improves accuracy of symbols obtained by demodulation. The process of log signal demodulation is described in this figure, but may include more or fewer operational steps based on conventional or non-inventive labor. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When a system or apparatus product in practice is executed, it may be executed sequentially or in parallel according to the method shown in the embodiments or the drawings. As shown in fig. 2, the method may include:
s210, determining a first time domain function and a power spectrum image corresponding to a signal to be demodulated in response to receiving the signal to be demodulated;
s220, determining target harmonic waves for demodulation according to the power spectrum image;
s230, processing a plurality of preset second time domain functions and first time domain functions according to target harmonic waves to obtain a plurality of filter functions and functions to be processed;
s240, demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements.
According to the method, firstly, a target frequency segment with the minimum interference degree is determined according to the power spectrum image of the signal to be demodulated. And constructing a plurality of filter functions based on the target harmonics corresponding to the target frequency band. And demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed corresponding to the signal to be demodulated, thereby improving the accuracy of the code element obtained by demodulation.
According to embodiments of the present disclosure, the signal to be demodulated comprises a signal acquired by an oil downhole instrument. After the signal to be demodulated is received, the signal to be demodulated is coded by utilizing a signal coding model to obtain a coded function, and the coded function is transmitted back to the electronic equipment on the well. The downhole electronics process the received signal to obtain a first time domain function. The signal encoding model may be any model that can encode a signal, for example, a manchester encoding model. And processing the frequency of the signal to be demodulated according to the received signal to be demodulated to determine a power spectrum image.
And integrating the power of different frequency segments to determine the energy value carried by the frequency segment. And determining the interference degree of the signal to be demodulated in the frequency band according to the energy value of each frequency band. The interference degree of each frequency segment is compared, and a target frequency segment for demodulation is determined. And then determining the target harmonic according to the target frequency band. For example, if the signal to be demodulated is determined to have the least interference degree in the frequency segments of 0Hz-40Hz, 40Hz-80Hz, 80Hz-120Hz and 120Hz-160Hz according to the energy, the 40Hz-80Hz is determined to be the target frequency segment. Since the fundamental wave is employed for demodulation in the field of signal demodulation. The frequency domain distribution is divided into a plurality of odd harmonic components, so that the target frequency range of 40Hz-80Hz corresponds to the third harmonic component, and the target harmonic is determined to be the third harmonic. Therefore, the method and the device realize that the target frequency segment with the least interference is determined according to the power spectrum image of the signal to be demodulated, and further determine the target harmonic wave which is most suitable for demodulation.
Before demodulating the signal to be demodulated, a time domain function for generating a filter function is preset. In this disclosure, this time domain function is referred to as a second time domain function. After the target harmonic is determined, processing the plurality of second time domain functions according to the order value of the target harmonic to obtain a plurality of filter functions, wherein the filter functions are used for constructing filters, and the plurality of filter functions are used for constructing a plurality of filters. And processing the first time domain function according to the order value of the target harmonic wave to obtain a function to be processed.
And inputting the function to be processed into a plurality of filters to obtain a plurality of objective functions. The numerical value at each time is respectively brought into the objective function, and a plurality of target values respectively corresponding to each time are obtained. And judging a plurality of target values at each moment to obtain a code element corresponding to each moment. And splicing the determined multiple code elements according to the corresponding time sequence to obtain a code element diagram for subsequent processing.
According to another embodiment of the present disclosure, the plurality of second time domain functions may further include:
wherein g 1 (t) and g 2 (T) is a second time domain function, T is time of day, and T is period, respectively.
According to another embodiment of the present disclosure, for higher accuracy of demodulation, processing the preset plurality of second time domain functions according to the target harmonic, to obtain a plurality of filter functions further includes: performing Fourier transform on each second time domain function to obtain a plurality of first Fourier series expansions; determining a plurality of first zeroing terms in each first fourier series expansion according to the target harmonics; and determining that coefficients of a plurality of first zeroing terms in each first Fourier series expansion are zero, so as to obtain a plurality of filter functions.
The determining of the first zeroing term may be, for example, determining, for each of the plurality of first fourier series expansions, that a term other than a term corresponding to the order value of the target harmonic is the first zeroing term.
For example, in the case where the subvalue of the target harmonic is determined to be 3 (i.e., 3 subharmonics), the term other than the third term (n=3) in the first fourier series expansion is determined to be the first zeroing term. The description will be continued taking the second time domain function in the equation (1) and the equation (2) as an example. Fourier transforms are performed on the formula (1) and the formula (2) respectively, so as to obtain a first fourier series expansion (formula (3) -formula (6) as an example of a derivation process) shown in the formula (6), and a first fourier series expansion (formula (7) -formula (10) as an example of a derivation process) shown in the formula (10).
Wherein,
thus, when t.epsilon.0, T,
wherein,
thus, when t.epsilon.0, T,
after determining a plurality of first fourier series expansions such as equation (6) and equation (10),n corresponding entries except n=3 are determined to be the first zeroed entry. Further, the coefficient of the first zeroing term is determined to be 0, and a plurality of filter functions are obtained as shown in the following formula (11)And +.>
It should be noted that, in the process of using fundamental demodulation, a second time domain function is often directly used to construct a filter, and the second time domain function is not processed and converted. The present disclosure determines corresponding higher harmonic components for each of a plurality of preset second time domain functions, and constructs a filter using the higher harmonic components. Therefore, the accuracy in the signal demodulation process is further improved.
According to another embodiment of the present disclosure, processing the first time domain function according to the target harmonic, obtaining the function to be processed further includes: performing Fourier transform on the first time domain function to obtain a second Fourier series expansion; determining a plurality of second zeroing terms in a second fourier series expansion according to the target harmonic; determining that coefficients of a plurality of second zeroing terms in the second Fourier series expansion are zero to obtain a function to be processed; and performing inverse Fourier transform on the pre-processing function to obtain the pre-processing function.
For example, the first time domain function is s (t), and the derivation process of the above formula (3) -formula (6) is performed for s (t) to perform fourier transformation, so as to obtain a second fourier series expansion corresponding to s (t). In the case where the subvalue of the target harmonic is determined to be 3, the term other than the third term (n=3) in the second fourier series expansion is determined to be the second zeroing term. The coefficients of the plurality of second zeroing terms are set to zero (as in the process of converting the above formula (6) into the formula (11) or the formula (10) into the formula (12)), resulting in the pre-processing function m (t). And (3) carrying out inverse Fourier transform on m (t) to obtain a function f (t) to be processed corresponding to m (t). The fourier transform and the inverse fourier transform are reciprocal transforms.
According to another embodiment of the present disclosure, demodulating the signal to be demodulated according to a plurality of filter functions and a function to be processed, obtaining a plurality of symbols further includes: performing convolution operation on the function to be processed and each filtering function to obtain a plurality of objective functions; determining a plurality of target values for each moment according to the objective function; and judging a plurality of target values at each moment to obtain a code element corresponding to each moment, and obtaining a plurality of code elements corresponding to the signal to be demodulated according to the code element corresponding to each moment.
For example, in the case where the order value of the target harmonic is determined to be 3 and the plurality of second time domain functions is determined to be two, the above-determined function f (t) to be processed and the plurality of filter functions areAnd->Respectively performing convolution operation to obtain multiple objective functions, such as y obtained by the following formulas (13) and (14) 1 (t) and y 1 (t)。
Wherein x is convolution operation, y 1 (t) and y 2 (t) are objective functions, respectively.
Will be associated with each momentThe corresponding n values are respectively substituted into the objective function y 1 (t) and y 1 In (t), two target values corresponding to each time (each n value) are obtained. A decision is made for the two target values, determining the symbol corresponding to each instant. Thereby determining a symbol map for subsequent processing.
According to another embodiment of the present disclosure, in the case where the plurality of target values is two, the decision includes the following formula (15), formula (16), formula (17), and formula (18).
Wherein n is a positive integer, T is a period, y 1 () And y 2 () Respectively as objective functions d n For symbols at the corresponding time instant, d n+1 Is the symbol at the next time of the corresponding time.
For each time instant (each n), two target values can be obtained, and the larger first target value of the two target values is determined. It is determined whether the first target value is greater than zero. Further, decisions as shown in the formulas (15), (16), (17) and (18) are made to determine the symbol (d) corresponding to the time (n) n ) And a symbol (d) corresponding to a time n+1 next to the time (d) n+1 )。
It should be noted that the decision point selected in the decision of the present disclosure is the middle point ((3/2) T) of one cycle to adapt to the filtering function in the present disclosure, instead of the boundary point (T) of one cycle.
According to another embodiment of the present disclosure, in order to improve the accuracy of the demodulated symbol, in the case where the plurality of target values is two, only the decision of the n odd or even terms is made at the time of the decision. That is, n is odd or n is even.
In the case where only n odd or even term decisions are made, only one corresponding symbol may be generated for each time instant. If the decision is made for all n, two corresponding symbols are generated for one time, thus making the generated symbols inaccurate.
Fig. 3A is a flowchart illustrating a method for demodulating a logging signal according to another embodiment of the present disclosure.
According to another embodiment of the present disclosure, as shown in fig. 3A, after receiving the signal 301 to be demodulated, the signal 301 to be demodulated is encoded by using a manchester encoding model, an encoded function is obtained, and the encoded function is transmitted back. And processing the received signal corresponding to the interfered coded function to obtain a first time domain function 303. For this signal to be demodulated 301, the energy value for each frequency segment is determined. And according to the energy value of each frequency segment, obtaining the target harmonic wave 302 after interference analysis.
And performing Fourier transform on the first time domain function 303 to obtain a second Fourier series expansion, determining that the term inconsistent with the term corresponding to the frequency of the target harmonic 302 is a second zeroing term, and setting the coefficients of the second zeroing term to be zero to obtain a function 330 to be processed. Similarly, the same operation is performed for each of the plurality of second time domain functions 310, resulting in a plurality of filter functions 320.
Convolution operation is performed on the to-be-processed function 330 and the plurality of filter functions 320 to obtain a plurality of objective functions 340. A target value corresponding to each time is determined by a plurality of objective functions 340. And then makes decisions for each target value to determine a plurality of symbols 350.
According to another embodiment of the present disclosure, in order to more accurately determine the target harmonic for demodulation, determining the target harmonic for demodulation from the power spectrum image further includes: according to the power spectrum image, determining energy values corresponding to a plurality of frequency segments to obtain a plurality of energy values; determining, for each frequency bin, at least one energy difference between the frequency bin and at least one target frequency bin, respectively; and determining the target harmonic according to the frequency segment corresponding to at least one energy difference value meeting the preset condition.
In the absence of interference, the ideal power spectrum image is shown in fig. 4. For this ideal power spectrum image, the energy of each frequency bin is determined, resulting in an energy distribution table for the different frequency bins as shown in table 1 below.
TABLE 1
As can be seen from table 1, the preset conditions were determined by setting the 3 rd order harmonic to be 10.9dBm less than the fundamental wave energy, the 5 th order harmonic to be 4.52dBm less than the 3 rd order harmonic energy, the 5 th order harmonic to be 15.42dBm less than the fundamental wave energy, the 7 th order harmonic to be 2.94dBm less than the 5 th order harmonic energy, the 7 th order harmonic to be 7.46dBm less than the 3 rd order harmonic, and the 7 th order harmonic to be 19.36dBm less than the fundamental wave energy. Because the preset conditions are determined based on the ideal power spectrum image, the rationality of the determined preset conditions is ensured, so that the target harmonic wave of the signal to be demodulated can be determined from a plurality of harmonic waves, and the accuracy of the code element obtained by demodulation is improved.
For each frequency bin, all frequency bins smaller than the frequency bin are determined as target frequency bins. For example, for a frequency band of 40Hz-80Hz (3 rd harmonic), the target frequency band is 0Hz-40Hz (fundamental); aiming at the frequency band of 80Hz-120Hz (5 th harmonic), the target frequency band of 0Hz-40Hz (fundamental wave) and 40Hz-80Hz (3 rd harmonic); for frequency bands of 120Hz-160Hz (7 th harmonic), the target frequency bands are 0Hz-40Hz (fundamental wave), 40Hz-80Hz (3 rd harmonic) and 80Hz-120Hz (5 th harmonic).
According to the frequency segment corresponding to the at least one energy difference value satisfying the preset condition, determining the target harmonic may include, for example: and determining a corresponding threshold according to the frequency segment and the target frequency segment, and determining that the frequency segment meets a preset condition and that the harmonic corresponding to the frequency segment is a target harmonic under the condition that the energy difference value of the frequency segment and each target frequency segment is larger than the corresponding threshold.
In the case of a frequency band of 40Hz-80Hz (3 rd harmonic) and a target frequency band of 0Hz-40Hz (fundamental wave), the corresponding threshold value is 10.9dBm.
Under the condition that the frequency band is 80Hz-120Hz (5 th harmonic) and the target frequency band is 40Hz-80Hz (3 rd harmonic), the corresponding threshold value is 4.52dBm; in the case of frequency bands 80Hz-120Hz (5 th harmonic) and target frequency bands 0Hz-40Hz (fundamental wave), the corresponding threshold is 15.42dBm.
Under the condition that the frequency band is 120Hz-160Hz (7 th harmonic) and the target frequency band is 80Hz-120Hz (5 th harmonic), the corresponding threshold value is 2.94dBm; under the condition that the frequency band is 120Hz-160Hz (7 th harmonic) and the target frequency band is 40Hz-80Hz (3 rd harmonic), the corresponding threshold value is 7.46dBm; in the case of frequency bands 120Hz-160Hz (7 th harmonic) and target frequency bands 0Hz-40Hz (fundamental wave), the corresponding threshold is 19.36dBm.
For example, for frequency bands 40Hz-80Hz (3 rd harmonics), the target frequency band is 0Hz-40Hz (fundamental). The energy value of the frequency band 40Hz-80Hz is determined to be 19dBm, the energy value of the target frequency band 0Hz-40Hz is determined to be 32dBm, and the energy difference is 13dBm. The threshold corresponding to the frequency band 40Hz-80Hz (3 rd harmonic) and the target frequency band 0Hz-40Hz (fundamental wave) is 10.9dBm. Comparing the energy difference value of 13dBm with a corresponding threshold value of 10.9dBm, and determining that the energy difference value is larger than the corresponding threshold value, and determining that the frequency range of 40Hz-80Hz (3 rd order harmonic) meets a preset condition, and determining that the 3 rd order harmonic is a target harmonic. It should be noted that the energy difference value must be a number equal to or greater than zero.
To illustrate the presence of multiple target frequency bins, for example, the target frequency bins are 0Hz-40Hz (fundamental) and 40Hz-80Hz (3 rd harmonic) for the frequency bins 80Hz-120Hz (5 th harmonic). Similarly, the energy value of the frequency band 80Hz-120Hz (5 th harmonic) is determined to be 14, the energy value of the target frequency band 40Hz-80Hz is determined to be 19dBm, the energy value of the target frequency band 0Hz-40Hz is determined to be 32dBm, the energy difference between the frequency band 80Hz-120Hz and the target frequency band 40Hz-80Hz is determined to be 18dBm, and the energy difference between the frequency band 80Hz-120Hz and the target frequency band 40Hz-80Hz is determined to be 5dBm. The threshold value of the frequency band 80Hz-120Hz (5 th harmonic) and the target frequency band 0Hz-40Hz (fundamental wave) is 15.42dBm, and the threshold value of the frequency band 80Hz-120Hz (5 th harmonic) and the target frequency band 40Hz-80Hz (fundamental wave) is 4.52dBm. I.e. 18dBm compared to the threshold of 15.42dBm and 5dBm compared to 4.52dBm. Based on the comparison, the energy difference values of the frequency band 80Hz-120Hz (5 th harmonic) and the target frequency band are larger than the corresponding threshold values, so that the frequency band 80Hz-120Hz (5 th harmonic) is determined to meet the preset condition, and the 5 th harmonic is determined to be the target harmonic. It should be noted that the target harmonic is not fundamental (1 st harmonic).
FIG. 3B is a flow chart illustrating a method of demodulating a log signal according to another embodiment of the present disclosure
According to another embodiment of the present disclosure, in the case that the plurality of second time domain functions is two, a flowchart of the logging signal demodulation method of the present embodiment may be as shown in fig. 3B. After the target harmonic 302 is determined, fourier transformation is performed on the second time domain function 311 to obtain a first fourier series expansion, a term inconsistent with a term corresponding to the number of times of the target harmonic 302 is determined to be a first zeroing term, and coefficients of the first zeroing terms are set to be zero to obtain a filter function 321. Similarly, the same operation is performed for the second time domain function 312, resulting in the filter function 322.
The objective function 341 is obtained by performing convolution operation on the function 330 to be processed and the filter function 321. The objective function 342 is obtained by convolving the function 330 to be processed with the filter function 322. Two target values corresponding to each time are determined by the objective function 341 and the objective function 342. And then, judging each target value to obtain a corresponding code element. For example, symbol 351 may be obtained for one time and symbol 352 may be obtained for another time.
Note that the number of symbols in this embodiment may be any number, and the two symbols in the above example are merely illustrative and not limiting in practice.
Fig. 5 is a schematic diagram illustrating demodulation of a logging signal according to another embodiment of the disclosure.
According to another embodiment of the disclosure, the original signal collected by the downhole instrument is shown in a signal map 501 to be demodulated, and is encoded by using a manchester encoding model to obtain an encoded function map 502, and due to interference of the underground environment, the electronic device on the well obtains an interfered signal shown in a first time domain function map 503. After receiving the first time domain function map 503 transmitted back by the downhole tool, determining a target harmonic according to the disclosure for the first time domain function map 503, determining a plurality of target functions according to a plurality of second time domain functions, the first time domain functions and the target harmonic, and determining an image obtained by corresponding target values for each moment, as shown in a harmonic component map 504. Further, the target value is decided to obtain a plurality of symbol patterns 505.
The embodiment of the disclosure also provides a logging signal demodulation apparatus, as shown in fig. 6, comprising,
a first determining unit 610, configured to determine a first time domain function and a power spectrum image corresponding to a signal to be demodulated in response to receiving the signal to be demodulated;
a second determining unit 620 for determining a target harmonic for demodulation from the power spectrum image;
the processing unit 630 is configured to process a plurality of second time domain functions and the first time domain functions according to the target harmonic, so as to obtain a plurality of filter functions and functions to be processed; and
and a third determining unit 640, configured to demodulate the signal to be demodulated according to the multiple filtering functions and the function to be processed, so as to obtain multiple symbols.
Since the principle of the device for solving the problem is similar to that of the method, the implementation of the device can be referred to the implementation of the method, and the repetition is omitted.
Fig. 7 is a schematic structural diagram of a computer device according to an embodiment of the disclosure, where an apparatus in the disclosure may be the computer device in the embodiment, and perform the method of the disclosure. The computer device 702 may include one or more processing devices 704, such as one or more Central Processing Units (CPUs), each of which may implement one or more hardware threads. The computer device 702 may also include any storage resources 706 for storing any kind of information, such as code, settings, data, etc. For example, and without limitation, storage resources 706 may include any one or more of the following combinations: any type of RAM, any type of ROM, flash memory devices, hard disks, optical disks, etc. More generally, any storage resource may store information using any technology. Further, any storage resource may provide volatile or non-volatile retention of information. Further, any storage resources may represent fixed or removable components of computer device 702. In one case, the computer device 702 can perform any of the operations of the associated instructions when the processing device 704 executes the associated instructions stored in any storage resource or combination of storage resources. The computer device 702 also includes one or more drive mechanisms 708, such as a hard disk drive mechanism, an optical disk drive mechanism, and the like, for interacting with any storage resources.
The computer device 702 may also include an input/output module 710 (I/O) for receiving various inputs (via an input device 712) and for providing various outputs (via an output device 714). One particular output mechanism may include a presentation device 716 and an associated Graphical User Interface (GUI) 718. In other embodiments, input/output module 710 (I/O), input device 712, and output device 714 may not be included as just one computer device in a network. The computer device 702 can also include one or more network interfaces 720 for exchanging data with other devices via one or more communication links 722. One or more communication buses 724 couple the above-described components together.
Communication link 722 may be implemented in any manner, for example, through a local area network, a wide area network (e.g., the internet), a point-to-point connection, etc., or any combination thereof. Communication link 722 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.
The disclosed embodiments also provide a computer readable storage medium storing a computer program which, when executed by a processor, implements the above method.
The disclosed embodiments also provide a computer program product comprising a computer program which, when executed by a processor, implements the above method.
It will be apparent to those skilled in the art that embodiments of the present disclosure may be provided as a method, system, or computer program product. Accordingly, the present disclosure may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present disclosure 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, etc.) having computer-usable program code embodied therein.
The present disclosure is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the disclosure. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing embodiments have been provided for the purpose of illustrating the general principles of the present disclosure, and are not meant to limit the scope of the disclosure, but to limit the scope of the disclosure.
Claims (10)
1. A method for demodulating a logging signal, comprising:
in response to receiving a signal to be demodulated, determining a first time domain function and a power spectrum image corresponding to the signal to be demodulated;
determining a target harmonic for demodulation according to the power spectrum image;
processing a plurality of preset second time domain functions and the first time domain functions according to the target harmonic wave to obtain a plurality of filter functions and functions to be processed; and
demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements.
2. The method of claim 1, wherein said determining a target harmonic for demodulation from said power spectral image comprises:
according to the power spectrum image, determining energy values corresponding to a plurality of frequency segments to obtain a plurality of energy values;
for each frequency segment, determining at least one energy difference between the frequency segment and at least one target frequency segment, respectively; and
and determining the target harmonic according to the frequency segment corresponding to the at least one energy difference value meeting a preset condition.
3. The method of claim 1, wherein processing a preset plurality of second time domain functions according to the target harmonic to obtain a plurality of filter functions includes:
performing Fourier transform on each second time domain function to obtain a plurality of first Fourier series expansions;
determining a plurality of first zeroing terms in each first fourier series expansion according to the target harmonic; and
and determining that coefficients of the first zeroing terms in each first Fourier series expansion are zero, and obtaining a plurality of filtering functions.
4. A method according to claim 1 or 3, wherein in case the plurality of second time domain functions is two, the plurality of second time domain functions comprises:
wherein the g 1 (t) and said g 2 (T) is a second time domain function, T is time of day, and T is period, respectively.
5. The method of claim 1, wherein processing the first time domain function according to the target harmonic to obtain a function to be processed comprises:
performing Fourier transform on the first time domain function to obtain a second Fourier series expansion;
determining a plurality of second zeroing terms in the second fourier series expansion according to the target harmonic;
determining that coefficients of the plurality of second zeroing terms in the second Fourier series expansion are zero to obtain a function to be processed; and
and carrying out inverse Fourier transform on the pre-processing function to obtain the pre-processing function.
6. The method of claim 1, wherein demodulating the signal to be demodulated according to the plurality of filter functions and the function to be processed to obtain a plurality of symbols comprises:
performing convolution operation on the function to be processed and each filtering function to obtain a plurality of objective functions;
determining a plurality of target values for each moment according to the objective function; and
making decisions for a plurality of target values at each time to obtain symbols corresponding to each time, and
and obtaining a plurality of code elements corresponding to the signal to be demodulated according to the code elements corresponding to each moment.
7. The method according to claim 6, wherein in case the plurality of target values is two, the deciding comprises:
d is then n ＝-1，d n+1 ＝-1
D is then n ＝1，d n+1 ＝1
D is then n ＝-1，d n+1 ＝-1
D is then n ＝-1，d n+1 ＝1
Wherein n is a positive integer, T is a period, and y is 1 () And said y 2 () Respectively as objective functions, d n For the symbol of the corresponding time instant, and d n+1 Is the symbol at the next time of the corresponding time.
8. A logging signal demodulation device is characterized by comprising,
a first determining unit, configured to determine a first time domain function and a power spectrum image corresponding to a signal to be demodulated in response to receiving the signal to be demodulated;
a second determining unit for determining a target harmonic for demodulation from the power spectrum image;
the processing unit is used for processing a plurality of preset second time domain functions and the first time domain functions according to the target harmonic wave to obtain a plurality of filter functions and functions to be processed; and
and the third determining unit is used for demodulating the signal to be demodulated according to the plurality of filtering functions and the function to be processed to obtain a plurality of code elements.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of the preceding claims 1-7 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium has stored thereon a computer program which, when executed by a processor, performs the method of any of the preceding claims 1-7.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211023979.5A CN117675091A (en) | 2022-08-24 | 2022-08-24 | Logging signal demodulation method and device, computer equipment and storage medium |
| PCT/CN2023/114755 WO2024041615A1 (en) | 2022-08-24 | 2023-08-24 | Harmonic decoding method and apparatus for underground electromagnetic-wave signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211023979.5A CN117675091A (en) | 2022-08-24 | 2022-08-24 | Logging signal demodulation method and device, computer equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117675091A true CN117675091A (en) | 2024-03-08 |
Family
ID=90012606
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211023979.5A Pending CN117675091A (en) | 2022-08-24 | 2022-08-24 | Logging signal demodulation method and device, computer equipment and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN117675091A (en) |
| WO (1) | WO2024041615A1 (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103338176B (en) * | 2013-06-13 | 2016-08-24 | 上海理工大学 | A kind of implementation method of waveform reuse modulator-demodulator |
| US9713099B2 (en) * | 2014-12-16 | 2017-07-18 | Netgear, Inc. | Systems and methods for cable and WLAN coexistence |
| CN112235222A (en) * | 2020-01-21 | 2021-01-15 | 南京新频点电子科技有限公司 | OFDM interference signal modulation device and modulation method |
| CN114609593B (en) * | 2022-03-09 | 2024-04-12 | 北京航空航天大学 | Harmonic radar based on FPGA and deep learning |
-
2022
- 2022-08-24 CN CN202211023979.5A patent/CN117675091A/en active Pending
-
2023
- 2023-08-24 WO PCT/CN2023/114755 patent/WO2024041615A1/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| WO2024041615A1 (en) | 2024-02-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110839016B (en) | Abnormal flow monitoring method, device, equipment and storage medium | |
| Olhede et al. | A generalized demodulation approach to time-frequency projections for multicomponent signals | |
| Frei et al. | Intrinsic time-scale decomposition: time–frequency–energy analysis and real-time filtering of non-stationary signals | |
| Borghesani | The envelope-based cyclic periodogram | |
| EP3621259A1 (en) | Method and device for fsk/gfsk demodulation | |
| CN110690931B (en) | Digital signal adaptive code rate estimation method and device based on multi-wavelet-base combination | |
| CN110084208B (en) | Calculation order tracking method capable of adaptively reducing noise and avoiding order aliasing | |
| CN103154932A (en) | Method for analyzing signals providing instantaneous frequencies and sliding fourier transforms, and device for analyzing signals | |
| CN111193627B (en) | Information processing method, device, equipment and storage medium | |
| CN115630277A (en) | Data signal processing method | |
| Daldal | A novel demodulation method for quadrate type modulations using a hybrid signal processing method | |
| CN117675091A (en) | Logging signal demodulation method and device, computer equipment and storage medium | |
| CN109783770A (en) | Frequency analysis method and device | |
| CN111611686A (en) | Detection method for communication signal time-frequency domain | |
| CN111125976A (en) | Automatic generation method of RTL model | |
| Chen et al. | Improvement of empirical mode decomposition based on correlation analysis | |
| CN101917239B (en) | Method and system for acquiring cyclic spectrum alpha section based on frequency domain smoothing | |
| CN115347978A (en) | Method, device and storage medium for identifying AIS frame header data | |
| Siohan et al. | Derivation of extended Gaussian functions based on the Zak transform | |
| WO2018142620A1 (en) | Packet format deduction device and packet format deduction program | |
| US20130124588A1 (en) | Encoding densely packed decimals | |
| CN113765838A (en) | DPSK signal demodulation method, device, equipment and storage medium | |
| CN108353045A (en) | A kind of signal processing system, method and device | |
| Wenbin et al. | Design and verification of MSK based on FPGA | |
| CN110782238A (en) | Abnormal resource transfer combination identification method, device, equipment and medium |
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 |