Ocean bottom node multi-component rotation quality control method and apparatus based on cross-correlation principle
By collecting multi-component seismic data at the seabed node, the polarity of the three-components X, Y, and Z are judged using the principle of cross-correlation, and a three-component polarity analysis diagram after rotation is formed for quality control, which solves the problem of difficult control of the rotation effect of the seabed node and improves the efficiency and accuracy of quality control.
Patent Information
- Application Number
- PCT/CN2024/116853
- Authority / Receiving Office
- WO · WO
- Patent Type
- Applications
- Current Assignee / Owner
- Priority Date
- 2023-12-04
- Filing Date
- 2024-09-04
- Publication Date
- 2025-06-12
AI Technical Summary
During the multi-component seismic data acquisition process at the seabed node, due to the uncertainty of construction conditions, it is difficult to ensure that the rotation effect of the three-component detector meets the design requirements, resulting in transverse wave energy on the Z component and longitudinal wave energy on the horizontal component, which affects subsequent wave field separation, imaging and inversion processing.
The multi-component rotation quality control method of subsea nodes based on the principle of cross-correlation is adopted. By obtaining the P component data collected by the subsea node and the rotated three-component data, the non-directional characteristics of the P component are used to calculate the correlation coefficients to judge the polarity of the X, Y, and Z components, thereby forming a rotating three-component polarity analysis diagram for quality control.
The need to control the rotation effect of multiple nodes simultaneously is achieved, the efficiency and accuracy of quality control is improved, the accuracy of the rotation effect is ensured, and the impact on subsequent processing is reduced.
Smart Images

