Data weak information processing method, system, equipment and medium

Abstract

The invention relates to the technical field of oil-gas exploration data processing, and discloses a data weak information processing method, system and device and a medium. The method comprises the following steps: optimizing time domain signal data of all frequencies in a target layer, and calculating full-band amplitude corresponding to each excitation frequency according to each optimized time domain signal data; judging whether each excitation frequency is greater than a preset frequency; if the excitation frequency is greater than the preset frequency, extracting a plurality of low-order harmonic amplitudes of the frequency according to the corresponding full-band amplitude to generate first target data; if the excitation frequency is smaller than or equal to the preset frequency, extracting a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-band amplitude to generate second target data; and performing constraint inversion based on all the first target data and all the second target data to determine an inversion profile of the target layer, and extracting weak information of the target layer according to the inversion profile. According to the scheme, the precision of data processing and interpretation of the controllable source electromagnetic method is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas exploration data processing, and in particular to a method, system, equipment and medium for processing weak data information. Background Art

[0002] Currently, controlled-source electromagnetic methods are widely used in oil and gas exploration. Conventional electromagnetic exploration techniques use an acquisition instrument at a fixed sampling rate to acquire discrete time-domain signals. These signals are then converted to the frequency domain using a Fourier transform. When using a zero-crossing square wave excitation, the time-frequency conversion of discrete time-domain signals can eliminate even harmonics, leaving only odd harmonics. Due to the fundamental properties of odd harmonics, the methods used in related techniques for handling these harmonics suffer from insufficient inversion information.

[0003] Therefore, it is urgent to propose a data weak information processing method to solve the above defects. Summary of the Invention

[0004] In view of this, the present invention proposes a data weak information processing method, system, computer device and computer-readable medium, which solves the problems of insufficient observation data during inversion calculation and low vertical resolution of inversion profile in the data interpretation process of electromagnetic exploration. It uses multiple harmonic information to realize frequency point encryption of the frequency range corresponding to the target layer, enriches the amount of observation data during inversion calculation, and improves the accuracy of controlled source electromagnetic data processing and interpretation.

[0005] Based on the above objectives, an embodiment of the present invention provides a method for processing weak data information, which specifically includes the following steps:

[0006] Optimizing the time domain signal data of all frequencies in the target layer, and calculating the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data;

[0007] Determining whether each of the excitation frequencies is greater than a preset frequency;

[0008] If the excitation frequency is greater than the preset frequency, determining the amplitudes of several low-order harmonics of the frequency according to the corresponding full-frequency amplitude to generate first target data;

[0009] If the excitation frequency is less than or equal to the preset frequency, determining a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-frequency amplitude to generate second target data;

[0010] Constrained inversion is performed based on all of the first target data and all of the second target data to determine an inversion section of the target layer, and weak information of the target layer is extracted according to the inversion section.

[0011] In some embodiments, the step of performing constrained inversion based on all of the first target data and all of the second target data to determine the inversion profile of the target layer includes:

[0012] sorting all the first target data and all the second target data according to the magnitude of the corresponding frequencies and storing them in a first database as amplitude data;

[0013] Optimizing each amplitude data in the first database, sequentially deriving the optimized amplitude data from the first database along the survey line and performing constrained inversion on the data to obtain a resistivity inversion profile and a polarizability inversion profile;

[0014] The resistivity inversion section and the polarizability inversion section are used as the inversion section of the target layer.

[0015] In some embodiments, the method for processing weak data information further includes:

[0016] Establishing a geoelectric model based on drilling data and well logging data, performing forward simulation based on the geoelectric model to calculate response results of the target layer at different frequencies, and determining a frequency response range based on all the response results;

[0017] The product of the lowest frequency corresponding to the frequency response range and the preset multiple is used as the preset frequency.

[0018] In some embodiments, the step of extracting weak information of the target layer according to the inversion profile includes:

[0019] The inversion section is subjected to layer flattening processing, and weak information of the target layer is extracted from the processed inversion section based on a residual gradient method and a local magnification method.

[0020] In some embodiments, the step of calculating the full-band amplitude corresponding to each excitation frequency based on the optimized time-domain signal data includes:

[0021] The optimized time domain signal data are respectively subjected to short-time Fourier transform calculation to obtain the full-band amplitude corresponding to each excitation frequency.

