Explosive source node data synthesis and quality control method and system based on shared memory

Abstract

The invention provides an explosive source node data synthesis and quality control method and system based on a shared memory. The method comprises the following steps: step 1, constructing a node data information table; step 2, constructing an observation system data table; step 3, creating a shared memory space; 4, cutting node data, and synthesizing a detection line record; 5, calculating to obtain a detection line record quality analysis result according to the detection line record; step 6, according to the detection line record, performing batch merging according to a field file number sequence to obtain all single shot records; and step 7, according to the synthesized single-shot record, calculating to obtain a single-shot quality analysis result. According to the explosive source node data synthesis and quality control method and system based on the shared memory, the problems that single-shot records are directly cut and synthesized one by one and quality control cannot be carried out in time in the past are solved, the cutting synthesis efficiency of the explosive source node acquisition data is greatly improved, and the quality of field seismic acquisition data can be monitored in time.

Description

Technical Field

[0001] This invention relates to the field of seismic acquisition node data processing technology, and in particular to a method and system for synthesizing and quality controlling explosive source node data based on shared memory. Background Technology

[0002] In my country's eastern exploration areas, with the continuous deepening of oil and gas exploration, complex surfaces, complex structures, and ultra-deep fine exploration have become the main battleground for oil and gas exploration, and its development increasingly relies on high-density exploration. High-density acquisition schemes are characterized by small area cells, wide azimuth, high density, and long arrangement, and field operations also require flexible observation system adaptability. However, due to the limitations of traditional seismograph track capacity and cable flexibility, cabled acquisition methods cannot meet the technical requirements of ultra-high shot density seismic acquisition. Nodal acquisition, on the other hand, allows nodal instruments to acquire data independently without cable or time constraints, possessing unlimited track capacity and flexible observation system adaptability, thus meeting the technical requirements of high-density seismic acquisition. However, nodal acquisition has the drawback of not being able to monitor the quality of individual seismic shots in real time. It requires waiting for all nodal data to be collected, segmented, and synthesized into complete individual shot data before copying it to a seismic data processing system (such as ProMAX, Omega, etc.) for manual quality inspection. This delay in the individual shot quality control process poses a serious risk to the seismic acquisition project.

[0003] Taking a high-density seismic acquisition project in an exploration area in eastern my country as an example, the seismic acquisition in the eastern exploration area of ​​my country is mainly based on explosive sources due to surface conditions. The observation system is a 6S42L336R, with the shot fired in the middle and receivers on both sides. The sampling interval is 1ms, and the seismic record length is 7s. According to the field construction organization mode, an average of 2,000 shots are fired per day, with 665 nodal instruments arranged in a single row. An average of 4 rows are retrieved per day, meaning that after 11 days of retrieval, the single-shot data acquired on the first day can be segmented and synthesized. The daily raw data volume collected from nodes is: 14404 (KB / hour) × 24 (hours) × 11 (days) × 665 (tracks) / 10243 = 2.36TB. The synthesized single-shot data volume is approximately 800 (MB / shot) × 2000 (shots) / 10242 = 1.53TB. Using the traditional method, directly cutting and synthesizing single-shot records takes approximately 39 seconds per shot. The time for cutting and synthesizing 2000 shot data is approximately 39 (seconds / shot) × 2000 (shots) = 21.7 hours. This means that for the explosive source-triggered node acquisition project in the eastern exploration area, it will take at least 12 days before the single-shot data collected on the first day can be seen. Therefore, it is impossible to check the quality of the data collected on the same day in a timely manner, posing a serious hidden danger to the seismic acquisition project. If quality problems are found in the previously collected data after the 12th day, the data needs to be rearranged, and new wells and blasts need to be drilled, affecting the construction period and causing incalculable economic losses.

[0004] Chinese patent application CN202210709041.2 discloses a method for stepwise synthesis of nodal instrument shot gather data, relating to the field of seismic exploration technology. The invention includes the following steps: collecting source point shot-receiver relationship files for synthesizing single-shot data to form a general source point shot-receiver relationship file; determining the start and end range of the general source point shot-receiver relationship file's file numbers based on the single-shot data file numbers corresponding to each source point in the general source point shot-receiver relationship file; and identifying the minimum and maximum station numbers among all receiver alignments based on the general source point shot-receiver relationship file to determine the alignment line station number range.

[0005] Chinese patent application CN202011170606.1 discloses a multi-host joint data acquisition system, comprising: an excitation host for receiving the working status of multiple microseismic acquisition hosts via a short message system; determining whether construction is permitted, and if so, sending construction instructions to the multiple microseismic acquisition hosts via the short message system; receiving node quality control data; acquiring microseismic data or single-shot data after construction begins; and sending the node quality control data, excitation result data and synthesis relation file, microseismic data, or single-shot data after construction ends; a microseismic acquisition host for sending the working status to the excitation host via the short message system; starting construction and acquiring microseismic data after receiving the construction instructions; and sending the microseismic data after construction ends; and a data processing unit for performing data separation and synthesis on the node quality control data, single-shot data, or microseismic data based on the excitation result data and synthesis relation file.

[0006] Chinese patent application CN202011473037.8 discloses a single-shot quality monitoring method for a wireless nodal instrument seismic data acquisition system, comprising the following steps: 1. The wireless nodal instrument stores data quality information during the acquisition process in the form of time nodes on an SD card; 2. The shotlog file is read to obtain the GPS time of the excitation pulse, and this is used as the zero time T0 of the data recorded by the wired device; 3. Each trace head of the Seg-D recording file is traversed; 4. The qclog file of the corresponding nodal station is found through the deployment results and observation system files; 5. The contents of the qclog file are scanned to find the most recent nodal station status information and daily inspection information before the T0 data, and each indicator is judged to be normal and recorded; 6. The synthesis is completed and the single-shot quality monitoring file is saved, otherwise, proceed to step 3.

[0007] The above-mentioned existing technologies are all quite different from the present invention and have failed to solve the technical problem we want to solve. Therefore, we have invented a new method and system for the synthesis and quality control of explosive source node data based on shared memory. Summary of the Invention

[0008] The purpose of this invention is to provide a method and system for explosive source node data synthesis and quality control based on shared memory, which solves the problems of direct cutting and synthesis of single-shot records without timely quality control in the past, greatly improves the efficiency of cutting and synthesizing explosive source node acquisition data, and can monitor the quality of field seismic acquisition data in a timely manner.

[0009] The objective of this invention can be achieved through the following technical measures: a method for synthesizing and quality controlling explosive source node data based on shared memory, which includes:

[0010] Step 1: Construct the node data information table;

[0011] Step 2: Construct the observation system data table;

[0012] Step 3, create a shared memory space;

[0013] Step 4: Cut the node data and synthesize the detector line record;

[0014] Step 5: Calculate the quality analysis results of the detector line records based on the detector line records;

[0015] Step 6: Based on the detector line records, merge all single-shot records in batches according to the field file number order;

[0016] Step 7: Calculate the single-shot quality analysis results based on the synthesized single-shot records.

[0017] The objective of this invention can also be achieved through the following technical measures:

[0018] In step 1, input the basic parameters of the work area, including sampling interval and seismic record length information; load the nodal instrument replacement files corresponding to the single shot collected on the day, and establish a nodal instrument replacement data table. The nodal instrument replacement data table includes: detector line number, detector station number, starting file number, equipment number before replacement, and equipment number after replacement; scan all the raw data collected by nodal instruments that were recovered on the day, and create a nodal data information table. The nodal data information table includes: detector line number, detector station number, nodal instrument equipment number, number of recovered nodal data files, number of abnormal nodal data files, nodal data storage path, and the starting field file number for replacing the nodal instrument.

[0019] Step 1, when creating the node data information table, includes the following steps:

[0020] Scan all node data downloaded that day to obtain the device number, detector line number, detector stake number and storage address of each node instrument;

[0021] Statistical analysis of the number of node data files collected by each node instrument;

[0022] Based on the sampling interval of the work area, the number of abnormal files is statistically analyzed, and the abnormal file paths and sizes are given. These data are not used in subsequent node data cutting and synthesis.

[0023] According to the node instrument replacement data sheet, mark the starting field file number for replacing the node instrument;

[0024] Construct a node data information table.

[0025] In step 2, load the shot-receiver relationship file and shot-point excitation file corresponding to the single shot collected on that day. The shot-receiver relationship file includes: field file number (ffid), shot point line number, shot point station number, receiver line number, receiver starting station number, and receiver ending station number. The shot-point excitation file includes: field file number (ffid) and GPS time. Based on the shot-receiver relationship file and shot-point excitation file, construct the observation system data table. The observation system data table includes: field file number, shot point line number, shot point station number, excitation time, number of receiver permutations, starting permutation number, ending permutation number, total number of receiver channels, number of missing node channels, data integrity, and node data retrieval rate.

[0026] Step 2, when constructing the observation system data table, includes the following steps:

[0027] Obtain the field file number, shot point line number, and shot point stake number from the shot checkpoint relationship file;

[0028] Obtain the field file number and firing time from the firing point firing file, match the field file number in the firing point firing file with the field file number in the firing point checkpoint relationship file, and establish a one-to-one correspondence between firing time and firing point line number and firing point station number;

[0029] The starting and ending receiver sequence numbers of each shot point in the shot-receiver relationship file, as well as the starting and ending receiver station numbers on each receiver sequence, are statistically analyzed to calculate the total number of receiver channels for each shot point, denoted as allTrace.

[0030] Match the field file number, receiver line number, receiver start station number, and receiver end station number in the shot-receiver relationship file with the receiver line number and receiver station number in the node data information table created in step 1. Count the missing receiver station numbers and calculate the number of missing node traces and node data recovery rate for each shot point. The number of missing node traces is denoted as noTrace, and the node data recovery rate is denoted as Q. The formula for calculating the node data recovery rate is:

[0031] Q=(allTrace-noTrace) / allTrace,

[0032] To ensure the quality of data collected in the field, the node data retrieval rate is required to be greater than 99%.

[0033] The system provides a judgment result on whether the node data is complete. The method for judging the integrity of node data is as follows: when the node data recovery rate is not equal to 100%, the node data of a single shot is determined to be incomplete. The system also compiles and outputs a list of missing geophone station numbers on each geophone line recorded by a single shot, so that construction personnel can locate the lost node instrument equipment in the field based on the geophone station numbers.

