Desert multi-work-area three-dimensional seismic data amplitude consistency processing method and device

Abstract

The invention discloses a desert multi-work-area three-dimensional seismic data amplitude consistency processing method and device, and relates to the technical field of geological exploration, and the technical scheme is characterized in that phase, frequency and energy consistency processing and noise reduction processing are carried out; wherein the noise reduction processing adopts a dip angle-oriented constructional filtering technology, and the filtering degree of the dip angle-oriented constructional filtering technology is controlled through a diffusion coefficient in a PM diffusion equation. According to the method, the influence of different earth surface conditions, acquisition modes, processing methods and acquisition years on data is fully considered, the problem that amplitude preservation and consistency of continuous data processed by a traditional post-stack consistency continuous processing means cannot be considered at the same time is solved, researchers are helped to more effectively improve the data difference between adjacent blocks, and the data processing efficiency is improved. A good data foundation is provided for work of contiguous horizon, fault interpretation, oil reservoir description, reservoir prediction and the like, and an important data supporting effect is also achieved for three-dimensional seismic data splicing work.

Description

Technical Field

[0001] The present invention relates to the field of geological prospecting technology, and more specifically, to a method and device for processing amplitude consistency of three-dimensional seismic data in multiple working areas of deserts. Background Art

[0002] As exploration and development continue to deepen, the study of regional overall reservoir and structural change characteristics plays an increasingly important role in reservoir analysis and understanding of oil and gas enrichment laws. Seismic data from a single work area can no longer meet the needs of research, and research needs to carry out the splicing of cross-work area 3D seismic data. At present, the consistency processing of seismic 3D data is divided into two methods: pre-stack and post-stack.

[0003] Among them, post-stack consistency joint processing usually considers the relevant information of a single block of stacked data, such as the previous processing flow, parameters, survey line direction, bin size, static correction, polarity, geodetic coordinate system, etc., determines the bin size and direction of each data body, and then considers the processing of the single block of stacked data itself, the matching processing between blocks, and the processing of the data body after the merger of multiple blocks. The advantage is that although it cannot replace the pre-stack joint processing, the superimposed joint data saves time, labor, and money, and can make full use of the merged data body for efficient joint structural interpretation, mapping, and prediction.

[0004] Traditional post-stack consistency processing usually stacks multiple seismic records to produce a stacked section. Then, the consistency is improved by applying the frequency domain weighted averaging method to the stacked section. The basic principle of the frequency domain weighted averaging method is to perform weighted averaging based on the spectral characteristics of the seismic data. Seismic signals at different frequencies have different energies and amplitudes, and weighted averaging can make the contribution of signals at different frequencies more balanced.

[0005] However, in the desert surface coverage area of ​​the Tarim Basin, due to the complex surface geological conditions, the surface dunes are undulating and have huge height differences, the surface elevation changes dramatically, and the thickness of the loose sand layer varies greatly. Compared with the diagenetic medium, the loose sand dunes on the surface have a stronger absorption and attenuation effect on the high-frequency components in the seismic waves, resulting in large differences in seismic excitation, receiving conditions and seismic observation systems. The energy, frequency and signal-to-noise ratio of the seismic data collected in the field are significantly different. In addition, due to the different acquisition methods, processing methods and acquisition years, the seismic data in each work area have great differences in phase, frequency, signal-to-noise ratio and energy. The contiguous data processed by the above-mentioned post-stack consistency contiguous processing method often have poor amplitude preservation.

[0006] Therefore, it is necessary to develop a method and device for amplitude consistency processing of multi-area 3D seismic data on the desert surface that takes amplitude preservation into account. Summary of the invention

[0007] In view of the shortcomings of the prior art, the purpose of the present invention is to provide a method and device for amplitude consistency processing of three-dimensional seismic data in multiple working areas in the desert, which fully considers the influence of surface conditions, acquisition methods, processing methods and acquisition years on the data, takes into account the amplitude preservation and consistency of post-stack consistency connection processing means, helps researchers to more effectively improve the data differences between adjacent blocks, provides a better data basis for the work of connection layer, fault interpretation, oil field description, reservoir prediction, etc., and also plays an important data support role in the splicing of three-dimensional seismic data.

[0008] To achieve the above object, the present invention provides the following technical solution: a method for processing amplitude consistency of three-dimensional seismic data in multiple working areas of deserts, comprising the following steps:

[0009] S1. Based on the calibration of synthetic seismic records of well logging in all work areas, phase consistency processing is performed to make the seismic polarity of the marker layer consistent and form a new data body.

[0010] S2. Perform spectrum statistical analysis on the new data body, determine the target work area based on the requirements of the main work area and target layer system, and determine the spectrum characteristic distribution of all work areas. Based on the target work area, perform frequency unification processing on other work areas from the perspectives of frequency waveform and energy.

[0011] S3. Use the dip-guided structural filtering technology to reduce the noise of seismic data in all work areas. Use the diffusion coefficient in the PM diffusion equation to control the degree of filtering by the dip-guided structural filtering technology, slow down the diffusion speed at the edge of the image, and retain the details at the edge of the image.

[0012] S4. Energy consistency processing is performed from two perspectives: energy balance between multi-working areas and energy balance within single-working areas.

[0013] Select a standard work area, and calculate the consistency coefficient of seismic data between each work area adjacent to the standard work area and the standard work area. Then, gradually expand outward from the standard work area to all work areas, and calculate the consistency coefficient of seismic data between every two adjacent work areas. For each work area except the standard work area, different weights are assigned according to the corresponding consistency coefficient of seismic data, and energy balance is performed for multiple work areas.

[0014] Given a time window, the time window is slid point by point, the amplitude is normalized within each time window, and the energy at different time window positions is normalized by scaling to achieve energy balance within a single working area.

[0015] S5, the result of comprehensive processing of output phase, frequency, noise reduction and energy.

[0016] The present invention is further configured as follows: the frequency uniformization process in S2 includes:

[0017] B1. With the target work area as the center, select superimposed seismic traces in the target work area and overlapping seismic traces in other work areas. Based on amplitude consistency correction and according to the reflection coefficient consistency principle, gradually expand outward to perform continuous sub-wave shaping on the seismic data of other work areas.

[0018] B2. Perform time-frequency decomposition on the reflection coefficient and seismic signal of the logging data to obtain the reflection coefficient spectrum and seismic spectrum.

[0019] B3. Gradually approach the seismic spectrum to the reflection coefficient spectrum to obtain a spectrum blueing operator that matches the two. Apply the spectrum blueing operator to all seismic spectra along the frequency direction to synthesize a new seismic spectrum. Convert the synthesized seismic spectrum back to the time domain signal, perform time-frequency spectrum reconstruction, and obtain the time domain representation of the seismic signal.

[0020] B4. Output new seismic signal.

[0021] The present invention is further configured as follows: B2 also includes taking the seismic spectrum of the target work area as a standard spectrum, and using the standard spectrum to perform spectrum equalization on the seismic spectra of other work areas.

[0022] The present invention is further configured as follows: in B2, time-frequency spectrum decomposition is performed by continuous wavelet transform.

[0023] The present invention is further configured as follows: in B2, the reflection coefficient spectrum and the seismic spectrum are also subjected to spectrum blueing processing.

[0024] The present invention is further configured as follows: in B3, inverse wavelet transform is used to perform time-frequency spectrum reconstruction.

[0025] The present invention is further configured such that: after the S2 frequency uniformization process, the main frequency difference between other work areas and the target work area is less than 3 Hz.

[0026] The present invention is further configured as follows: the dip-oriented structural filtering technology in S3 uses the dip and azimuth of the actual stratum for constraints, and controls the response of the filter in different dip ranges and different azimuth ranges.

[0027] The present invention is further controlled as follows: in S4, the method for obtaining the consistency coefficient of seismic data between every two adjacent work areas is to respectively calculate the amplitude envelopes of the seismic traces of the two work areas, align the amplitude envelopes extracted from the two work areas on the time axis, so that the time points corresponding to the envelope data of the same event in different work areas are consistent, and then divide the envelope values ​​of the same time points at the overlapping positions in the adjacent work areas to obtain the consistency coefficient between the two adjacent work areas.

[0028] The present invention is further controlled as follows: the calculation of the seismic trace amplitude envelope adopts the Hilbert transform method.

[0029] The present invention is further controlled as follows: in S4, the method for performing normalization processing within each time window is to calculate the absolute average value, median value or root mean square value of the amplitude within the time window as the amplitude reference value, and use the inverse of the amplitude reference value as the weight coefficient of the time window to adjust the amplitude within the time window.