[0022] In some embodiments, the step of optimizing the time domain signal data of all frequencies in the target layer includes:

[0023] The time domain signal data are respectively subjected to a DC drift removal process, a 50 Hz power frequency filtering process, and a multi-cycle superposition process to obtain optimized time domain signal data.

[0024] In some embodiments, the step of optimizing each amplitude data in the first database includes:

[0025] Preprocessing is performed on each of the amplitude data to optimize each of the amplitude data, wherein the preprocessing includes noise suppression, device coefficient normalization, and static displacement correction.

[0026] Another aspect of the present invention provides a weak data information processing system, comprising:

[0027] a calculation unit, configured to optimize the time domain signal data of all frequencies in the target layer, and calculate the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data;

[0028] a judging unit, configured to judge whether each of the excitation frequencies is greater than a preset frequency;

[0029] a first generating unit, configured to determine the amplitudes of several low-order harmonics of the frequency according to the corresponding full-frequency amplitude to generate first target data if the excitation frequency is greater than the preset frequency;

[0030] a second generating unit, configured to determine, if the excitation frequency is less than or equal to the preset frequency, a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-frequency amplitude to generate second target data;

[0031] An extraction unit is configured to perform constrained inversion based on all of the first target data and all of the second target data to determine an inversion profile of the target layer, and extract weak information of the target layer according to the inversion profile.

[0032] According to another aspect of the embodiments of the present invention, a computer device is provided, comprising: at least one processor; and a memory, wherein the memory stores a computer program that can be run on the processor, and the computer program implements the steps of the above method when executed by the processor.

[0033] According to another aspect of the embodiments of the present invention, a computer-readable storage medium is provided, which stores a computer program that implements the above method steps when executed by a processor.

[0034] The present invention has at least the following beneficial technical effects: The weak information processing method of the present invention, by combining the characteristics of strong low-frequency energy and weak high-frequency energy of the background field during excitation of a controlled source electromagnetic emission system, extracts the corresponding higher harmonics of the low-frequency portion in addition to the high-frequency portion, effectively obtaining multiple harmonic information of the target frequency band, thereby achieving frequency point encryption of the frequency band corresponding to the target layer. Furthermore, through fine modeling, constrained inversion, and weak information extraction, the resistivity and polarizability of the target layer are calculated, thereby increasing the profile interpretation accuracy from the 100-meter level to the 10-meter level. Furthermore, the weak information processing method of the present invention has the advantages of fast calculation speed, ease of implementation, and high efficiency.

[0035] In addition, the present invention also provides a data weak information processing system, a computer device and a computer-readable storage medium, which can also achieve the above-mentioned technical effects and will not be repeated here. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other embodiments can be obtained based on these drawings without paying any creative work.

[0037] Figure 1 A block diagram of an embodiment of a method for processing weak data information provided by the present invention;

[0038] Figure 2 A schematic diagram of an embodiment of the frequency domain amplitude curve comparison provided by the present invention;

[0039] Figure 3 A schematic diagram of an embodiment of the processing and interpretation cross-section of each sub-layer of a certain line dragon-1 section provided by the present invention;

[0040] Figure 4 A schematic diagram of an embodiment of a weak data information processing system provided by the present invention;

[0041] Figure 5 A schematic structural diagram of an embodiment of a computer device provided by the present invention;

[0042] Figure 6 This is a schematic structural diagram of an embodiment of a computer-readable storage medium provided by the present invention. DETAILED DESCRIPTION

[0043] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the embodiments of the present invention are further described in detail below in conjunction with specific embodiments and with reference to the accompanying drawings.

[0044] It should be noted that all expressions using "first" and "second" in the embodiments of the present invention are for distinguishing two non-identical entities with the same name or non-identical parameters. It can be seen that "first" and "second" are only for the convenience of expression and should not be understood as limitations on the embodiments of the present invention. Subsequent embodiments will not explain this one by one.

[0045] When controlled-source electromagnetic methods employ zero-passing square wave excitation in oil and gas exploration, even harmonics are eliminated during the time-frequency conversion of discrete time-domain signals, leaving only odd harmonics. Because odd harmonics become weaker as their multiples of the fundamental wave increase, related techniques for processing odd harmonics have focused solely on the fundamental, third, and fifth harmonics, ignoring higher harmonics. This results in insufficient observational data for inversion calculations, compromising the vertical resolution of the inversion profile.

