Dual-threaded seismic acquisition data weighted surface element analysis method

Abstract

This invention discloses a dual-threaded weighted pixel analysis method for seismic acquisition data, comprising two aspects: a method for calculating the effective coverage count and a programming implementation method. The method for calculating the effective coverage count first calculates the maximum and minimum shot-receiver spacing based on the observation system parameters. Using the average of the maximum and minimum shot-receiver spacings as a base (or an offset base specified interactively), the effective coverage count of the pixels is calculated according to a weighted algorithm. The programming implementation method involves a dual-threaded approach to achieve comprehensive pixel analysis of massive seismic acquisition data. One thread reads the shot-receiver point and relationship data; the other thread calculates the conventional coverage count, shot-receiver spacing, azimuth angle, and weighted coverage count. This method effectively extends conventional pixel analysis methods. This invention aims to achieve comprehensive pixel analysis of massive seismic acquisition data by focusing on the effective coverage count and utilizing an efficient and practical weighted pixel analysis algorithm.

Description

Technical Field

[0001] This invention relates to the field of seismic data acquisition technology, specifically to a weighted surface element analysis method for dual-thread seismic acquisition data. Background Technology

[0002] With the continuous upgrading of seismic acquisition hardware, the number of receiver channels for seismic data acquisition is increasing. Single shot gathers exceed ten thousand channels, shot spacing is becoming larger, and coverage times are increasing. As shot spacing increases, the contribution of data acquired through vibration to the overlay data of seismic elements decreases. Conventional element attribute calculations do not distinguish between shot spacing; a single shot-receiver pair is used to calculate the coverage times. While seismic data processing includes modules that weight seismic data based on shot spacing, this technology is not yet applied in the element attribute analysis software of the acquisition and observation system. Summary of the Invention

[0003] To address the shortcomings of existing technologies, the present invention aims to provide a weighted surface element analysis method for dual-threaded seismic acquisition data. This method calculates the weighted coverage times of the shot spacing while simultaneously calculating the conventional surface element attributes, thereby obtaining comprehensive surface element analysis data from the seismic acquisition data for subsequent comparative analysis of surface element attributes.

[0004] To achieve the above objectives, this invention employs the following technical solution: a dual-threaded seismic acquisition data weighted surface element analysis method, which expands the method for calculating seismic acquisition surface elements. In its programming implementation, it adopts a dual-threaded method driven by observation system relation patch data, comprising the following steps:

[0005] S01: Read data from the seismic acquisition and observation system, obtain the observation system parameters, and establish the correspondence between shot receiver points;

[0006] S02: Obtain the maximum and minimum gun spacing;

[0007] S03: Set a reasonable weighting base for shot spacing. Two methods are designed: one is to calculate a default value P0 = (Offset_max + Offset_min) / 2 based on the obtained maximum and minimum shot spacing; the other is to interactively input the weighting base. The weighting base for shot spacing is obtained through these two methods and used to calculate the number of effective coverages.

[0008] S04: Start the first thread to read a relational slice of data from the seismic acquisition and observation system data file; set the valid data signal or the data reading completion signal, and then wait for the calculation results of the second thread;

[0009] S05: Start the second thread. After receiving the data signal from the first thread, initialize the surface point data and start the calculation.

[0010] S06: Calculate the surface attribute data of the surface element points, including: shot spacing, azimuth angle, coverage count, and effective coverage count;

[0011] The calculation method for effective coverage count is as follows: When the current gun spacing Offset_current is less than or equal to the weighted base, the coverage count is 1; when the current gun spacing is greater than the weighted base, the attenuation coefficient of the coverage count is first calculated using the formula: C1 = (Offset_current - P0) / (Offset_max - P0). The value of C1 is compared, with a maximum of 0.9. If it is greater than 0.9, it is set to 0.9. Effective coverage count: F p =1–C1.

[0012] S07: To ensure program efficiency, when the element attribute data reaches 8Mb, save the current data volume to a binary file and clear the cache;

[0013] S08: Loop calculation. When all the seismic acquisition data has been read, the first thread outputs an end signal, the second thread saves the current data, and the program ends.

[0014] Preferably, in step S01, the data files of the seismic acquisition and observation system are browsed through to establish the shot-receiver relationship.

[0015] Preferably, in step S02, the global maximum and minimum shot spacing are obtained based on the shot-receiver relationship corresponding to the observation system.

[0016] Preferably, in step S03, a reasonable weighting base for the gun spacing is set, and the calculated weighting base for the gun spacing is reasonably adjusted based on the maximum and minimum gun spacing.

[0017] Preferably, step S03, which involves calculating the effective coverage count, calculates both the effective coverage count of a surface element and the regular surface element attributes.