[0030] The present invention is further controlled as follows: S4 also includes quality control of energy balance in multiple working intervals: given a threshold value, it is determined whether the energy difference of the overlapping areas of two adjacent working intervals after energy balance in multiple working intervals is less than the threshold value; when the energy differences of the overlapping areas of a working area and all its adjacent working areas are less than the threshold value, the energy balance in the multiple working intervals is effective in the working area, and energy balance processing in a single working area is performed on the working area; when the energy difference of the overlapping area between a working area and a certain adjacent working area is greater than or equal to the threshold value, energy balance processing in a single working area is performed on the working area, or amplitude consistency processing of three-dimensional seismic data is abandoned for the working area, or parameters and data of the working area are optimized until energy balance in the multiple working intervals is effective in the working area, and then energy balance processing in a single working area is performed on the working area.

[0031] The present invention is further configured as follows: the threshold value used for the energy balance quality control in the multiplexer interval in S4 is 5%.

[0032] The present invention is further configured that the phase consistency processing comprises the following steps:

[0033] A1. Convert all independent tilted measurement networks into a unified large measurement network through network rotation transformation.

[0034] A2. According to the changes in the oil group interface and formation thickness, the reflection characteristics of the marker layer are determined based on the well point calibration.

[0035] A3. Based on the same seismic phase tracking principle, the consistency of time thickness and depth thickness in the vertical direction of the well point is achieved.

[0036] A4. Through phase migration, the original seismic volume data is spatially resampled and interpolated to form new seismic data.

[0037] The present invention is further configured such that: the similarity coefficient between the new seismic data wellside trace and the well point synthetic trace in A4 is greater than 98%.

[0038] The present invention is further configured as follows: the time window size in S4 is 500ms-1000ms.

[0039] The present invention also provides a device for processing amplitude consistency of 3D seismic data in multiple desert working areas, which includes a storage medium on which a computer program is stored, and a processor for implementing the above-mentioned method for processing amplitude consistency of 3D seismic data in multiple desert working areas when executing the computer program.

[0040] In summary, compared with the prior art, the present invention has the following beneficial effects: the present invention fully considers the impact of surface conditions, acquisition methods, processing methods and acquisition years on data, and makes up for the problem that the amplitude preservation and consistency of the contiguous data after processing by the traditional post-stack consistency contiguous processing method cannot be taken into account at the same time, and realizes a more comprehensive consideration of the phase, frequency, signal-to-noise ratio and energy of three-dimensional seismic data, ensures the relative relationship of the effective information of the original data within the time window, and helps researchers to more effectively improve the data differences between adjacent blocks, providing a better data basis for the work of contiguous layers, fault interpretation, oil reservoir description, reservoir prediction, etc., and also plays an important data support role in the splicing of three-dimensional seismic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Figure 1 It is a flow chart of a method for processing amplitude consistency of desert multi-area 3D seismic data in an embodiment;

[0042] Figure 2 This is the seismic profile before S1 phase processing;

[0043] Figure 3 This is the seismic profile after S1 phase processing;

[0044] Figure 4 It is a seismic profile diagram and spectrum statistical analysis diagram with a main frequency of 30hz;

[0045] Figure 5 It is a seismic profile diagram and spectrum statistical analysis diagram with a main frequency of 18hz;

[0046] Figure 6 This is the signal-to-noise ratio analysis diagram before S3 noise reduction processing;

[0047] Figure 7 This is the signal-to-noise ratio analysis diagram after S3 noise reduction processing;

[0048] Figure 8 a is the original seismic data, b is the seismic profile after energy balance in the single-work area using a 300ms time window, c is the seismic profile after energy balance in the single-work area using a 500ms time window; d is the seismic profile after energy balance in the single-work area using a 700ms time window; e is the seismic profile after energy balance in the single-work area using a 1000ms time window;

[0049] Fig. 9 is the original seismic data volume before processing;

[0050] Fig.10 This is the spliced ​​result output by S5 in this embodiment. DETAILED DESCRIPTION

[0051] The technical solution of the present invention will be clearly described below in conjunction with the accompanying drawings. Obviously, the described embodiments are not all embodiments of the present invention, and all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of the invention.