[0046] In response to the above-mentioned problem of insufficient inversion information, the present invention proposes a method for processing weak data information, which performs frequency boosting processing on the target layer, fully utilizes multiple harmonic information, and realizes frequency point encryption of the frequency range corresponding to the target layer, thereby increasing the amount of data information to enrich the amount of observation data during inversion calculation, and improving the accuracy of controlled source electromagnetic data processing and interpretation.

[0047] Based on the above purpose, the first aspect of the embodiment of the present invention provides an embodiment of a method for processing weak data information. Figure 1 As shown, the data weak information processing method includes the following steps:

[0048] Step S100: Optimizing the time domain signal data of all frequencies in the target layer, and calculating the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data;

[0049] Step S200, determining whether each excitation frequency is greater than a preset frequency;

[0050] Step S300 , if the excitation frequency is greater than the preset frequency, determining the amplitudes of several low-order harmonics of the frequency according to the corresponding full-band amplitude to generate first target data;

[0051] Step S400 , if the excitation frequency is less than or equal to the preset frequency, determining a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-band amplitude to generate second target data;

[0052] Step S500 : performing constrained inversion based on all the first target data and all the second target data to determine an inversion profile of the target layer, and extracting weak information of the target layer according to the inversion profile.

[0053] In some embodiments, the time domain signal data of each frequency in the target layer are subjected to DC drift processing, industrial frequency 50Hz filtering processing and multi-cycle superposition processing respectively, which can extract useful information from the time domain signal data of each frequency, remove noise and interference, and improve the quality and reliability of the time domain signal data, so as to optimize the time domain signal data of all frequencies in the target layer.

[0054] In some embodiments, a window function is selected to perform short-time Fourier transform calculation on the optimized time domain signal data to obtain the amplitude of the full frequency band corresponding to the excitation frequency. Commonly used window functions include rectangular window, Hamming window and Hanning window. In the short-time Fourier transform calculation, the full-band amplitude can be obtained by calculating the spectrum of the signal in each time window and obtaining the amplitude of each frequency point. If the frequency in the target layer is greater than a preset multiple of the lowest frequency of the corresponding frequency response range, for example, greater than one-fifth of the lowest frequency, all low-order harmonic amplitudes of the frequency are extracted from the full-band amplitude corresponding to the frequency. The low-order harmonic amplitudes include the fundamental amplitude A0, the 3rd harmonic amplitude A3 and the 5th harmonic amplitude A5. If the frequency in the target layer is less than or equal to a preset multiple of the lowest frequency in the corresponding frequency response range, for example, less than or equal to one-fifth of the lowest frequency, then in addition to extracting all low-order harmonic amplitudes of that frequency from the full-band amplitude corresponding to that frequency, it is also necessary to extract a preset number of higher-order harmonic amplitudes of that frequency. These higher-order harmonic amplitudes include the 7th harmonic amplitude A7, the 9th harmonic amplitude A9, the 11th harmonic amplitude A11, and the nth harmonic amplitude A2n+1, where n is a preset value set according to actual research needs and is not limited by the present invention. After compensating the harmonic amplitude data for all frequencies in the target layer, the harmonic amplitude data for all frequencies are merged together and sorted by frequency to obtain frequency-domain amplitude data for all measurement points. The obtained amplitude data are then aggregated to form a first database. The amplitude data in the first database is then preprocessed, including noise suppression, device coefficient normalization, and static displacement correction. All preprocessed amplitude data are then aggregated to form a second database. Observation system information and preprocessed amplitude data are derived from the second database by survey line. A horizon constraint model is established using drilling and seismic data, and constrained inversion is performed on the derived amplitude data to obtain resistivity and polarizability inversion profiles. Layer flattening is performed on these resistivity and polarizability inversion profiles to obtain the target layer inversion profile. The residual gradient method and local amplification techniques are then used to extract weak information from the target layer, thereby improving the resolution of the inversion profile interpretation.

