Seismic data imaging method and device, electronic equipment and storage medium

Abstract

The embodiment of the invention provides a seismic data imaging method and device, electronic equipment and a storage medium. The method comprises the steps of determining frequency-division seismic data according to full-frequency seismic data of a target area, obtaining a well parameter matching factor of the target area, obtaining a to-be-output imaging result according to the well parameter matching factor and the frequency-division seismic data, and outputting the to-be-output imaging result when the to-be-output imaging result meets an output standard. And determining the to-be-output imaging result as a seismic frequency conversion imaging result corresponding to the full-frequency seismic data. According to the method provided by the scheme, the full-frequency seismic data of the target area is processed to obtain the frequency-division seismic data, and then imaging is performed, so that full-information migration imaging of specific favorable frequency bands and frequencies can be realized, possible false frequencies or technical risks and technical traps caused by filtering are eliminated, and the imaging accuracy is improved. The technical problem that the imaging result is inaccurate due to the fact that the seismic data of the target stratum cannot be effectively utilized in the prior art is solved, and the well fixing success rate is increased.

Description

Technical Field

[0001] The present invention relates to the technical field of petroleum geophysical exploration, and in particular to a seismic data imaging method, device, electronic equipment and storage medium. Background Art

[0002] In the field of petroleum geophysical exploration technology, the well location needs to be determined based on the imaging results of seismic data. Therefore, the accuracy of the seismic data imaging results will directly affect the technical risks of subsequent exploration.

[0003] In practical applications, seismic data has a wide bandwidth, and both low-frequency and high-frequency energies are very strong, while the seismic reflection band of the real target stratum is narrower, with lower energy, and has a complex relationship with nearby small faults and small contacts. These geological targets are not clear and intuitive enough in conventional seismic data, and their ranges cannot be accurately determined and quantified. The current imaging method based on improving the resolution of seismic data is to first superimpose the imaging results, and then filter and optimize the imaging results, so that most of the effective information has been submerged and overlapped, and cannot be effectively extracted and effectively utilized, which brings great technical pitfalls to subsequent geological exploration and well location implementation.

[0004] Therefore, the current seismic data imaging technology has the technical problem of being unable to effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results. Summary of the invention

[0005] In order to alleviate the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, the embodiments of the present invention provide a seismic data imaging method, device, electronic device and storage medium.

[0006] In a first aspect, an embodiment of the present invention provides a seismic data imaging method, comprising:

[0007] Determine the frequency-divided seismic data based on the full-frequency seismic data of the target area;

[0008] Obtaining the well parameter matching factor of the target area;

[0009] Obtaining an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data;

[0010] When the imaging result to be output meets the output standard, the imaging result to be output is determined as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data.

[0011] In some implementations, determining the fractional-frequency seismic data based on the full-frequency seismic data of the target area includes:

[0012] Obtaining a first frequency division frequency band range to be optimized;

[0013] Based on the first frequency division frequency band range and the full-frequency seismic data, a frequency division stacking difference is obtained;

[0014] Optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range;

[0015] According to the second frequency division frequency band range, the full-frequency seismic data is filtered to obtain the frequency division seismic data.

[0016] In some implementations, obtaining the frequency division stacking difference based on the first frequency division frequency band range and the full-frequency seismic data includes:

[0017] According to the first frequency division frequency range, filtering the full-frequency seismic data to obtain reference frequency division seismic data;

[0018] Performing stacking imaging processing on the reference frequency-divided seismic data to obtain frequency-divided stacking data;

[0019] Performing stack imaging processing on the full-frequency seismic data to obtain full-frequency stack data;

[0020] According to the first frequency division frequency band range, filtering the full-frequency superposition data to obtain superposition frequency division data;

[0021] The frequency division and superposition difference is determined according to the frequency division and superposition data and the superposition frequency division data.

[0022] In some implementations, optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range includes:

[0023] Determining an optimization direction of the first frequency division frequency band range according to the positive and negative attributes of the frequency division superposition difference;

[0024] Determining an optimized step length of the first frequency division frequency band range according to the absolute value of the frequency division and superposition difference;

[0025] According to the optimization direction and the optimization step, the first frequency division frequency band range is optimized to obtain the second frequency division frequency band range.

[0026] In some implementations, after optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range, the method further includes:

[0027] When the imaging result to be output does not meet the output standard, the first frequency division frequency band range is reset, and the frequency division seismic data is re-determined based on the full-frequency seismic data.