[0052] like Figure 1 FIG. 1 is a flow chart of a method for processing amplitude consistency of desert multi-area 3D seismic data in a preferred embodiment of the present invention, which includes the following steps:

[0053] S1: Based on the calibration of synthetic seismic records of well logging in all work areas, phase unification processing is performed to make the seismic polarity of the marker layer consistent and form a new data body. The schematic diagrams before and after S1 processing in this embodiment are respectively Figure 2 , Figure 3 , Figure 2 , Figure 3 The left side of the figure is other working areas, and the right side is the reference working area for phase migration. In this embodiment, the phase consistency processing includes the following steps:

[0054] A1. Convert all independent tilted measurement networks into a unified large measurement network through network rotation transformation.

[0055] A2. According to the changes in the oil layer interface and the formation thickness, the reflection characteristics of the marker layer are determined based on the well point calibration. That is, according to the changes in the oil layer and the formation thickness, the reflection characteristics of the marker layer are determined at the wellhead.

[0056] A3. Based on the same seismic phase tracking principle, the consistency of time thickness and depth thickness in the vertical direction of the well point is achieved, ensuring that the geological thickness above the wellhead remains consistent on the time axis and depth axis.

[0057] A4. Perform phase shift on the original seismic volume data, perform spatial resampling and interpolation processing according to the large measurement network, and form new seismic data. The similarity coefficient between the new seismic data wellside trace and the well point synthetic trace is greater than 98%.

[0058] S2, perform spectrum statistical analysis on the new data body, Figure 4 The frequency spectrum of the seismic profile with a main frequency of 30 Hz was analyzed. The analysis results are located in Figure 4 Lower left corner; Figure 5 The frequency spectrum of the seismic profile with a main frequency of 18 Hz was analyzed. The analysis results are located in Figure 5The lower right corner. Then, the target work area is determined in combination with the requirements of the main work area and the target layer system, and the frequency spectrum characteristic distribution of all work areas is determined. Based on the target work area, the frequency of other work areas is uniformly processed from the perspectives of frequency waveform and energy to ensure that the frequency of other work areas has a high consistency with the frequency of the target work area, that is, the main frequency difference between other work areas and the target work area is less than the set value, which is 3 Hz in this embodiment.

[0059] The target work area can be a work area with high data reliability and similar or close geological features to other work areas, or a work area with a greater possibility of well deployment based on experience. The frequency uniformization process includes the following steps:

[0060] B1. With the target work area as the center, select superimposed seismic traces in the target work area, select overlapping seismic traces in other work areas, and gradually expand outward to perform contiguous sub-wave shaping on the seismic data of other work areas based on amplitude consistency correction and the reflection coefficient consistency principle. Data consistency processing is achieved for other work areas around the target work area, except for amplitude differences caused by instrument differences, recording conditions and other factors, and ensure that the seismic records of different work areas at the same stratigraphic position reflect consistent underground structural changes, thereby improving the reliability of underground structure interpretation.

[0061] B2. Perform time-frequency decomposition of the reflection coefficient and seismic signal of the logging data through continuous wavelet transform to obtain the reflection coefficient spectrum and seismic spectrum, and perform spectrum blueing on the reflection coefficient spectrum and seismic spectrum. Then, the seismic spectrum of the target work area is used as the standard spectrum, and the seismic spectrum of other work areas is spectrally balanced using the standard spectrum, so that the seismic spectra of different work areas tend to be consistent in overall form, and the seismic response between each work area is more unified, which is convenient for comparative analysis.

[0062] B3. For each work area, the seismic spectrum is gradually approached to the reflection coefficient spectrum to obtain a spectrum blueing operator that matches the two. The spectrum blueing operator is applied to all seismic spectra along the frequency direction to achieve high-frequency migration, frequency reduction, and frequency increase processing, synthesize a new seismic spectrum, enhance the information related to the reflection coefficient spectrum in the seismic spectrum, and enhance the seismic spectrum's ability to interpret underground geological structures. The synthesized seismic spectrum is converted back to the time domain signal, and the time-frequency spectrum is reconstructed using the inverse wavelet transform to obtain the time domain representation of the seismic signal.

[0063] B4. Output new seismic signal.