[0055] In some embodiments, a geoelectric model is established based on drilling and well logging data. The geoelectric model is a mathematical model that describes the electrical characteristics of the formation. After the geoelectric model is established, a forward simulation is performed to predict the electromagnetic response characteristics of the target layer at different frequencies. Based on the results of the forward simulation, the frequency response range corresponding to the target layer is analyzed. The determination of the frequency response range is helpful in providing guidance for subsequent exploration work. For example, the frequency response range can be from frequency f1 to frequency f2 (f1 <f2)。

[0056] The weak data information processing method of the present invention combines the characteristics of strong low-frequency energy and weak high-frequency energy of the background field during excitation of a controlled source electromagnetic emission system. Compared with the high-frequency part, the corresponding higher harmonics of the low-frequency part are extracted separately, effectively obtaining multiple harmonic information of the target frequency band, thereby achieving frequency point encryption of the frequency band corresponding to the target layer. Through detailed modeling, constrained inversion and weak information extraction, the resistivity and polarizability of the target layer are calculated, and the profile interpretation accuracy is increased from hundreds of meters to tens of meters. At the same time, the weak data information processing method of the present invention has the advantages of fast calculation speed, ease of implementation, and high efficiency.

[0057] In some embodiments, the step of performing constrained inversion based on all the first target data and all the second target data to determine the inversion profile of the target layer includes sorting all the first target data and all the second target data according to the size of their corresponding frequencies and storing them as amplitude data in the first database; optimizing each amplitude data in the first database, exporting the optimized amplitude data from the first database in sequence according to the survey line and performing constrained inversion on them to obtain a resistivity inversion profile and a polarizability inversion profile; and using the resistivity inversion profile and the polarizability inversion profile as the inversion profile of the target layer.

[0058] In some embodiments, the data weak information processing method of the present invention also includes: establishing a geoelectric model based on drilling data and well logging data, performing forward simulation based on the geoelectric model to calculate the response results of the target layer at different frequencies, and determining the frequency response range based on all response results; and taking the product of the lowest frequency corresponding to the frequency response range and a preset multiple as the preset frequency.

[0059] In some embodiments, the step of extracting weak information of the target layer according to the inversion profile includes: performing layer flattening processing on the inversion profile, and extracting the weak information of the target layer in the processed inversion profile based on the residual gradient method and the local magnification method.

[0060] In some embodiments, the step of calculating the full-band amplitude corresponding to each frequency based on the optimized time domain signal data includes: performing short-time Fourier transform calculations on each optimized time domain signal data to obtain the full-band amplitude corresponding to each frequency.

[0061] In some embodiments, the step of optimizing the time domain signal data of all frequencies in the target layer includes: performing DC drift removal processing, industrial frequency 50Hz filtering processing, and multi-cycle superposition processing on the time domain signal data to obtain optimized time domain signal data.

[0062] In some embodiments, the step of optimizing each amplitude data in the first database includes: preprocessing each amplitude data to optimize each amplitude data, wherein the preprocessing includes noise suppression, device coefficient normalization, and static displacement correction.

[0063] In some embodiments, the processing process of the method for processing weak data information is described below:

[0064] (1) Based on the drilling and logging data, a geoelectric model is established, and the frequency response range of the target layer is determined to be f1 to f2 (f1 <f2);

[0065] (2) Extract the time domain signal data D1 of any frequency F, perform DC drift removal, power frequency 50 Hz filtering, and multi-cycle superposition processing to obtain data D2;

[0066] (3) Select a window function to perform short-time Fourier transform calculation on the data D2 to obtain the full-band amplitude A;

[0067] (4) If the frequency F is greater than or equal to 1 / 5f1, extract the fundamental wave amplitude A0, the third harmonic amplitude A3, and the fifth harmonic amplitude A5 of the frequency F from the full-band amplitude A, and compensate for the above harmonic amplitudes to form data D3;

[0068] (5) If the frequency F is less than 1 / 5f1, in addition to extracting the fundamental wave amplitude A0, the 3rd harmonic amplitude A3, and the 5th harmonic amplitude A5 of the frequency F from the full-band amplitude A, it is also necessary to extract the 7th harmonic amplitude A7, the 9th harmonic amplitude A9, the 11th harmonic amplitude A11 ... A2n+1 of the frequency F from the full-band amplitude A, where n is a preset value, and at the same time compensate the above harmonic amplitudes to form data D4;