[0028] In some implementations, obtaining the imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data includes:

[0029] Determining imaging velocity according to the frequency-divided seismic data and the full-frequency seismic data;

[0030] According to the imaging speed and the well parameter matching factor, the frequency-divided seismic data is imaged to obtain the imaging result to be output.

[0031] In some embodiments, obtaining the well parameter matching factor of the target area includes:

[0032] Acquire well data curves of the target area;

[0033] The well data curve is fitted according to the seismic attribute parameters to determine the well parameter matching factor.

[0034] In some embodiments, after obtaining the imaging result to be output, the method further includes:

[0035] Obtaining a calibration result of the well data curve;

[0036] Comparing and analyzing the imaging result to be output with the calibration result to obtain a calibration result matching value;

[0037] When the calibration result coincidence value is greater than a preset value, determining that the imaging result to be output meets the output standard;

[0038] When the calibration result matching value is not greater than a preset value, it is determined that the imaging result to be output does not meet the output standard.

[0039] In a second aspect, an embodiment of the present invention provides a seismic data imaging device, comprising:

[0040] A first determination module is used to determine the fractional frequency seismic data according to the full frequency seismic data of the target area;

[0041] An acquisition module, used for acquiring the well parameter matching factor of the target area;

[0042] A second determination module is used to obtain an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data;

[0043] The third determination module is used to determine the imaging result to be output as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data when the imaging result to be output meets the output standard.

[0044] In a third aspect, an embodiment of the present invention provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the method described in the first aspect is implemented.

[0045] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by one or more processors, the method described in the first aspect is implemented.

[0046] Compared with the prior art, one or more embodiments of the present invention can at least bring the following beneficial effects:

[0047] The embodiment of the present invention provides a method, device, electronic device and storage medium for imaging seismic data; the method includes: determining frequency-divided seismic data based on full-frequency seismic data of the target area, obtaining the well parameter matching factor of the target area, obtaining the imaging result to be output based on the well parameter matching factor and the frequency-divided seismic data, and when the imaging result to be output meets the output standard, determining the imaging result to be output as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data. In the method provided by this scheme, the full-frequency seismic data of the target area is first processed to obtain the frequency-divided seismic data, and then imaging is performed, which can realize full-information offset imaging of specific favorable frequency bands, eliminate the technical risks and technical traps caused by possible false frequencies or filtering, alleviate the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, and improve the success rate of well determination. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.

[0049] Figure 1 A schematic diagram of a first flow chart of a seismic data imaging method provided in an embodiment of the present application;

[0050] Figure 2 A second flow chart of the seismic data imaging method provided in the embodiment of the present application;

[0051] Figure 3 A schematic diagram of a frequency scanning spectrum provided in an embodiment of the present application;

[0052] Figure 4 A schematic diagram of a well parameter curve fitted to well data provided in an embodiment of the present application;

[0053] Figure 5 A schematic diagram showing a cross-sectional comparison between the frequency conversion offset imaging and conventional imaging results provided in an embodiment of the present application;

[0054] Figure 6 A schematic diagram of the structure of a seismic data imaging device provided in an embodiment of the present application.

[0055] In the drawings, the same reference numerals are used for the same components, and the drawings are not drawn to scale. DETAILED DESCRIPTION

[0056] The following will describe the implementation methods of the present application in detail with reference to the accompanying drawings and embodiments, so that the implementation process of how the present application applies technical means to solve technical problems and achieve corresponding technical effects can be fully understood and implemented accordingly. The embodiments of the present application and the various features in the embodiments can be combined with each other without conflict, and the technical solutions formed are all within the scope of protection of the present application.

[0057] Meanwhile, in the following description, many specific details are set forth for the purpose of explanation to provide a thorough understanding of the embodiments of the present invention. However, it is obvious to those skilled in the art that the present invention can be implemented without the specific details or the specific manner described herein.

[0058] Embodiment 1

[0059] Figure 1 For a first flow chart of the seismic data imaging method provided in the embodiment of the present application, please refer to Figure 1 , the seismic data imaging method provided in this embodiment includes:

[0060] Step S110: Determine the fractional-frequency seismic data according to the full-frequency seismic data of the target area.