Figure CN2024116853_12062025_PF_FP_ABST
Abstract
Description
Multi-component rotation quality control method and device for seabed nodes based on cross-correlation principle Technical Field
[0001] The present invention relates to a quality control technology for the multi-component rotation effect of a seabed node, and specifically to a multi-component rotation quality control method for a seabed node based on the cross-correlation principle, a multi-component rotation quality control device for a seabed node based on the cross-correlation principle, an electronic device, and a machine-readable storage medium. Background Art
[0002] With the continuous development of geophysical exploration technology, the use of seabed nodes for acquisition has become the main trend in marine seismic exploration. Placing seabed nodes, i.e., detectors, directly on the seabed not only makes the observation system layout more flexible and can better meet the exploration needs of all-round, high-density, and large offset, but also can receive shear waves and longitudinal waves, realize four-component acquisition, and make the seismic records contain richer information, providing a data basis for converted wave imaging, fracture determination, reservoir inversion, etc.
[0003] When acquiring multi-component seismic data on the seafloor, especially seismic data, the uncertainties in construction conditions prevent direct control of the placement direction of seafloor nodes. This makes it difficult to ensure that the Z component of each three-component geophone is placed vertically and the X component is placed parallel to the survey line in accordance with the acquisition design. Therefore, it is necessary to use the three-component geophone directional information recorded by the inclinometer, gyroscope, or compass during data acquisition to determine the true placement direction of the three-component geophone and complete three-component rotation correction. If the data is not rotated properly, shear wave energy will be present in the Z component and longitudinal wave energy in the horizontal component, which will seriously affect subsequent wavefield separation, imaging, and inversion processing. Therefore, quality control of the three-component directional rotation effect is a key foundational step.
[0004] Currently, a commonly used quality control method involves performing single-node, multi-component RMS amplitude quality control on raw continuous gather data. This method relies on comprehensively presenting the energy value of each node at each shot to achieve quality control of the rotation effect. However, this method has limitations. First, due to the current large volume of node data, reviewing the RMS amplitude values of all nodes is very time-consuming, and the efficiency of quality control needs to be improved. Second, it cannot meet the requirements for simultaneous rotation quality control of multiple nodes.
[0005] Summary of the Invention
[0006] The purpose of the embodiments of the present invention is to provide a multi-component rotation quality control method, device, equipment and storage medium for seabed nodes. The method utilizes the non-directional physical characteristics of the P component received by the seabed node. First, the polarity characteristics of the P component are judged according to the amplitude value of the P component. Then, the principle of cross-correlation is used to perform cross-correlation operations on the other three components with the P component respectively. The polarities of the X, Y, and Z components are judged according to the obtained correlation coefficients. Finally, the three-component polarity analysis diagrams after rotation of all nodes are subjected to quality control comparison to meet the needs of simultaneous rotation effect quality control of multiple nodes, thereby improving the efficiency and accuracy of quality control.
[0007] In order to achieve the above objectives, the present invention provides a method for multi-component rotation quality control of a seabed node based on the cross-correlation principle in a first aspect, the method comprising:
[0008] Obtain the P component data collected by the seabed node and the rotated X, Y, and Z component data, and separate the P, X, Y, and Z four-component data;
[0009] Preprocessing the separated four-component data to obtain preprocessed four-component data;
[0010] Obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data;
[0011] Calculate the correlation coefficients between the P component water detection and the rotated X, Y, and Z components land detection respectively;
[0012] Based on the polarity quality control diagram of the P component and the correlation coefficient between the water test of the P component and the land test of the rotated X, Y, Z components, a polarity analysis diagram of the rotated X, Y, Z components is formed;
[0013] The rotation effect of the rotated X, Y, and Z components is quality controlled based on the polarity analysis diagram of the rotated X, Y, and Z components.
[0014] In the embodiment of the present application, a polarity analysis diagram of the rotated X, Y, and Z components is formed based on the polarity quality control diagram of the P component and the correlation coefficient between the water detection of the P component and the land detection of the rotated X, Y, and Z components, including:
[0015] The polarity of each component is determined based on the correlation coefficient between the P component water test and the rotated X, Y, and Z components land test:
[0016] If the correlation coefficient is greater than or equal to the first preset value, it is determined that the P component is positively correlated with the corresponding component, and the polarity of the corresponding component is consistent with the polarity of the P component;
[0017] If the correlation coefficient is less than or equal to the second preset value, it is determined that the P component is negatively correlated with the corresponding component, and the polarity of the corresponding component is opposite to that of the P component;
[0018] If the correlation coefficient is greater than the second preset value and less than the first preset value, it is determined that the P component has no correlation with the corresponding component;
[0019] A polarity analysis diagram of the corresponding component is formed according to the polarity of each component.
[0020] According to the above technical means, the polarity of the X, Y, and Z component data can be determined based on the correlation coefficient between the P component and the rotated X, Y, and Z component data, thereby determining the polarity analysis diagram of the corresponding component according to the polarity of each component, providing a data basis for quality control.
[0021] In the embodiment of the present application, the rotation effect of the rotated X, Y, and Z three-component data is quality-controlled based on the polarity analysis diagram of the rotated X, Y, and Z three-component data, including:
[0022] Traverse each node in the polarity analysis diagram of each component to determine whether the polarity of each node is consistent. If there is an abnormal polarity, it is determined that the corresponding component is not rotated in place or rotated incorrectly.
[0023] According to the above technical means, it is possible to determine whether there is an abnormality in the corresponding polarity based on the X, Y, and Z three-component polarity analysis diagram, and to quickly perform rotation quality control.
[0024] In the embodiment of the present application, the separated four-component data is preprocessed to obtain preprocessed four-component data, including:
[0025] Perform offset limiting processing on the separated component data to obtain the component data after offset limiting;
[0026] Linear correction is performed on each component data after the restricted offset processing to obtain the four-component data after linear correction.
[0027] According to the above technical means, the time difference between the travel time of the reflected waves at different offsets and the travel time of the reflected waves obtained at zero offset (self-excited) can be eliminated through linear correction.
[0028] In the embodiment of the present application, obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data includes:
[0029] Reduce the offset based on the shot point coordinates and seabed node coordinates in the four-component data;
[0030] Calculate the average amplitude value of the P component according to the P component data in the four-component data after line calibration;
[0031] According to the physical characteristics of the P component and the average amplitude value of the P component corresponding to each node, a P component polarity analysis diagram of each node is formed;
[0032] A comprehensive polarity analysis diagram of the P component is formed according to the P component polarity analysis diagrams of different nodes.
[0033] According to the above technical means, the shot points detected simultaneously by each seabed node can be separated by reducing the shot distance. Based on the calculated average amplitude value, the physical characteristics of the P component can be determined more accurately, and ultimately an accurate polar quality control map of the P component can be formed.
[0034] In the embodiment of the present application, the offset is reduced according to the shot point coordinates and the seabed node coordinates in the four-component data, including:
[0035] The shot point coordinates are calculated according to the following formula: S x -R x =(S x -R x ) / D l ×D s ; S y -R y =(S y -R y ) / D l ×D s ;
[0036] Among them, S x is the X coordinate of the shot point, S y is the Y coordinate of the shot point, R x is the X coordinate of the seabed node, R y is the Y coordinate of the seabed node, D l is the preset offset distance, D s is the distance between the gunpoints.
[0037] According to the above technical means, the shot offset is reduced with the corresponding seabed node as the center, so as to separate the shot points detected simultaneously by each seabed node.
[0038] In the embodiment of the present application, the average amplitude value of the P component is calculated based on the P component data in the four-component data after line calibration, including:
[0039] The time window of the first arrival of the P component below the preset time is extracted using the seismic trace attributes;
[0040] The average amplitude of the P component in each time window is calculated based on the P component data.
[0041] According to the above technical means, the amplitude value of the first arrival of the P component can be extracted.
[0042] In the embodiment of the present application, a P component polarity analysis diagram of each node is formed based on the physical characteristics of the P component and the average amplitude value of the P component corresponding to each node, including:
[0043] Determine whether the average amplitude value of each P component corresponding to the node is positive. If so, the corresponding P component is determined to meet the physical characteristics of the P component negative-to-positive polarity, and a polarity analysis diagram of the first color is formed. Otherwise, the corresponding P component is determined not to meet the physical characteristics of the P component negative-to-positive polarity, and a polarity analysis diagram of the second color is formed.
[0044] A polarity analysis graph corresponding to the node is formed according to the polarity analysis graphs of all P components corresponding to the node.
[0045] According to the above technical means, a P component polarity quality control diagram of each shot line can be formed according to the average amplitude value of the P component, and different polarities can be distinguished by using different colors.
[0046] In an embodiment of the present application, the method further includes:
[0047] Obtain the original four-component data collected by the seabed node;
[0048] Calculate the correlation coefficient between the P component water detection and the X, Y, Z three-component land detection before rotation based on the original four-component data;
[0049] Based on the polarity quality control diagram of the P component and the correlation coefficient between the water test of the P component and the land test of the X, Y, Z components before rotation, a polarity analysis diagram of the X, Y, Z components before rotation is formed;
[0050] The rotation effects of the X, Y, Z components before and after rotation are quality controlled based on the polarity analysis diagrams of the X, Y, Z components before and after rotation.
[0051] According to the above technical means, quality control comparison can be performed based on the three-component polarity analysis diagrams before and after rotation to improve the accuracy of quality control.
[0052] In the embodiment of the present application, the rotation effect of the X, Y, Z three-component data before and after rotation is quality-controlled based on the polarity analysis diagram of the X, Y, Z three-component data before and after rotation, including:
[0053] Compare the corresponding relationship between the polarity analysis diagram of the X, Y, and Z components of each seabed node before rotation and the polarity analysis diagram of the X, Y, and Z components after rotation to see if they are consistent. If there is any inconsistency, it is determined that the rotation angle of the corresponding node is abnormal.
[0054] According to the above technical means, the three-component polarity analysis diagram before and after rotation is used for quality control comparison, which can more accurately correspond to the seabed nodes with abnormalities during the quality control process, thereby improving the accuracy.
[0055] A second aspect of the present application provides a multi-component rotation quality control device for a seabed node based on the cross-correlation principle, the multi-component rotation quality control device for a seabed node based on the cross-correlation principle comprising:
[0056] The data acquisition unit is used to obtain the P component data collected by the seabed node and the rotated X, Y, and Z three-component data, and separate the P, X, Y, and Z four-component data;
[0057] A data preprocessing unit is used to preprocess the separated four-component data to obtain preprocessed four-component data and separate the P, X, Y, and Z four-component data;
[0058] A P component polarity analysis diagram generating unit is used to obtain a comprehensive polarity analysis diagram of the P component based on the P component data in the preprocessed four-component data;
[0059] A correlation coefficient calculation unit is used to calculate the correlation coefficients of the P component water detection and the rotated X, Y, and Z components land detection respectively;
[0060] An X, Y, Z three-component polarity analysis diagram generating unit is used to generate a polarity analysis diagram of the rotated X, Y, Z three-component based on the polarity quality control diagram of the P component and the correlation coefficient between the P component water test and the rotated X, Y, Z three-component land test;
[0061] The quality control judgment unit is used to quality control the rotation effect of the X, Y, and Z three-component data after rotation based on the polarity analysis diagram of the X, Y, and Z three-component data.
[0062] According to the above technical means, the non-directional physical characteristics of the P component received by the seabed node are utilized. First, the polarity characteristics of the P component are judged according to the amplitude value of the P component. Then, the principle of cross-correlation is used to perform cross-correlation operations on the other three components with the P component respectively. The polarities of the X, Y, and Z components are judged according to the obtained correlation coefficients. Finally, the three-component polarity analysis diagrams after rotation of all nodes are compared for quality control to meet the needs of simultaneous rotation effect quality control of multiple nodes, thereby improving the efficiency and accuracy of quality control.
[0063] The third aspect of the present application provides an electronic device, comprising: one or more processors; a memory; a module with multiple applications installed; and one or more programs, wherein the one or more programs are stored in the memory, and when the one or more programs are executed by the processor, the electronic device executes the multi-component rotation quality control method of the seabed node based on the cross-correlation principle.
[0064] A fourth aspect of the present application provides a machine-readable storage medium having stored thereon instructions for enabling a machine to execute the multi-component rotation quality control method for seabed nodes based on the cross-correlation principle.
[0065] Through the above technical solution, the principle of cross-correlation is utilized to perform cross-correlation operations on the X, Y, and Z components with the P component respectively. The polarities of the three components are determined based on the obtained correlation coefficients. The polarity comparison diagram of the three components of all nodes before and after rotation is used to quality control the rotation effect, which meets the needs of simultaneous rotation effect quality control of multiple nodes and improves the efficiency and accuracy of quality control.
[0066] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS
[0067] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following detailed description, they are used to explain the embodiments of the present invention, but do not constitute a limitation of the embodiments of the present invention. In the accompanying drawings:
[0068] FIG1 is a flow chart of a method for multi-component rotation quality control of a seabed node based on the cross-correlation principle provided by one embodiment of the present invention;
[0069] FIG2 is a schematic diagram showing the effect of limiting the offset distance processing provided by an embodiment of the present invention;
[0070] FIG3 is a schematic diagram of the effect of linear correction provided by an embodiment of the present invention;
[0071] FIG4 is a schematic diagram showing the effect of reducing the offset provided by one embodiment of the present invention;
[0072] FIG5 is a schematic diagram of the cross-correlation principle provided by one embodiment of the present invention;
[0073] FIG6 is a comparison diagram of different shot line polarities of a four-component common detection point gather provided by one embodiment of the present invention;
[0074] FIG7 is a four-component polar quality control diagram after rotation provided by an embodiment of the present invention;
[0075] FIG8 is a flow chart of a method for multi-component rotation quality control of a seabed node based on the cross-correlation principle provided by another embodiment of the present invention;
[0076] FIG9 is a block diagram of a device for multi-component rotation quality control of a seabed node based on the cross-correlation principle provided in one embodiment of the present invention. DETAILED DESCRIPTION
[0077] The following describes the specific embodiments of the present invention in detail with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the present invention and are not intended to limit the present invention.
[0078] The seabed node is a four-component acquisition device. The P component of the water detector is a pressure detector that receives the pressure wave field. The land detector is a velocity detector, which includes three components, X, Y, and Z, and receives the particle motion velocity wave field. Under horizontal layered medium conditions, the seismic response of conventional single-component seismic records is independent of the shot detection orientation. However, the multi-component node instrument receives a vector wave field. The polarity and amplitude of the seismic response recorded by the two horizontal components, X and Y, change with the shot detection orientation. Therefore, the true direction of the three-component geophone placement is extremely important information for multi-component data processing. Due to the different physical mechanisms of pressure detectors and velocity detectors, they have different responses to the received seismic wave fields. The signal received by the P-component pressure geophone is a scalar and has nothing to do with the direction. When the seismic wave field near it expands, a positive pulse is generated, and under compression, a negative pulse is generated. The purpose of the present invention is to make full use of the non-directional physical characteristics of the P component. First, the polarity characteristics of the P component are judged according to its initial amplitude value. Then, the other three components are cross-correlated with the P component according to the principle of cross-correlation. The polarities of the X, Y, and Z components are judged according to the obtained correlation coefficients. Finally, the three-component polarity analysis diagrams after all nodes are rotated are compared for quality control to meet the needs of simultaneous rotation effect quality control of multiple nodes, thereby improving the efficiency and accuracy of quality control.
[0079] FIG1 is a flow chart of a method for multi-component rotation quality control of a seabed node based on the cross-correlation principle provided by an embodiment of the present invention. As shown in FIG1 , the method for multi-component rotation quality control of a seabed node based on the cross-correlation principle includes:
[0080] S1: Obtain the P component data collected by the seabed node and the rotated X, Y, and Z component data, and separate the P, X, Y, and Z four-component data.
[0081] In this embodiment of the present application, the four-component data collected by the seafloor node includes at least the collected P, X, Y, and Z four-component data, the coordinates of the seafloor node, and the coordinates of each shot point. In this embodiment, only the collected P component data, the coordinates of the seafloor node, and the coordinates of each shot point are required. The rotated X, Y, and Z three-component data can be directly obtained from a system that implements three-component rotation correction, or the rotation-corrected three-component data can be input.
[0082] In the embodiment of the present application, the seabed node is a four-component acquisition, the water detection P component is a pressure detector, and the land detection is a velocity detector, including three components of X, Y, and Z. The four components are separated by component type.
[0083] S2: Preprocessing the separated four-component data to obtain preprocessed four-component data.
[0084] In the embodiment of the present application, the four-component data is preprocessed to obtain corrected four-component data, including:
[0085] S201: Perform offset-limiting processing on the separated data components to obtain offset-limited data for each component. In the embodiment of the present application, the offset-limiting processing primarily limits the distance between the seabed node and the shot point. Data exceeding a preset offset distance is eliminated through the offset-limiting processing. The preset offset distance is set based on demand. In one embodiment of the present application, the preset offset distance is 300 meters. The effect of the offset-limiting processing is shown in Figure 2.
[0086] S202: Linear correction is performed on each component data after constrained offset processing to obtain four-component data after linear correction. Linear correction, also known as linear moveout correction (LMO), eliminates the travel time difference between the reflection waves at different offsets and the reflection waves at zero offset (self-excited). The effect of linear correction is shown in Figure 3.
[0087] According to the above technical means, by limiting the offset distance processing, data exceeding the preset offset distance can be eliminated, eliminating the deviation caused by distance. Through linear correction, the time difference between the travel time of the reflected waves at different shot offsets and the travel time of the reflected waves obtained at zero shot offset (self-excited spontaneous) can be eliminated.
[0088] S3: Obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data.
[0089] In the embodiment of the present application, obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data includes:
[0090] S301: reducing the offset according to the shot point coordinates and the seabed node coordinates in the four-component data.
[0091] In the embodiment of the present application, the offset is reduced according to the shot point coordinates and the seabed node coordinates in the four-component data, including:
[0092] The shot point coordinates are calculated according to the following formula: S x -R x =(S x -R x ) / D l ×D s ; S y -R y =(S y -R y ) / D l ×D s ;
[0093] Among them, S xis the X coordinate of the shot point, S y is the Y coordinate of the shot point, R x is the X coordinate of the seabed node, R y is the Y coordinate of the seabed node, D l is the preset offset distance, D s is the distance between shot points. The effect of reducing the offset is shown in Figure 4.
[0094] According to the above technical means, the shot offsets are reduced to the corresponding seabed node as the center, so as to separate the shot points detected simultaneously by each seabed node.
[0095] S302: Calculate the average amplitude value of the P component according to the P component data in the four-component data after line calibration.
[0096] In the embodiment of the present application, the average amplitude value of the P component is calculated as follows:
[0097] The seismic trace attributes are used to extract a time window of a preset time below the first arrival of the P component; in one embodiment, the preset time is 10 ms.
[0098] The average amplitude of the P component in each time window is calculated based on the P component data.
[0099] According to the above technical means, the average amplitude value of each shot at each node is calculated and extracted, and the extracted average amplitude value can form a statistical text file of the average amplitude value of the P component.
[0100] S303: Based on the physical characteristics of the P component and the average amplitude value of the P component corresponding to each node, a P component polarity analysis diagram of each node is formed. Specifically:
[0101] Determine whether the average amplitude value of each P component corresponding to the node is positive. If so, the corresponding P component is determined to meet the physical characteristics of the P component negative-to-positive polarity, and a polarity analysis diagram of the first color is formed. Otherwise, the corresponding P component is determined not to meet the physical characteristics of the P component negative-to-positive polarity, and a polarity analysis diagram of the second color is formed.
[0102] A polarity analysis graph corresponding to the node is formed according to the polarity analysis graphs of all P components corresponding to the node.
[0103] According to the above technical means, a P component polarity quality control diagram of each shot line can be formed according to the average amplitude value of the P component, and different polarities can be distinguished by using different colors.
[0104] S304: forming a comprehensive polarity analysis diagram of the P component according to the P component polarity analysis diagrams of different nodes.
[0105] According to the above technical means, the shot points detected simultaneously by each seabed node can be separated by reducing the shot distance. Based on the calculated average amplitude value, the physical characteristics of the P component can be determined more accurately, and ultimately an accurate polar quality control map of the P component can be formed.
[0106] S4: Calculate the correlation coefficients between the P component water detection and the rotated X, Y, and Z components land detection respectively.
[0107] In the embodiment of the present application, the cross-correlation principle is adopted by calculating the correlation coefficient, that is, in signal processing, cross-correlation is used to represent the similarity between two signals f(x) and g(x). It is a function of the two signals directly relative to time, sometimes also called "sliding dot product", and the principle is shown in Figure 5.
[0108] Correlation coefficient: When the correlation coefficient is closer to +1 or -1, it indicates a positive correlation (+1) or negative correlation (-1) between the arrays; when the correlation coefficient is close to 0, it indicates no or weak correlation. Semblance: The correlation coefficient between two input data within a specified time window. The calculation formula of the correlation coefficient is:
[0109] Right now
[0110] Where Cov(X,Y) represents the covariance, and D(x) and D(Y) represent the variance.
[0111] In the embodiment of the present application, a time window range of 10ms below the first arrival is selected, and the P component is cross-correlated with the X, Y, and Z components respectively. According to the above formula, the correlation coefficients of PX, PY, and PZ are obtained.
[0112] S5: A polarity analysis diagram of the rotated X, Y, and Z components is formed based on the polarity quality control diagram of the P component, the water inspection of the P component, and the land inspection of the rotated X, Y, and Z components.
[0113] In the embodiment of the present application, a polarity analysis diagram of the rotated X, Y, and Z components is formed based on the polarity quality control diagram of the P component and the correlation coefficient between the water detection of the P component and the land detection of the rotated X, Y, and Z components, including:
[0114] The polarity of each component is determined based on the correlation coefficient between the P component water test and the rotated X, Y, and Z components land test:
[0115] If the correlation coefficient is greater than or equal to the first preset value, it is determined that the P component is positively correlated with the corresponding component, and the polarity of the corresponding component is consistent with the polarity of the P component;
[0116] If the correlation coefficient is less than or equal to the second preset value, it is determined that the P component is negatively correlated with the corresponding component, and the polarity of the corresponding component is opposite to that of the P component;
[0117] If the correlation coefficient is greater than the second preset value and less than the first preset value, it is determined that the P component has no correlation with the corresponding component;
[0118] A polarity analysis diagram of the corresponding component is formed according to the polarity of each component.
[0119] In practical applications, different colors are used to represent different polarities. For example, in the aforementioned scheme, the first color is used to represent the negative jump to positive polarity of the P component, and the second color is used to represent the positive jump to negative polarity with the opposite polarity of the P component. Therefore, when the correlation coefficient is greater than or equal to the first preset value, it is determined that the P component is positively correlated with the corresponding component, that is, the polarity color of the corresponding node is consistent with the polarity color of the P component; when the correlation coefficient is less than or equal to the second preset value, it is determined that the P component is negatively correlated with the corresponding component, that is, the polarity color of the corresponding node is opposite to the polarity color of the P component; that is, when the P component is the first color, the corresponding node uses the second color, and when the P component is the second color, the corresponding node uses the first color; if the correlation coefficient is greater than the second preset value and less than the first preset value, it is determined that the P component has no correlation with the corresponding component, and the third color is introduced for representation.
[0120] According to the above technical means, the polarity of the X, Y, and Z component data can be determined based on the correlation coefficient between the P component and the rotated X, Y, and Z component data, thereby determining the polarity analysis diagram of the corresponding component according to the polarity of each component, providing a data basis for quality control.
[0121] S6: Quality control the rotation effect of the rotated X, Y, and Z three-component data based on the polarity analysis diagram of the rotated X, Y, and Z three-component data.
[0122] In the embodiment of the present application, the rotation effect of the rotated X, Y, and Z three-component data is quality-controlled based on the polarity analysis diagram of the rotated X, Y, and Z three-component data, including:
[0123] Traverse each node in the polarity analysis diagram of each component to determine whether the polarity of each node is consistent. If there is an abnormal polarity, it is determined that the corresponding component has an abnormal rotation angle, that is, the rotation is not in place or the rotation is wrong.
[0124] The above method can determine whether the corresponding component has an abnormal rotation angle. In order to more quickly determine the abnormal nodes, during the verification process, a small number of inconsistent nodes are prioritized for verification.
[0125] According to the above technical means, it is possible to determine whether there is an abnormality in the corresponding polarity based on the polarity analysis diagram of the X, Y, and Z components, and to quickly perform polarity quality control.
[0126] As shown in Figure 6, the polarity of the P component is negative-to-positive. By extracting the amplitude value, the polarity quality control of the P component is shown in Figure 7. The polarities of the three components X, Y, and Z can be determined based on the correlation coefficients of PX, PY, and PZ. As shown in Figure 7, the Z component is positive-to-negative, and the correlation coefficient between the P and Z components is -1, indicating that the polarity of the P component and the polarity of the Z component are negatively correlated. The polarity quality control chart is shown in Figure 7. As shown in the X component polarity quality control chart in Figure 7, the correlation coefficient of the upper left portion of each seafloor node in the X component is -1, and the correlation coefficient of the lower right portion is 1, indicating that the polarity of the upper left portion is positive-to-negative, while the polarity of the lower right portion is negative-to-positive. As shown in the Y component polarity quality control chart in Figure 7, the correlation coefficient of the upper right portion of the Y component is -1, and the correlation coefficient of the lower left portion is 1, indicating that the polarity of the upper right portion is positive-to-negative, while the polarity of the lower left portion is negative-to-positive.
[0127] Due to the difference in receiving direction and incident angle, as shown in Figure 7, the polarity distribution of the X component of each seabed node shows symmetrical opposite polarity along the vertical survey line (CrossLine), and the polarity distribution of the Y component shows symmetrical opposite polarity along the survey line (InLine). The node can intuitively control the accuracy of rotation based on the polarity QC charts of the X, Y, and Z components. As shown in Figure 7, the X component polarity QC chart and the Y component polarity QC chart can be used to judge whether the polarity of the rotated node is symmetrical and opposite along the CrossLine direction or the InLine direction. At the circled node in the figure, there are deviations in the polarity of both the X and Y components, which can be judged as abnormal.
[0128] Through the above technical solution, the principle of cross-correlation is utilized to perform cross-correlation operations on the X, Y, and Z components with the P component respectively. The polarities of the three components are determined based on the obtained correlation coefficients. The polarity comparison diagram of the three components of all nodes before and after rotation is used to quality control the rotation effect, which meets the needs of simultaneous rotation effect quality control of multiple nodes and improves the efficiency and accuracy of quality control.
[0129] FIG8 is a flow chart of a method for quality control of multi-component rotation of a seabed node based on the cross-correlation principle according to an embodiment of the present invention. As shown in FIG8 , the method for quality control of multi-component rotation of a seabed node based on the cross-correlation principle includes:
[0130] S1: Obtain the original four-component data collected by the seabed node and the rotated X, Y, and Z three-component data, and separate the P, X, Y, and Z four components.
[0131] S2: Preprocessing the separated four-component data to obtain preprocessed four-component data.
[0132] S3: Obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data.
[0133] S4: Calculate the correlation coefficient between the P component water detection and the X, Y, Z three-component land detection after rotation and the correlation coefficient between the P component water detection and the X, Y, Z three-component land detection before rotation based on the four-component data.
[0134] S5: forming a polarity analysis diagram of the rotated X, Y, and Z components based on the polarity quality control diagram of the P component and the correlation coefficient between the water test of the P component and the land test of the rotated X, Y, and Z components;
[0135] S6: forming a polarity analysis diagram of the X, Y, and Z components before rotation based on the polarity quality control diagram of the P component and the correlation coefficient between the water detection of the P component and the land detection of the X, Y, and Z components before rotation;
[0136] S7: Quality control the rotation effect of the X, Y, Z three-component data before and after rotation based on the polarity analysis diagram of the X, Y, Z three-component data before and after rotation.
[0137] According to the above technical means, quality control comparison can be performed based on the three-component polarity analysis diagrams before and after rotation to improve the accuracy of quality control.
[0138] In the embodiment of the present application, the rotation effect of the X, Y, Z three-component data before and after rotation is quality-controlled based on the polarity analysis diagram of the X, Y, Z three-component data before and after rotation, including:
[0139] Compare the corresponding relationship between the polarity analysis diagram of the X, Y, and Z components of each seabed node before rotation and the polarity analysis diagram of the X, Y, and Z components after rotation to see if they are consistent. If there is any inconsistency, it is determined that the rotation angle of the corresponding node is abnormal.
[0140] According to the above technical means, the three-component polarity analysis diagram before and after rotation is used for quality control comparison, which can more accurately correspond to the seabed nodes with abnormalities during the quality control process, thereby improving the accuracy.
[0141] A second aspect of the present application provides a multi-component rotation quality control device for a seabed node based on the cross-correlation principle, as shown in FIG9 . The multi-component rotation quality control device for a seabed node based on the cross-correlation principle includes:
[0142] The data acquisition unit is used to obtain the P component data collected by the seabed node and the rotated X, Y, and Z three-component data, and separate the P, X, Y, and Z four-component data;
[0143] A data preprocessing unit, configured to preprocess the separated four-component data to obtain preprocessed four-component data;
[0144] A P component polarity analysis diagram generating unit is used to obtain a comprehensive polarity analysis diagram of the P component based on the P component data in the preprocessed four-component data;
[0145] A correlation coefficient calculation unit is used to calculate the correlation coefficients of the P component water detection and the rotated X, Y, and Z components land detection respectively;
[0146] An X, Y, Z three-component polarity analysis diagram generating unit is used to generate a polarity analysis diagram of the rotated X, Y, Z three-component based on the polarity quality control diagram of the P component and the correlation coefficient between the P component water test and the rotated X, Y, Z three-component land test;
[0147] The quality control judgment unit is used to quality control the rotation effect of the X, Y, and Z three-component data after rotation based on the polarity analysis diagram of the X, Y, and Z three-component data.
[0148] According to the above technical means, the non-directional physical characteristics of the P component received by the seabed node are utilized. First, the polarity characteristics of the P component are judged according to the amplitude value of the P component. Then, the principle of cross-correlation is used to perform cross-correlation operations on the other three components with the P component respectively. The polarities of the X, Y, and Z components are judged according to the obtained correlation coefficients. Finally, the three-component polarity analysis diagrams after rotation of all nodes are compared for quality control to meet the needs of simultaneous rotation effect quality control of multiple nodes, thereby improving the efficiency and accuracy of quality control.
[0149] The third aspect of the present application provides an electronic device, comprising: one or more processors; a memory; a module with multiple applications installed; and one or more programs, wherein the one or more programs are stored in the memory, and when the one or more programs are executed by the processor, the electronic device executes the multi-component rotation quality control method of the seabed node based on the cross-correlation principle.
[0150] A fourth aspect of the present application provides a machine-readable storage medium having stored thereon instructions for enabling a machine to execute the multi-component rotation quality control method for seabed nodes based on the cross-correlation principle.
[0151] Through the above technical solution, the principle of cross-correlation is utilized to perform cross-correlation operations on the X, Y, and Z components with the P component respectively. The polarities of the three components are determined based on the obtained correlation coefficients. The polarity comparison diagram of the three components of all nodes before and after rotation is used to quality control the rotation effect, which meets the needs of simultaneous rotation effect quality control of multiple nodes and improves the efficiency and accuracy of quality control.
[0152] Those skilled in the art will appreciate that all or part of the steps in the methods of the aforementioned embodiments can be accomplished by instructing the relevant hardware through a program, which is stored in a storage medium and includes a number of instructions for causing a single-chip microcomputer, chip, or processor to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.
[0153] The above describes in detail the optional embodiments of the present invention in conjunction with the accompanying drawings. However, the embodiments of the present invention are not limited to the specific details in the above embodiments. Within the technical concept of the embodiments of the present invention, a variety of simple modifications can be made to the technical solutions of the embodiments of the present invention, and these simple modifications all fall within the scope of protection of the embodiments of the present invention. It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner unless there is any contradiction. In order to avoid unnecessary repetition, the embodiments of the present invention will no longer describe the various possible combinations separately.
[0154] In addition, the various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the embodiments of the present invention, they should also be regarded as the contents disclosed in the embodiments of the present invention.
Claims
1. A method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle, characterized in that: The method comprises: Obtain the P component data collected by the seabed node and the rotated X, Y, and Z three-component data, and separate the P, X, Y, and Z four-component data; Preprocessing the separated four-component data to obtain preprocessed four-component data; Obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data; Calculate the correlation coefficients of the P component water detection and the rotated X, Y, and Z components land detection respectively; According to the polarity quality control diagram of the P component and the correlation coefficient between the water inspection of the P component and the land inspection of the rotated X, Y, Z components, a polarity analysis diagram of the rotated X, Y, Z components is formed; The rotation effect of the X, Y, Z three-component data after rotation is quality controlled based on the polarity analysis diagram of the X, Y, Z three-component data after rotation.
2. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 1 is characterized in that: According to the polarity quality control diagram of the P component and the correlation coefficient between the water inspection of the P component and the land inspection of the rotated X, Y, Z components, a polarity analysis diagram of the rotated X, Y, Z components is formed, including: The polarity of each component is determined based on the correlation coefficient between the P component water test and the rotated X, Y, and Z components land test: If the correlation coefficient is greater than or equal to the first preset value, it is determined that the P component is positively correlated with the corresponding component, and the polarity of the corresponding component is consistent with the polarity of the P component; If the correlation coefficient is less than or equal to the second preset value, it is determined that the P component is negatively correlated with the corresponding component, and the polarity of the corresponding component is opposite to that of the P component; If the correlation coefficient is greater than the second preset value and less than the first preset value, it is determined that the P component has no correlation with the corresponding component; A polarity analysis diagram of the corresponding component is formed according to the polarity of each component.
3. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 1 is characterized in that: The rotation effect of the X, Y, and Z components after rotation is quality controlled based on the polarity analysis diagram of the X, Y, and Z components after rotation, including: Traverse each node in the polarity analysis diagram of each component to determine whether the polarity of each node is consistent. If there is an abnormal polarity, it is determined that the corresponding component has an abnormal rotation angle.
4. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 1 is characterized in that: The separated four-component data is preprocessed to obtain preprocessed four-component data, including: The data of each component after separation are processed by limiting the offset distance to obtain the data of each component after limiting the offset; Linear correction is performed on each component data after the restricted offset processing to obtain the four-component data after linear correction.
5. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 4 is characterized in that: According to the P component data in the preprocessed four-component data, a comprehensive polarity analysis diagram of the P component is obtained, including: Reduce the offset according to the shot point coordinates and seabed node coordinates in the four-component data; Calculate the average amplitude value of the P component according to the P component data in the four-component data after line calibration; According to the physical characteristics of the P component and the average amplitude value of the P component corresponding to each node, a P component polarity analysis diagram of each node is formed; A comprehensive polarity analysis diagram of the P component is formed according to the polarity analysis diagrams of the P components of different nodes.
6. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 5 is characterized in that: Reduce the offset based on the shot point coordinates and seabed node coordinates in the four-component data, including: The shot point coordinates are calculated according to the following formula: S x -R x =(S x -R x ) / D l ×D s ; S y -R y =(S y -R y ) / D l ×D s ; Among them, S x is the X coordinate of the shot point, S y is the Y coordinate of the shot point, R x is the X coordinate of the seafloor node, R y is the Y coordinate of the seafloor node, D l is the preset offset distance, D s is the distance between the gunpoints.
7. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 5 is characterized in that: The average amplitude value of the P component is calculated based on the P component data in the four-component data after line calibration, including: The time window of the first arrival of the P component below the preset time is extracted using the seismic trace attributes; The average amplitude value of the P component in each time window is calculated based on the P component data.
8. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 5 is characterized in that: According to the physical characteristics of the P component and the average amplitude value of the P component corresponding to each node, a P component polarity analysis diagram of each node is formed, including: Determine whether the average amplitude value of each P component corresponding to the node is positive. If so, determine that the corresponding P component meets the physical characteristics of the P component negative jump positive polarity, and form a polarity analysis diagram of the first color; otherwise, determine that the corresponding P component does not meet the physical characteristics of the P component negative jump positive polarity, and form a polarity analysis diagram of the second color; According to the polarity analysis diagrams of all P components corresponding to the node, a polarity analysis diagram corresponding to the node is formed.
9. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 1 is characterized in that: The method further comprises: Obtain the original four-component data collected by the seabed nodes; The correlation coefficient between the P component water detection and the X, Y, and Z components land detection before rotation is calculated based on the original four-component data; Based on the polarity quality control diagram of the P component and the correlation coefficient between the water inspection of the P component and the land inspection of the X, Y, Z components before rotation, a polarity analysis diagram of the X, Y, Z components before rotation is formed; The rotation effect of the X, Y, Z three-component data before and after the rotation is quality controlled based on the polarity analysis diagram of the X, Y, Z three-component data before and after the rotation.
10. The method for multi-component rotation quality control of seabed nodes based on the cross-correlation principle according to claim 9, characterized in that: The rotation effect of the X, Y, Z three-component data before and after rotation is quality controlled based on the polarity analysis diagram of the X, Y, Z three-component data before and after rotation, including: Compare the corresponding relationship between the polarity analysis diagram of the X, Y, and Z components of each seabed node before rotation and the polarity analysis diagram of the X, Y, and Z components after rotation to see if they are consistent. If there is any inconsistency, it is determined that the rotation angle of the corresponding node is abnormal.
11. A multi-component rotation quality control device for seabed nodes based on the cross-correlation principle, characterized in that: The multi-component rotation quality control device for seabed nodes based on the cross-correlation principle comprises: A data acquisition unit is used to acquire the P component data collected by the seabed node and the rotated X, Y, and Z three-component data, and separate the P, X, Y, and Z four-component data; A data preprocessing unit, used for preprocessing the separated four-component data to obtain preprocessed four-component data; A P component polarity analysis diagram generating unit, used for obtaining a comprehensive polarity analysis diagram of the P component according to the P component data in the preprocessed four-component data; A correlation coefficient calculation unit is used to calculate the correlation coefficients of the P component water detection and the rotated X, Y, and Z three-component land detection respectively; The X, Y, Z three-component polarity analysis diagram generating unit is used to form the polarity analysis diagram of the rotated X, Y, Z three-component according to the polarity analysis diagram of the P component and the correlation coefficient between the water detection of the P component and the land detection of the rotated X, Y, Z three-component; The quality control judgment unit is used to quality control the rotation effect of the X, Y, and Z three-component data after rotation according to the polarity analysis diagram of the X, Y, and Z three-component data.
12. An electronic device, characterized in that: include: one or more processors; Memory; a module with multiple applications installed; And one or more programs, wherein the one or more programs are stored in the memory, and when the one or more programs are executed by the processor, the electronic device executes the multi-component rotation quality control method of the seabed node based on the cross-correlation principle as described in any one of claims 1-10.
13. A machine-readable storage medium having instructions stored thereon, the instructions being used to enable a machine to execute the multi-component rotation quality control method for seafloor nodes based on the cross-correlation principle as described in any one of claims 1 to 10.
Citation Information
Patent Citations
Three-component directional correction method of submarine node seismograph
CN110687606A
Combined type ocean bottom seismograph attitude automatic correction device and method
CN111257940A
Underwater attitude orientation method for seabed node detector
CN112147695A
Method and device for determining horizontal azimuth angle of seabed node
CN112394415A
Seabed node secondary positioning method and device
CN112415595A