[0069] (6) Repeat steps (2) to (5) until the amplitudes of all frequencies of the target layer are calculated, merge the amplitude data of all frequencies together, and sort the amplitude data according to the corresponding frequency size to generate data D5;

[0070] (7) The frequency domain amplitude data D5 of all the measuring points obtained above are collected into a database M1, and then the database M1 is preprocessed by noise suppression, device coefficient normalization, and static displacement correction to obtain a database M2;

[0071] (8) Derive observation system information and amplitude data D6 from database M2 according to survey line to prepare for inversion calculation;

[0072] (9) Using drilling and seismic data to establish a horizon constraint model, appropriately reduce the vertical interval when segmenting the target layer, and use the data D6 derived in step (8) to perform constrained inversion.

[0073] Obtain the resistivity and polarizability inversion profile P1;

[0074] (10) After layer flattening processing on the polarizability inversion profile P1, the target layer inversion profile P2 is obtained. Then, the residual gradient method and local amplification technology are used to extract the weak information of the target layer from the target layer inversion profile P2, which ultimately improves the resolution of the profile interpretation.

[0075] In some embodiments, Figure 2 A schematic diagram of an embodiment of the frequency domain amplitude curve comparison provided by the present invention. Figure 3 This is a schematic diagram of an embodiment of the processing and interpretation section of each sub-layer of the Long-1 section of a certain line provided by the present invention. Low-resistivity shale reservoirs are common in the Qiongzhusi Formation and Wufeng Formation-Longmaxi Formation in the Sichuan Basin. Their lateral distribution is unclear, which seriously restricts the selection of favorable areas for future development. At the end of 2022, a controlled source electromagnetic exploration was carried out on the shale layer of the Wufeng-Longmaxi Formation in the Baozang-Yunjin syncline in Sichuan. 2 The acquisition, processing and interpretation of time-frequency electromagnetic data from 13 survey lines totaling 300 km was completed. In one example, using controlled source electromagnetic data from a certain area in Sichuan as an example, the specific process of processing and interpreting the weak information processing method of the present invention is as follows:

[0076] 1. According to the present invention, the above step (1) is performed to determine that the frequency response range corresponding to the Wufeng-Longmaxi Formation shale target layer is 1-10 Hz.

[0077] 2. For frequencies greater than or equal to 0.2 Hz, perform steps (2), (3), and (4). For frequencies less than 0.2 Hz (e.g., 0.1 Hz, 0.122 Hz, 0.138 Hz, 0.156 Hz, 0.177 Hz, and 0.199 Hz), perform steps (2), (3), and (5) to additionally extract the higher harmonic amplitudes A7, A9, A11, A13, and A15 of the frequency.

[0078] 3. Combine the amplitude data of all frequencies and sort them by frequency to generate data D5, as shown in Figure 2 As shown, the method of the present invention increases the number of frequency points in the target frequency band from 50 to 80.

[0079] 4. Perform steps (7) to (10) on data D5 to obtain Figure 3 The resistivity variation pattern of each sublayer in the Long-1 Member is shown.

[0080] The weak data information processing method of the present invention combines the characteristics of strong low-frequency energy and weak high-frequency energy of the background field during excitation of a controlled source electromagnetic emission system. Compared with the high-frequency part, the corresponding higher harmonics of the low-frequency part are extracted separately, effectively obtaining multiple harmonic information of the target frequency band, thereby achieving frequency point encryption of the frequency band corresponding to the target layer. Through detailed modeling, constrained inversion and weak information extraction, the resistivity and polarizability of the target layer are calculated, and the profile interpretation accuracy is increased from hundreds of meters to tens of meters. At the same time, the weak data information processing method of the present invention has the advantages of fast calculation speed, ease of implementation, and high efficiency.

[0081] It should be noted that the above figures are merely illustrative of the processes included in the method according to exemplary embodiments of the present invention and are not intended to be limiting. It is readily understood that the processes illustrated in the above figures do not indicate or limit the temporal order of these processes. Furthermore, it is readily understood that these processes may be executed synchronously or asynchronously, for example, in multiple modules.