[0061] In some embodiments, determining the divided-frequency seismic data based on the full-frequency seismic data of the target area includes: obtaining a first divided-frequency band range to be optimized; obtaining a divided-frequency superposition difference based on the first divided-frequency band range and the full-frequency seismic data; optimizing the first divided-frequency band range according to the divided-frequency superposition difference to obtain an optimized second divided-frequency band range; and filtering the full-frequency seismic data according to the second divided-frequency band range to obtain the divided-frequency seismic data.

[0062] In the present application, the first frequency division frequency band range refers to the initial low frequency, middle frequency and high frequency band range.

[0063] In some embodiments, the frequency division and superposition difference is obtained based on the first frequency division frequency band range and the full-frequency seismic data, including: filtering the full-frequency seismic data according to the first frequency division frequency band range to obtain reference frequency division seismic data; performing superposition imaging processing on the reference frequency division seismic data to obtain frequency division and superposition data; performing superposition imaging processing on the full-frequency seismic data to obtain full-frequency superposition data; filtering the full-frequency superposition data according to the first frequency division frequency band range to obtain superposition frequency division data; and determining the frequency division and superposition difference according to the frequency division and superposition frequency division data.

[0064] Specifically, this embodiment performs spectrum analysis on the target layer segment demarcation time window of the full-frequency seismic data of the target area. According to the spectrum characteristics, the initial low-frequency, medium-frequency, and high-frequency segment ranges are first divided (i.e., the first frequency division frequency band range to be optimized is obtained); secondly, the full-frequency seismic data is divided into low-frequency, medium-frequency, and high-frequency seismic data by filtering according to this frequency band range (i.e., the full-frequency seismic data is filtered to obtain reference frequency division seismic data); then, the low-frequency, medium-frequency, and high-frequency seismic data divided by filtering are stacked and imaged to obtain frequency division stacking data; then, the full-frequency data without frequency division in the seismic data is stacked and imaged to obtain full-frequency stacking data; then, the full-frequency stacking data is filtered according to the initial low-frequency, medium-frequency, and high-frequency segment division ranges to obtain the corresponding low-frequency, medium-frequency, and high-frequency stacking frequency division data; finally, the frequency division stacking data and the stacked frequency division data are subtracted to obtain the frequency division stacking difference.

[0065] In some embodiments, the first frequency division frequency band range is optimized according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range, including: determining the optimization direction of the first frequency division frequency band range according to the positive and negative attributes of the frequency division superposition difference; determining the optimization step of the first frequency division frequency band range according to the absolute value of the frequency division superposition difference; and optimizing the first frequency division frequency band range according to the optimization direction and the optimization step to obtain the second frequency division frequency band range.

[0066] Specifically, this embodiment obtains the optimization direction of the first frequency division frequency band range according to the positive and negative attributes of the frequency division superposition difference, that is, increase or decrease; obtains the optimization step size of the first frequency division frequency band range according to the absolute value of the frequency division superposition difference; that is, the frequency division superposition difference of the low frequency, medium frequency and high frequency is quantitatively analyzed, and according to the analysis results, the initial frequency band range (that is, the first frequency division frequency band range) is optimized, and it is iterated repeatedly until the frequency division superposition difference is 0, and the final start and end frequency band ranges of the low frequency, medium frequency and high frequency are determined, that is, the second frequency division frequency band range is obtained.

[0067] In some embodiments, after optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain the optimized second frequency division frequency band range, it also includes: when the imaging result to be output does not meet the output standard, resetting the first frequency division frequency band range, and re-determining the frequency division seismic data based on the full-frequency seismic data.

[0068] Step S120: Obtain the well parameter matching factor of the target area.

[0069] In some embodiments, obtaining the well parameter matching factor of the target area includes: obtaining a well data curve of the target area; fitting the well data curve according to seismic attribute parameters to determine the well parameter matching factor.

[0070] In the present application, the seismic attribute parameters may include waveform, energy, amplitude, frequency, etc. of the seismic attributes.

[0071] Specifically, this embodiment quantitatively analyzes the relevant well data curves such as measurement, recording, and drilling that can be collected in the seismic work area (i.e., the target area), optimizes the curves by calibrating the artificial synthetic records, and performs nonlinear fitting calculations based on the waveform, energy, amplitude, frequency, etc. of the seismic attributes to obtain the well parameter fitting matching factor (i.e., the well parameter matching factor) suitable for this work area.

[0072] Step S130: obtaining an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data.

[0073] In some embodiments, obtaining the imaging result to be output based on the well parameter matching factor and the frequency-divided seismic data includes: determining an imaging speed based on the frequency-divided seismic data and the full-frequency seismic data; and imaging the frequency-divided seismic data based on the imaging speed and the well parameter matching factor to obtain the imaging result to be output.

