Method, device and equipment for determining measurement rule of three-component detector

Abstract

The invention discloses a measurement rule determination method, device and equipment for a three-component detector, and belongs to the technical field of seismic exploration. The method comprises the following steps: acquiring first multi-component data; determining a first time shifting radial component of each of the plurality of shot points; determining a first time shifting radial component model channel of each of a plurality of receiving points in the plurality of shot points; determining respective first correlation coefficients of the plurality of shot points; determining respective second time shifting radial components of the plurality of shot points; determining a second time-shifting radial component model channel of each of a plurality of receiving points in the plurality of shot points; determining respective second correlation coefficients of the plurality of shot points; and determining a measurement rule of the three-component detector according to the respective first correlation coefficients of the plurality of shot points and the respective second correlation coefficients of the plurality of shot points. According to the method, the left and right hand rules followed by the three-component geophone can be accurately determined by comparing the difference between the time-shifted radial component model channels and the time-shifted radial component model channels under the left hand rule and the right hand rule.

Description

Technical Field

[0001] This application relates to the field of seismic exploration technology, and in particular to a method, apparatus and equipment for determining the measurement rules of a three-component geophone. Background Technology

[0002] With the deepening of oil and gas exploration, multi-wave seismic exploration has received increasing attention. Among these methods, three-component geophones are widely used to acquire multi-component data. This multi-component data is used to record various types of wavefields, including P-waves and S-waves. Specifically, S-waves are projected onto the horizontal components in the X and Y directions for processing.

[0003] When processing the horizontal components of multi-component data, in order to extract transverse wave information, the horizontal components in the X and Y directions need to be rotated to the radial and tangential directions to obtain the radial and tangential components. There is a prerequisite for this rotation process: the sign of the rotation angle must be consistent with the measurement rules followed by the three-component detector. In other words, it is necessary to determine whether the three-component detector used follows a left-handed or right-handed system of reference to ensure the correctness of the rotation operation.

[0004] However, in practical multi-component data processing, the left-hand or right-hand system information of the three-component detector used during acquisition is often missing or inconsistent with the actual situation. This lack of information or inconsistency brings great difficulties to subsequent multi-component data processing. Therefore, further research is needed on the measurement rules of the three-component detector for the acquired multi-component data. Summary of the Invention

[0005] This application provides a method, apparatus, and device. The technical solutions provided by this application are as follows:

[0006] According to one aspect of the embodiments of this application, a method for determining the measurement rules of a three-component detector is provided, the method comprising:

[0007] First multi-component data is acquired. The first multi-component data is determined based on a three-component geophone according to measurement rules. The measurement rules are used to instruct the coordinate system criteria for the three-component geophone to record the multi-component data of seismic waves. The measurement rules include left-hand rule and right-hand rule. The first multi-component data includes the horizontal components of multiple shot-receiver points. The shot-receiver points are determined based on the shot point and the receiver point. The shot point is used to emit seismic waves, and the receiver point is used to receive the seismic waves emitted by the shot point.

[0008] Based on the horizontal components of the plurality of shot receivers, a first time-shifted radial component is determined for each of the plurality of shot receivers. The first time-shifted radial component is determined by time-shifting the first radial component in the time domain. The first radial component is determined based on the horizontal component under the right-hand rule.

[0009] Based on the first time-shift radial components of each of the plurality of shot receivers, a first time-shift radial component model trace for each of the plurality of shot receivers is determined. The first time-shift radial component model trace for each receiver is determined based on the plurality of first time-shift radial components included in the receiver.

[0010] Based on the first time-shift radial component of each of the plurality of shot receivers and the first time-shift radial component model trace of each of the plurality of receivers, a first correlation coefficient is determined for each of the plurality of shot receivers. The first correlation coefficient refers to the correlation coefficient between the first time-shift radial component and the first time-shift radial component model trace.

[0011] Based on the horizontal components of each of the plurality of shot receivers, a second time-shifted radial component is determined for each of the plurality of shot receivers. The second time-shifted radial component is determined by time-shifting the second radial component in the time domain. The second radial component is determined based on the horizontal component under the left-hand rule.

[0012] Based on the second time-shift radial components of each of the plurality of shot receivers, a second time-shift radial component model trace for each of the plurality of shot receivers is determined. The second time-shift radial component model trace for each receiver is determined based on the plurality of second time-shift radial components included in the receiver.

[0013] Based on the second time-shift radial component of each of the plurality of shot receivers and the second time-shift radial component model trace of each of the plurality of receivers, a second correlation coefficient is determined for each of the plurality of shot receivers. The second correlation coefficient refers to the correlation coefficient between the second time-shift radial component and the second time-shift radial component model trace.

[0014] The measurement rules of the three-component detector are determined based on the first correlation coefficient of each of the plurality of shot detectors and the second correlation coefficient of each of the plurality of shot detectors.

[0015] In some embodiments, determining the first time-shift radial component model trace of each of the plurality of receiver points among the plurality of shot receiver points based on the first time-shift radial component of each of the plurality of shot receiver points includes:

[0016] For each of the plurality of shot receivers, determine the plurality of first time-shift radial components included in the receiving point of the shot receiver;

[0017] The sum of the multiple first time-shift radial components included in the receiving point is determined as the first time-shift radial component model channel of the receiving point;

[0018] The step of determining the second time-shift radial component model trace of each of the multiple shot receivers based on the second time-shift radial component of each of the multiple shot receivers includes:

[0019] For each of the plurality of shot receivers, determine the plurality of second time-shift radial components included in the receiving point of the shot receiver;

[0020] The sum of the multiple second time-shift radial components included in the receiving point is determined as the second time-shift radial component model channel of the receiving point.

[0021] According to one aspect of the embodiments of this application, a measurement rule determination device for a three-component detector is provided, the device comprising:

[0022] The acquisition module is used to acquire first multi-component data, which is determined based on a three-component geophone according to measurement rules. The measurement rules are used to instruct the coordinate system criteria for the three-component geophone to record the multi-component data of seismic waves. The measurement rules include left-hand rule and right-hand rule. The first multi-component data includes the horizontal components of multiple shot-receiver points. The shot-receiver points are determined based on the shot point and the receiver point. The shot point is used to emit seismic waves, and the receiver point is used to receive the seismic waves emitted by the shot point.

[0023] The first determining module is used to determine the first time-shift radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers. The first time-shift radial component is determined by time-shifting the first radial component in the time domain. The first radial component is determined based on the horizontal component under the right-hand rule.

[0024] The first determining module is further configured to determine the first time-shift radial component model trace of each of the plurality of gun receivers based on the first time-shift radial component of each of the plurality of gun receivers, wherein the first time-shift radial component model trace of the receiving point is determined based on the plurality of first time-shift radial components included in the receiving point;

[0025] The first determining module is further configured to determine a first correlation coefficient for each of the plurality of shot receivers based on the first time-shift radial component of each of the plurality of shot receivers and the first time-shift radial component model trace of each of the plurality of receivers. The first correlation coefficient refers to the correlation coefficient between the first time-shift radial component and the first time-shift radial component model trace.