[0082] It should be understood that, although the above is described in a certain order, these steps are not necessarily performed in sequence according to the above order. Unless there is clear explanation in this article, the execution of these steps does not have strict order restriction, and these steps can be performed in other orders. Moreover, a part of the steps of the present embodiment may include multiple steps or multiple stages, and these steps or stages are not necessarily performed at the same time, but can be performed at different times, and the execution order of these steps or stages is not necessarily carried out in sequence, but can be performed in turn or alternately with at least a portion of the steps or stages in other steps or other steps.

[0083] Based on the same inventive concept, according to another aspect of the present invention, Figure 4 As shown, a data weak information processing system is also provided, and the data weak information processing system specifically includes:

[0084] The calculation unit 110 is used to optimize the time domain signal data of all frequencies in the target layer, and calculate the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data;

[0085] The judging unit 120 is used to judge whether each excitation frequency is greater than a preset frequency;

[0086] The first generating unit 130 is configured to determine the amplitudes of several low-order harmonics of the frequency according to the corresponding full-band amplitudes to generate first target data if the excitation frequency is greater than a preset frequency;

[0087] The second generating unit 140 is configured to determine the amplitudes of several low-order harmonics and several high-order harmonics of the frequency according to the corresponding full-band amplitudes to generate second target data if the excitation frequency is less than or equal to the preset frequency;

[0088] The extraction unit 150 is configured to perform constrained inversion based on all the first target data and all the second target data to determine an inversion profile of the target layer, and extract weak information of the target layer according to the inversion profile.

[0089] The weak data information processing system of the present invention combines the characteristics of strong low-frequency energy and weak high-frequency energy of the background field during excitation of a controlled source electromagnetic emission system. Compared with the high-frequency part, it extracts the corresponding higher harmonics of the low-frequency part, effectively obtaining multiple harmonic information of the target frequency band, thereby achieving frequency point encryption of the frequency band corresponding to the target layer. Through detailed modeling, constrained inversion and weak information extraction, the resistivity and polarizability of the target layer are calculated, which increases the profile interpretation accuracy from hundreds of meters to tens of meters. At the same time, the weak data information processing method of the present invention has the advantages of fast calculation speed, ease of implementation, and high efficiency.

[0090] Based on the same inventive concept, according to another aspect of the present invention, Figure 5 As shown, an embodiment of the present invention further provides a computer device 30, which includes a processor 310 and a memory 320. The memory 320 stores a computer program 321 that can be run on the processor. When the processor 310 executes the program, the steps of the above method are performed.

[0091] Based on the same inventive concept, according to another aspect of the present invention, Figure 6 As shown, an embodiment of the present invention further provides a computer-readable storage medium 40 , which stores a computer program 410 for executing the above method when executed by a processor.

[0092] Finally, it should be noted that those skilled in the art will understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The program can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium of the program can be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM). The above-mentioned computer program embodiments can achieve the same or similar effects as any of the corresponding aforementioned method embodiments.

[0093] It will also be appreciated by those skilled in the art that the various exemplary logic blocks, modules, circuits and algorithmic steps described in conjunction with the disclosure herein can be implemented as electronic hardware, computer software or a combination of the two. In order to clearly illustrate this interchangeability of hardware and software, a general description has been given of the functions of various schematic components, blocks, modules, circuits and steps. Whether this function is implemented as software or hardware depends on specific applications and the design constraints imposed on the entire system. Those skilled in the art can implement the function in various ways for each specific application, but this implementation decision should not be interpreted as causing a departure from the disclosed scope of the embodiments of the present invention.

[0094] The above are exemplary embodiments disclosed in the present invention, but it should be noted that various changes and modifications can be made without departing from the scope of the disclosure of the embodiments of the present invention as defined in the claims. The functions, steps and / or actions of the method claims according to the disclosed embodiments described herein do not need to be performed in any particular order. The serial numbers of the embodiments disclosed in the above embodiments of the present invention are for description only and do not represent the advantages and disadvantages of the embodiments. In addition, although the elements disclosed in the embodiments of the present invention can be described or required in individual form, they can also be understood as multiple unless expressly limited to the singular.

[0095] It should be understood that, as used herein, the singular forms "a" and "an" are intended to include the plural forms as well, unless the context clearly supports an exception. It should also be understood that, as used herein, "and / or" is intended to include any and all possible combinations of one or more of the associated listed items.