[0074] Specifically, this embodiment will obtain accurate imaging velocity based on the frequency-divided seismic data and full-frequency seismic data in the low, medium and high frequency bands determined in step S110; then the imaging velocity and the well parameter matching factor fitted in step S120 will be used as input data for pre-stack imaging, and pre-stack frequency conversion imaging processing will be performed to obtain the imaging results to be output.

[0075] In some embodiments, after obtaining the imaging result to be output, it also includes: obtaining a calibration result of the well data curve; comparing and analyzing the imaging result to be output with the calibration result to obtain a calibration result matching value; when the calibration result matching value is greater than a preset value, determining that the imaging result to be output meets the output standard; when the calibration result matching value is not greater than the preset value, determining that the imaging result to be output does not meet the output standard.

[0076] Specifically, after obtaining the imaging result to be output, the calibration result of the well data curve is obtained, and the imaging result to be output is compared and analyzed with the calibration result to obtain the calibration result matching value, so as to determine whether to output the imaging result to be output; if the calibration result is completely consistent (that is, when the calibration result matching value is greater than the preset value), it means that the imaging result to be output meets the output standard, and the imaging data is output; if the calibration result does not meet the requirement (that is, when the calibration result matching value is not greater than the preset value), it means that the imaging result to be output does not meet the output standard, and then return to step S110, re-set the initial low-frequency, medium-frequency, and high-frequency bands, and then iterate to determine that the imaging result to be output that is completely consistent with the calibration result can be finally obtained.

[0077] Step S140: when the imaging result to be output meets the output standard, the imaging result to be output is determined as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data.

[0078] In this embodiment, following step S130, when the imaging result to be output meets the output standard, the imaging result to be output is determined as a seismic frequency conversion imaging result corresponding to the full-frequency seismic data, and the seismic frequency conversion imaging result is output.

[0079] Through the above-mentioned seismic data imaging method, full information offset imaging of specific favorable frequency bands can be achieved, eliminating the technical risks and technical traps caused by possible false frequencies or filtering, alleviating the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, and improving the success rate of well positioning.

[0080] Figure 2 For a second flow chart of the seismic data imaging method provided in the embodiment of the present application, please refer to Figure 2 , the seismic data imaging method provided in this embodiment includes:

[0081] Step S210: Acquire seismic data.

[0082] In this embodiment, full-frequency seismic data of the target area is acquired.

[0083] Step S220: Frequency scanning and feature analysis of seismic data.

[0084] In this embodiment, spectrum analysis is performed on the target layer segment of the full-frequency seismic data of the target area during the time window demarcation. According to the spectrum characteristics, first, the initial low-frequency, medium-frequency, and high-frequency segment ranges are divided (i.e., corresponding to the first frequency division frequency band range in the above text); second, the full-frequency seismic data is divided into low-frequency, medium-frequency, and high-frequency seismic data by filtering according to this frequency band range, and the low-frequency, medium-frequency, and high-frequency seismic data divided by filtering are stacked and imaged to obtain frequency division stacked data; third, the full-frequency data without frequency division in the seismic data is stacked and imaged to obtain full-frequency stacked data; fourth, the full-frequency stacked data is filtered out according to the initial frequency band division to obtain the corresponding low, medium, and high-frequency stacked profiles (corresponding to the stacked frequency division data in the above text) ; Fifth, subtract the frequency division superposition data from the superposition frequency division data to obtain the frequency division superposition difference; Sixth, quantitatively analyze the three frequency division superposition differences of low frequency, medium frequency and high frequency, that is, obtain the optimization direction of the initial frequency band range according to the positive and negative attributes of the frequency division superposition difference, that is, increase or decrease; obtain the optimization step of the initial frequency band range according to the absolute value of the frequency division superposition difference; Seventh, according to the analysis results, return to the first step to optimize the initial frequency band range, iterate repeatedly until the frequency division superposition difference is 0, and determine the final low frequency, medium frequency and high frequency start and end frequency band ranges (corresponding to the second frequency division frequency band range mentioned above).

[0085] Step S230: well data fitting.

[0086] In this embodiment, the relevant well data curves such as measurement, recording and drilling that can be collected in the seismic work area are quantitatively analyzed, and after calibration through artificial synthetic records, the curves are optimized using mature reservoir software, and nonlinear fitting calculations are performed in combination with the waveform energy, amplitude and frequency of seismic attributes to obtain the well parameter fitting matching factor suitable for this work area (i.e., corresponding to the well parameter matching factor mentioned above).