[0026] The second determining module is used to determine the second time-shift radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers. The second time-shift radial component is determined by time-shifting the second radial component in the time domain. The second radial component is determined based on the horizontal component under the left-hand rule.

[0027] The second determining module is further configured to determine the second time-shift radial component model trace of each of the plurality of gun receivers based on the second time-shift radial component of each of the plurality of gun receivers, wherein the second time-shift radial component model trace of the receiving point is determined based on the plurality of second time-shift radial components included in the receiving point;

[0028] The second determining module is further configured to determine the second correlation coefficient of each of the plurality of shot receivers based on the second time-shift radial component of each of the plurality of shot receivers and the second time-shift radial component model trace of each of the plurality of receivers. The second correlation coefficient refers to the correlation coefficient between the second time-shift radial component and the second time-shift radial component model trace.

[0029] The third determining module is used to determine the measurement rules of the three-component detector based on the first correlation coefficient of each of the plurality of shot detectors and the second correlation coefficient of each of the plurality of shot detectors.

[0030] According to one aspect of the embodiments of this application, a computer device is provided, the computer device including a processor and a memory, the memory storing a computer program, the computer program being loaded and executed by the processor to implement the above-described method for determining the measurement rules of a three-component detector.

[0031] According to one aspect of the embodiments of this application, a computer-readable storage medium is provided, wherein a computer program is stored in the storage medium, the computer program being loaded and executed by a processor to implement the above-described method for determining the measurement rules of a three-component detector.

[0032] According to one aspect of the embodiments of this application, a computer program product is provided, the computer program product including a computer program stored in a computer-readable storage medium, and a processor reading from the computer-readable storage medium and executing the computer program to implement the above-described method for determining the measurement rules of a three-component detector.

[0033] The technical solutions provided in this application have at least the following beneficial effects:

[0034] Based on the different ways in which the horizontal component of the three-component detector rotates to the radial component under the left-hand rule and the right-hand rule, the left-hand and right-hand rules followed by the three-component detector can be accurately determined by comparing the differences between the time-shifted radial component and the time-shifted radial component model trace under the left-hand rule and the right-hand rule. Attached Figure Description

[0035] Figure 1 This is a schematic diagram of the implementation environment of a solution provided in one embodiment of this application;

[0036] Figure 2 This is a flowchart of a method for determining the measurement rules of a three-component detector according to an embodiment of this application;

[0037] Figure 3 This is a schematic diagram of the left-handed kinship rule (left figure) and the right-handed kinship rule (right figure) provided in one embodiment of this application;

[0038] Figure 4 This is a schematic diagram of the horizontal components of a three-component seismic record with right-handed and left-handed systems provided in one embodiment of this application;

[0039] Figure 5 This is a schematic diagram of the R and T components of a three-component seismic record in a right-handed and left-handed system provided in one embodiment of this application;

[0040] Figure 6 This is a block diagram of a measurement rule determination device for a three-component detector provided in one embodiment of this application;

[0041] Figure 7 This is a structural block diagram of a computer device provided in one embodiment of this application. Detailed Implementation

[0042] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0043] Please refer to Figure 1 The diagram illustrates an implementation environment for a solution provided in one embodiment of this application. This implementation environment may include: a terminal device 10 and a server 20.

[0044] Terminal devices 10 include, but are not limited to, PCs (Personal Computers), host computers of cloud computing technology platforms, mobile phones, tablet computers, smart voice interaction devices, game consoles, wearable devices, multimedia playback devices, vehicle terminals, smart home appliances, AR (Augmented Reality) devices, VR (Virtual Reality) devices, and other electronic devices.

[0045] Server 20 can be an independent physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN (Content Delivery Network), and big data and artificial intelligence platforms, but it is not limited to these.

[0046] Terminal device 10 and server 20 can communicate with each other via a network. This network can be a wired network or a wireless network.

[0047] Please refer to Figure 2 This document illustrates a flowchart of a method for determining measurement rules for a three-component detector according to an embodiment of this application. The execution entity for each step of this method can be a computer device; for example, the computer device could be... Figure 1 Terminal device 10 in the middle can also be Figure 1 The method may include at least one of the following steps 210 to 280.

[0048] Step 210: Acquire the first multi-component data. The first multi-component data is determined based on the three-component geophone according to the measurement rules. The measurement rules are used to instruct the coordinate system criteria for the three-component geophone to record the multi-component data of seismic waves. The measurement rules include the left-hand rule and the right-hand rule. The first multi-component data includes the horizontal components of multiple shot-receiver points. The shot-receiver points are determined based on the shot point and the receiver point. The shot point is used to emit seismic waves, and the receiver point is used to receive the seismic waves emitted by the shot point.

[0049] A three-component geophone is a device capable of simultaneously measuring the vibration of seismic waves in three vertical directions. It typically includes sensors in two directions on a horizontal plane (x-axis and y-axis) and a vertical sensor (z-axis) to acquire multi-component data of the seismic waves. This multi-component data can be used to analyze the various components of the seismic waves and for further geological exploration analysis.

[0050] In seismic exploration, multi-component data refers to data that simultaneously records the vibration of seismic waves in multiple directions (which can be three orthogonal directions: x-axis, y-axis, and z-axis). This data comprehensively reflects the propagation characteristics of seismic waves and information about underground structures. The horizontal component is a part of the multi-component data, specifically referring to the vibration data in the x-axis and y-axis directions. It records the motion of seismic waves on the horizontal plane and is used to analyze the shear waves of seismic waves.

[0051] In some embodiments, the measurement rules for a three-component geophone specify how the geophone records seismic wave vibration data in three-dimensional space. The measurement rules include left-handed and right-handed measurement rules. For the left-handed rule, in a spatial rectangular coordinate system, such as... Figure 2 As shown, three mutually perpendicular number axes are drawn through a fixed point O. These axes all have O as their origin and generally share the same unit of length. These three axes are called the x-axis (horizontal axis), y-axis (vertical axis), and z-axis (vertical axis), collectively known as coordinate axes. The x-axis and y-axis are usually placed on a horizontal plane, while the z-axis is a vertical line. If your left thumb points in the positive direction of the x-axis, your index finger in the positive direction of the y-axis, and your middle finger in the positive direction of the z-axis, then this coordinate system is called a left-handed rectangular coordinate system. Similarly, if your right thumb points in the positive direction of the x-axis, your index finger in the positive direction of the y-axis, and your middle finger in the positive direction of the z-axis, then this coordinate system is called a right-handed rectangular coordinate system. For the left-hand rule, if you grasp the axis of rotation with your left hand, and your thumb points in the positive direction of the axis of rotation, the rotation angle is positive if the direction of rotation is the same as the direction of your bent four fingers; otherwise, it is negative. For the right-hand rule, hold the axis of rotation with your right hand, with your thumb pointing in the positive direction of the axis of rotation. If the direction of rotation is the same as the direction of your bent four fingers, the rotation angle is positive; otherwise, it is negative.