[0096] Those skilled in the art should understand that the discussion of any of the above embodiments is merely illustrative and is not intended to imply that the scope of the disclosure of the embodiments of the present invention (including the claims) is limited to these examples. Within the spirit of the embodiments of the present invention, the technical features of the above embodiments or different embodiments may be combined, and there are many other variations of different aspects of the above embodiments of the present invention, which are not provided in detail for the sake of simplicity. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present invention should be included in the scope of protection of the embodiments of the present invention.

Claims

1. A method for processing weak data information, characterized in that: include: Optimizing the time domain signal data of all frequencies in the target layer, and calculating the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data; Determining whether each of the excitation frequencies is greater than a preset frequency; If the excitation frequency is greater than the preset frequency, determining the amplitudes of several low-order harmonics of the frequency according to the corresponding full-frequency amplitude to generate first target data; If the excitation frequency is less than or equal to the preset frequency, determining a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-frequency amplitude to generate second target data; Constrained inversion is performed based on all of the first target data and all of the second target data to determine an inversion section of the target layer, and weak information of the target layer is extracted according to the inversion section.

2. The method for processing weak data information according to claim 1, characterized in that: The step of performing constrained inversion based on all the first target data and all the second target data to determine the inversion profile of the target layer includes: sorting all the first target data and all the second target data according to the magnitude of the corresponding frequencies and storing them in a first database as amplitude data; Optimizing each amplitude data in the first database, sequentially deriving the optimized amplitude data from the first database along the survey line and performing constrained inversion on the data to obtain a resistivity inversion profile and a polarizability inversion profile; The resistivity inversion section and the polarizability inversion section are used as the inversion section of the target layer.

3. The method for processing weak data information according to claim 1, characterized in that: Also includes: Establishing a geoelectric model based on drilling data and well logging data, performing forward simulation based on the geoelectric model to calculate response results of the target layer at different frequencies, and determining a frequency response range based on all the response results; The product of the lowest frequency corresponding to the frequency response range and the preset multiple is used as the preset frequency.

4. The method for processing weak data information according to claim 1, characterized in that: The step of extracting weak information of the target layer according to the inversion profile includes: The inversion section is subjected to layer flattening processing, and weak information of the target layer is extracted from the processed inversion section based on a residual gradient method and a local magnification method.

5. The method for processing weak data information according to claim 1, characterized in that: The step of calculating the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data includes: The optimized time domain signal data are respectively subjected to short-time Fourier transform calculation to obtain the full-band amplitude corresponding to each excitation frequency.

6. The method for processing weak data information according to claim 1, characterized in that: The step of optimizing the time domain signal data of all frequencies in the target layer includes: The time domain signal data are respectively subjected to a DC drift removal process, a 50 Hz power frequency filtering process, and a multi-cycle superposition process to obtain optimized time domain signal data.

7. The method for processing weak data information according to claim 2, characterized in that: The step of optimizing each amplitude data in the first database includes: Preprocessing is performed on each of the amplitude data to optimize each of the amplitude data, wherein the preprocessing includes noise suppression, device coefficient normalization, and static displacement correction.

8. A data weak information processing system, characterized in that: include: a calculation unit, configured to optimize the time domain signal data of all frequencies in the target layer, and calculate the full-band amplitude corresponding to each excitation frequency based on the optimized time domain signal data; a judging unit, configured to judge whether each of the excitation frequencies is greater than a preset frequency; a first generating unit, configured to determine the amplitudes of several low-order harmonics of the frequency according to the corresponding full-frequency amplitude to generate first target data if the excitation frequency is greater than the preset frequency; a second generating unit, configured to determine, if the excitation frequency is less than or equal to the preset frequency, a plurality of low-order harmonic amplitudes and a plurality of high-order harmonic amplitudes of the frequency according to the corresponding full-frequency amplitude to generate second target data; An extraction unit is configured to perform constrained inversion based on all of the first target data and all of the second target data to determine an inversion profile of the target layer, and extract weak information of the target layer according to the inversion profile.

9. A computer device comprising: at least one processor; as well as A memory storing a computer program that can be run on the processor, wherein the processor executes the steps of the method according to any one of claims 1 to 7 when executing the program.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are performed.

Images