[0018] Preferably, the dual-threaded data-driven programming implementation method in steps S04 and S05 can accurately and efficiently calculate the regular surface attributes and effective coverage times of surface points.

[0019] The present invention has the following beneficial effects:

[0020] 1. Based on the calculation of conventional surface element attributes (shot spacing, azimuth, coverage times), this invention implements an analysis method for effective coverage times, which can be compared and analyzed with conventional coverage times to analyze the strength of coverage times from the perspective of actual reception effect;

[0021] 2. This invention employs a dual-threaded coding implementation to improve program execution efficiency;

[0022] 3. This invention is designed to save the current data volume to a binary file when the calculated element attribute data reaches 8Mb, thereby reducing I / O operations and improving program running efficiency. Attached Figure Description

[0023] The present invention will now be described in detail with reference to the accompanying drawings and specific embodiments;

[0024] Figure 1 This invention provides a flowchart for reading SPS files of seismic acquisition data to obtain the maximum and minimum shot spacing.

[0025] Figure 2 This is a flowchart of the weighted surface element analysis of dual-thread seismic acquisition data according to the present invention;

[0026] Figure 3 This serves as a template for the relationship slices of the observation system in this invention;

[0027] Figure 4 This is a schematic diagram showing the weighted base of the gun spacing, the current gun spacing, and the maximum gun spacing of the present invention;

[0028] Figure 5 This is a schematic diagram of the gun spacing and azimuth angle of a surface element attribute of the present invention. Detailed Implementation

[0029] To make the technical means, creative features, objectives and effects of this invention easier to understand, the invention will be further described below in conjunction with specific embodiments.

[0030] Reference Figure 1 The specific implementation method adopts the following technical solution: The main purpose of the dual-thread seismic acquisition data weighted surface element analysis method is to increase the means of surface element analysis of the observation system through weighted surface element analysis, so as to reflect the influence of vibration on the superimposed surface element data with distance; in terms of programming implementation method, the computational efficiency is improved by driving the relationship slice data of the dual-thread observation system.

[0031] It includes the following steps:

[0032] S01 reads seismic acquisition data:

[0033] At the start of the program, the initial maximum and minimum shot spacing are set to 0. The seismic acquisition observation system data (SPS file) is read to obtain the observation system parameters and establish the correspondence between shot and receiver points. During the file reading process, the offset of the shot and receiver pair is calculated.

[0034] like Figure 1 As shown, three files in SPS format from the observation system (sensor point, shot-sensor, and shot-sensor relationship) are read respectively, and the shot-shot distance is calculated based on the correspondence between the shot point and the sensor point.

[0035] S02 calculates the maximum and minimum gun spacing:

[0036] During the process of reading the file by S01, the maximum and minimum gun spacings are reset by comparing the calculated gun spacings with the current maximum and minimum gun spacings. After the file reading is completed, the maximum (Offset_max) and minimum gun spacing (Offset_min) are obtained.

[0037] S03 calculates the weighted base for gun spacing:

[0038] Setting a reasonable weighting base for shot spacing is achieved through two methods. One method is to calculate a default value P0 = (Offset_max + Offset_min) / 2 based on the obtained maximum and minimum shot spacing. The other method is to interactively input the weighting base. The weighting base for shot spacing is obtained through these two methods. In the programming implementation, the interactive input is to modify and confirm the method 1 according to the actual situation, and the confirmed weighting base is used to calculate the number of effective coverages.

[0039] S04's first thread reads the observation system relation slice data:

[0040] Start the first thread to restart data reading. Read a data patch of observation system relationships (including shot points, receiver points, and relationships) from the seismic acquisition observation system data file and put it into the corresponding buffer; if the read is successful, set the observation system relationship patch data driver signal to 1 (valid data signal); if the file read is complete, set the observation system relationship patch data driver signal to 2 (file read complete signal); if the read fails, set the observation system relationship patch data driver signal to 0 (read failure signal).

[0041] Figure 3 This is a template for the relationship slice of the observation system of this invention.

[0042] The second thread in S05 retrieves observation system relation slice data from the cache:

[0043] Start the second thread and wait for the data drive signal from the first thread. Upon receiving the first thread's data signal 1 (valid data signal), initialize the surface point data, retrieve the observation system relation slice data information from the cache, and begin calculation.

[0044] S06 Calculates surface element attribute data:

[0045] Calculate the surface element attribute data of the surface element points, including: shot spacing, azimuth angle, coverage count, and effective coverage count;

[0046] Calculation methods: shot spacing (distance from the shot point to the receiver point), azimuth (angle between the shot point and the receiver point), and number of coverages (one shot-receiver pair accumulates one coverage).