[0034] Step 3 includes:

[0035] Step 31: Obtain the server memory capacity, denoted as RAM. To fully utilize server memory resources, improve the computational efficiency of cutting and synthesizing detector line records, and ensure operational stability, the default maximum allowed memory usage rate is 80%.

[0036] Step 32: Select the detector lines to be cut and synthesized, and count the number of detector lines to be cut and synthesized, as well as the number of nodal devices on the longest detector line. The specific method is as follows:

[0037] Select the detector line number to be cut and synthesized, and denote the number of detector lines to be cut and synthesized as LN;

[0038] Obtain the starting and ending station numbers of each geophone line, calculate the required number of geophones on the line, and denot it as AN[LN], where AN[LN] = {an1, an2, ..., an...} LN}, where an i Let i be the number of detector points on the i-th detector line, where i = 1, 2, ..., LN;

[0039] Statistics i The maximum value of (i = 1, 2, ..., LN) is denoted as an. max ;

[0040] Step 33: Obtain the sampling interval of the work area, denoted as Dt, and automatically calculate the data file size of a single node, denoted as Bin;

[0041] Step 34: Select the field file number to be cut and synthesized, and obtain the maximum number of blasts per hour, specifically:

[0042] The field file numbers are grouped according to the GPS time of the firing point. Single shots fired within an hour are grouped together, and the number of shots fired in each hour is counted.

[0043] The maximum number of explosions per hour is denoted as sn. max ;

[0044] Step 35, calculate the number of sampling points for a single shot record, specifically:

[0045] Obtain the sampling interval and seismic record length for the work area. The sampling interval is denoted as Dt and the seismic record length is denoted as SeisLen. Calculate the number of sampling points per shot record, denoted as Samp. The calculation method is: Samp = SeisLen ÷ Dt.

[0046] Step 36: Calculate the maximum number of detector lines that the server can cut simultaneously, that is, calculate the maximum number of parallel threads that can be created. Specifically:

[0047] First, assume that the server can cut the maximum number of synthesized detector lines at the same time as CL, that is, the maximum number of parallel threads that can be opened is CL;

[0048] Let N1 be the memory space required for the original node data when cutting and synthesizing the longest detector line. The calculation method is: N1 = an max ×Bin×2, where 2 represents creating two shared memory spaces for each detector station number, i.e. caching 2 hours of node data;

[0049] The buffer space required for the synthesized detector line record during the hour with the most firings, when cutting and synthesizing the longest detector line, is denoted as N2. The calculation method is as follows:

[0050] N2=an max ×sn max ×samp×4, where 4 indicates that each sample of the detector line data is stored using 4 bytes;

[0051] The maximum number of parallel threads that the server can create is denoted as CL, where CL is a positive integer. The calculation method is as follows:

[0052]

[0053] CL≤LN

[0054] To make the above system of inequalities hold, take the largest positive integer value of CL;

[0055] Step 37: Based on the CL value, create a shared memory space, denoted as shareM, where shareM = CL × (N1 + N2), to maximize the use of server memory resources and improve the efficiency of detector line data cutting and synthesis.

[0056] Step 4 includes:

[0057] Step 41, extract node data into shared memory, specifically:

[0058] Based on the GPS time of single shot firing, extract the node data for 2 hours corresponding to all detector point station numbers on the detector line in sequence, and load them into the two shared memory spaces created in step 3.

[0059] Create two memory pointers: the processing pointer points to the first memory block, and the storage pointer points to the second memory block.

[0060] Step 42: Cut the node data and merge it into detector line records until all detector line records of all single shots fired in the first hour are synthesized.

[0061] Step 43: Repeat steps 41 and 42 until all detector line records acquired that day have been cut and the shared memory space has been released.

[0062] In step 5, based on the cut and synthesized detector line data, calculate the total number of channels, the number of zero-value channels, the number of abnormal channels, the nodal data recovery rate, and the nodal data effectiveness rate for each detector line record, and output the detector line record quality analysis results, specifically including:

[0063] Batch loading of cut and synthesized detector line data;

[0064] To count the total number of channels recorded for each geophone line, specifically: obtain the starting geophone station number and the ending geophone station number for each geophone line, and calculate the total number of channels for the geophone line. The calculation formula is: Total number of channels = Ending geophone station number - Starting geophone station number + 1;

[0065] The calculation of single-channel energy is specifically as follows: calculate the sum of the amplitude values ​​of all sampling points in each seismic trace on the detector line record, denoted as A. k The calculation formula is: Where, k = 1, 2, ..., RN, RN is the total number of traces recorded by the detector line, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points in the seismic record;

[0066] The number of zero-value channels is counted, specifically: the number of channels recorded on detector line A. k The number of seismic traces with a value of 0, k = 1, 2, ..., RN, where RN is the total number of traces recorded by the detector lines;

[0067] The number of abnormal traces is counted. The method for calculating abnormal traces is as follows: Take the current seismic trace, assuming it is the k-th trace, with 10 traces on each side. The energy of a single trace is denoted as {A}. k-10 A k-9 A k-1 A K+1 ,...,A K+！0} Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the energy of the seismic channel was abnormally high. It is assumed that the energy of the seismic traces is abnormally small. The number of abnormal traces is the sum of the number of seismic traces with large single-trace energy anomalies and the number of seismic traces with small single-trace energy anomalies. The threshold value for judging single-trace energy anomalies can be appropriately adjusted according to the actual situation of the work area.

[0068] The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels

[0069] The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels

[0070] Output the quality analysis results of the detector line.

[0071] In step 6, all receiver line records corresponding to each shot point are obtained in order of field file number, arranged in order of receiver line number, and batch-merged to obtain all single-shot records, specifically including:

[0072] Select the field file number to be merged;

[0073] Input the number of concurrent threads for merging single-shot records;

[0074] Based on the observation system data table, the detector line numbers contained in each field file number are obtained in order of field file number, arranged in order of detector line number, and single-shot records are synthesized in batches.

[0075] In step 7, based on the synthesized single-shot records, the single-shot energy, frequency, bandwidth, and signal-to-noise ratio are calculated. The node data recovery rate and node data effectiveness are statistically analyzed, and the single-shot quality analysis results are output, specifically:

[0076] The energy, frequency, bandwidth, and signal-to-noise ratio of a single shot were calculated based on the entire time window using a general geophysical algorithm. Specifically: the energy analysis of a single shot used the root mean square amplitude method; the signal-to-noise ratio analysis used the spectrum estimation method; the dominant frequency was calculated using Fourier transform, converting the single-channel data from the time domain to the frequency domain, and the frequency with the highest amplitude was recorded as the dominant frequency; the bandwidth was calculated as the difference between the maximum and minimum dominant frequencies under 22dB attenuation.

[0077] The calculation of single-channel energy is specifically as follows: calculate the sum of the amplitude values ​​of all sampling points in each seismic trace on a single shot record, denoted as B. k The calculation formula is: Where, k = 1, 2, ..., M, M is the total number of traces recorded by a single shot, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points in the seismic record;

[0078] The number of zero-value channels is counted, specifically: the number of single-shot records at point B. kThe number of seismic traces with a value of 0, k = 1, 2, ..., M, where M is the total number of traces recorded by a single shot;

[0079] The number of abnormal traces is counted. The calculation method for abnormal traces is as follows: Take the current seismic trace, assuming it is the k-th trace, with 10 traces on each side. The energy of a single trace is denoted as {B}. k-10 B k-9 B k-1 B K+1 ,...,B K+！0} Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the energy of the seismic channel was abnormally high. It is assumed that the energy of the seismic traces is abnormally small, and the number of abnormal traces is the sum of the number of seismic traces with abnormally large energy and the number of seismic traces with abnormally small energy. The threshold value for judging energy anomalies can be adjusted appropriately according to the actual situation of the work area.

[0080] The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels;

[0081] The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels;

[0082] Output single-shot recording quality analysis results, specifically including: single-shot energy, main frequency, bandwidth, signal-to-noise ratio, number of zero-value channels, number of abnormal channels, node data recovery rate, and node data effectiveness rate;

[0083] Based on the analysis of data from the test sites in the work area, an evaluation standard for seismic acquisition data was established, specifically as follows:

[0084] Analyze the single-shot energy at the test points in the work area, determine the minimum single-shot energy under different activation well depths and activation charge amounts, and consider the single-shot quality to be unqualified if it is lower than the minimum energy value;

[0085] Analyze the frequency of the test data from the work area and give the main frequency range and bandwidth range. When the main frequency and bandwidth values ​​are abnormal, they can be used as reference factors to determine the quality of a single shot is unqualified, but not as the deciding factor. The single shot record can be opened for manual evaluation.

[0086] Analyze the signal-to-noise ratio of the test data in the work area, and give the minimum signal-to-noise ratio for the work area. If the signal-to-noise ratio is lower than the minimum, the quality of a single shot is considered unqualified.

[0087] Analyze the node data recovery rate; if it falls below the set first threshold, the quality of a single gun is considered unqualified.

[0088] The efficiency of node data analysis is determined. If the efficiency falls below the set second threshold, the quality of a single gun is considered unqualified.

[0089] The objective of this invention can also be achieved through the following technical measures: a shared memory-based explosive source node data synthesis and quality control system, comprising:

[0090] The node data information table construction module constructs the node data information table.

[0091] The observation system data table construction module constructs the observation system data table collected on the same day.

[0092] The shared memory space creation module allocates shared memory space for detector line record cutting and synthesis, improving the efficiency of node data cutting and synthesis.

[0093] The detector line recording cutting and synthesis module cuts node data and synthesizes detector line records.

[0094] The detector line recording quality analysis module analyzes the quality of the synthesized detector line recordings, which is to monitor the quality of the data acquired from the field earthquake on the same day.

[0095] The single-shot record batch merging module combines detector data into single-shot records in batches.

[0096] The single-shot record quality analysis module analyzes the quality of the synthesized single-shot records.

[0097] The objective of this invention can also be achieved through the following technical measures:

[0098] The module for constructing the node's data information table specifically includes:

[0099] The work area information input unit allows users to input the sampling interval and seismic record length information for the construction work area.

[0100] The node instrument replacement file loading unit creates a node instrument replacement data table.

[0101] The node data scanning unit extracts the device number, detector line number, and detector station number of each node instrument, and records the data storage address of each node instrument.

[0102] The node data analysis unit counts the total number of node data files collected that day and analyzes the number of abnormal files based on work area information.

[0103] The node data information table construction unit establishes a one-to-one correspondence between the device number, detector line number, detector stake number and storage address of each node instrument, and marks the starting file number of the replaced node instrument.

[0104] The data table construction module for this observation system includes:

[0105] The blasting checkpoint relationship file loading unit retrieves the blasting checkpoint relationships corresponding to a single blast collected in the construction area on that day.

[0106] The firing point firing file loading unit acquires the firing time of a single shot collected in the construction area on that day;

[0107] The observation system data table construction unit establishes the relationship between the nodal instrument equipment and the observation system, counts the number of missing nodal traces for each shot point, calculates the nodal data recovery rate, and provides the judgment result on whether the nodal data is complete.

[0108] The shared memory space creation module includes:

[0109] The server memory capacity acquisition unit retrieves the current memory capacity of the server.

[0110] The detector line and field file number selection unit allows you to select the detector line that needs to be cut and synthesized, and the field file number that needs to be cut and synthesized on the detector line.

[0111] The node data file size calculation unit calculates the size of a single node data file in the work area.

[0112] Field file number grouping and statistical unit: Sort the field file numbers according to the GPS triggering time, and group them by hour. Count the number of blasts in each hour and obtain the maximum number of blasts in a single hour.

[0113] Maximum number of parallel threads calculation unit, which calculates the maximum number of synthesized detector lines that the server can cut simultaneously;

[0114] Create shared memory units, creating two shared memory spaces for each detector station number. Each shared memory space can cache one hour of node data.

[0115] The detector line recording cut synthesis module includes:

[0116] The node data extraction unit extracts the node data to be cut and synthesized into data blocks in the shared memory space;

[0117] The node data cutting unit cuts the node data according to the shot point excitation GPS time and the number of seismic record sampling points to form single-channel seismic data.

[0118] The receiver line data synthesis unit arranges the segmented single-channel seismic data according to the receiver point station number order and outputs the receiver line record.

[0119] The detector line recording quality analysis module includes:

[0120] Detector line record loading unit, batch loading detector line records to be quality controlled;

[0121] Total Channel Count Unit: Calculates the total number of channels for each detector line;

[0122] A single-channel energy calculation unit calculates the sum of amplitude values ​​at all sampling points for each seismic trace;

[0123] Zero-value trace statistics unit: counts the number of seismic traces with a single trace energy value of 0;

[0124] Abnormal channel statistics unit: calculates and counts the total number of abnormal low-energy channels and abnormal high-energy channels;

[0125] The node data recovery rate calculation unit calculates the node data recovery rate of the detector line records;

[0126] The node data efficiency calculation unit calculates the node data efficiency recorded by the detector line.

[0127] The single-shot record batch merging module includes:

[0128] The field file number selection unit to be synthesized obtains the number of single guns to be synthesized and the field file number;

[0129] The number of parallel threads for merging single-shot records is input to determine the number of parallel threads for batch single-shot record synthesis.

[0130] The automatic batch merging unit for single-shot data obtains the detector line number and detector point number from the detector line record, rearranges them in order, merges them to obtain the single-shot record, and outputs it.

[0131] The single-shot record quality analysis module includes:

[0132] Single-shot record loading unit, batch loading of single-shot records to be quality controlled;

[0133] Total Track Count Unit: Calculates the total track count for each firing point;

[0134] The single-channel energy calculation unit calculates the sum of the amplitude values ​​of all sampling points for each seismic trace;

[0135] The single-shot energy calculation unit calculates the sum of all single-track energy recorded by a single shot;

[0136] The single-gun main frequency calculation unit calculates the main frequency value recorded for a single gun.

[0137] Single-shot bandwidth calculation unit, calculates the bandwidth of a single shot record;

[0138] The single-shot signal-to-noise ratio calculation unit calculates the signal-to-noise ratio of a single shot record;

[0139] Zero-value trace statistics unit: counts the number of seismic traces with a single trace energy value of 0.

[0140] Abnormal channel statistics unit: calculates and counts the total number of abnormal low-energy channels and abnormal high-energy channels;

[0141] The node data recovery rate calculation unit calculates the node data recovery rate of a single shot record;

[0142] The node data efficiency calculation unit calculates the node data efficiency of a single shot record.

[0143] The single-shot record quality analysis result output unit outputs the single-shot record quality analysis results.

[0144] The present invention discloses a method and system for synthesizing and quality controlling nodal data of explosive seismic sources based on shared memory. The method includes: scanning all raw data collected by nodal instruments on the same day (hereinafter referred to as nodal data), loading a nodal instrument equipment replacement list, and establishing a nodal data information table; loading the shot-detector relationship file and shot-point excitation file corresponding to each shot collected on the same day, and establishing an observation system data table; selecting the detector line to be cut and synthesized and the field file number (i.e., FFID) to be cut and synthesized on that detector line, and creating a shared memory space; cutting the nodal data and synthesizing detector line records; calculating and statistically analyzing the nodal data recovery rate and nodal data effectiveness rate of the detector line records, and outputting the quality analysis results of the detector line records; automatically merging the detector line records into single-shot records according to the field file number order; calculating the energy, frequency, bandwidth, and signal-to-noise ratio of the single-shot records, and statistically analyzing the nodal data recovery rate and nodal data effectiveness rate, and outputting the quality analysis results of the single-shot records. This invention can segment and synthesize detector line data from the day's collected nodal data within 1.5 hours, and automatically analyze the nodal data recovery rate and nodal data effectiveness of the detector line records, thus monitoring the quality of the day's field acquisition data in a timely manner. Simultaneously, it can merge detector line data from 2000 shots performed that day into single-shot data within 30 minutes, automatically calculating the energy, frequency, bandwidth, and signal-to-noise ratio of the synthesized single-shot data, counting the number of zero-value channels and abnormal channels, calculating the nodal data recovery rate and nodal data effectiveness of the single-shot records, and evaluating the quality of the single-shot seismic acquisition. In other words, it can complete the segmentation and synthesis of 2000 shot records daily within 2 hours with automatic quality control, significantly improving the efficiency of segmentation and synthesis of explosive source nodal acquisition data and ensuring the quality of seismic acquisition data. Attached Figure Description

[0145] Figure 1 This is a flowchart of the method for synthesizing and quality controlling explosive source node data based on shared memory according to the present invention;

[0146] Figure 2 This is a structural diagram of the explosive source node data synthesis and quality control system based on shared memory according to the present invention;

[0147] Figure 3 is a schematic diagram of the node instrument replacement data table and node data information created in a specific embodiment of the present invention;

[0148] Figure 4 is a schematic representation of the shot-receiver relationship file, shot-point excitation file, and created observation system data designed in a specific example of the present invention.

[0149] Figure 5 This is a schematic diagram illustrating the result of creating a shared memory space in a specific embodiment of the present invention;

[0150] Figure 6 This is a schematic diagram of a detector line record in a specific embodiment of the present invention, where node data is cut and synthesized.

[0151] Figure 7 This is a schematic diagram of the quality analysis results of the detector line recording in a specific embodiment of the present invention;

[0152] Figure 8 This is a schematic diagram illustrating the selection of field file numbers to be merged and the number of parallel computing threads in a specific embodiment of the present invention.

[0153] Figure 9 This is a schematic diagram of a single-shot record obtained by merging in a specific embodiment of the present invention;

[0154] Figure 10 This is a schematic diagram of the single-shot recording quality analysis results in a specific embodiment of the present invention;

[0155] Figure 11 This is a schematic diagram of the list of node instruments manually compiled when directly synthesizing single-shot records shot by shot in the traditional method. Detailed Implementation

[0156] It should be noted that the following detailed descriptions are exemplary and intended to provide further illustration of the invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0157] It should be noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the exemplary embodiments of the present invention. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Furthermore, it should be understood that when the terms "comprising" and / or "including" are used in this specification, they indicate the presence of features, steps, operations, and / or combinations thereof.

[0158] The method for synthesizing and quality control of explosive source node data based on shared memory in this invention mainly includes:

[0159] Step 1: Construct the node data information table;

[0160] Input the basic parameters of the work area, load the node instrument equipment replacement file, scan all the raw data collected by the node instruments that were recovered on the same day, and create a node data information table;

[0161] Step 2: Construct the observation system data table;

[0162] Load the shot-detector relationship file and shot-point excitation file corresponding to the single shot collected on the same day, and create an observation system data table;

[0163] Step 3: Select the detector line to be cut and synthesized and the field file number to be cut and synthesized on the detector line, and create a shared memory space;

[0164] Select the detector line to be cut and synthesized and the field file number to be cut and synthesized on the detector line, and create a shared memory space for each detector point station number that can cache node data for 2 hours.

[0165] Step 4: Cut the node data and synthesize the detector line record;

[0166] Step 5: Batch load the detector line records, automatically calculate the number of zero-value channels and the number of abnormal channels, calculate the node data recovery rate and node data effectiveness, and output the detector line record quality analysis results;

[0167] Step 6: Based on the detector line records, automatically batch merge all single shot records in the order of field file number to obtain all detector line records corresponding to each shot point.

[0168] Step 7: Batch load the synthesized single-shot records, automatically calculate the single-shot energy, frequency, bandwidth and signal-to-noise ratio, statistically analyze the node data recovery rate and node data effectiveness, and output the single-shot quality analysis results.