[0087] Step S240: frequency conversion offset imaging.

[0088] In this embodiment, the divided-frequency seismic data and full-frequency seismic data obtained by frequency division according to the low, medium and high frequency bands determined in step S220 are used to obtain an accurate imaging velocity; the imaging velocity and the well parameter matching factor fitted in step S230 are then used as input data for pre-stack imaging, and pre-stack frequency conversion imaging is performed to obtain pre-stack frequency conversion imaging results (i.e., corresponding to the above "determining the imaging velocity according to the divided-frequency seismic data and the full-frequency seismic data; imaging the divided-frequency seismic data according to the imaging velocity and the well parameter matching factor to obtain the imaging result to be output").

[0089] Step S250: iterative quality control.

[0090] In this embodiment, the prestack frequency conversion imaging result obtained in step S240 (corresponding to the imaging result to be output in the above text) is subjected to manual intervention result quality control. The obtained prestack frequency conversion imaging result and the calibration result of the well data are compared and analyzed to obtain the calibration result matching value to determine whether to output the final imaging data; if the calibration result is completely consistent, the final imaging data (the final imaging data is the prestack frequency conversion imaging result) is output, and step S260 is executed; if the calibration result does not meet the requirements (i.e., it is not completely consistent), return to step S220, re-set the initial frequency band, and then iterate to determine that the prestack frequency conversion imaging result with completely consistent calibration can be finally obtained.

[0091] Step S260: outputting final imaging data.

[0092] In this embodiment, following step S250, if the output pre-stack frequency conversion imaging result is completely consistent with the well calibration result, the pre-stack frequency conversion imaging result is output as the final imaging data.

[0093] In summary, the method provided in this embodiment first processes the full-frequency seismic data of the target area to obtain frequency-divided seismic data, and then performs imaging, which can achieve full-information offset imaging of frequencies in specific favorable frequency bands, eliminate technical risks and technical traps caused by possible false frequencies or filtering, alleviate the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, and improve the success rate of well positioning.

[0094] Based on the above records, this application is now described in detail in combination with specific scenarios.