[0047] like Figure 4 The diagram shown illustrates the weighted base of the gun spacing, the current gun spacing, and the maximum gun spacing of this invention.

[0048] The calculation method for effective coverage count is as follows: When the current gun spacing Offset_current is less than or equal to the weighted base, the coverage count is 1; when the current gun spacing is greater than the weighted base, the attenuation coefficient of the coverage count is first calculated using the formula: C1 = (Offset_current - P0) / (Offset_max - P0). The value of C1 is compared, with a maximum of 0.9. If it is greater than 0.9, it is set to 0.9. Effective coverage count: F p =1–C1.

[0049] Figure 5 This is a schematic diagram showing the calculated offset distance and azimuth angle of a surface element.

[0050] S07 outputs the calculation results in a centralized manner:

[0051] To ensure program efficiency, when the element attribute data reaches 8Mb, the current data volume is saved to a binary file and the cache is cleared.

[0052] S08 Loop Calculation and Saving:

[0053] The calculation is performed in a loop. When all the seismic acquisition data has been read, the first thread outputs an end signal, the second thread saves the current data, and the program ends. The saved binary file contains the shot spacing, azimuth, coverage count, effective coverage count, and corresponding index information of the surface elements.

[0054] The method of this invention achieves effective coverage times F by setting a weighted base for shot spacing. p By employing a dual-threaded, cached approach, the system fully leverages the multi-core CPU and the efficient I / O capabilities of memory cache to achieve rapid computation. The program design utilizes a data-driven approach based on observation system relation slices to enable efficient thread coordination.

[0055] The foregoing has shown and described the basic principles, main features, and advantages of the present invention. Those skilled in the art should understand that the present invention is not limited to the above embodiments. The embodiments and descriptions in the specification are merely illustrative of the principles of the invention. Various changes and modifications can be made to the invention without departing from its spirit and scope, and all such changes and modifications fall within the scope of the present invention as claimed. The scope of protection of this invention is defined by the appended claims and their equivalents.

Claims

1. A weighted surface element analysis method for dual-thread seismic acquisition data, characterized in that, The method for calculating seismic acquisition elements has been expanded. In terms of programming implementation, a dual-threaded approach driven by observation system relation patch data is adopted, including the following steps: (S01): Read the data from the seismic acquisition and observation system, obtain the observation system parameters, and establish the correspondence between the shot receivers; (S02): Obtain the maximum and minimum gun spacing; (S03): Set a reasonable weighting base for the gun spacing. Two methods are designed: one is to calculate a default value P0 = (Offset_max + Offset_min) / 2 based on the obtained maximum and minimum gun spacing; the other is to interactively input the weighting base. The weighting base for the gun spacing is obtained through these two methods and used to calculate the number of effective coverages. (S04): Start the first thread to read a relational data patch from the seismic acquisition observation system data file; set the valid data signal or the data reading completion signal, and then wait for the calculation result of the second thread; (S05): Start the second thread. After receiving the data signal from the first thread, initialize the surface point data and start the calculation. (S06): Calculate the surface attribute data of the surface element points, including: shot spacing, azimuth angle, number of coverages, and number of effective coverages; The calculation method for effective coverage count is as follows: When the current gun spacing Offset_current is less than or equal to the weighted base, the coverage count is 1; when the current gun spacing is greater than the weighted base, the attenuation coefficient of the coverage count is first calculated using the formula: C1 = (Offset_current - P0) / (Offset_max - P0). The value of C1 is compared, with a maximum of 0.

9. If it is greater than 0.9, it is set to 0.

9. Effective coverage count: F p =1–C1; (S07): To ensure the efficiency of the program, when the element attribute data reaches 8Mb, save the current data volume to a binary file and clear the cache; (S08): Loop calculation. When all the seismic acquisition data has been read, the first thread outputs an end signal, the second thread saves the current data, and the program ends.

2. The weighted surface element analysis method for dual-thread seismic acquisition data according to claim 1, characterized in that, In step (S01), the data files of the seismic acquisition and observation system are browsed through to establish the shot-receiver relationship.

3. The weighted surface element analysis method for dual-thread seismic acquisition data according to claim 1, characterized in that, In step (S02), the global maximum and minimum shot spacing are obtained based on the shot-receiver relationship corresponding to the observation system.

4. The weighted surface element analysis method for dual-thread seismic acquisition data according to claim 1, characterized in that, In step (S03), a reasonable weighting base for the gun spacing is set, and the calculated weighting base for the gun spacing is adjusted reasonably based on the maximum and minimum gun spacing.

5. The weighted surface element analysis method for dual-thread seismic acquisition data according to claim 1, characterized in that, The method for calculating the effective coverage count in step (S03) is used to calculate the effective coverage count of surface elements and the regular surface element attributes.

Images