[0064] S3. Analyze the signal-to-noise ratio of seismic data from all work areas in the study area, such as Figure 6As shown; then the dip-guided structural filtering technology is used to reduce the noise of seismic data in all work areas, and the diffusion coefficient in the PM diffusion equation is used to control the degree of filtering by the dip-guided structural filtering technology, so that the diffusion speed at the edge of the image is slower to retain more details at the edge of the image, as shown Figure 7 As shown in the figure, while protecting the effective information of seismic data, noise reduction is performed to improve the signal-to-noise ratio.

[0065] The dip-guided structural filtering technology in this embodiment uses the dip and azimuth of the actual stratum for constraints, and controls the response of the filter in different dip ranges and different azimuth ranges. The dip constraint can highlight or suppress the stratum characteristics at a specific dip, reduce noise and interference, and improve the interpretation ability and clarity of seismic images. The azimuth constraint further constrains the directional characteristics of a specific stratum, realizes the enhancement or suppression of the stratum characteristics in a specific direction, helps to extract structural information with a specific direction, and helps to interpret the geometric morphology and directional characteristics of underground structures.

[0066] S4. Energy consistency processing is performed from two perspectives: energy balance between multi-working areas and energy balance within single-working areas.

[0067] Energy balance of multiple work areas: Select a standard work area, and calculate the consistency coefficient of seismic data of each work area adjacent to the standard work area and the standard work area. Then, gradually expand outward from the standard work area to all work areas, and calculate the consistency coefficient of seismic data between every two adjacent work areas. For each work area except the standard work area, different weights are assigned according to the corresponding consistency coefficient of seismic data, and energy balance of multiple work areas is performed.

[0068] Specifically, the consistency coefficient of seismic data between each two adjacent work areas is obtained by calculating the amplitude envelope of the seismic traces of the two work areas respectively, aligning the amplitude envelopes extracted from the two work areas on the time axis, so that the time points corresponding to the envelope data of the same event in different work areas are consistent, and then dividing the envelope values ​​at the same time point at the overlapping position in the adjacent work areas to obtain the consistency coefficient between the two adjacent work areas. The division of the envelope values ​​can obtain the amplitude difference of the seismic waveform, more accurately reflect the consistency of the seismic waveform, and reduce the influence of high-frequency noise in the seismic data.

[0069] In this embodiment, the calculation of the amplitude envelope of the seismic trace adopts the Hilbert transform method. The Hilbert transform is a complex domain transform that combines the real and imaginary parts of the signal into a complex signal and extracts its envelope. The Hilbert transform can effectively extract high-frequency and low-frequency amplitude changes by processing the signal in the full frequency range, which helps to capture subtle amplitude changes.

[0070] Before energy balancing in a single work area, quality control of the energy balancing effect of multiple work areas must be performed to determine whether the energy balancing of multiple work areas is effective in different work areas. The quality control means are: given a threshold value, it is determined whether the energy difference of the overlapping area of ​​two adjacent work areas after energy balancing of multiple work areas is less than the threshold value; when the energy difference of the overlapping area of ​​a work area and all its adjacent work areas is less than the threshold value, the energy balancing of multiple work areas is effective in the work area, and the work area is subjected to energy balancing in a single work area; when the energy difference of the overlapping area of ​​a work area and a work area adjacent to it is greater than or equal to the threshold value, the work area is subjected to energy balancing in a single work area, or the amplitude consistency processing of three-dimensional seismic data is abandoned for the work area, or the parameters and data of the work area are optimized until the energy balancing of multiple work areas is effective in the work area, and then the work area is subjected to energy balancing in a single work area. In this embodiment, the threshold value used for quality control of energy balancing in multiple work areas is 5%.

[0071] Energy balance in a simplex zone: Given a time window, slide the time window point by point, normalize the amplitude within each time window, and normalize the energy of different time windows by scaling the energy at different time window positions to achieve energy balance in a simplex zone. Figure 8 As shown in b, c, d, and e in FIG, this embodiment uses time windows of different sizes to perform energy balancing in a simplex area. It is finally determined that in the area where this embodiment is applied, the time window that can be used is 500ms-1000ms. If the time window is less than 500ms, the amplitude preservation becomes poor, and if the time window is greater than 1000ms, the energy balancing effect is too poor and it is difficult to meet the energy consistency requirement. The time window used for the best energy balancing effect in a simplex area is 700ms.