[0169] The main advantages of this invention are: (1) It creates a node data information table. First, it automatically extracts the equipment number, detector line number, and detector station number of each node instrument from the node data, and establishes the correspondence between the node instrument equipment number and the detector line number and detector station number, thus changing the previous need to manually create a node instrument equipment list based on the field node instrument deployment record (e.g., Figure 11The system employs several methods to improve efficiency and reduce errors in manual data processing. First, the node data information table records the storage address of the original data for each node instrument, enabling rapid location of node data during subsequent node data segmentation and improving extraction efficiency. Second, it matches the replacement node instrument with an initial field file number, allowing for quick location of the initial field file number for the replacement node instrument's detector station number during subsequent data segmentation and synthesis, further enhancing node data extraction efficiency.

[0170] (2) An observation system data table was created, and the relationship between field file number, shot point excitation time and receiver arrangement was established. The efficiency of database retrieval is much higher than that of frequently opening and closing shot point relationship files, shot point excitation files, nodal instrument equipment lists and nodal instrument equipment replacement files in the past, and the efficiency of frequently searching and matching data among the four files is improved, thereby improving the efficiency of nodal data cutting and synthesis.

[0171] (3) In the step of creating shared memory space, firstly, two shared memory spaces are allocated for each receiver station number. Each memory space can cache one hour of node data (the I-Nodal node instrument generates one node data file per hour during seismic acquisition, and 24 node data files are generated in 24 hours of acquisition). This is because single-shot records may have cross-hour phenomena, that is, the first half of the single-shot data is in the first hour node data file, and the second half of the single-shot data is in the second hour node data file. Therefore, this invention is designed to cache two hours of node data. When the single-shot cross-hour phenomenon occurs in the first hour, the node data of the second hour is already in the buffer memory, avoiding frequent reading and releasing of node data when synthesizing cross-hour single shots, thereby improving the node data cutting efficiency; secondly, for the two shared memory spaces allocated for each receiver station number, two pointers are used to record the starting position of the two memory blocks respectively. The first pointer records the address of the memory block to be processed (referred to as the processing pointer), and the second pointer records the address of the memory block of the next hour (referred to as the storage pointer). During initial loading, the first hour node data is loaded into the storage pointer. Data is loaded into the first memory block, and node data for the second hour is loaded into the second memory block. The processing pointer points to the first memory block, and the storage pointer points to the second memory block. When all single shots received by this detector line are cut and synthesized within the first hour, the processing pointer points to the second memory block, and the storage pointer points to the first memory block. The node data for the third hour is automatically loaded into the first memory block, and so on, until all detector line records are merged. The shared memory space is then released to avoid frequent memory allocation and release, improve the efficiency of node data extraction and release, and thus improve the efficiency of node data cutting and synthesis. Third, based on the detector line to be cut and synthesized and the field file number to be cut and synthesized on that detector line, the maximum number of shots fired and the number of receiving channels of the longest detector line within an hour are counted. A memory space is allocated to cache detector line records. All cut and synthesized detector line records share this memory space. This space can ensure that the longest receiving arrangement can be temporarily stored within the hour with the largest number of shots fired, avoiding frequent memory allocation and release due to different detector line lengths and different number of shots fired per hour, thus improving the efficiency of node data cutting and synthesis.

[0172] (4) During the process of cutting and synthesizing detector line data, the maximum number of detector lines that can be cut and synthesized simultaneously is automatically calculated based on the number of detector lines to be cut, the number of detector points on the longest detector line, the maximum number of shots fired in one hour, the length of the seismic record, the sampling interval, and the server memory size. That is, the maximum number of threads that can be opened at the same time, so as to improve the utilization rate of computer memory resources and thus improve the efficiency of node data cutting. Second, during the process of node data extraction and release, the field file numbers to be cut and synthesized are grouped by hour according to the GPS time of single shot excitation, and the node data of 2 hours are extracted sequentially. After all the single shots received by the detector line in the first hour are cut and synthesized, the single shot data of the second hour is cut and synthesized. This method can avoid frequent reading and release of node data of different hours, thereby improving the efficiency of node data cutting and synthesis.

[0173] (5) The invention method can cut and synthesize the node data collected on the same day into detector line records within 1.5 hours, and automatically calculate the node data recovery rate and node data effectiveness, and monitor the quality of field data collected on the same day in a timely manner. This changes the previous technical status quo that could not monitor the quality of field data collected on the same day, and ensures the quality of collected data.

[0174] (6) In the process of merging single shot records, firstly, the detector line number of each detector line record is directly read and automatically merged according to the detector line number order to obtain the single shot record. This process does not require loading the observation system, which improves the merging efficiency of single shot data; secondly, multi-threaded parallel processing technology is adopted to generate single shot data in batches, which improves the data merging efficiency.

[0175] (7) During the quality control of single-shot records, the energy, frequency, bandwidth and signal-to-noise ratio of a single shot can be automatically analyzed in batches, and the node data recovery rate and node data effectiveness can be statistically analyzed to give the single-shot quality analysis results. This changes the previous technical status quo of having to transfer the synthesized single-shot data to the seismic data processing system and manually analyze the quality of each shot, thus improving the efficiency of data acquisition quality control.

[0176] (8) The method of the present invention greatly improves the efficiency of cutting and synthesizing node acquisition data. Starting from the 11th day of seismic acquisition, it only takes 2 hours per day to cut and synthesize the 2,000 shot records produced on that day (first, it takes 1.5 hours to cut and synthesize the node data collected on that day into detector line data, and then it takes another 0.5 hours to merge all the detector line data cut and synthesized in the previous 11 days into shot data); while the traditional method of cutting and synthesizing shot records directly requires 21.7 hours, which can easily cause data backlog in actual production. Usually, after the field acquisition is completed, it still takes 2-3 days for shot data synthesis, which delays the construction period.

[0177] The following are several specific embodiments of the application of the present invention.

[0178] Example 1:

[0179] Figure 1 This is a flowchart of the method for synthesizing and quality controlling explosive source node data based on shared memory, as described in this invention. Figure 1 As shown, a method and system for synthesizing and quality controlling explosive source node data based on shared memory includes:

[0180] Step 101: Establish a node data information table, which specifically includes:

[0181] Input the basic parameters of the work area, including sampling interval and seismic record length information;

[0182] Download all raw data from I-Nodal nodal instruments recovered from the field that day;

[0183] Load the nodal instrument replacement files corresponding to the single shot collected on the same day, and create a nodal instrument replacement data table. The nodal instrument replacement data table includes: detector line number, detector stake number, starting file number, equipment number before replacement, and equipment number after replacement. Figure 3-1 As shown;

[0184] Establish a node data information table, which includes: detector line number, detector stake number, node instrument equipment number, number of recovered node data files, number of abnormal node data files, node data storage path, and the starting field file number for changing the node instrument.

[0185] The specific method for constructing the node data information table is as follows:

[0186] Scan all node data downloaded that day to obtain the device number, detector line number, detector stake number and storage address of each node instrument;

[0187] Statistical analysis of the number of node data files collected by each node instrument;

[0188] Based on the sampling interval of the work area, the number of abnormal files is statistically analyzed (according to the working principle of Sinopec's proprietary I-Nodal node instrument, when the node instrument is powered on, one node data file is generated per hour. Under the condition of fixed seismic data sampling interval, the size of a single node data file is a fixed value, such as: when the sampling interval is 2ms, the size of the node data in one hour is 7202KB, and when the sampling interval is 1ms, the size of the node data in one hour is 14404KB), and the abnormal file path and abnormal file size are given. These data do not participate in the subsequent node data cutting and synthesis.

[0189] According to the node instrument replacement data table, mark the starting field file number for replacing the node instrument;

[0190] Construct a node data information table, such as Figure 3-2 As shown, the node instrument with detector line number 1201 and detector point number 5325 has two abnormal files. These two abnormal files will not participate in the subsequent node data cutting and synthesis. The node instrument with detector line number 1201 and detector point number 5312 is marked with the starting file number 858.

[0191] Step 102: Construct a data table of the observation system collected on the same day in the construction area, specifically including:

[0192] Load the shot-detector relationship file corresponding to the single shot collected on the same day. The shot-detector relationship file includes: field file number (i.e., ffid), shot point line number, shot point stake number, detector point line number, detector point start stake number, detector point end stake number, etc. Figure 4-1 As shown;

[0193] Load the firing point files corresponding to the single shot collected that day. The firing point files include: field file number (i.e., fffid), GPS time, etc. Figure 4-2 As shown;

[0194] Based on the aforementioned shot-receiver relationship file and shot-point firing file, an observation system data table is constructed. The observation system data table includes: field file number, shot-point line number, shot-point stake number, firing time, number of receiver permutations, starting permutation number, ending permutation number, total number of receiver channels, number of missing node channels, data integrity, and node data retrieval rate. Figure 4-3 As shown;

[0195] The specific method for constructing the observation system data table is as follows:

[0196] Obtain the field file number, shot point line number, and shot point stake number from the shot checkpoint relationship file;

[0197] Obtain the field file number and firing time from the firing point firing file, match the field file number in the firing point firing file with the field file number in the firing point checkpoint relationship file, and establish a one-to-one correspondence between firing time and firing point line number and firing point station number;

[0198] The starting and ending receiver sequence numbers of each shot point in the shot-receiver relationship file, as well as the starting and ending receiver station numbers on each receiver sequence, are statistically analyzed to calculate the total number of receiver channels for each shot point, denoted as allTrace.

[0199] Match the field file number, receiver line number, receiver start station number, and receiver end station number in the shot-receiver relationship file with the receiver line number and receiver station number in the node data information table created in step 101. Count the missing receiver station numbers and calculate the number of missing node traces and node data recovery rate for each shot point. The number of missing node traces is denoted as noTrace, and the node data recovery rate is denoted as Q.

[0200] The node data recovery rate is calculated as Q = (allTrace - noTrace) / allTrace. To ensure the quality of data collected in the field, the node data recovery rate is required to be greater than 99%.

[0201] The method for judging the integrity of node data is as follows: when the node data recovery rate is not equal to 100%, the system automatically determines that the node data of the single shot is incomplete, and the system automatically counts and outputs a list of missing geophone station numbers on each geophone line recorded by the single shot. This makes it easier for construction personnel to locate the lost node instrument equipment in the field based on the geophone station numbers, thereby improving the quality of the collected data.

[0202] Step 103: Select the receiver line to be cut and synthesized, and the field file number to be cut and synthesized on that receiver line. Create a shared memory space. Create two shared memory spaces for each receiver station number. Each shared memory space can cache one hour of node data, specifically including:

[0203] The first step is to automatically obtain the server's memory capacity, denoted as RAM. To fully utilize the server's memory resources, improve the computational efficiency of cutting and synthesizing detector line records, and ensure operational stability, the default maximum allowed memory usage rate is 80% (this value can be manually adjusted as needed, but the system limits the maximum memory usage rate to no more than 95%, otherwise it may cause system instability). Assuming the server's memory capacity is 128GB, the system defaults to allowing a maximum shared memory space of 102.4GB, calculated as follows: 0.8 × RAM = 0.8 × 128 = 102.4GB.

[0204] The second step is to select the detector lines that need to be cut and synthesized, and to count the number of detector lines that need to be cut and synthesized, as well as the number of nodal devices on the longest detector line. The specific method is as follows:

[0205] Select the detector lines to be cut and synthesized. Suppose we select 1141, 1146, 1151, 1156, 1161, 1166, 1171 and 1176, a total of 8 detector lines. Then we denote the number of detector lines to be cut and synthesized as LN, LN = 8.

[0206] Obtain the starting and ending station numbers of each geophone line, calculate the required number of geophones on that line, and denot it as AN[LN], where AN[LN] = {an1, an2, ..., an...} LN}, where an i Let i be the number of detector points on the i-th detector line, where i = 1, 2, ..., LN;

[0207] Statistics i The maximum value of (i = 1, 2, ..., LN) is denoted as an. max ;

[0208] In this embodiment of the invention, AN[8]={665,665,665,665,665,665,665,665,665}, then an max =665;

[0209] The third step is to obtain the sampling interval of the work area, denoted as Dt, and automatically calculate the data file size of a single node, denoted as Bin;

[0210] In this embodiment of the invention, Dt = 1, then Bin = 14404KB (According to the acquisition principle of the I-Nodal node instrument, under the condition of fixed seismic data sampling interval, the size of a single node data file is a fixed value. If the sampling interval is 1ms, then the node data size for one hour is 14404KB; if the sampling interval is 2ms, then the node data size for one hour is 7202KB).

[0211] The fourth step is to select the field file number to be cut and synthesized, and obtain the maximum number of blasts per hour, specifically:

[0212] The field file numbers are grouped according to the GPS time of the firing point. Single shots fired within an hour are grouped together, and the number of shots fired in each hour is counted.

[0213] The maximum number of explosions per hour is denoted as sn. max ;

[0214] Assuming the maximum number of shots fired per hour on the first day is 210, then sn max =210;

[0215] The fifth step is to calculate the number of sampling points for a single shot record, specifically:

[0216] Obtain the sampling interval and seismic record length for the work area. The sampling interval is denoted as Dt and the seismic record length is denoted as SeisLen. Calculate the number of sampling points per shot record, denoted as Samp. The calculation method is: Samp = SeisLen ÷ Dt.

[0217] In this embodiment of the invention, Dt = 1, SeisLen = 7000ms, and Samp = 7000 is calculated;

[0218] The sixth step is to calculate the maximum number of detector lines that the server can cut simultaneously, that is, to calculate the maximum number of parallel threads that can be created. Specifically:

[0219] First, assume that the server can cut the maximum number of synthesized detector lines at the same time as CL, that is, the maximum number of parallel threads that can be opened is CL;

[0220] Let N1 be the memory space required for the original node data when cutting and synthesizing the longest detector line. The calculation method is: N1 = an max ×Bin×2, where “2” indicates that two shared memory spaces are created for each detector station number, that is, to cache 2 hours of node data;

[0221] In this embodiment of the invention:

[0222] The buffer space required for the synthesized detector line record during the hour with the most firings, when cutting and synthesizing the longest detector line, is denoted as N2. The calculation method is as follows:

[0223] N2=an max ×sn max ×samp×4, where “4” indicates that each sample of the detector line data is stored using 4 bytes;

[0224] In this embodiment of the invention:

[0225]

[0226] The maximum number of parallel threads that the server can create is denoted as CL, where CL is a positive integer. The calculation method is as follows:

[0227]

[0228] CL≤LN

[0229] To make the above system of inequalities hold, take the largest positive integer value of CL;

[0230] In this embodiment of the invention:

[0231]

[0232] CL≤8

[0233] To make the above system of inequalities hold, we calculate CL=4, which means that 4 threads can be opened simultaneously. Each thread can independently cut and synthesize 1 detector line at the same time. The other 4 detector lines are queued in the queue. The thread that completes the cutting of the first detector line data first will start cutting the fifth detector line data, and so on, until all detector line data is cut.

[0234] Step 7: Based on the CL value, create a shared memory space, denoted as shareM, where shareM = CL × (N1 + N2), to maximize the use of server memory resources and improve the efficiency of detector line data cutting and synthesis.

[0235] In this embodiment of the invention:

[0236] shareM=CL×(N1+N2)=4×(18.27+3.65)=87.68GB;

[0237] like Figure 5 As shown, 1141, 1146, 1151, 1156, 1161, 1166, 1171, and 1176 are the selected detector line numbers to be cut and synthesized. 1-1838 is the range of field file numbers to be cut and synthesized on these 8 detector lines (because the image cannot fully display the data, only a portion is shown, and the field file number range is 1-20). This method automatically calculates the number of lines to be synthesized simultaneously to be 4, that is, the number of parallel threads opened (CL) is 4, and the shared memory space created is 87.68GB.

[0238] Step 104: Cut the node data and synthesize the detector line record, specifically including:

[0239] Within each thread, the synthesized detector line record is cut and synthesized independently, without interference between threads. The following example illustrates the process of cutting and synthesizing the detector line record using a single thread. Assume the first thread to start cuts and synthesizes detector line 1141, specifically as follows:

[0240] The first step is to extract node data into shared memory, specifically:

[0241] Based on the GPS time of single shot firing, extract the node data for 2 hours corresponding to all detector point station numbers on the 1141 detector line in sequence and load them into the two shared memory spaces created in step 103.

[0242] Assuming the earliest GPS trigger time is 2022 / 02 / 13 13:35:34.895, then determine whether the node data of each node on the 1141 detector line for the two hours of 13:00 and 14:00 on February 13, 2022 has been loaded into the shared memory space created in step 103;

[0243] If so, then abandon the loading process;

[0244] If not, then load the node data for the two hours of 13:00 and 14:00 on February 13, 2022 into the two shared memory spaces created in step 103, where the node data for the first hour is loaded into the first memory block and the data for the second hour is loaded into the second memory block;

[0245] Create two memory pointers: the processing pointer points to the first memory block, and the storage pointer points to the second memory block.

[0246] The second step is to cut the node data and merge it into detector line records, specifically:

[0247] Obtain the sampling interval and seismic record length for the work area, where the sampling interval Dt = 1 and the seismic record length SeisLen = 7000. Calculate the number of sampling points per shot record as Samp = SeisLen ÷ Dt = 7000.

[0248] Based on the GPS excitation time at the shot point, Samp sample points are cut from the data at each node on the 1141 receiver line, and this is recorded as the single-channel seismic record of that shot point on the 1141 receiver line, denoted as X. j [Samp], where j is the detector point number on detector line 1141, j = 1, 2, ..., RN, RN is the number of detector points on detector line 1141, RN = 665, and Samp is the number of sampling points for single shot recording, Samp = 7000.

[0249] X j [Samp] is arranged in order of receiver station number, merged into a receiver line record, and output;

[0250] Repeat step two until all the 1141 detector records from all the single guns excited in the first hour have been synthesized.

[0251] The third step is to repeat steps one and two until all 1141 receiver records from each individual shot acquired that day have been cut and the shared memory space has been released. Specifically:

[0252] Once all the 1141 detector records from all the firing points fired within the one hour of 13:00 on February 13, 2022 have been cut and synthesized, the processing pointer will point to the second memory block, the storage pointer will point to the first memory block, and the node data of the third hour (i.e., the data at 15:00 on February 13, 2022) will be automatically loaded into the first memory block.

[0253] This process continues until all the shot points fired on the 1141 detector line that day have been cut and synthesized, releasing the created shared memory space.

[0254] like Figure 6 As shown, this is a schematic diagram of a detector line record synthesized from node data. The field file number is 43, and the detector point station number range is 5291 to 5955.

[0255] Step 105: Load the above-mentioned cut and synthesized detector line data, automatically calculate the total number of channels, zero-value channels, abnormal channels, nodal data retrieval rate, and nodal data effectiveness for each detector line record, and output the detector line record quality analysis results, specifically including:

[0256] Batch load the above-mentioned cut and synthesized detector line data;

[0257] To count the total number of channels recorded for each geophone line, specifically: obtain the starting geophone station number and the ending geophone station number for each geophone line, and calculate the total number of channels for that geophone line. The calculation formula is: Total number of channels = Ending geophone station number - Starting geophone station number + 1;

[0258] The energy of a single channel is calculated as follows: the sum of the amplitude values ​​of all sampling points in each seismic channel recorded by the detector line is calculated, denoted as A. k The calculation formula is: Where, k = 1, 2, ..., RN, RN is the total number of channels recorded by the detector line, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points of the seismic record. In this embodiment of the invention: on line 1141, RN = 665, sample = 7000;

[0259] To count the number of zero-value channels, specifically: count the number of channels A on the detector line record. k The number of seismic traces with a value of 0, k = 1, 2, ..., RN, where RN is the total number of traces recorded by the detector line;

[0260] The number of abnormal traces is counted. The method for calculating abnormal traces is as follows: take 10 traces to the left and 10 to the right of the current seismic trace (assuming the k-th trace), and denote the energy of a single trace as {A}. k-10 A k-9 A k-1 A K+1 ,...,A K+！0} Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the seismic trace had abnormally high energy. The seismic trace is considered to have an abnormally low energy level. The number of abnormal traces is the sum of the number of traces with abnormally high single-channel energy and the number of traces with abnormally low single-channel energy. The threshold value for judging single-channel energy anomalies can be adjusted appropriately according to the actual situation of the work area. According to Sinopec's enterprise standards, in the node acquisition project of the eastern exploration area, this threshold value is [1 / 10, 30], that is, the energy value A of the seismic trace. kBelow average energy value One-tenth of the average energy value, or higher than the average energy value 30 times that of the previous year, which was deemed abnormal.