[0052] A shot-receiver point consists of a shot point and a receiver point, also known as a geophone point. The shot point is used to emit seismic waves during seismic exploration. Seismic waves are typically generated at the shot point by an explosion or a seismic source (such as a seismograph) and propagate underground. These seismic waves will propagate through underground structures and reflect back to the surface. The receiver point is used to receive the seismic waves reflected from underground. The receiver point is typically equipped with a geophone (i.e., a three-component geophone) to record the data of the seismic waves returning from underground. This data is used to analyze underground structures and geological features.

[0053] In some embodiments, after acquiring the first multi-component data, a data preprocessing operation is further performed on the first multi-component data. Specifically, the x', y', and z' component data in the acquired coordinate system are rotated to the north, east, and vertical component data in the geodetic coordinate system, as shown in Formula 1 below:

[0054]

[0055] in, c θ =cos(θ), s θ =sin(θ), c ψ =cos(ψ), s ψ =sin(ψ), c Φ =cos(Φ), s Φ=sin(Φ). q, f, and y are three Euler angles describing the placement direction of the three-component detector, where q is the detector tilt angle, f is the detector roll angle, y is the detector azimuth angle, and t is the time axis.

[0056] The above method, by rotating the x', y', and z' component data in the collected coordinate system to the north, east, and vertical component data in the geodetic coordinate system, aims to transform the data into a standard, unified coordinate system for more accurate analysis and interpretation.

[0057] Step 220: Determine the first time-shift radial component of each of the multiple shot receivers based on their respective horizontal components. The first time-shift radial component is determined by time-shifting the first radial component in the time domain. The first radial component is determined based on the horizontal component under the right-hand rule.

[0058] In some embodiments, there is a one-to-one correspondence between shot receiver points and horizontal components. That is, one shot receiver point corresponds to one horizontal component, and different shot receiver points correspond to different horizontal components.

[0059] In some embodiments, step 220 above can be implemented as steps 221 to 223.

[0060] Step 221: Determine the first radial component of each of the multiple shot receivers based on their respective horizontal components using the right-hand rotation formula.

[0061] In some embodiments, the step of determining the first radial component of each of the plurality of shot receiver points may include: obtaining the shot point coordinates and receiver coordinates of each of the plurality of shot receiver points; determining the shot receiver direction angle of each of the plurality of shot receiver points based on the shot point coordinates and receiver coordinates of each of the plurality of shot receiver points; and determining the first radial component of each of the plurality of shot receiver points based on the horizontal component of each of the plurality of shot receiver points and the shot receiver direction angle of each of the plurality of shot receiver points using a right-hand rotation formula.

[0062] Shot coordinates refer to the coordinates of the source point, which can be represented as (Xsi, Ysi, Zsi). Here, Xsi, Ysi, and Zsi represent the x, y, and z coordinates of the shot point in three-dimensional space, respectively. Receiver coordinates refer to the coordinates of the seismic wave receiver, which can be represented as (Xgj, Ygj, Zgj). Here, Xgj, Ygj, and Zgj represent the x, y, and z coordinates of the receiver point in three-dimensional space, respectively.

[0063] The shot-receiver azimuth angle is the angle between the direction from the shot point to the receiver point and a reference direction (usually true north). It can be determined by calculating the coordinates of the shot point and the receiver point. The formula for calculating the shot-receiver azimuth angle is shown in Formula 2 below:

[0064]

[0065] Where i is the shot vantage point sequence number, i = 1, 2, ..., N, where N is the total number of shot vantage points, and j is the receiver vantage point sequence number, j = 1, 2, ..., M, where M is the total number of receiver vantage points. α ij This is the firing and receiving azimuth angle corresponding to the i-th firing point and the j-th receiving point. It can be understood that there is a one-to-one correspondence between firing and receiving points and firing and receiving azimuth angles. One firing and receiving point corresponds to one firing and receiving azimuth angle. Different firing and receiving points correspond to different firing and receiving azimuth angles.

[0066] In some embodiments, the right-hand rotation formula is used based on the horizontal components x of each of the multiple shot receiver points. ij y ij and the firing direction angle α ij The first radial component after the horizontal component is rotated is obtained. and the first tangential component The right-hand rotation formula is shown in Formula 3 below:

[0067]

[0068] Step 222: Obtain the first arrival time of each of the multiple shot receivers. The first arrival time refers to the earliest time when the receiver in the shot receiver receives the seismic wave emitted by the shot.

[0069] It is understandable that there is a one-to-one correspondence between shot-receiver points and first arrival times. That is, one shot-receiver point corresponds to one first arrival time, and different shot-receiver points correspond to different first arrival times. For example, in a seismic exploration project, three shot-receiver points are set up: Shot-receiver point 1 includes shot point A1 and receiver point B1; Shot-receiver point 2 includes shot point A2 and receiver point B2; and Shot-receiver point 3 includes shot point A3 and receiver point B3. The first arrival times corresponding to the three shot-receiver points are S1, S2, and S3, respectively. The first arrival time of each shot-receiver point is recorded as follows: The first arrival time of shot-receiver point 1 is S1, which indicates the earliest time that receiver point B1 receives the seismic wave emitted from shot point A1. The first arrival time of shot-receiver point 2 is S2, which indicates the earliest time that receiver point B2 receives the seismic wave emitted from shot point A2. The first arrival time of shot-receiver point 3 is S3, which indicates the earliest time that receiver point B3 receives the seismic wave emitted from shot point A3. The first arrival time of each shot receiver is recorded independently. This one-to-one correspondence allows for the analysis of the first arrival time of each shot receiver separately, in order to determine the distance between different shot receivers and receivers, as well as the propagation of seismic waves.

[0070] In some embodiments, the first arrival time S corresponding to each receiver point can be picked from the z-component seismic record. ijIn this method, i is the shot point sequence number, i = 1, 2, ..., N, where N is the total number of shot points, and j is the receiver sequence number, j = 1, 2, ..., M, where M is the total number of receiver points. In seismic exploration, the z-component typically represents seismic waves in the vertical direction, which is one of the strongest directions for seismic waves, and the main energy of seismic waves usually propagates along this vertical direction. By picking the first arrival time from the z-component, the main signals and propagation characteristics of seismic waves can be captured more accurately.

[0071] Step 223: Determine the first time-shift radial component of each of the multiple shot receivers based on their respective first radial components and their respective arrival times.

[0072] In some embodiments, based on the initial arrival time S of the shot receiver... ij For each first radial component Perform a time shift to obtain the time-shifted radial component, which is then used as the first time-shifted radial component. Wherein, the time shift can be S ij The time shift direction is the direction in which time decreases. That is, it is necessary to determine the initial arrival time S. ij The first radial component Adjustments are made accordingly on the time axis. The time shift direction is the direction of decreasing time, indicating that the first radial component needs to be shifted to the left of the time axis to compensate for the difference in the first arrival time of the seismic waves. For example, for receiver point j, the first arrival time of the seismic wave emitted from shot point i is S. ij .but

[0073] The above method, by performing a time-shift operation on the first radial component, can more accurately represent the actual propagation of seismic waves, remove the influence of the first arrival time, and facilitate further analysis.