[0072] The method for normalization within each time window is to calculate the absolute average value, median value or root mean square value of the amplitude within the time window as the amplitude reference value, and use the inverse of the amplitude reference value as the weight coefficient of the time window to perform adaptive gain control (AGC control) on the amplitude within the time window. Among them, the absolute average value helps to reflect the overall amplitude level of the signal, and has relatively little impact on signal spikes and outliers. It is suitable for situations where the signal fluctuates smoothly and there are no obvious outliers. The median value can better resist the interference of noise and pulses. When there are obvious outliers in the signal, the median value can be used as the amplitude reference value to improve the stability and reliability of the algorithm. The root mean square has a more balanced response to high and low frequencies. When the signal needs to comprehensively consider the signal energy distribution, the root mean square can be used as the amplitude reference value to provide more comprehensive information.

[0073] S5, output phase, frequency, noise reduction, energy comprehensive processing results, such as Fig.10 As shown. Fig. 9 The original seismic data volume not processed by this embodiment is Fig.10In comparison, it can be clearly seen that the present embodiment effectively improves the data difference between adjacent blocks.

[0074] The present invention also provides a device for processing amplitude consistency of three-dimensional seismic data in multiple desert working areas, which includes a storage medium on which a computer program is stored, and a processor for implementing the above-mentioned method for processing amplitude consistency of three-dimensional seismic data in multiple desert working areas when executing the computer program.

[0075] In summary, this embodiment fully considers the impact of surface conditions, acquisition methods, processing methods and acquisition years on data, and makes up for the problem that the amplitude preservation and consistency of the contiguous data after processing by the traditional post-stack consistency contiguous processing method cannot be taken into account at the same time, and realizes a more comprehensive consideration of the phase, frequency, signal-to-noise ratio and energy of the three-dimensional seismic data, ensuring the relative relationship of the effective information of the original data in the time window, especially the relative relationship of the effective information in the 700ms time window, helping researchers to more effectively improve the data differences between adjacent blocks, providing a better data basis for the work of contiguous layers, fault interpretation, reservoir description, reservoir prediction, etc., and also plays an important data support role in the splicing of three-dimensional seismic data.

[0076] The above description is only the preferred embodiment of the present invention and is not intended to be exclusive. The present invention may be modified and varied in various ways. Any limitation to the present invention made within the spirit and principle of the present invention, modification, equivalent substitution, improvement, etc. by a person skilled in the art shall be included in the protection scope of the present invention.

Claims

1. A method for processing amplitude consistency of 3D seismic data in multiple desert areas, characterized by: The steps include: S1. Based on the calibration of synthetic seismic records of well logging in all work areas, phase consistency processing is performed to make the seismic polarity of the marker layer consistent and form a new data body; S2. Perform spectrum statistical analysis on the new data body to determine the target work area and the spectrum characteristic distribution of all work areas; perform frequency unification processing on other work areas from the perspectives of frequency waveform and energy based on the target work area; S3. Use the dip-guided structural filtering technology to reduce the noise of seismic data in all work areas. Use the diffusion coefficient in the PM diffusion equation to control the degree of filtering by the dip-guided structural filtering technology, slow down the diffusion speed at the edge of the image, and retain the details at the edge of the image; S4, energy consistency processing is performed from two perspectives: energy balance between multiple working areas and energy balance within simple working areas; Select a standard work area, calculate the consistency coefficient of seismic data between each work area adjacent to the standard work area and the standard work area respectively; then gradually expand outward to all work areas with the standard work area as the center, calculate the consistency coefficient of seismic data between every two adjacent work areas; for each work area except the standard work area, assign different weights according to the corresponding consistency coefficient of seismic data, and perform energy balance among multiple work areas; Given a time window, the time window is slid point by point, the amplitude is normalized within each time window, and the energy at different time window positions is normalized by scaling to achieve energy balance within a single working area; S5, the result of comprehensive processing of output phase, frequency, noise reduction and energy.

2. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: The frequency uniformization process in S2 includes: B1. Taking the target work area as the center, select the superimposed seismic traces in the target work area, select the overlapping seismic traces in other work areas, and gradually expand outward to perform contiguous sub-wave shaping on the seismic data of other work areas based on amplitude consistency correction and the reflection coefficient consistency principle; B2. Performing time-frequency spectrum decomposition on the reflection coefficient and seismic signal of the logging data simultaneously to obtain a reflection coefficient spectrum and a seismic spectrum; B3. Gradually approach the seismic spectrum to the reflection coefficient spectrum to obtain a spectrum blueing operator that matches the two; apply the spectrum blueing operator to all seismic spectra along the frequency direction to synthesize a new seismic spectrum; convert the synthesized seismic spectrum back to the time domain signal, perform time-frequency spectrum reconstruction, and obtain the time domain representation of the seismic signal; B4. Output new seismic signal.

3. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 2, characterized in that: B2 also includes taking the seismic spectrum of the target work area as a standard spectrum and using the standard spectrum to perform spectrum equalization on seismic spectra of other work areas.

4. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 2, characterized in that: In B2, time-frequency spectrum decomposition is performed by continuous wavelet transform.

5. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 2, characterized in that: In B2, the reflection coefficient spectrum and seismic spectrum are also processed by spectrum blueing.

6. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 2, characterized in that: In B3, inverse wavelet transform is used for time-frequency spectrum reconstruction.

7. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: After S2 frequency consistency processing, the main frequency difference between other work areas and the target work area is less than 3hz.

8. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: The dip-guided structural filtering technology in S3 uses the dip and azimuth of the actual stratum as constraints to control the response of the filter in different dip ranges and different azimuth ranges.

9. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: In S4, the consistency coefficient of seismic data between every two adjacent work areas is obtained by respectively calculating the amplitude envelopes of the seismic traces of the two work areas, aligning the amplitude envelopes extracted from the two work areas on the time axis, so that the time points corresponding to the envelope data of the same event in different work areas are consistent, and then dividing the envelope values ​​at the same time points at the overlapping positions in the adjacent work areas to obtain the consistency coefficient between the two adjacent work areas.

10. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 9, characterized in that: The calculation of the seismic trace amplitude envelope adopts the Hilbert transform method.

11. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: In S4, the method for performing normalization processing within each time window is to calculate the absolute average value, median value or root mean square value of the amplitude within the time window as the amplitude reference value, and use the reciprocal of the amplitude reference value as the weight coefficient of the time window to adjust the amplitude within the time window.

12. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: S4 also includes quality control of energy balance in multiple working intervals: given a threshold value, it is determined whether the energy difference of the overlapping areas of two adjacent working intervals after energy balance in multiple working intervals is less than the threshold value; when the energy differences of the overlapping areas of a working area and all its adjacent working areas are less than the threshold value, the energy balance in multiple working intervals is effective in the working area, and energy balance processing in a single working area is performed on the working area; when the energy difference of the overlapping area of ​​a working area and a certain adjacent working area is greater than or equal to the threshold value, energy balance in multiple working intervals is not effective in the working area, and energy balance processing in a single working area is performed on the working area, or amplitude consistency processing of three-dimensional seismic data is abandoned for the working area, or parameters and data of the working area are optimized until energy balance in multiple working intervals is effective in the working area, and then energy balance processing in a single working area is performed on the working area.

13. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 12, characterized in that: The threshold value used for energy balance quality control in multiplexed intervals in S4 is 5%.

14. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: Phase consistency processing includes the following steps: A1. All independent tilt measurement networks are transformed into a unified large measurement network through network rotation transformation; A2. According to the changes in the oil group interface and formation thickness, the reflection characteristics of the marker layer are determined based on the well point calibration; A3. Based on the same seismic phase tracking principle, the consistency of time thickness and depth thickness in the vertical direction of the well point is achieved; A4. Through phase migration, the original seismic volume data is spatially resampled and interpolated to form new seismic data.

15. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 14, characterized in that: The similarity coefficient between the new seismic data well bypass trace and the well point synthetic trace in A4 is greater than 98%.

16. The method for processing amplitude consistency of desert multi-area 3D seismic data according to claim 1, characterized in that: The time window size in S4 is 500ms-1000ms.

17. Amplitude consistency processing device for 3D seismic data in multiple desert areas, characterized by: It includes a storage medium on which a computer program is stored; and it also includes a processor, which is used to implement the method for processing amplitude consistency of desert multi-area three-dimensional seismic data as claimed in any one of claims 1 to 16 when executing the computer program.

Images