[0095] A 3D seismic work area was randomly selected for testing using the seismic data imaging method of the present application. The work area is a mountainous area with large terrain changes, and the altitude of the work area is between 600-800m. The work area adopts 3D seismic acquisition; the design area element is 25m×25m, the design coverage times are 18 (vertical)×12 (horizontal)=216 times, and the horizontal-to-vertical ratio is 0.66. After pre-stack high-precision preprocessing of the data, the quality control line of the data is subjected to frequency scanning and feature analysis (such as Figure 3 As shown, Figure 3 The frequency scanning spectrum diagram provided in the embodiment of the present application is shown in Figure 1. Seismic data of different dominant frequency bands are split according to the analysis results. The well parameter curves of the main strata and target strata of the work area are fitted based on the drilling and logging data results collected from the work area (such as Figure 4 As shown, Figure 4Schematic diagram of the well data fitting well parameter curve provided in the embodiment of the present application). The seismic data split into different dominant frequency bands and the fitted well parameter data are used as input items, and pre-stack frequency conversion imaging processing is performed to obtain pre-stack frequency conversion imaging results; thereafter, data quality control is performed on the obtained pre-stack frequency conversion imaging results. If the accuracy of the data imaging results (i.e., the degree of consistency of the calibration results mentioned above) does not meet the requirements, the characteristic analysis of the returned frequency scanning results is performed to re-divide the data segments, and data optimization and frequency conversion imaging are continuously performed; if the quality control results reach the designed accuracy, the results are output and applied.

[0096] Figure 5 Schematic diagram of cross-sectional comparison between frequency conversion offset imaging and conventional imaging results provided in the embodiment of the present application; Figure 5 As shown, the prestack frequency conversion imaging results obtained by the seismic data imaging method of the present application show that the application effect is relatively ideal. The seismic data imaging method of the present application is not only effective for target processing areas, but also very effective for large desert areas and piedmont target imaging. It has a better effect on re-exploitation of old data and better profile application effect; seismic profiles are also better improved than other imaging methods. The drilling success rate is increased by about 5-8%, which verifies the feasibility of the seismic data imaging method of the present application.

[0097] In summary, the method provided in this embodiment first processes the full-frequency seismic data of the target area to obtain frequency-divided seismic data, and then performs imaging, which can achieve full-information offset imaging of frequencies in specific favorable frequency bands, eliminate technical risks and technical traps caused by possible false frequencies or filtering, alleviate the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, and improve the success rate of well positioning.

[0098] Embodiment 2

[0099] Figure 6 For a structural diagram of a seismic data imaging device provided in an embodiment of the present application, please refer to Figure 6 , the seismic data imaging device provided in this embodiment includes:

[0100] A first determination module 610 is used to determine the fractional frequency seismic data according to the full frequency seismic data of the target area;

[0101] An acquisition module 620 is used to acquire the well parameter matching factor of the target area;

[0102] A second determination module 630 is used to obtain an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data;

[0103] The third determination module 640 is used to determine the imaging result to be output as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data when the imaging result to be output meets the output standard.

[0104] In some embodiments, the first determination module 610 is further used to obtain a first frequency division frequency band range to be optimized; obtain a frequency division superposition difference based on the first frequency division frequency band range and the full-frequency seismic data; optimize the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range; and filter the full-frequency seismic data according to the second frequency division frequency band range to obtain the frequency division seismic data.

[0105] In some embodiments, the first determination module 610 is further used to filter the full-frequency seismic data according to the first frequency division frequency band range to obtain reference frequency division seismic data; perform stacking imaging processing on the reference frequency division seismic data to obtain frequency division stacking data; perform stacking imaging processing on the full-frequency seismic data to obtain full-frequency stacking data; filter the full-frequency stacking data according to the first frequency division frequency band range to obtain stacked frequency division data; and determine the frequency division stacking difference according to the frequency division stacking data and the stacked frequency division data.

[0106] In some embodiments, the first determination module 610 is further used to determine an optimization direction of the first frequency division frequency band range according to the positive and negative attributes of the frequency division superposition difference; determine an optimization step of the first frequency division frequency band range according to the absolute value of the frequency division superposition difference; and optimize the first frequency division frequency band range according to the optimization direction and the optimization step to obtain the second frequency division frequency band range.

[0107] In some embodiments, the second determination module 630 is further used to determine an imaging velocity based on the fractionated-frequency seismic data and the full-frequency seismic data; and to perform imaging processing on the fractionated-frequency seismic data based on the imaging velocity and the well parameter matching factor to obtain the imaging result to be output.

[0108] The specific implementation examples based on the above-mentioned modules for executing the method have been described in detail in the first embodiment and will not be repeated here.

[0109] Embodiment 3

[0110] This embodiment provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the seismic data imaging method described in Embodiment 1 is implemented. It is understood that the electronic device may also include an input / output (I / O) interface and a communication component.

[0111] The processor is used to execute all or part of the steps in the seismic data imaging method in Embodiment 1. The memory is used to store various types of data, which may include instructions of any application or method in the terminal device, and data related to the application.

[0112] The processor can be an application specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is used to execute the seismic data imaging method in the above-mentioned embodiment 1.

[0113] The memory can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk or optical disk.

[0114] Embodiment 4

[0115] This embodiment also provides a computer-readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card-type memory (for example, an SD or DX memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a disk, an optical disk, a server, an App application store, etc., wherein a computer program is stored on the storage medium, and when the computer program is executed by a processor, the following method steps can be implemented:

[0116] Determine the frequency-divided seismic data based on the full-frequency seismic data of the target area;

[0117] Obtaining the well parameter matching factor of the target area;

[0118] Obtaining an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data;

[0119] When the imaging result to be output meets the output standard, the imaging result to be output is determined as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data.

[0120] The specific implementation process of the above method steps can be found in Example 1, and this embodiment will not be repeated here.

[0121] In summary, the embodiments of the present invention provide a method, device, electronic device and storage medium for imaging seismic data; the method includes: determining frequency-divided seismic data based on the full-frequency seismic data of the target area, obtaining the well parameter matching factor of the target area, obtaining the imaging result to be output based on the well parameter matching factor and the frequency-divided seismic data, and when the imaging result to be output meets the output standard, determining the imaging result to be output as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data. In the method provided by this scheme, the full-frequency seismic data of the target area is first processed to obtain the frequency-divided seismic data, and then imaging is performed, which can realize full-information offset imaging of specific favorable frequency bands, eliminate the technical risks and technical traps caused by possible false frequencies or filtering, alleviate the technical problem that the current technology cannot effectively utilize the seismic data of the target formation, resulting in inaccurate imaging results, and improve the success rate of well determination.

[0122] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0123] Although the implementation methods disclosed in this application are as above, the contents described are only implementation methods adopted for facilitating the understanding of this application, and are not intended to limit this application. Any technician in the technical field to which this application belongs can make any modifications and changes in the form and details of implementation without departing from the spirit and scope disclosed in this application, but the scope of patent protection of this application shall still be based on the scope defined in the attached claims.

Claims

1. A seismic data imaging method, characterized in that: include: Determine the frequency-divided seismic data based on the full-frequency seismic data of the target area; Obtaining the well parameter matching factor of the target area; Obtaining an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data; When the imaging result to be output meets the output standard, the imaging result to be output is determined as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data.

2. The seismic data imaging method according to claim 1, characterized in that: Determining the frequency-divided seismic data according to the full-frequency seismic data of the target area includes: Obtaining the first frequency division frequency band range to be optimized; Based on the first frequency division frequency band range and the full-frequency seismic data, a frequency division stacking difference is obtained; Optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range; According to the second frequency division frequency band range, the full-frequency seismic data is filtered to obtain the frequency division seismic data.

3. The seismic data imaging method according to claim 2, characterized in that: The obtaining of the frequency division stacking difference based on the first frequency division frequency band range and the full-frequency seismic data comprises: According to the first frequency division frequency range, filtering the full-frequency seismic data to obtain reference frequency division seismic data; Performing stacking imaging processing on the reference frequency-divided seismic data to obtain frequency-divided stacking data; Performing stack imaging processing on the full-frequency seismic data to obtain full-frequency stack data; According to the first frequency division frequency band range, filtering the full-frequency superposition data to obtain superposition frequency division data; The frequency division and superposition difference is determined according to the frequency division and superposition data and the superposition and division data.

4. The seismic data imaging method according to claim 2, characterized in that: The step of optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain an optimized second frequency division frequency band range includes: Determining an optimization direction of the first frequency division frequency band range according to the positive and negative attributes of the frequency division superposition difference; Determining an optimized step length of the first frequency division frequency band range according to the absolute value of the frequency division superposition difference; According to the optimization direction and the optimization step, the first frequency division frequency band range is optimized to obtain the second frequency division frequency band range.

5. The seismic data imaging method according to claim 2, characterized in that: After optimizing the first frequency division frequency band range according to the frequency division superposition difference to obtain the optimized second frequency division frequency band range, the method further includes: When the imaging result to be output does not meet the output standard, the first frequency division frequency band range is reset, and the frequency division seismic data is re-determined based on the full-frequency seismic data.

6. The seismic data imaging method according to claim 1, characterized in that: The step of obtaining the imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data includes: Determining imaging velocity according to the frequency-divided seismic data and the full-frequency seismic data; According to the imaging speed and the well parameter matching factor, the frequency-divided seismic data is imaged to obtain the imaging result to be output.

7. The seismic data imaging method according to any one of claims 1 to 6, characterized in that: The step of obtaining the well parameter matching factor of the target area includes: Acquire well data curves of the target area; The well data curve is fitted according to the seismic attribute parameters to determine the well parameter matching factor.

8. The seismic data imaging method according to claim 7, characterized in that: After obtaining the imaging result to be output, the method further includes: Obtaining a calibration result of the well data curve; Comparing and analyzing the imaging result to be output with the calibration result to obtain a calibration result matching value; When the calibration result coincidence value is greater than a preset value, determining that the imaging result to be output meets the output standard; When the calibration result matching value is not greater than a preset value, it is determined that the imaging result to be output does not meet the output standard.

9. A seismic data imaging device, characterized in that: include: A first determination module is used to determine the fractional frequency seismic data according to the full frequency seismic data of the target area; An acquisition module, used for acquiring the well parameter matching factor of the target area; A second determination module is used to obtain an imaging result to be output according to the well parameter matching factor and the frequency-divided seismic data; The third determination module is used to determine the imaging result to be output as the seismic frequency conversion imaging result corresponding to the full-frequency seismic data when the imaging result to be output meets the output standard.

10. An electronic device, characterized in that: The method comprises a memory and a processor; a computer program is stored in the memory, and when the computer program is executed by the processor, the method according to any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by one or more processors, the method according to any one of claims 1 to 8 is implemented.

Images