[0074] Step 230: Based on the first time-shift radial components of each of the multiple shot receivers, determine the first time-shift radial component model traces of each of the multiple shot receivers. The first time-shift radial component model traces of the receiving points are determined based on the multiple first time-shift radial components included in the receiving points.

[0075] A radial component model trace is a comprehensive model representing the radial component of seismic waves at a receiver point, obtained in seismic exploration by summing or other combinations of the first time-shifted radial components from multiple shot-receiver points. The radial component model trace is used to represent the radial response characteristics of seismic waves emitted from multiple shot points at different receiver locations.

[0076] In some embodiments, for each of a plurality of shot receiver points, a plurality of first time-shift radial components included in the receiver point of the shot receiver point are determined; the sum of the plurality of first time-shift radial components included in the receiver point is determined as the first time-shift radial component model trace of the receiver point. The calculation formula for the first time-shift radial component model trace can be as shown in Formula 4 below:

[0077]

[0078] in, Let N represent the first time-shift radial component model trace of the j-th receiving point, where N is the total number of shot points.

[0079] Step 240: Based on the first time-shift radial component of each of the multiple shot receivers and the first time-shift radial component model trace of each of the multiple receivers, determine the first correlation coefficient of each of the multiple shot receivers. The first correlation coefficient refers to the correlation coefficient between the first time-shift radial component and the first time-shift radial component model trace.

[0080] The correlation coefficient between the first time-shifted radial component and its model trace is used to assess the degree of matching between the seismic wave responses at different shot-receiver points and the receiver point. By calculating the correlation coefficient, the similarity between the first time-shifted radial component obtained from different shot-receiver points and the integrated model trace at the receiver point can be determined. This helps to understand whether the seismic wave responses at the shot-receiver points are consistent at the receiver point, thereby assessing the quality of the first time-shifted radial component data and the accuracy of the model. It is understandable that a higher correlation (i.e., a correlation coefficient close to 1) indicates a relatively consistent response pattern between the radial component calculated based on the right-hand rule and the model trace at the receiver point. This suggests that the seismic wave signal from the shot-receiver point can better interpret the model trace at the receiver point, meaning that the first radial component calculated based on the right-hand rotation formula conforms to the actual radial component determination rules. Conversely, a lower correlation (i.e., a correlation coefficient close to 0) indicates an inconsistent response pattern between the radial component calculated based on the right-hand rule and the model trace at the receiver point, meaning that the first radial component calculated based on the right-hand rotation formula does not conform to the actual radial component determination rules.

[0081] In some embodiments, for each of a plurality of shot receivers, a first time-shift radial component of the shot receiver within the first arrival window is obtained, wherein the start and end times of the first arrival window are determined based on the first arrival time of the shot receiver, which is the earliest time point at which the receiver at the shot receiver receives the seismic wave emitted by the shot; a model trace of the first time-shift radial component of the shot receiver within the first arrival window is obtained; and a first correlation coefficient of each of the plurality of shot receivers is determined based on the first time-shift radial component of each of the plurality of shot receivers within the first arrival window and the model trace of the first time-shift radial component of each of the plurality of shot receivers within the first arrival window.

[0082] The first arrival window is a time interval based on the first arrival time, used to capture the response of seismic waves after they arrive at the receiving point. The start and end times of this time window can be determined based on the first arrival time of the shot-receiver point. In some embodiments, the start time of the first arrival window for the shot-receiver point can be less than or equal to the first arrival time of that point, and the end time of the first arrival window can be greater than or equal to the first arrival time of that point. For example, if the first arrival time of the shot-receiver point is S1, then the first arrival window can be [S2, S3], where S1 can be greater than or equal to S2, and S3 can be less than or equal to S3.

[0083] In some embodiments, the length of the first arrival window is not limited in this application. First correlation coefficient The calculation formula can be shown in Formula 5 below:

[0084]

[0085] Where k is the sequence number of the seismic data sample point within the first arrival window, and L is the number of sample points within the length of the first arrival window. This application does not limit the value of L. This is the expected value of the first time-shift radial component. This is the expected value of the first time-shift radial component model channel.

[0086] In some embodiments, the expected value of the first time-shifted radial component can be the average of a plurality of first time-shifted radial components, and the expected value of the first time-shifted radial component model channel can be the average of a plurality of first time-shifted radial component model channels.

[0087] In some embodiments, the formula for calculating the correlation coefficient described above may also be based on other methods, such as the Spearman rank correlation coefficient, which is not limited in this application.

[0088] Step 250: Determine the second time-shifted radial component of each of the multiple shot receivers based on their respective horizontal components. The second time-shifted radial component is determined by time-shifting the second radial component in the time domain. The second radial component is determined based on the horizontal component under the left-hand rule.

[0089] Similarly, step 250 above can be implemented as steps 251 to 253.

[0090] Step 251: Determine the second radial component of each of the multiple shot receivers based on their respective horizontal components using the left-hand rotation formula.

[0091] In some embodiments, the step of determining the second radial component of each of the plurality of shot-receiver points may include: obtaining the shot point coordinates and receiver coordinates of each of the plurality of shot-receiver points; determining the shot-receiver heading angle of each of the plurality of shot-receiver points based on the shot point coordinates and receiver coordinates of each of the plurality of shot-receiver points; and determining the second radial component of each of the plurality of shot-receiver points based on the horizontal component of each of the plurality of shot-receiver points and the shot-receiver heading angle of each of the plurality of shot-receiver points using a right-hand rotation formula. Figure 4 As shown, the sub Figure 1 The horizontal component (geodetic coordinate system) and sub-component of the three-component seismic record acquired by a right-handed three-component geophone are... Figure 2 This represents the horizontal component (geodetic coordinate system) of a three-component seismic record acquired by a left-handed three-component geophone. The horizontal axis represents the x and y coordinates, and the vertical axis represents time. For example... Figure 5 As shown, the sub Figure 1 The horizontal component, R component, and T component, and sub-component of the right-handed three-component detector are rotated after recording. Figure 2 The left-handed three-component detector records the horizontal component, followed by the rotated R and T components. The horizontal axis represents the R and T coordinates, and the vertical axis represents time. Please refer to step 221 for the above related content.

[0092] In some embodiments, the left-handed rotation formula is used based on the horizontal components x of each of the multiple shot receiver points. ij y ij and the firing direction angle α ij The second radial component is obtained after rotating the horizontal component. Second tangential component The left-hand rotation formula is shown in Formula 6 below:

[0093]

[0094] Step 252: Obtain the first arrival times for each of the multiple shot receivers. The first arrival time refers to the earliest time when the receiver at a shot receiver receives the seismic wave emitted by the shot. For related content, please refer to step 222 above.

[0095] Step 253: Determine the second time-shift radial component of each of the multiple shot receivers based on their respective second radial components and their respective first arrival times.

