CN115236748A - Frequency domain combined wave normalization expression method and device and storage medium - Google Patents
Frequency domain combined wave normalization expression method and device and storage medium Download PDFInfo
- Publication number
- CN115236748A CN115236748A CN202210777269.5A CN202210777269A CN115236748A CN 115236748 A CN115236748 A CN 115236748A CN 202210777269 A CN202210777269 A CN 202210777269A CN 115236748 A CN115236748 A CN 115236748A
- Authority
- CN
- China
- Prior art keywords
- frequency
- oct
- determining
- combined wave
- combined
- 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
- 238000010606 normalization Methods 0.000 title claims abstract description 37
- 238000010183 spectrum analysis Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
A frequency domain combined wave normalization expression method, a device and a storage medium. The invention discloses a frequency domain combined wave normalization expression method which comprises the steps of obtaining the depth range of an area to be explored; determining a fundamental frequency according to the depth range of the region to be explored; determining frequency density according to the target scale of the region to be explored; determining a combined signal according to the fundamental frequency and the frequency density; expressing the formula F = { a1 × F from the combined signal base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230, determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal; realizing the description of all combined signals through parameter matrixes of parameters an, fbase, oct and bn of the combined wave; and sending the described combined wave to the area to be explored so as to solve the technical problems of inconvenient signal application and poor expansibility in a fixed naming mode in the prior art.
Description
Technical Field
The invention relates to the technical field of exploration, in particular to a frequency domain combined wave normalization expression method, a frequency domain combined wave normalization expression device and a storage medium.
Background
At present, square wave emission is easy to realize on hardware, and the same peak-to-peak value has stronger energy than that of a sine wave, so that most frequency domain exploration equipment adopts the square wave as an excitation signal.
Conventional frequency domain surveys use frequency conversion, i.e., a frequency-swept method of transmitting at a frequency, where the waveform description is characterized by a popular waveform frequency, such as 1000Hz.
At present, a pseudo-random 7-frequency wave (a combined wave with a relatively low frequency) is mainly used as a signal for transmitting and receiving, the signal is divided into 12 frequency groups, and because 7 main frequency signals (so called 7-frequency waves) exist in the signal of each frequency group, a single frequency cannot well describe the signal at this time, so that the 12 frequency groups are named as 0-11 frequency groups respectively by adopting a naming mode, for example, the 7 frequency groups have main frequencies of 1, 2, 4, 8, 16, 32 and 64Hz. The nomenclature at this time has no relation to the waveform characteristics.
With the requirements of further research and exploration on pseudo-random signals, the development and application of high-order pseudo-random signals (a combined wave with a large frequency) are promoted, at the moment, a frequency description mode only suitable for single-frequency signals cannot meet the requirements, a corresponding table needs to be established in a fixed naming mode, and the pseudo-random signals are inconvenient to apply and poor in expansibility.
Disclosure of Invention
The invention provides a frequency domain combined wave normalization expression method, which aims to solve the technical problems that a fixed naming mode is inconvenient to apply to combined waves and has poor expansibility in the prior art.
In order to achieve the above object, the present invention provides a frequency domain combined wave normalization expression method, where the frequency domain combined wave normalization expression method includes:
acquiring the depth range of a region to be explored;
determining a fundamental frequency according to the depth range of the region to be explored;
determining frequency density according to the target scale of the region to be explored;
determining a combined signal according to the fundamental frequency and the frequency density;
according to the combined wave formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal;
realizing the description of all combined signals through parameter matrixes of parameters an, fbase, oct and bn of the combined wave;
and sending the described combined wave to the area to be explored.
Optionally, the method further comprises, according to the formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining parameters an, fbase, oct and bn of the combined wave comprises:
performing spectral analysis on the combined signal to determine useful frequencies and fundamental frequencies in the signal;
determining order coefficients, octave coefficients and octaves corresponding to the number of the useful frequencies according to the useful frequencies and the fundamental frequency;
and determining a combined wave according to the order coefficient, the octave coefficient, the fundamental frequency and the octave.
Optionally, the step of performing a spectral analysis on the combined wave signal to determine a useful frequency and a fundamental frequency in the combined wave signal comprises:
performing spectrum analysis on the combined wave signal to obtain a spectrogram;
when the frequency amplitude in the spectrogram is larger than a preset threshold value, determining that the current frequency is a useful frequency;
the lowest of all the useful frequencies is the fundamental frequency.
Optionally, the step of determining order coefficients, octave coefficients and octaves corresponding to the number of useful frequencies according to the useful frequencies and the fundamental frequency includes:
carrying out ratio processing on each useful frequency and the fundamental frequency to obtain a plurality of frequency coefficients;
grouping a plurality of the frequency coefficients to obtain order coefficients;
determining a common ratio according to each frequency coefficient;
defining the common ratio as an octave;
determining the number of the frequency coefficients of each group after grouping;
and determining octave coefficients according to the number of the frequency coefficients of each group.
Optionally, the step of grouping a plurality of the frequency coefficients to obtain order coefficients includes:
performing prime factor decomposition on a plurality of the frequency coefficients to obtain at least one prime factor;
and determining an order coefficient according to the quality factor.
Optionally, the step of determining the octave coefficients according to the number of each group includes:
and performing difference operation on the number of the frequency coefficients of each group to obtain the octave coefficients.
In order to achieve the above object, the present invention further provides a frequency domain combined wave normalization expression device, where the frequency domain combined wave normalization expression method includes:
the parameter acquisition module is used for acquiring the depth range of the area to be explored;
the calculation module is used for determining a fundamental frequency according to the depth range of the region to be explored; determining frequency density according to the target scale of the region to be explored; determining a combined signal according to the fundamental frequency and the frequency density;
a control module that, according to a combined wave formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining order coefficients an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal to realize the description of all the combined signals through a parameter matrix of parameters an, fbase, oct and bn of the combined wave;
and the sending module is used for sending the described combined wave to the area to be explored.
In order to achieve the above object, the present invention further proposes a storage medium storing a computer program, which, when executed by a processor, causes the processor to execute the frequency-domain combined wave normalization expression method as described above.
In order to achieve the above object, the present invention further provides a frequency domain combined wave normalization expression device, which includes a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor is caused to execute the steps of the frequency domain combined wave normalization expression method.
The depth range of an area to be explored is obtained; determining a fundamental frequency according to the depth range of the region to be explored; determining frequency density according to the target scale of the region to be explored; determining a combined signal according to the fundamental frequency and the frequency density; according to the formula of combined wave F = { a 1. F base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230, determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal; and realizing the description of all combined signals through the parameter matrixes of the parameters an, fbase, oct and bn of the combined wave. By reestablishing a set of frequency combination description mode, the expression form description of the combination wave can be conveniently adjusted according to the detection requirement without adopting a fixed naming mode for description, and the technical problems that the fixed naming mode combined signal in the prior art is inconvenient to apply and poor in expansibility are solved.
Drawings
The invention is further described below with reference to the accompanying drawings and examples;
fig. 1 is a schematic flow chart of a frequency domain combined wave normalization expression method in an embodiment.
FIG. 2 is a before-after comparison graph of spectral analysis converting a time-amplitude time domain graph to a frequency-amplitude spectrogram in one embodiment.
Fig. 3 is a schematic flowchart of a normalized expression method of a frequency domain combined wave in an embodiment.
Fig. 4 is a schematic block diagram of a frequency domain combined wave normalization expression method in an embodiment.
Detailed Description
Reference will now be made in detail to the present preferred embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to like elements throughout.
The method aims to solve the technical problems that a fixed naming mode combined signal is inconvenient to apply and poor in expansibility in the prior art. The invention provides a frequency domain combined wave normalization expression method.
It should be noted that the frequency domain combined wave normalization expression method provided by the invention is mainly applied to a detailed exploration stage, that is, a geological exploration drawing with depth and range distribution can be generated according to exploration data at the moment.
In an embodiment, as shown in fig. 1, the method for normalizing and expressing a frequency domain combined wave includes:
s1, acquiring a depth range of a region to be explored;
when the exploration is carried out, the depth range needing to be explored can be input by a user, wherein the depth range can be the depth range needing to be further explored and obtained after the initial detection, and can also be the depth range of a region to be explored, which needs to be explored and is directly determined.
S2, determining a fundamental frequency according to the depth range of the region to be explored;
wherein, the exploration depth of the exploration signals with different frequencies is different, therefore, the fundamental frequency needs to be determined according to the depth. It should be noted that the fundamental frequency at this time can also be determined in other manners. And not exclusively. Or may be entered directly by the user.
S3, determining frequency density according to the target scale of the area to be explored;
the target scale of the area to be explored actually refers to the longitudinal extension of a target body of the area to be explored, generally, the smaller the scale of the target body of the area to be explored is, the higher the requirement on the longitudinal resolution capability is, the denser the required frequency density is, and the frequency density at this time refers to the frequency interval between each main frequency of the high-order pseudorandom signal for measurement.
S4, determining a combined signal according to the fundamental frequency and the frequency density;
in this case, the combined signal may be generated according to the above parameters by a wide-area electromagnetic method or an excitation method. The combined signal may be a 7-frequency wave, which needs to be named for each frequency group when in use, and the main frequency may include multiple frequencies, such as 1, 2, 4, 8, 16, 32, and 64Hz, so that no unification is achieved, and the more frequency groups, the more the description is disordered. The combined wave signal may also be encoded pseudo-randomly at a higher order, such as a3 rd order 39-frequency wave, which is also complicated by the random encoding.
S5, according to a combined wave formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined wave;
s6, describing all combined signals through parameter matrixes of the parameters an, fbase, oct and bn of the combined wave;
because the existing frequency group naming modes (such as F7-0 to F7-11 respectively represent 12 frequency groups of 7 frequency waves from low frequency to high frequency) cannot meet the definition of working frequency combination, and meanwhile, the construction flexibility is not provided. Thus. In the embodiment of the present invention, a unique combined wave can be determined according to the fundamental frequency and the frequency density, and the combined wave can express that a high-order pseudo-random signal with multiple frequencies exists, which is not described by the conventional 7-frequency wave name. Thus, different waveforms can be characterized by the formula, and the definition is simple and clear and is convenient for the control of the transmitting instrument.
And S7, sending the described combined wave to the area to be explored.
In the above embodiment, when the combined signal is used for exploration, the expression form description of the combined wave can be conveniently adjusted according to the detection requirement by reestablishing a set of frequency combination description mode of the combined signal, without adopting a fixed naming mode for description, thereby solving the technical problems of inconvenient application and poor expansibility of the fixed naming mode combined signal in the prior art.
The working principle of the present scheme is explained by an embodiment, the conventional frequency group naming method is that F7-0 to F7-11 respectively represent 12 frequency groups of 7 frequency waves from low frequency to high frequency, and the sequential pseudo-random driving control signal is composed of the conventional single-frequency, 3 frequency, 5 frequency, 7 frequency and high-order 29, 39, 51, 81 combined waves. For example, the lowest frequency of 0.1Hz for wide-area electromagnetic exploration in an exploration area can meet the requirement on the exploration depth, and according to the working mode of high, medium and low frequency groups, the lowest frequency of the medium frequency group is 0.65Hz and cannot meet the requirement at present, but the lowest frequency of the low frequency group is 0.0097Hz, and in order to meet the requirement on the lowest frequency, the low frequency group has to be selected, so that the observation time is long, the superposition times in the same observation time are few, and the noise is difficult to suppress better. After the combination wave is adopted, the required lowest frequency is determined according to the exploration depth requirement, and the fundamental frequency of a 7-frequency wave lowest frequency band group or a high-order pseudo-random wave in the naming description is set to be the required lowest frequency, so that the transmission and the reception of signals in the required frequency band can be flexibly realized, the construction efficiency is improved, or under the same construction time length, the frequency density is determined as required, and the description mode can realize more overlapping times and improve the exploration data quality. In practical application, after the exploration depth and the frequency density are determined, if the combined wave determined according to the wide-area electromagnetic method or the induced polarization method is determined as the 7-frequency wave lowest frequency band group or the high-order pseudo-random wave, the combined wave can be directly expressed through a combined wave formula to realize the application on the computer level, so that the unified expression of the traditional pseudo-random wave, the single-frequency combined wave and the high-order combined wave is realized.
Optionally, the step of determining a fundamental frequency from the depth range of the area to be surveyed comprises:
acquiring background resistivity and exploration depth of an exploration area;
determining a fundamental frequency according to an empirical formula, the background resistivity, and the depth of investigation.
It should be noted that the lowest frequency is generally determined according to the empirical formula, the background resistivity and the exploration depth, and the fundamental frequency is set to be smaller than the lowest frequency when the fundamental frequency is actually determined, so as to further ensure that the exploration signal at this time can certainly detect the depth range of the area to be explored.
Optionally, the formula for determining the fundamental frequency from the empirical formula, the background resistivity and the depth of investigation is as follows:
where H is the exploration depth, ρ is the background resistivity of the exploration area, and f is the exploration frequency.
An empirical formula for frequency domain electromagnetic sounding survey depth is shown above, where H is the survey depth, ρ is the background resistivity of the survey area, and f is the survey frequency. Background resistivity rho can be known by collecting past geophysical prospecting, geological or well logging data of an exploration area, the required lowest frequency f is calculated according to the required exploration depth, and in order to ensure the exploration depth, fbase which is slightly lower than the lowest frequency f is set in the process of finally applying to construction operation.
Alternatively, the formula F = { a1 × F base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230one initial expression form is determined according to the following method, in practical practice, for example, the dominant frequency of one intermediate frequency group of the conventional 7-frequency waves in wide-area electromagnetic prospecting is {12 48 16 64}, the dominant frequency of one pseudo-random 3-frequency wave for frequency-domain excitation is { 14 16}, the frequency set is an equal ratio series, and therefore, the frequency set is described by defining octave Oct and octave coefficient b, namely, the parameters of the 7-frequency waves are Oct =2, b =6; the parameters of the 3-frequency wave are Oct =4, b =2;
to meet the requirements of different depths of exploration, the frequency range needs to be adjusted for Oct (0:b) The formula being increased by a factor, i.e.Fundamental frequency f base Then the expression is f base *Oct (0:b) E.g. the above-mentioned 7-frequency wave, by a factor f base =128, the frequency combination is {128 256 512 1024 2048 8192}, and is more suitable for shallow exploration.
The combined wave is composed of a plurality of frequency sets with different orders, so that one group of wave forms comprises more frequencies, and the requirements of frequency density and frequency bandwidth are met. Order coefficient a is introduced for defining pseudo-random combined wave, single-frequency combined wave and high-order combined wave with different orders, and f is introduced due to the fact that the order coefficients are inconsistent base *Oct (0:b) Rewritten as f base *Oct (0:bn) Finally, combining the frequencies of each order to form a set to obtain:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) …}。
f represents a frequency combination of the main frequencies, an represents an order coefficient, bn represents an octave coefficient, fbase represents a fundamental frequency, and Oct represents an octave, e.g., a 7-frequency group of a conventional 7-frequency wave, which in this application may be determined by the parameters a1=1,fbase =1,b1=6,oct =2, and the main frequency thereof is {12 48 16 32 }, and when transmission is performed, only the relevant a1=1,fbase =1,b1=6,oct =2 parameter may be transmitted, and then the relevant parameter can be interpreted according to the determination rule of F; and as another example, 27-frequency waves with parameters a = [1 3], b = [9 8 ], oct =2, with the primary frequency of {8 16 24 32 40 48 64 96 128 160 192 256 320 384 512 768 1024 1280 1536 2048 2560 3072 306 5120 6144}.
The lowest frequency can be deduced according to the background resistivity and an empirical formula, the combination of a and b is used for defining the frequency density, the two are used for determining the longitudinal resolution capability of exploration, and the Oct parameter mainly adopts 2 or 4 at present and is mainly used for being compatible with an electroexcitation method and an electromagnetic method. The product of a1 and fbase is equal to the lowest frequency.
Alternatively, as shown with reference to fig. 3, said equation F = { a1 × F base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determiningThe steps of the parameters an, fbase, oct, and bn of the combined wave include:
s51, performing spectrum analysis on the combined signal to determine a useful frequency and a fundamental frequency in the combined signal;
the spectral analysis in this case can be implemented by means of software and tools.
S52, determining order coefficients, octave coefficients and octaves corresponding to the number of the useful frequencies according to the useful frequencies and the fundamental frequencies;
s53, determining a combined wave according to the order coefficient, the octave coefficient, the fundamental frequency and the octave.
The invention can adjust the expression form of the combined wave conveniently according to the detection requirement by reestablishing a set of frequency combination mode, can perfectly describe the combined wave and single frequency in various forms, can be highly adapted to various application occasions without carrying out frequency description of each detection signal, avoids that a user forgets various frequency group names and expression forms, directly expresses in the form of the combined wave, and is convenient for the direct use of the user. The technical problems that the frequency description mode of the single-frequency signal in the prior art cannot meet the requirement and a corresponding table needs to be established in a fixed naming mode are solved.
Optionally, the step of performing spectral analysis on the combined wave signal to determine a useful frequency and a fundamental frequency in the combined wave signal comprises:
carrying out spectrum analysis on the combined wave signal to obtain a spectrogram;
taking the original signal as a high-order pseudo-random encoded signal as an example, referring to fig. 2, the spectral analysis converts a time-amplitude time domain diagram into a frequency-amplitude spectrogram.
When the frequency amplitude in the spectrogram is larger than a preset threshold value, determining that the current frequency is a useful frequency;
the preset threshold value at this time can be set to any frequency amplitude value, and is set according to actual needs, based on fig. 2. The preset threshold value at this time may be set to 0.12, so that useful frequencies as in table 1 below may be screened out:
| 3072 | 2560 | 2048 | 1536 | 1280 | 1024 | 768 | 640 | 512 | 384 | 320 | 256 | 192 |
| 160 | 128 | 96 | 80 | 64 | 48 | 40 | 32 | 24 | 20 | 16 | 12 | 10 |
| 8 | 6 | 5 | 4 | 3 | 2.5 | 2 | 1.5 | 1.25 | 1 | 0.75 | 0.5 | 0.25 |
TABLE 1
The lowest of the useful frequencies is the fundamental frequency.
Based on the above table, the lowest frequency is therefore 0.25, i.e. the fundamental frequency fbase =0.25.
Optionally, the step of determining order coefficients, octave coefficients and octaves corresponding to the number of useful frequencies according to the useful frequencies and the fundamental frequency includes:
carrying out ratio processing on each useful frequency and the fundamental frequency to obtain a plurality of frequency coefficients;
useful frequencies/fundamental frequency = frequency coefficient, taking the data in table 1 above as an example, the useful frequencies are all divided by 0.25 to obtain table 2 below:
TABLE 2
Grouping a plurality of the frequency coefficients to obtain order coefficients;
the grouping at this time can be selected according to actual conditions, and in this embodiment, in order to fully describe the combined wave, decomposition is generally performed by a prime factor decomposition rule; determining a common ratio according to each frequency coefficient;
determining the number of the frequency coefficients of each group after grouping;
and determining octave coefficients according to the number of the frequency coefficients of each group.
Through the process, a plurality of parameters of the expression form of the combined wave can be determined, and the method has the characteristics of uniqueness and simplicity, so that the finally determined combined wave can completely describe the initial combined wave signal without adopting complicated Chinese nomenclature, and great convenience is brought to the use of the combined wave.
Optionally, the step of grouping a plurality of the frequency coefficients to obtain order coefficients includes:
performing prime factor decomposition on a plurality of the frequency coefficients to obtain at least one prime factor;
performing quality factor decomposition on all frequency coefficients to obtain a quality factor combination [ 52 ],
and determining an order coefficient according to the quality factor.
The frequency coefficients which can be divided by 5 are extracted [10240 5120 2560 1280 640 320 160 80 20 10 5], the remaining frequency coefficients which can be divided by 3 are extracted [12288 6144 3072 1536 768 384 192 48 12 6 ], all the remaining frequency coefficients except 1 can be divided by 2, and are [ 814096 2048 1024 512 256 128 32 16 84 ], the order coefficient a1=5, the order coefficient a2=3, and the order coefficient a3=2.
According to the above grouping, the common ratio of the frequency coefficient groups is 2, and therefore, the common ratio can be defined as an octave; i.e., oct =2.
Optionally, the step of determining the octave coefficients according to the number of each group includes:
and performing difference operation on the number of the frequency coefficients of each group to obtain the octave coefficients.
Based on the data in table 2, the number of frequencies in each group minus the difference value 1 is obtained as b, where b1=11, b2=12, and b3=12. Of course the difference may be selected as desired.
And restoring by substituting the parameters into a combined wave formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) …}
F＝{5*0.25*2 (0:11) ,3*0.25*2 (0:12) ,2*0.25*2 (0:12) get the following deformation:
F＝{5*0.25*2 (0:11) ,3*0.25*2 (0:12) ,1*0.25*2 (1:13) }
when the fundamental frequency coefficient is 1, it can be classified as 1 x 0.25 x 2 (0:12) Within a frequency group, get
F＝{5*0.25*2 (0:11) ,3*0.25*2 (0:12) ,1*0.25*2 (0:13) }
Therefore, the restoration can find that the combined wave formula provided by the application can perfectly describe the combined waves and single frequencies in various forms, can be highly adapted to various application occasions, avoids the user from forgetting various frequency group names and expression forms, directly expresses the combined waves in the form of the combined waves, and is convenient for the user to directly use the combined waves.
Optionally, the combined wave formula is:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) …}
wherein F represents the frequency combination of the main frequencies, an represents the order coefficient, bn represents the octave coefficient, F base Representing the fundamental frequency and Oct representing the octave.
In order to solve the above problem, referring to fig. 4, the present invention further provides a frequency domain combined wave normalization expression device, where the frequency domain combined wave normalization expression device includes:
the parameter acquisition module 10 is used for acquiring the depth range of the area to be explored;
the calculation module 20 determines a fundamental frequency according to the depth range of the region to be explored; determining frequency density according to the target scale of the region to be explored; determining a combined signal according to the fundamental frequency and the frequency density;
a control module 30 for controlling the operation of the apparatus according to a combined wave formula F = { a1= F base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining order coefficients an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal to realize the description of all the combined signals through a parameter matrix of parameters an, fbase, oct and bn of the combined wave;
and the sending module 40 is used for sending the described combined wave to the area to be explored.
In the above embodiment, when the combined signal is used for exploration, the expression form description of the combined wave can be conveniently adjusted according to the detection requirement by reestablishing a set of signal combination description mode without adopting a fixed naming mode description, thereby solving the technical problems of inconvenient signal application and poor expansibility in the fixed naming mode in the prior art.
In order to solve the above problem, the present invention also proposes a storage medium, and the computer program, when executed by a processor, causes the processor to execute the frequency domain combined wave normalization expression method as described above.
It should be noted that, since the storage medium of the present application includes all the steps of the above frequency domain combined wave normalization expression method, the storage medium may also implement all the schemes of the frequency domain combined wave normalization expression method, and has the same beneficial effects, and details are not described herein again.
And executing a frequency domain combined wave normalization expression method in the method embodiment. The above-described embodiments of the apparatus are merely illustrative, wherein the units illustrated as separate components may or may not be physically separate, i.e. may be located in one place, or may also be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art will appreciate that all or some of the steps, systems, and methods disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as is well known to those skilled in the art. Computer storage 15 storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by a computer. In addition, as is well known to those of ordinary skill in the art, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism 20 and includes any information delivery media
In order to solve the above problem, the present invention further provides a frequency domain combined wave normalization expression device, which includes a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor is caused to execute the steps of the frequency domain combined wave normalization expression method.
It should be noted that, since the frequency domain combined wave normalization expression apparatus of the present application includes all the steps of the above frequency domain combined wave normalization expression method, the frequency domain combined wave normalization expression apparatus can also implement all the schemes of the frequency domain combined wave normalization expression method, and has the same beneficial effects, and details are not repeated herein.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
Claims (9)
1. A frequency domain combined wave normalization expression method is characterized by comprising the following steps:
acquiring a depth range of a region to be explored;
determining a fundamental frequency according to the depth range of the region to be explored;
determining frequency density according to the target scale of the region to be explored;
determining a combined signal according to the fundamental frequency and the frequency density;
expressing the formula according to the combined signal:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) …}
determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal;
realizing the description of all combined signals through parameter matrixes of the parameters an, fbase, oct and bn of the combined wave;
and sending the described combined wave to the area to be explored.
2. The method according to claim 1, wherein the expression is expressed according to the formula F = { a1 x F = base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230determining the parameters an, fbase, oct and bn of the combined wave comprises:
performing spectral analysis on the combined signal to determine a useful frequency and a fundamental frequency in the combined signal;
determining order coefficients, octave coefficients and octaves corresponding to the number of the useful frequencies according to the useful frequencies and the fundamental frequency;
and determining a combined wave according to the order coefficient, the octave coefficient, the fundamental frequency and the octave.
3. The frequency-domain combined wave normalization expression method according to claim 2, wherein the step of performing spectral analysis on the combined wave signal to determine a useful frequency and a fundamental frequency in the combined wave signal comprises:
performing spectrum analysis on the combined wave signal to obtain a spectrogram;
when the frequency amplitude in the spectrogram is larger than a preset threshold value, determining that the current frequency is a useful frequency;
the lowest of all the useful frequencies is the fundamental frequency.
4. The method according to claim 3, wherein the step of determining the order coefficients, octave coefficients and octaves corresponding to the number of useful frequencies from the useful frequencies and the fundamental frequency comprises:
carrying out ratio processing on each useful frequency and the fundamental frequency to obtain a plurality of frequency coefficients;
grouping a plurality of the frequency coefficients to obtain order coefficients;
determining a common ratio according to each frequency coefficient;
defining the common ratio as an octave;
determining the number of the frequency coefficients of each group after grouping;
and determining octave coefficients according to the number of the frequency coefficients of each group.
5. The method as claimed in claim 4, wherein the step of grouping a plurality of the frequency coefficients to obtain order coefficients comprises:
performing prime factor decomposition on a plurality of the frequency coefficients to obtain at least one prime factor;
and determining an order coefficient according to the quality factor.
6. The method for normalized expression of a frequency-domain combined wave according to claim 5, wherein the step of determining the octave coefficients according to the number of each group comprises:
and performing difference operation on the number of the frequency coefficients of each group to obtain the octave coefficients.
7. A frequency domain combined wave normalization expression device, characterized in that the frequency domain combined wave normalization expression device comprises:
the parameter acquisition module is used for acquiring the depth range of the area to be explored;
the calculation module is used for determining a fundamental frequency according to the depth range of the region to be explored; determining frequency density according to the target scale of the region to be explored; determining a combined signal according to the fundamental frequency and the frequency density;
a control module that, according to a combined wave formula:
F＝{a1*f base *Oct (0:b1) ,a2*f base *Oct (0:b2) ,a3*f base *Oct (0:b3) 8230, determining an order coefficient an, a fundamental frequency fbase, an octave Oct and an octave coefficient bn corresponding to the combined signal to realize the description of all the combined signals through a parameter matrix of parameters an, fbase, oct and bn of the combined wave;
and the sending module is used for sending the described combined wave to the area to be explored.
8. A storage medium storing a computer program which, when executed by a processor, causes the processor to execute the frequency-domain combined wave normalization expression method according to any one of claims 1 to 6.
9. A frequency-domain combined wave normalization expression apparatus, characterized by comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of the frequency-domain combined wave normalization expression method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210777269.5A CN115236748B (en) | 2022-06-30 | 2022-06-30 | Frequency domain combination wave normalization expression method and device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210777269.5A CN115236748B (en) | 2022-06-30 | 2022-06-30 | Frequency domain combination wave normalization expression method and device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115236748A true CN115236748A (en) | 2022-10-25 |
| CN115236748B CN115236748B (en) | 2023-10-27 |
Family
ID=83671010
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210777269.5A Active CN115236748B (en) | 2022-06-30 | 2022-06-30 | Frequency domain combination wave normalization expression method and device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115236748B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117169880A (en) * | 2023-11-03 | 2023-12-05 | 深圳市吉奥地球科技有限公司 | Pseudo-random signal transmitting method, receiving method and system |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009158052A1 (en) * | 2008-06-24 | 2009-12-30 | Exxonmobil Upstream Research Company | Method for electromagnetic prospecting waveform design |
| CN108427145A (en) * | 2018-01-26 | 2018-08-21 | 吉林大学 | Air-ground frequency domain electromagnetic methods controllable frequency source detection signal pulse duration modulation method |
| CN111239839A (en) * | 2020-02-10 | 2020-06-05 | 山东大学 | Frequency spectrum density calculation method and device for frequency domain electromagnetic method |
| CN111505723A (en) * | 2020-04-27 | 2020-08-07 | 山东大学 | 2n sequence pseudorandom electromagnetic exploration signal generation method and system |
| CN111505722A (en) * | 2020-04-27 | 2020-08-07 | 山东大学 | Logarithmic non-uniform pseudorandom electromagnetic exploration signal generation method and system |
| CN111522064A (en) * | 2020-04-27 | 2020-08-11 | 山东大学 | High-order pseudorandom electromagnetic exploration signal generation method and system |
| CN112859177A (en) * | 2020-12-31 | 2021-05-28 | 中南大学 | Control method and system for automatic acquisition by wide-area electromagnetic method |
-
2022
- 2022-06-30 CN CN202210777269.5A patent/CN115236748B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009158052A1 (en) * | 2008-06-24 | 2009-12-30 | Exxonmobil Upstream Research Company | Method for electromagnetic prospecting waveform design |
| CN108427145A (en) * | 2018-01-26 | 2018-08-21 | 吉林大学 | Air-ground frequency domain electromagnetic methods controllable frequency source detection signal pulse duration modulation method |
| CN111239839A (en) * | 2020-02-10 | 2020-06-05 | 山东大学 | Frequency spectrum density calculation method and device for frequency domain electromagnetic method |
| CN111505723A (en) * | 2020-04-27 | 2020-08-07 | 山东大学 | 2n sequence pseudorandom electromagnetic exploration signal generation method and system |
| CN111505722A (en) * | 2020-04-27 | 2020-08-07 | 山东大学 | Logarithmic non-uniform pseudorandom electromagnetic exploration signal generation method and system |
| CN111522064A (en) * | 2020-04-27 | 2020-08-11 | 山东大学 | High-order pseudorandom electromagnetic exploration signal generation method and system |
| CN112859177A (en) * | 2020-12-31 | 2021-05-28 | 中南大学 | Control method and system for automatic acquisition by wide-area electromagnetic method |
Non-Patent Citations (1)
| Title |
|---|
| 张友山,席振铢,王鹤,王文,穆建宏: "复合波相对相位频率谱异常区分法", no. 05, pages 80 - 84 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117169880A (en) * | 2023-11-03 | 2023-12-05 | 深圳市吉奥地球科技有限公司 | Pseudo-random signal transmitting method, receiving method and system |
| CN117169880B (en) * | 2023-11-03 | 2024-02-09 | 深圳市吉奥地球科技有限公司 | Pseudo-random signal transmitting method, receiving method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115236748B (en) | 2023-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100966904B1 (en) | Apparatus, Methods and Recording medium for imaging a subsurface using the waveform inversion in the Laplace-Fourier domain | |
| Chen et al. | Detecting a known near-surface target through application of frequency-dependent traveltime tomography and full-waveform inversion to P-and SH-wave seismic refraction data | |
| CN102112894B (en) | Estimation of soil properties using waveforms of seismic surface waves | |
| RU2309436C2 (en) | Data processing method | |
| CN105467444B (en) | A kind of elastic wave full waveform inversion method and device | |
| US7679991B2 (en) | Processing of seismic data acquired using twin over/under streamers | |
| US7379385B2 (en) | Processing of seismic data acquired using over/under streamers and/or over/under sources | |
| US20190277993A1 (en) | Integrated method for estimation of seismic wavelets and synthesis of seismic records in depth domain | |
| US20220308249A1 (en) | Method and system for generating logarithmic non-uniform pseudo-random electromagnetic exploration signal | |
| US20090276159A1 (en) | Modeling and filtering coherent noise in seismic surveying | |
| CN111239839B (en) | Frequency spectrum density calculation method and device for frequency domain electromagnetic method | |
| US8139440B2 (en) | Spectral conditioning for surface seismic data | |
| CN115236748A (en) | Frequency domain combined wave normalization expression method and device and storage medium | |
| CN109738950B (en) | The noisy-type data primary wave inversion method of domain inverting is focused based on sparse 3 D | |
| Padhi et al. | 2-D ocean temperature and salinity images from pre-stack seismic waveform inversion methods: An example from the South China Sea | |
| CN106934183A (en) | Dispersion curve determines method and apparatus, and p-and s-wave velocity determines method and apparatus | |
| Montejo et al. | An empirical relationship between Fourier and response spectra using spectrum-compatible times series | |
| Aimar et al. | Novel techniques for in-situ estimation of shear-wave velocity and damping ratio through MASW testing Part I: A beamforming procedure for extracting Rayleigh-wave phase velocity and phase attenuation | |
| Oth et al. | Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling | |
| Rix | Near-surface site characterization using surface waves | |
| GB2599849A (en) | Low frequency acquisition with towed streamers | |
| CN115453629A (en) | Waveform correction method and device based on time-frequency domain adaptive filtering | |
| Luo et al. | Improvement of Frequency–Bessel Phase‐Velocity Spectra of Multicomponent Cross‐Correlation Functions from Seismic Ambient Noise | |
| Mahvelati et al. | Effects of near-surface spatial variability of soil stiffness on surface wave dispersion | |
| CN116088049B (en) | Least square inverse time migration seismic imaging method and device based on wavelet transformation |
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 |