[0261] The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels

[0262] The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels

[0263] According to Sinopec's corporate standards, the node data recovery rate of the node acquisition project in the eastern exploration area should not be less than 99%, and the node data validity rate should not be less than 96%. If any detector line records found to be lower than these two standards are found during inspections, the field should take timely measures to rectify the situation and improve the quality of the acquired data.

[0264] Output the detector line quality analysis results;

[0265] like Figure 7 As shown in the diagram, this is a schematic diagram of the quality analysis results of the detector line records. Since this is the first day of data collection in the field, the nodal data recovery rate and nodal data effectiveness are both low. For unqualified detector line records, it is necessary to find the lost nodal instrument equipment in the field, or re-download the missing nodal data and redo the detector line data cutting and synthesis.

[0266] Step 106: Obtain all receiver line records corresponding to each shot point in order of field file number, arrange them in order of receiver line number, and automatically batch merge them to obtain all single-shot records, specifically including:

[0267] Select the field file number to be merged;

[0268] Input the number of concurrent threads for merging single-shot records. Specifically, the default number of concurrent threads for merging single-shot records is 99 (this value is determined based on work experience, with the single-shot record data size being approximately 1GB. This number of threads can be dynamically adjusted according to the actual size of the single-shot records in the work area and the server's hard drive read / write speed).

[0269] Based on the observation system data table, the detector line numbers contained in each field file number are automatically obtained in the order of the field file numbers, arranged in the order of the detector line numbers, and single shot records are automatically synthesized in batches.

[0270] like Figure 8 As shown, this diagram illustrates the selection of the field file number to be merged and the number of parallel computing threads. The system defaults to a number of 99 parallel threads.

[0271] like Figure 9As shown, the final merged single-shot record diagram (only 3 arrangements are shown in the diagram) has a field file number of 43, detector line numbers of 1141, 1146, and 1151, and detector point chainage range of 5291 to 5955.

[0272] Step 107: Load the synthesized single-shot records, automatically calculate the single-shot energy, frequency, bandwidth, and signal-to-noise ratio, statistically analyze the node data recovery rate and node data effectiveness, and output the single-shot quality analysis results, specifically:

[0273] The calculation of single-shot energy, frequency, bandwidth, and signal-to-noise ratio (SNR) is based on the entire time window and is obtained using a general geophysical algorithm. This is not considered a technological innovation. Specifically: single-shot energy analysis uses the root mean square amplitude method; SNR analysis uses the spectrum estimation method; the dominant frequency is calculated using Fourier transform, converting single-channel data from the time domain to the frequency domain, and the frequency with the highest amplitude is recorded as the dominant frequency; the bandwidth is calculated as the difference between the maximum and minimum dominant frequencies under 22dB attenuation.

[0274] The energy of a single trace is calculated as follows: the sum of the amplitude values ​​of all sampling points in each seismic trace of the single shot record is calculated, denoted as B. k The calculation formula is: Where, k = 1, 2, ..., M, M is the total number of traces recorded by a single shot, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points in the seismic record;

[0275] To count the number of zero-value channels, specifically: count the number of channels B recorded for this single shot. k The number of seismic traces with a value of 0, k = 1, 2, ..., M, where M is the total number of traces recorded by a single shot;

[0276] The number of abnormal traces is counted. The method for calculating abnormal traces is as follows: take 10 traces to the left and 10 to the right of the current seismic trace (assuming the k-th trace), and denote the energy of a single trace as {B}. k-10 B k-9 B k-1 B K+1 ,...,B K+！0} Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the seismic trace had abnormally high energy. The seismic trace is considered to have an abnormally low energy level. The number of abnormal traces is the sum of the number of traces with abnormally high energy levels and the number of traces with abnormally low energy levels. The threshold value for judging energy anomalies can be adjusted appropriately according to the actual situation of the work area. According to Sinopec's enterprise standards, in the node acquisition project of the eastern exploration area, this threshold value is [1 / 10, 30], which is the energy value B of the seismic trace. kBelow average energy value One-tenth of the average energy value, or higher than the average energy value 30 times that of the previous year, which was deemed abnormal.