[0096] In some embodiments, based on the initial arrival time S of the shot receiver... ij For each second radial component Perform a time shift to obtain the time-shifted radial component, which is then used as the second time-shifted radial component. Wherein, the time shift can be S ij The time shift direction is the direction in which time decreases. That is, it is necessary to determine the initial arrival time S. ij The second radial component Adjustments are made accordingly on the time axis. The time shift direction is the direction of decreasing time, indicating that the second radial component needs to be shifted to the left of the time axis to compensate for the difference in the first arrival time of the seismic waves. For example, for receiver point j, the first arrival time of the seismic wave emitted from shot point i is S. ij .but

[0097] Step 260: Based on the second time-shift radial components of each of the multiple shot receivers, determine the second time-shift radial component model traces of each of the multiple shot receivers. The second time-shift radial component model traces of the receiving points are determined based on the multiple second time-shift radial components included in the receiving points.

[0098] In some embodiments, for each of a plurality of shot receiver points, a plurality of second time-shift radial components included in the receiver point of the shot receiver point are determined; the sum of the plurality of second time-shift radial components included in the receiver point is determined as the second time-shift radial component model trace of the receiver point. The calculation formula for the second time-shift radial component model trace can be as shown in Formula 7 below:

[0099]

[0100] in, Let N represent the second time-shift radial component model trace of the j-th receiving point, where N is the total number of shot points.

[0101] Step 270: Based on the second time-shift radial component of each of the multiple shot receivers and the second time-shift radial component model trace of each of the multiple receivers, determine the second correlation coefficient of each of the multiple shot receivers. The second correlation coefficient refers to the correlation coefficient between the second time-shift radial component and the second time-shift radial component model trace.

[0102] The correlation coefficient between the second time-shifted radial component and its model trace is used to assess the degree of matching between the seismic wave responses at different shot-receiver points and the receiver point. By calculating the correlation coefficient, the similarity between the second time-shifted radial component obtained from different shot-receiver points and the integrated model trace at the receiver point can be determined. This helps to understand whether the seismic wave responses at the shot-receiver points are consistent at the receiver point, thereby assessing the quality of the second time-shifted radial component data and the accuracy of the model. Understandably, a higher correlation (i.e., a correlation coefficient close to 1) indicates a relatively consistent response pattern between the radial component of the shot-receiver point calculated based on a left-handed system and the model trace at the receiver point. This suggests that the seismic wave signal from the shot-receiver point can better interpret the model trace at the receiver point, meaning that the second radial component calculated based on the left-handed rotation formula conforms to the actual radial component determination rules. Conversely, a lower correlation (i.e., a correlation coefficient close to 0) indicates an inconsistent response pattern between the radial component of the shot-receiver point calculated based on a left-handed system and the model trace at the receiver point, meaning that the second radial component calculated based on the left-handed rotation formula does not conform to the actual radial component determination rules.

[0103] In some embodiments, for each of a plurality of shot receivers, a second time-shifted radial component of the shot receiver within the first arrival window is obtained. The start and end times of the first arrival window are determined based on the first arrival time of the shot receiver, which is the earliest time point at which the receiver at the shot receiver receives the seismic wave emitted by the shot. A model trace of the second time-shifted radial component of the shot receiver within the first arrival window is obtained. Based on the second time-shifted radial components of each of the plurality of shot receivers within the first arrival window, and the model trace of the second time-shifted radial component of each of the plurality of shot receivers within the first arrival window, a second correlation coefficient for each of the plurality of shot receivers is determined.

[0104] Second correlation coefficient The calculation formula can be shown in Formula 8 below:

[0105]

[0106] Where k is the sequence number of the seismic data sample point within the first arrival window, and L is the number of sample points within the length of the first arrival window. This application does not limit the value of L. This is the expected value of the second time-shift radial component. This is the expected value of the second time-shifted radial component model channel.

[0107] In some embodiments, the expected value of the second time-shifted radial component can be the average of a plurality of second time-shifted radial components, and the expected value of the second time-shifted radial component model channel can be the average of a plurality of second time-shifted radial component model channels.

[0108] Step 280: Determine the measurement rules for the three-component detector based on the first correlation coefficients of each of the multiple shot detectors and the second correlation coefficients of each of the multiple shot detectors.

[0109] In some embodiments, the first correlation coefficient and the second correlation coefficient are determined based on right-handed and left-handed coordinate system rules, respectively. By comparing these two correlation coefficients, it can be determined which coordinate system criterion is more suitable for determining the radial component. Choosing a criterion with a larger correlation coefficient makes the measurement rules more consistent with the actual three-component detector data acquisition rules, which is beneficial to the accuracy of subsequent multi-component data processing.

[0110] In some embodiments, the measurement rules for the three-component detector can be determined based on the relationship between the first correlation coefficient and the second correlation coefficient. In some embodiments, when the first correlation coefficient is greater than the second correlation coefficient, since the first correlation coefficient is calculated based on the right-hand rule, the measurement rules for the three-component detector are right-hand rules. When the second correlation coefficient is greater than the second correlation coefficient, since the second correlation coefficient is calculated based on the left-hand rule, the measurement rules for the three-component detector are left-hand rules.

[0111] Table 1 below shows the correlation coefficients of different nodes calculated based on right-handed and left-handed rules. Please refer to Table 1; for example, in row 1, the correlation coefficient C... 1 (That is, the first correlation coefficient) is greater than the correlation coefficient C. 2 (That is, the second correlation coefficient). Since the first correlation coefficient is calculated based on the right-hand rule, the estimated three-component detector type (that is, the measurement rule of the three-component detector) follows the right-hand rule and is consistent with the actual measurement rule of the three-component detector. The method proposed in this application ensures that the actual three-component detector type matches the estimated three-component detector type, thus providing an accurate estimation result of the measurement rule of the three-component detector.

[0112] Table 1