[0277] The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels;

[0278] The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels;

[0279] Output single-shot recording quality analysis results, specifically including: single-shot energy, main frequency, bandwidth, signal-to-noise ratio, number of zero-value channels, number of abnormal channels, node data recovery rate, and node data effectiveness rate;

[0280] Based on the analysis of data from the test sites in the work area, an evaluation standard for seismic acquisition data was established, specifically as follows:

[0281] The energy of a single shot at the test site in the work area is analyzed to determine the minimum energy of a single shot under different well depths and charge amounts. If the energy is lower than the minimum energy value, the single shot is considered to be unqualified. According to the data analysis of the test site in the work area of ​​this invention, a single shot with an energy of less than 30,000 is considered to be unqualified.

[0282] The frequency of the test data from the work area assessment points is analyzed to give the main frequency range and bandwidth range. When the main frequency and bandwidth values ​​are abnormal, they can be used as reference factors to determine whether a single shot is unqualified, but not as the deciding factor. The single shot record can be opened for manual evaluation. The analysis of the test data from the work area assessment points in this embodiment of the invention considers a main frequency range of 30-60Hz and a bandwidth range of 30-50Hz to be a qualified single shot.

[0283] The signal-to-noise ratio of the test data of the work area is analyzed, and the minimum signal-to-noise ratio of the work area is given. If the signal-to-noise ratio is lower than the minimum signal-to-noise ratio, the quality of a single shot is considered unqualified. In this embodiment of the invention, the minimum signal-to-noise ratio of a single shot in the work area is 1.0.

[0284] Analyzing the node data recovery rate, according to Sinopec's enterprise standards, the node data recovery rate of the node acquisition project in the eastern exploration area shall not be lower than 99%. When it is lower than this threshold, the quality of a single shot is considered unqualified.

[0285] The effectiveness rate of node data analysis is determined according to Sinopec's enterprise standards. The effectiveness rate of node data acquisition projects in the eastern exploration area must not be lower than 96%. If it is lower than this threshold, the quality of a single shot is considered unqualified.

[0286] like Figure 10As shown in the figure, this is a schematic diagram of the single-shot recording quality analysis results. The diagram shows the single-shot quality analysis results of field files 53 to 60. Among them, the single shots with field files 57 and 59 have energy below 30,000 and signal-to-noise ratio below 1.0, and are classified as unqualified shots. The remaining single shots are qualified shots.

[0287] Example 2:

[0288] Figure 2 This is a structural diagram of the explosive source node data synthesis and quality control system based on shared memory, as described in this invention. Figure 2 As shown, a shared memory-based system for synthesizing and quality controlling explosive source node data includes:

[0289] Node data information table construction module 201 is used to construct the node data information table;

[0290] The observation system data table construction module 202 is used to construct the observation system data table collected on the same day;

[0291] The shared memory space creation module 203 is used to allocate shared memory space for detector line record cutting and synthesis, thereby improving the efficiency of node data cutting and synthesis.

[0292] The detector line record cutting and synthesis module 204 is used to cut node data and synthesize detector line records.

[0293] The detector line recording quality analysis module 205 is used to analyze the quality of the synthesized detector line recordings, that is, to monitor the quality of the data acquired from the field earthquake on the same day.

[0294] The single-shot record batch merging module 206 is used to batch synthesize the detector line data into single-shot records;

[0295] The single-shot record quality analysis module 207 is used to analyze the quality of synthesized single-shot records.

[0296] The node data information table construction module 201 specifically includes:

[0297] The work area information input unit is used to input the sampling interval and seismic record length information of the construction work area;

[0298] The node instrument replacement file loading unit is used to create a node instrument replacement data table.

[0299] The node data scanning unit is used to automatically extract the equipment number, detector line number, and detector station number of each node instrument, and record the data storage address of each node instrument.

[0300] The node data analysis unit is used to count the total number of node data files collected on the same day and analyze the number of abnormal files based on work area information.

[0301] The node data information table construction unit is used to establish a one-to-one correspondence between the device number, detector line number, detector station number and storage address of each node instrument, and to mark the starting file number of the replaced node instrument.

[0302] The observation system data table construction module 202 specifically includes:

[0303] The blasting checkpoint relationship file loading unit is used to obtain the blasting checkpoint relationship corresponding to a single blast collected in the construction area on that day.

[0304] The firing point firing file loading unit is used to obtain the single firing time collected in the construction area on the same day;

[0305] The observation system data table construction unit is used to establish the relationship between the nodal instrument equipment and the observation system, count the number of missing nodal traces for each shot point, calculate the nodal data recovery rate, and provide a judgment result on whether the nodal data is complete.

[0306] The shared memory space creation module 203 specifically includes:

[0307] The server memory capacity acquisition unit is used to obtain the current memory capacity of the server.

[0308] The detector line and field file number selection unit is used to select the detector line that needs to be cut and synthesized, and the field file number that needs to be cut and synthesized on the detector line.

[0309] A node data file size calculation unit is used to calculate the size of a single node data file in the work area;

[0310] The field file number grouping and statistics unit is used to sort the field file numbers according to the GPS triggering time, group them by hour, count the number of blasts in each hour, and obtain the maximum number of blasts in a single hour.

[0311] The maximum number of parallel threads calculation unit is used to calculate the maximum number of synthesized detector lines that the server can cut simultaneously.

[0312] Create shared memory units to create two shared memory spaces for each detector station number. Each shared memory space can cache one hour of node data.

[0313] The detector line recording, cutting, and synthesis module 204 specifically includes:

[0314] The node data extraction unit is used to extract the node data to be cut and synthesized into data blocks in the shared memory space;

[0315] The node data cutting unit is used to cut node data into single-channel seismic data based on the shot point excitation GPS time and the number of seismic record sampling points.

[0316] The detector line data synthesis unit is used to arrange the cut single-channel seismic data according to the detector point station number order and output the detector line record.

[0317] The detector line recording quality analysis module 205 specifically includes:

[0318] Detector line record loading unit, used for batch loading of detector line records to be quality controlled;

[0319] Total Channel Count Unit, used to calculate the total number of channels for each detector line;

[0320] A single-channel energy calculation unit is used to calculate the sum of the amplitude values ​​of all sampling points for each seismic trace;

[0321] Zero-value trace statistics unit, used to count the number of seismic traces with a single trace energy value of 0;

[0322] The abnormal channel statistics unit is used to calculate and count the total number of abnormal low-energy channels and abnormal high-energy channels.

[0323] The node data recovery rate calculation unit is used to calculate the node data recovery rate of the detector line record.

[0324] The node data efficiency calculation unit is used to calculate the node data efficiency of the detector line record.

[0325] The single-shot record batch merging module 206 specifically includes:

[0326] The field file number selection unit is used to obtain the number of single guns to be synthesized and the field file number;

[0327] The input unit for the number of parallel threads to start merging single-shot records is used to allocate the number of parallel threads for batch synthesis of single-shot records.

[0328] The automatic batch merging unit for single-shot data is used to obtain the detector line number and detector point number from the detector line record, rearrange them in order, merge them to obtain the single-shot record, and output it.

[0329] The single-shot recording quality analysis module 207 specifically includes:

[0330] The single-shot record loading unit is used for batch loading of single-shot records to be quality controlled;

[0331] The total number of shots is a counting unit used to calculate the total number of shots for each firing point.

[0332] A single-channel energy calculation unit is used to calculate the sum of the amplitude values ​​of all sampling points for each seismic trace;

[0333] The single-shot energy calculation unit is used to calculate the sum of all single-track energy recorded by a single shot;

[0334] The single-gun main frequency calculation unit is used to calculate the main frequency value recorded by a single gun.

[0335] Single-shot bandwidth calculation unit, used to calculate the bandwidth of a single shot record;

[0336] A single-shot signal-to-noise ratio calculation unit is used to calculate the signal-to-noise ratio of a single shot record;

[0337] Zero-value trace statistics unit, used to count the number of seismic traces with a single trace energy value of 0;

[0338] The abnormal channel statistics unit is used to calculate and count the total number of abnormal low-energy channels and abnormal high-energy channels.

[0339] The node data recovery rate calculation unit is used to calculate the node data recovery rate of a single shot record.

[0340] The node data efficiency calculation unit is used to calculate the node data efficiency of a single shot record.

[0341] The single-shot record quality analysis result output unit is used to output the single-shot record quality analysis results.

[0342] The method and system for synthesizing and quality control of explosive source node data based on shared memory of the present invention can replace the previous method of reading the shot-detector relationship file, shot-point excitation file, node instrument equipment list and node instrument equipment replacement file, and directly cutting and synthesizing single shot records according to the field file number order, which can greatly improve the efficiency of node data cutting and synthesis and monitor the quality of node data acquisition in a timely manner.

[0343] Finally, it should be noted that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

[0344] Except for the technical features described in the specification, all other technologies are known to those skilled in the art.

Claims

1. A method for synthesizing and quality controlling explosive source node data based on shared memory, characterized in that, The method for synthesizing and quality control of explosive source node data based on shared memory includes: Step 1: Construct the node data information table; Step 2: Construct the observation system data table; Step 3, create a shared memory space; Step 4: Cut the node data and synthesize the detector line record; Step 5: Calculate the quality analysis results of the detector line records based on the detector line records; Step 6: Based on the detector line records, merge all single-shot records in batches according to the field file number order; Step 7: Calculate the single-shot quality analysis results based on the synthesized single-shot records.

2. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, In step 1, input the basic parameters of the work area, including the sampling interval and seismic record length information; Load the nodal instrument replacement files corresponding to the single shot collected on the same day, and establish a nodal instrument replacement data table. The nodal instrument replacement data table includes: detector line number, detector station number, starting file number, equipment number before replacement, and equipment number after replacement. Scan all the raw data collected by the nodal instruments that were recovered on the same day, and create a nodal data information table. The nodal data information table includes: detector line number, detector station number, nodal instrument equipment number, number of recovered nodal data files, number of abnormal nodal data files, nodal data storage path, and the starting field file number for replacing the nodal instrument.

3. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 2, characterized in that, Step 1, when creating the node data information table, includes the following steps: Scan all node data downloaded that day to obtain the device number, detector line number, detector stake number and storage address of each node instrument; Statistical analysis of the number of node data files collected by each node instrument; Based on the sampling interval of the work area, the number of abnormal files is statistically analyzed, and the abnormal file paths and sizes are given. These data are not used in subsequent node data cutting and synthesis. According to the node instrument replacement data sheet, mark the starting field file number for replacing the node instrument; Construct a node data information table.

4. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, In step 2, load the shot-receiver relationship file and shot-point excitation file corresponding to the single shot collected on that day. The shot-receiver relationship file includes: field file number (ffid), shot point line number, shot point station number, receiver line number, receiver starting station number, and receiver ending station number. The shot-point excitation file includes: field file number (ffid) and GPS time. Based on the shot-receiver relationship file and shot-point excitation file, construct the observation system data table. The observation system data table includes: field file number, shot point line number, shot point station number, excitation time, number of receiver permutations, starting permutation number, ending permutation number, total number of receiver channels, number of missing node channels, data integrity, and node data retrieval rate.

5. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 4, characterized in that, Step 2, when constructing the observation system data table, includes the following steps: Obtain the field file number, shot point line number, and shot point stake number from the shot checkpoint relationship file; Obtain the field file number and firing time from the firing point firing file, match the field file number in the firing point firing file with the field file number in the firing point checkpoint relationship file, and establish a one-to-one correspondence between firing time and firing point line number and firing point station number; The starting and ending receiver sequence numbers of each shot point in the shot-receiver relationship file, as well as the starting and ending receiver station numbers on each receiver sequence, are statistically analyzed to calculate the total number of receiver channels for each shot point, denoted as allTrace. Match the field file number, receiver line number, receiver start station number, and receiver end station number in the shot-receiver relationship file with the receiver line number and receiver station number in the node data information table created in step 1. Count the missing receiver station numbers and calculate the number of missing node traces and node data recovery rate for each shot point. The number of missing node traces is denoted as noTrace, and the node data recovery rate is denoted as Q. The formula for calculating the node data recovery rate is: Q=(allTrace-noTrace) / allTrace, To ensure the quality of data collected in the field, the node data retrieval rate is required to be greater than 99%. The system provides a judgment result on whether the node data is complete. The method for judging the integrity of node data is as follows: when the node data recovery rate is not equal to 100%, the node data of a single shot is determined to be incomplete. The system also compiles and outputs a list of missing geophone station numbers on each geophone line recorded by a single shot, so that construction personnel can locate the lost node instrument equipment in the field based on the geophone station numbers.

6. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, Step 3 includes: Step 31: Obtain the server memory capacity, denoted as RAM. To fully utilize server memory resources, improve the computational efficiency of cutting and synthesizing detector line records, and ensure operational stability, the default maximum allowed memory usage rate is 80%. Step 32: Select the detector lines to be cut and synthesized, and count the number of detector lines to be cut and synthesized, as well as the number of nodal devices on the longest detector line. The specific method is as follows: Select the detector line number to be cut and synthesized, and denote the number of detector lines to be cut and synthesized as LN; Obtain the starting and ending station numbers of each geophone line, calculate the required number of geophones on the line, and denot it as AN[LN], where AN[LN] = {an1, an2, ..., an...} LN }, where an i Let i be the number of detector points on the i-th detector line, where i = 1, 2, ..., LN; Statistics i The maximum value of (i = 1, 2, ..., LN) is denoted as an. max ; Step 33: Obtain the sampling interval of the work area, denoted as Dt, and automatically calculate the data file size of a single node, denoted as Bin; Step 34: Select the field file number to be cut and synthesized, and obtain the maximum number of blasts per hour, specifically: The field file numbers are grouped according to the GPS time of the firing point. Single shots fired within an hour are grouped together, and the number of shots fired in each hour is counted. The maximum number of explosions per hour is denoted as sn. max ; Step 35, calculate the number of sampling points for a single shot record, specifically: Obtain the sampling interval and seismic record length for the work area. The sampling interval is denoted as Dt and the seismic record length is denoted as SeisLen. Calculate the number of sampling points per shot record, denoted as Samp. The calculation method is: Samp = SeisLen ÷ Dt. Step 36: Calculate the maximum number of detector lines that the server can cut simultaneously, that is, calculate the maximum number of parallel threads that can be created. Specifically: First, assume that the server can cut the maximum number of synthesized detector lines at the same time as CL, that is, the maximum number of parallel threads that can be opened is CL; Let N1 be the memory space required for the original node data when cutting and synthesizing the longest detector line. The calculation method is: N1 = an max ×Bin×2, where 2 represents creating two shared memory spaces for each detector station number, i.e. caching 2 hours of node data; The buffer space required for the synthesized detector line record during the hour with the most firings, when cutting and synthesizing the longest detector line, is denoted as N2. The calculation method is as follows: N2=an max ×sn max ×samp×4, where 4 indicates that each sample of the detector line data is stored using 4 bytes; The maximum number of parallel threads that the server can create is denoted as CL, where CL is a positive integer. The calculation method is as follows: CL≤LN To make the above system of inequalities hold, take the largest positive integer value of CL; Step 37: Based on the CL value, create a shared memory space, denoted as shareM, where shareM = CL × (N1 + N2), to maximize the use of server memory resources and improve the efficiency of detector line data cutting and synthesis.

7. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 6, characterized in that, Step 4 includes: Step 41, extract node data into shared memory, specifically: Based on the GPS time of single shot firing, extract the node data for 2 hours corresponding to all detector point station numbers on the detector line in sequence, and load them into the two shared memory spaces created in step 3. Create two memory pointers: the processing pointer points to the first memory block, and the storage pointer points to the second memory block. Step 42: Cut the node data and merge it into detector line records until all detector line records of all single shots fired in the first hour are synthesized. Step 43: Repeat steps 41 and 42 until all detector line records acquired that day have been cut and the shared memory space has been released.

8. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, In step 5, based on the cut and synthesized detector line data, calculate the total number of channels, the number of zero-value channels, the number of abnormal channels, the nodal data recovery rate, and the nodal data effectiveness rate for each detector line record, and output the detector line record quality analysis results, specifically including: Batch loading of cut and synthesized detector line data; To count the total number of channels recorded for each geophone line, specifically: obtain the starting geophone station number and the ending geophone station number for each geophone line, and calculate the total number of channels for the geophone line. The calculation formula is: Total number of channels = Ending geophone station number - Starting geophone station number + 1; The calculation of single-channel energy is specifically as follows: calculate the sum of the amplitude values ​​of all sampling points in each seismic trace on the detector line record, denoted as A. k The calculation formula is: Where, k = 1, 2, ..., RN, RN is the total number of traces recorded by the detector line, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points in the seismic record; The number of zero-value channels is counted, specifically: the number of channels recorded on detector line A. k The number of seismic traces with a value of 0, k = 1, 2, ..., RN, where RN is the total number of traces recorded by the detector lines; The number of abnormal traces is counted. The method for calculating abnormal traces is as follows: Take the current seismic trace, assuming it is the k-th trace, with 10 traces on each side. The energy of a single trace is denoted as {A}. k-10 A k-9 A k-1 A K+1 ,...,A K+！0 } Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the energy of the seismic channel was abnormally high. It is assumed that the energy of the seismic traces is abnormally small. The number of abnormal traces is the sum of the number of seismic traces with large single-trace energy anomalies and the number of seismic traces with small single-trace energy anomalies. The threshold value for judging single-trace energy anomalies can be appropriately adjusted according to the actual situation of the work area. The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels Output the quality analysis results of the detector line.

9. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, In step 6, all receiver line records corresponding to each shot point are obtained in order of field file number, arranged in order of receiver line number, and batch-merged to obtain all single-shot records, specifically including: Select the field file number to be merged; Input the number of concurrent threads for merging single-shot records; Based on the observation system data table, the detector line numbers contained in each field file number are obtained in order of field file number, arranged in order of detector line number, and single-shot records are synthesized in batches.

10. The method for synthesizing and quality control of explosive source node data based on shared memory according to claim 1, characterized in that, In step 7, based on the synthesized single-shot records, the single-shot energy, frequency, bandwidth, and signal-to-noise ratio are calculated. The node data recovery rate and node data effectiveness are statistically analyzed, and the single-shot quality analysis results are output, specifically: The energy, frequency, bandwidth, and signal-to-noise ratio of a single shot were calculated based on the entire time window using a general geophysical algorithm. Specifically: the energy analysis of a single shot used the root mean square amplitude method; the signal-to-noise ratio analysis used the spectrum estimation method; the dominant frequency was calculated using Fourier transform, converting the single-channel data from the time domain to the frequency domain, and the frequency with the highest amplitude was recorded as the dominant frequency; the bandwidth was calculated as the difference between the maximum and minimum dominant frequencies under 22dB attenuation. The calculation of single-channel energy is specifically as follows: calculate the sum of the amplitude values ​​of all sampling points in each seismic trace on a single shot record, denoted as B. k The calculation formula is: Where, k = 1, 2, ..., M, M is the total number of traces recorded by a single shot, X is the amplitude value, i = 1, 2, ..., sample, sample is the number of sampling points in the seismic record; The number of zero-value channels is counted, specifically: the number of single-shot records at point B. k The number of seismic traces with a value of 0, k = 1, 2, ..., M, where M is the total number of traces recorded by a single shot; The number of abnormal traces is counted. The calculation method for abnormal traces is as follows: Take the current seismic trace, assuming it is the k-th trace, with 10 traces on each side. The energy of a single trace is denoted as {B}. k-10 B k-9 B k-1 B K+1 ,...,B K+！0 } Calculate the average energy of these 20 seismic traces, denoted as . The calculation method is as follows: when At that time, it was believed that the energy of the seismic channel was abnormally high. It is assumed that the energy of the seismic traces is abnormally small, and the number of abnormal traces is the sum of the number of seismic traces with abnormally large energy and the number of seismic traces with abnormally small energy. The threshold value for judging energy anomalies can be adjusted appropriately according to the actual situation of the work area. The node data recovery rate is calculated as: (Total number of channels - Number of channels with zero values) / Total number of channels; The efficiency of node data is calculated as follows: (Total number of channels - Number of channels with zero values ​​- Number of abnormal channels) / Total number of channels; Output single-shot recording quality analysis results, specifically including: single-shot energy, main frequency, bandwidth, signal-to-noise ratio, number of zero-value channels, number of abnormal channels, node data recovery rate, and node data effectiveness rate; Based on the analysis of data from the test sites in the work area, an evaluation standard for seismic acquisition data was established, specifically as follows: Analyze the single-shot energy at the test points in the work area, determine the minimum single-shot energy under different activation well depths and activation charge amounts, and consider the single-shot quality to be unqualified if it is lower than the minimum energy value; Analyze the frequency of the test data from the work area and give the main frequency range and bandwidth range. When the main frequency and bandwidth values ​​are abnormal, they can be used as reference factors to determine the quality of a single shot is unqualified, but not as the deciding factor. The single shot record can be opened for manual evaluation. Analyze the signal-to-noise ratio of the test data in the work area and give the minimum signal-to-noise ratio for the work area. If the signal-to-noise ratio is lower than the minimum, the quality of a single shot is considered unqualified. Analyze the node data recovery rate; if it falls below the set first threshold, the quality of a single gun is considered unqualified. The efficiency of node data analysis is determined. If the efficiency falls below the set second threshold, the quality of a single gun is considered unqualified.

11. A shared-memory-based system for synthesizing and quality controlling explosive source node data, characterized in that, This shared-memory-based explosive source node data synthesis and quality control system includes: The node data information table construction module constructs the node data information table. The observation system data table construction module constructs the observation system data table collected on the same day. The shared memory space creation module allocates shared memory space for detector line record cutting and synthesis, improving the efficiency of node data cutting and synthesis. The detector line recording cutting and synthesis module cuts node data and synthesizes detector line records. The detector line recording quality analysis module analyzes the quality of the synthesized detector line recordings, which is to monitor the quality of the data acquired from the field earthquake on the same day. The single-shot record batch merging module combines detector data into single-shot records in batches. The single-shot record quality analysis module analyzes the quality of the synthesized single-shot records.

12. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The module for constructing the node's data information table specifically includes: The work area information input unit allows users to input the sampling interval and seismic record length information for the construction work area. The node instrument replacement file loading unit creates a node instrument replacement data table. The node data scanning unit extracts the device number, detector line number, and detector station number of each node instrument, and records the data storage address of each node instrument. The node data analysis unit counts the total number of node data files collected that day and analyzes the number of abnormal files based on work area information. The node data information table construction unit establishes a one-to-one correspondence between the device number, detector line number, detector station number and storage address of each node instrument, and marks the starting file number of the replaced node instrument.

13. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The data table construction module for this observation system includes: The blasting checkpoint relationship file loading unit retrieves the blasting checkpoint relationships corresponding to a single blast collected in the construction area on that day. The firing point firing file loading unit acquires the firing time of a single shot collected in the construction area on that day. The observation system data table construction unit establishes the relationship between the nodal instrument equipment and the observation system, counts the number of missing nodal traces for each shot point, calculates the nodal data recovery rate, and provides the judgment result on whether the nodal data is complete.

14. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The shared memory space creation module includes: The server memory capacity acquisition unit retrieves the current memory capacity of the server. The detector line and field file number selection unit allows you to select the detector line that needs to be cut and synthesized, and the field file number that needs to be cut and synthesized on the detector line. The node data file size calculation unit calculates the size of a single node data file in the work area. Field file number grouping and statistical unit: Sort the field file numbers according to the GPS triggering time, and group them by hour. Count the number of blasts in each hour and obtain the maximum number of blasts in a single hour. Maximum number of parallel threads calculation unit, which calculates the maximum number of synthesized detector lines that the server can cut simultaneously; Create shared memory units, creating two shared memory spaces for each detector station number. Each shared memory space can cache one hour of node data.

15. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The detector line recording cut synthesis module includes: The node data extraction unit extracts the node data to be cut and synthesized into data blocks in the shared memory space; The node data cutting unit cuts the node data according to the shot point excitation GPS time and the number of seismic record sampling points to form single-channel seismic data. The receiver line data synthesis unit arranges the segmented single-channel seismic data according to the receiver point station number order and outputs the receiver line record.

16. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The detector line recording quality analysis module includes: Detector line record loading unit, batch loading detector line records to be quality controlled; Total Channel Count Unit: Calculates the total number of channels for each detector line; The single-channel energy calculation unit calculates the sum of the amplitude values ​​of all sampling points for each seismic trace; Zero-value trace statistics unit: counts the number of seismic traces with a single trace energy value of 0. Abnormal channel statistics unit: calculates and counts the total number of abnormal low-energy channels and abnormal high-energy channels; The node data recovery rate calculation unit calculates the node data recovery rate of the detector line records; The node data efficiency calculation unit calculates the node data efficiency recorded by the detector line.

17. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The single-shot record batch merging module includes: The field file number selection unit to be synthesized obtains the number of single guns to be synthesized and the field file number; The number of parallel threads for merging single-shot records is input to determine the number of parallel threads for batch synthesis of single-shot records. The automatic batch merging unit for single-shot data obtains the detector line number and detector point number from the detector line record, rearranges them in order, merges them to obtain the single-shot record, and outputs it.

18. The explosive source node data synthesis and quality control system based on shared memory according to claim 11, characterized in that, The single-shot record quality analysis module includes: Single-shot record loading unit, batch loading of single-shot records to be quality controlled; Total Track Count Unit: Calculates the total track count for each firing point; The single-channel energy calculation unit calculates the sum of the amplitude values ​​of all sampling points for each seismic trace; The single-shot energy calculation unit calculates the sum of all single-track energy recorded by a single shot; The single-gun main frequency calculation unit calculates the main frequency value recorded for a single gun. Single-shot bandwidth calculation unit, calculates the bandwidth of a single shot record; The single-shot signal-to-noise ratio calculation unit calculates the signal-to-noise ratio of a single shot record; Zero-value trace statistics unit: counts the number of seismic traces with a single trace energy value of 0. Abnormal channel statistics unit: calculates and counts the total number of abnormal low-energy channels and abnormal high-energy channels; The node data recovery rate calculation unit calculates the node data recovery rate of a single shot record; The node data efficiency calculation unit calculates the node data efficiency of a single shot record. The single-shot record quality analysis result output unit outputs the single-shot record quality analysis results.

Images