[0113] Serial Number node station number Detector type <![CDATA[Coefficient of correlation C 1 > <![CDATA[Coefficient of correlation C 2 > Estimated detector type 1 2005 Right hand 0.95 0.03 Right hand 2 2006 Left-handed 0.04 0.97 Left-handed 3 2010 Right hand 0.96 0.04 Right hand 4 2011 Left-handed 0.06 0.96 Left-handed 5 2028 Right hand 0.75 0.21 Right hand 6 2029 Left-handed 0.13 0.95 Left-handed

[0114] The technical solution provided in this application, based on the different rotation methods of the horizontal component to the radial component collected by the three-component detector under the left-hand rule and the right-hand rule, can accurately determine the left-hand and right-hand criteria followed by the three-component detector by comparing the time-shifted radial component and the time-shifted radial component model trace under the left-hand rule and the right-hand rule.

[0115] The following describes the specific implementation method for determining the measurement fraction of the three-component detector. Step 280 above can be implemented as the following steps 281 to 283.

[0116] Step 281: Determine the sum of the first correlation coefficients of multiple shot receivers based on their individual first correlation coefficients. As shown in Formula 9 below:

[0117]

[0118] in, This is the sum of the first correlation coefficients. The sum of the first correlation coefficients refers to summing the first correlation coefficients of all shot-receiver points within the first arrival time window to obtain an overall correlation assessment value. This sum of the first correlation coefficients represents the comprehensive situation of the correlation between the first time-shift radial component and the first time-shift radial component model traces among all shot-receiver points under the right-hand rule.

[0119] Step 282: Determine the sum of the second correlation coefficients of the multiple shot receivers based on their respective second correlation coefficients. As shown in Formula 10 below:

[0120]

[0121] in, This is the sum of the second correlation coefficients. The sum of the second correlation coefficients refers to summing the second correlation coefficients of all shot-receiver points within the first arrival time window to obtain an overall correlation assessment value. This sum of the second correlation coefficients represents the comprehensive situation of the correlation between the second time-shift radial component and the second time-shift radial component model traces among all shot-receiver points under the left-handed criterion.

[0122] Step 283: Determine the measurement rules for the three-component detector based on the sum of the first correlation coefficient and the sum of the second correlation coefficient.

[0123] In some embodiments, if the sum of the first correlation coefficients is greater than the sum of the second correlation coefficients, the measurement rule of the three-component detector is determined to be the right-hand rule; if the sum of the first correlation coefficients is less than the sum of the second correlation coefficients, the measurement rule of the three-component detector is determined to be the left-hand rule.

[0124] In some embodiments, if the sum of the first correlation coefficients is greater than the sum of the second correlation coefficients, it indicates that the radial component obtained by rotation under the right-hand rule matches the actual observed data better than the radial component obtained under the left-hand rule. Therefore, the measurement rule for the three-component detector is determined to be the right-hand rule. If the sum of the first correlation coefficients is less than the sum of the second correlation coefficients, it indicates that the radial component obtained by rotation under the left-hand rule matches the actual observed data better than the radial component obtained under the right-hand rule. Therefore, the measurement rule for the three-component detector is determined to be the left-hand rule.

[0125] In some embodiments, the average value of the first correlation coefficients of the multiple shot receivers is determined based on the first correlation coefficients of each of the multiple shot receivers. Based on the second correlation coefficients of each of the multiple shot receivers, determine the average value of the second correlation coefficients of the multiple shot receivers. Based on the average value of the first correlation coefficient and the average value of the second correlation coefficient Determine the measurement rules for the three-component detector. In some embodiments, if the average value of the first correlation coefficient is... The average value of the second correlation coefficient Therefore, the measurement rule for the three-component detector is determined to be the right-hand rule. If the average value of the first correlation coefficient is less than the average value of the second correlation coefficient, then the measurement rule for the three-component detector is determined to be the left-hand rule.

[0126] The method described above, by comprehensively calculating the correlation coefficients of multiple shot-receiver points, can more fully evaluate which coordinate system rule (right-handed or left-handed) more accurately describes the seismic data overall. This aggregation method can reduce random errors at individual shot-receiver points and improve the reliability of overall data processing.

[0127] The following are embodiments of the apparatus described in this application, which can be used to execute the embodiments of the method described in this application. For details not disclosed in the apparatus embodiments of this application, please refer to the embodiments of the method described in this application.

[0128] Please refer to Figure 6 This diagram illustrates a block diagram of a digital signature generation apparatus according to an embodiment of this application. The apparatus has the function of implementing the aforementioned digital signature generation method; this function can be implemented in hardware or by hardware executing corresponding software. The apparatus can be a computer device or can be installed within a computer device. The apparatus 600 may include: an acquisition module 610, a first determination module 620, a second determination module 630, and a third determination module 640.

[0129] The acquisition module 610 is used to acquire first multi-component data. The first multi-component data is determined based on a three-component geophone according to measurement rules. The measurement rules are used to instruct the coordinate system criteria for the three-component geophone to record the multi-component data of seismic waves. The measurement rules include left-hand rule and right-hand rule. The first multi-component data includes the horizontal components of multiple shot-receiver points. The shot-receiver points are determined based on the shot point and the receiver point. The shot point is used to emit seismic waves, and the receiver point is used to receive the seismic waves emitted by the shot point.

[0130] The first determining module 620 is used to determine the first time-shift radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers. The first time-shift radial component is determined by time-shifting the first radial component in the time domain. The first radial component is determined based on the horizontal component under the right-hand rule.

[0131] The first determining module 620 is further configured to determine the first time-shift radial component model trace of each of the plurality of gun receivers based on the first time-shift radial component of each of the plurality of gun receivers, wherein the first time-shift radial component model trace of the receiving point is determined based on the plurality of first time-shift radial components included in the receiving point.

[0132] The first determining module 620 is further configured to determine the first correlation coefficient of each of the plurality of shot receivers based on the first time-shift radial component of each of the plurality of shot receivers and the first time-shift radial component model trace of each of the plurality of receivers, wherein the first correlation coefficient refers to the correlation coefficient between the first time-shift radial component and the first time-shift radial component model trace.

[0133] The second determining module 630 is further configured to determine the second time-shift radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers. The second time-shift radial component is determined by time-shifting the second radial component in the time domain. The second radial component is determined based on the horizontal component under the left-hand rule.

[0134] The second determining module 630 is further configured to determine the second time-shift radial component model trace of each of the plurality of gun receivers based on the second time-shift radial component of each of the plurality of gun receivers, wherein the second time-shift radial component model trace of the receiving point is determined based on the plurality of second time-shift radial components included in the receiving point.

[0135] The second determining module 630 is further configured to determine the second correlation coefficient of each of the plurality of shot receivers based on the second time-shift radial component of each of the plurality of shot receivers and the second time-shift radial component model trace of each of the plurality of receivers, wherein the second correlation coefficient refers to the correlation coefficient between the second time-shift radial component and the second time-shift radial component model trace.

[0136] The third determining module 640 is used to determine the measurement rules of the three-component detector based on the first correlation coefficient of each of the plurality of shot detectors and the second correlation coefficient of each of the plurality of shot detectors.

[0137] In some embodiments, the third determining module 640 is configured to determine the sum of the first correlation coefficients of the plurality of shot detectors based on the first correlation coefficients of each of the plurality of shot detectors; determine the sum of the second correlation coefficients of the plurality of shot detectors based on the second correlation coefficients of each of the plurality of shot detectors; and determine the measurement rules of the three-component detector based on the sum of the first correlation coefficients and the sum of the second correlation coefficients.

[0138] In some embodiments, the third determining module 640 is used to determine the measurement rule of the three-component detector as the right-hand rule if the sum of the first correlation coefficients is greater than the sum of the second correlation coefficients; and to determine the measurement rule of the three-component detector as the left-hand rule if the sum of the first correlation coefficients is less than the sum of the second correlation coefficients.

[0139] In some embodiments, the first determining module 620 is configured to, for each of the plurality of shot receivers, obtain the first time-shift radial component of the shot receiver within the first arrival window, wherein the start and end times of the first arrival window are determined based on the first arrival time of the shot receiver, which refers to the earliest time point at which the receiver at the shot receiver receives the seismic wave emitted by the shot; obtain the model trace of the first time-shift radial component of the shot receiver within the first arrival window; and, based on the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, and the first time-shift radial component of each of the plurality of shot receivers, obtain the model trace of the first time-shift radial component of the shot receiver within the first arrival window; and, based on the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, obtain the model trace of the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, and the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, obtain the model trace of the first time-shift radial component of the shot receiver within the first arrival window; and obtain the model trace of the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, and the first time-shift radial component of each of the plurality of shot receivers within the first arrival window, obtain the model trace of the first time-shift radial component of the first time-shift radial component of each of the plurality of shot receivers within the first arrival window; obtain the model trace of the first time-shift radial component of each of the plurality of shot receivers within the first arrival window; and ... The first correlation coefficient of each of the plurality of shot receivers is determined by the first time-shift radial component model trace within the first arrival window; the second determining module 630 is used to obtain the second time-shift radial component of the shot receiver within the first arrival window for each of the plurality of shot receivers; obtain the second time-shift radial component model trace of the shot receiver within the first arrival window; and determine the second correlation coefficient of each of the plurality of shot receivers based on the second time-shift radial component of each of the plurality of shot receivers within the first arrival window and the second time-shift radial component model trace of each of the plurality of shot receivers within the first arrival window.

[0140] In some embodiments, the first determining module 620 is configured to determine the first radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers using a right-handed rotation formula; obtain the first arrival time of each of the plurality of shot receivers, wherein the first arrival time refers to the earliest time point at which the receiving point in the shot receiver receives the seismic wave emitted by the shot point; and determine the first time-shifted radial component of each of the plurality of shot receivers based on the first radial component of each of the plurality of shot receivers and the first arrival time of each of the plurality of shot receivers. The second determining module 630 is configured to determine the second radial component of each of the plurality of shot receivers based on the horizontal component of each of the plurality of shot receivers using a left-handed rotation formula; and determine the second time-shifted radial component of each of the plurality of shot receivers based on the second radial component of each of the plurality of shot receivers and the first arrival time of each of the plurality of shot receivers.

[0141] In some embodiments, the first determining module 620 is used to obtain the shot point coordinates and receiver coordinates of each of the plurality of shot and receiver points; determine the shot and receiver direction angle of each of the plurality of shot and receiver points based on the shot point coordinates and receiver coordinates of each of the plurality of shot and receiver points; and determine the first radial component of each of the plurality of shot and receiver points based on the horizontal component and the shot and receiver direction angle of each of the plurality of shot and receiver points using the right-hand rotation formula; the second determining module 630 is used to determine the second radial component of each of the plurality of shot and receiver points based on the horizontal component and the shot and receiver direction angle of each of the plurality of shot and receiver points using the left-hand rotation formula.

[0142] It should be noted that the apparatus provided in the above embodiments is only illustrated by the division of the above functional modules when implementing its functions. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the apparatus and method embodiments provided in the above embodiments belong to the same concept, and the specific implementation process can be found in the method embodiments, which will not be repeated here.

[0143] Please refer to Figure 7 It shows a structural block diagram of a computer device 700 provided in one embodiment of this application.

[0144] Typically, computer device 700 includes a processor 710 and a memory 720.

[0145] Processor 710 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. Processor 710 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field Programmable Gate Array), and PLA (Programmable Logic Array). Processor 710 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, processor 710 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the screen. In some embodiments, processor 710 may also include an AI processor for handling computational operations related to machine learning.

[0146] The memory 720 may include one or more computer-readable storage media, which may be non-transitory. The memory 720 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory devices. In some embodiments, the non-transitory computer-readable storage media in the memory 720 is used to store a computer program configured to be executed by one or more processors to implement the measurement rule determination method for the three-component detector described above.

[0147] Those skilled in the art will understand that Figure 7 The structure shown does not constitute a limitation on the computer device 700, and may include more or fewer components than shown, or combine certain components, or use different component arrangements.

[0148] In some embodiments, a computer-readable storage medium is also provided, wherein a computer program is stored therein, the computer program being loaded and executed by a processor to implement the measurement rule determination method for the three-component detector described above.

[0149] Optionally, the computer-readable storage medium may include: ROM (Read-Only Memory), RAM (Random-Access Memory), SSD (Solid State Drives), or optical disc, etc. The random access memory may include ReRAM (Resistance Random Access Memory) and DRAM (Dynamic Random Access Memory).

[0150] In some embodiments, a computer program product is also provided, the computer program product including a computer program stored in a computer-readable storage medium, the processor reading from the computer-readable storage medium and executing the computer program to implement the above-described method for determining the measurement rules of the three-component detector.

[0151] It should be understood that "multiple" as used herein refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. Furthermore, the step numbers described herein are merely illustrative of one possible execution order. In some other embodiments, the steps may not be executed in numerical order, such as two steps with different numbers being executed simultaneously, or two steps with different numbers being executed in the reverse order of the illustration. This application does not limit this.

[0152] The above description is merely an optional embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for determining a measurement rule of a three-component geophone, characterized in that, The method comprises: acquiring first multi-component data, the first multi-component data being determined based on a three-component geophone according to a measurement rule, the measurement rule being used to indicate a coordinate system criterion of the three-component geophone recording multi-component data of seismic waves, the measurement rule comprising a left-hand system rule and a right-hand system rule, the first multi-component data comprising a plurality of respective horizontal components of shot-receiver points, the shot-receiver points being determined based on a shot point and a receiving point, the shot point being used to emit seismic waves, and the receiving point being used to receive the seismic waves emitted by the shot point; determining, according to the respective horizontal components of the plurality of shot-receiver points, a first time-shifted radial component of each of the plurality of shot-receiver points, the first time-shifted radial component being determined by time-shifting a first radial component in a time domain, the first radial component being determined based on the horizontal component under the right-hand system rule; determining, according to the respective first time-shifted radial components of the plurality of shot-receiver points, a first time-shifted radial component model trace of each of a plurality of receiving points in the plurality of shot-receiver points, the first time-shifted radial component model trace of the receiving point being determined based on a plurality of first time-shifted radial components included by the receiving point; determining, according to the respective first time-shifted radial components of the plurality of shot-receiver points and the respective first time-shifted radial component model traces of the plurality of receiving points, a first correlation coefficient of each of the plurality of shot-receiver points, the first correlation coefficient being a correlation coefficient between the first time-shifted radial component and the first time-shifted radial component model trace; determining, according to the respective horizontal components of the plurality of shot-receiver points, a second time-shifted radial component of each of the plurality of shot-receiver points, the second time-shifted radial component being determined by time-shifting a second radial component in a time domain, the second radial component being determined based on the horizontal component under the left-hand system rule; determining, according to the respective second time-shifted radial components of the plurality of shot-receiver points, a second time-shifted radial component model trace of each of a plurality of receiving points in the plurality of shot-receiver points, the second time-shifted radial component model trace of the receiving point being determined based on a plurality of second time-shifted radial components included by the receiving point; determining, according to the respective second time-shifted radial components of the plurality of shot-receiver points and the respective second time-shifted radial component model traces of the plurality of receiving points, a second correlation coefficient of each of the plurality of shot-receiver points, the second correlation coefficient being a correlation coefficient between the second time-shifted radial component and the second time-shifted radial component model trace; determining, according to the respective first correlation coefficients of the plurality of shot-receiver points and the respective second correlation coefficients of the plurality of shot-receiver points, a measurement rule of the three-component geophone.

2. The method of claim 1, wherein, The determining, according to the respective first correlation coefficients of the plurality of shot-receiver points and the respective second correlation coefficients of the plurality of shot-receiver points, a measurement rule of the three-component geophone, comprises: determining a sum of the first correlation coefficients of the plurality of shot-receiver points according to the respective first correlation coefficients of the plurality of shot-receiver points; determining a sum of the second correlation coefficients of the plurality of shot-receiver points according to the respective second correlation coefficients of the plurality of shot-receiver points; determining the measurement rule of the three-component geophone according to the sum of the first correlation coefficients and the sum of the second correlation coefficients.

3. The method of claim 2, wherein, The method comprises the following steps: If the sum of the first correlation coefficients is greater than the sum of the second correlation coefficients, the measurement rule of the three-component geophone is determined as the right-hand rule; If the sum of the first correlation coefficients is less than the sum of the second correlation coefficients, the measurement rule of the three-component geophone is determined as the left-hand rule.

4. The method of claim 1, wherein, The method comprises the following steps: For each of the plurality of shot-receiver points, the first time-shifted radial component of the shot-receiver point in a first arrival time window is obtained, the start time and the end time of the first arrival time window being determined based on the first arrival time of the shot-receiver point, the first arrival time being the earliest time point at which the receiver in the shot-receiver point receives the seismic wave emitted by the shot point; The first time-shifted radial component model trace of the shot-receiver point in the first arrival time window is obtained; The first correlation coefficient of each of the plurality of shot-receiver points is determined based on the first time-shifted radial component of each of the plurality of shot-receiver points in the first arrival time window and the first time-shifted radial component model trace of each of the plurality of shot-receiver points in the first arrival time window. The method comprises the following steps: For each of the plurality of shot-receiver points, the second time-shifted radial component of the shot-receiver point in the first arrival time window is obtained; The second time-shifted radial component model trace of the shot-receiver point in the first arrival time window is obtained; The second correlation coefficient of each of the plurality of shot-receiver points is determined based on the second time-shifted radial component of each of the plurality of shot-receiver points in the first arrival time window and the second time-shifted radial component model trace of each of the plurality of shot-receiver points in the first arrival time window.

5. The method of claim 1, wherein, The method comprises the following steps: The first time-shifted radial component of each of the plurality of shot-receiver points is determined based on the horizontal component of each of the plurality of shot-receiver points by using the right-hand rule rotation formula; The first arrival time of each of the plurality of shot-receiver points is obtained, the first arrival time being the earliest time point at which the receiver in the shot-receiver point receives the seismic wave emitted by the shot point; The first time-shifted radial component of each of the plurality of shot-receiver points is determined based on the first radial component of each of the plurality of shot-receiver points and the first arrival time of each of the plurality of shot-receiver points. The method comprises the following steps: The second time-shifted radial component of each of the plurality of shot-receiver points is determined based on the horizontal component of each of the plurality of shot-receiver points by using the left-hand rule rotation formula; The second time-shifted radial component of each of the plurality of shot-receiver points is determined based on the second radial component of each of the plurality of shot-receiver points and the first arrival time of each of the plurality of shot-receiver points.

6. The method of claim 5, wherein, The determining the first radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points by the right-hand system rotation formula comprises: obtaining the shot point coordinate and the receiving point coordinate of each of the plurality of shot-receiver points; determining the shot-receiver direction angle of each of the plurality of shot-receiver points according to the shot point coordinate and the receiving point coordinate of each of the plurality of shot-receiver points; determining the first radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points and the shot-receiver direction angle of each of the plurality of shot-receiver points by the right-hand system rotation formula; The determining the second radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points by the left-hand system rotation formula comprises: determining the second radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points and the shot-receiver direction angle of each of the plurality of shot-receiver points by the left-hand system rotation formula.

7. A measurement rule determination apparatus of a three-component geophone, characterized by comprising: a three-component geophone measurement rule determination unit that determines a measurement rule of the three-component geophone based on a type of the three-component geophone. The device comprises: an obtaining module, configured to obtain first multi-component data, the first multi-component data being determined based on a three-component geophone according to a measurement rule, the measurement rule being used to indicate a coordinate system criterion of the three-component geophone recording seismic wave multi-component data, the measurement rule comprising a left-hand system rule and a right-hand system rule, the first multi-component data comprising a horizontal component of each of a plurality of shot-receiver points, the shot-receiver point being determined based on a shot point and a receiving point, the shot point being used to emit a seismic wave, and the receiving point being used to receive the seismic wave emitted by the shot point; a first determining module, configured to determine a first time-shift radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points, the first time-shift radial component being determined by time-shifting a first radial component in a time domain, the first radial component being determined based on the horizontal component under the right-hand system rule; the first determining module is further configured to determine a first time-shift radial component model trace of each of a plurality of receiving points in the plurality of shot-receiver points according to the first time-shift radial component of each of the plurality of shot-receiver points, the first time-shift radial component model trace of the receiving point being determined based on a plurality of first time-shift radial components included by the receiving point; the first determining module is further configured to determine a first correlation coefficient of each of the plurality of shot-receiver points according to the first time-shift radial component of each of the plurality of shot-receiver points and the first time-shift radial component model trace of each of the plurality of receiving points, the first correlation coefficient being a correlation coefficient between the first time-shift radial component and the first time-shift radial component model trace; a second determining module, configured to determine a second time-shift radial component of each of the plurality of shot-receiver points according to the horizontal component of each of the plurality of shot-receiver points, the second time-shift radial component being determined by time-shifting a second radial component in a time domain, the second radial component being determined based on the horizontal component under the left-hand system rule; The second determining module is further configured to determine, according to the second time-shift radial components of the plurality of shot-receiver points, second time-shift radial component model traces of a plurality of receiver points in the plurality of shot-receiver points, the second time-shift radial component model trace of the receiver point being determined based on the plurality of second time-shift radial components included in the receiver point; The second determining module is further configured to determine, according to the second time-shift radial components of the plurality of shot-receiver points and the second time-shift radial component model traces of the plurality of receiver points, second correlation coefficients of the plurality of shot-receiver points, the second correlation coefficient being a correlation coefficient between the second time-shift radial component and the second time-shift radial component model trace; The third determining module is configured to determine, according to the first correlation coefficients of the plurality of shot-receiver points and the second correlation coefficients of the plurality of shot-receiver points, a measurement rule of the three-component geophone.

8. A computer device, comprising: The computer device comprises a processor and a memory, and the memory stores a computer program, which is loaded and executed by the processor to implement the method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, which is loaded and executed by the processor to implement the method according to any one of claims 1 to 6.

10. A computer program product, characterised in that, The computer program product comprises a computer program, which is loaded and executed by the processor to implement the method according to any one of claims 1 to 6.

Images