Seismic data correction methods, devices and computer equipment

Abstract

This application discloses a seismic data correction method, apparatus, and computer equipment, belonging to the field of geophysical exploration technology. The method uses a first measurement angle parameter, multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters to determine a first measurement intensity parameter on a first horizontal component, a second measurement intensity parameter on a second horizontal component, and a third measurement intensity parameter on a first vertical component. Based on the first, second, and third measurement intensity parameters corresponding to the multiple measurement intensity parameters, a first correction matrix is ​​determined. This first correction matrix can then be used to correct the actual first seismic data to obtain theoretical seismic data, thereby improving the accuracy of the seismic data.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration technology. In particular, it relates to a seismic data correction method, apparatus, and computer equipment. Background Technology

[0002] Currently, in seismic exploration, seismic waves are primarily generated at the excitation point, and seismic data are acquired at the receiving point using geophones. For a single excitation point, it can be generated by sources in multiple directions. These sources can be sources along the survey line (equivalent to X-sources), sources perpendicular to the survey line and in the same plane as the survey line (equivalent to Y-sources), and sources perpendicular to the plane of the survey line (equivalent to Z-sources). For a single receiving point, a three-component geophone can be used to record the three-dimensional particle vibrations caused by different sources. The three components of this geophone are the X, Y, and Z components. Therefore, if an excitation point is generated three times sequentially by X-sources, Y-sources, and Z-sources, and each time the data is received by a three-component geophone, nine components of seismic data can be obtained, resulting in nine-component seismic data.

[0003] However, in actual seismic exploration, due to objective reasons such as topographic relief, multiple seismic sources may fail to be activated according to the preset activation direction, resulting in a certain deviation between the actual activation direction and the preset activation direction. In this case, the actual nine-component seismic data deviates from the theoretical nine-component seismic data, leading to poor accuracy of the seismic data. Therefore, it is necessary to correct the actual nine-component seismic data. Summary of the Invention

[0004] This application provides a seismic data correction method, apparatus, and computer equipment, which can improve the accuracy of seismic data. The specific technical solution is as follows:

[0005] On the one hand, embodiments of this application provide a seismic data correction method, the method comprising:

[0006] Acquire the first earthquake data obtained from multiple seismic sources;

[0007] Obtain the first measurement angle parameter corresponding to the first seismic data, whereby the first measurement angle parameter is used to represent the angle between the survey line direction corresponding to the first seismic data and the first horizontal component;

[0008] Multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters corresponding to the multiple earthquake sources are obtained. The second measurement angle parameter is used to represent the angle between the excitation direction of the earthquake source and the first vertical component. The third measurement angle parameter is used to represent the angle between the excitation direction of the earthquake source and the first horizontal component. The measurement intensity parameter is used to represent the magnitude of the force excited by the earthquake source. One earthquake source corresponds to one second measurement angle parameter, one third measurement angle parameter, and one measurement intensity parameter.

[0009] Based on the first measurement angle parameter, the plurality of second measurement angle parameters, and the plurality of third measurement angle parameters, a first measurement intensity parameter on the first horizontal component and a second measurement intensity parameter on the second horizontal component are determined for the plurality of measurement intensity parameters, respectively. The first measurement intensity parameter and the second measurement intensity parameter are used to represent the measurement intensity of the measurement intensity parameter on the first horizontal component and the measurement intensity on the second horizontal component, respectively.

[0010] Based on the plurality of second measurement angle parameters, a third measurement intensity parameter is determined on the first vertical component of the plurality of measurement intensity parameters respectively. The first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other. The third measurement intensity parameter is used to represent the measurement intensity of the measurement intensity parameter on the first vertical component. One measurement intensity parameter corresponds to the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter respectively.

[0011] Based on the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters, a first correction matrix is ​​determined;

[0012] Based on the first correction matrix, the first seismic data is corrected to obtain the second seismic data.

[0013] In one possible implementation, the process of determining the first measurement intensity parameter of the measurement intensity parameter corresponding to the earthquake source on the first horizontal component, based on the first measurement angle parameter, the second measurement angle parameter corresponding to a seismic source, and the third measurement angle parameter, includes:

[0014] Determine the sine value of the second measured angle parameter to obtain the first parameter value;

[0015] The difference between the first and third measured angle parameters is determined to obtain the fourth angle parameter;

[0016] Determine the sine value of the fourth angle parameter to obtain the value of the second parameter;

[0017] The first measurement intensity parameter is determined based on the first parameter value, the second parameter value, and the measurement intensity parameter.

[0018] In another possible implementation, determining the first measurement intensity parameter based on the first parameter value, the second parameter value, and the measurement intensity parameter includes:

[0019] The first measurement intensity parameter is obtained by multiplying the measured intensity parameter, the first parameter value, and the second parameter value.

[0020] In another possible implementation, the process of determining the second measurement intensity parameter of the measurement intensity parameter corresponding to the earthquake source on the second horizontal component, based on the first measurement angle parameter, the second measurement angle parameter corresponding to the earthquake source, and the third measurement angle parameter, includes:

[0021] Determine the cosine value of the fourth angle parameter to obtain the value of the third parameter;

[0022] The second measurement intensity parameter is determined based on the third parameter value, the first parameter value, and the measurement intensity parameter.

[0023] In another possible implementation, the process of determining the third measurement intensity parameter of the measurement intensity parameter corresponding to the earthquake source on the first vertical component, based on the second measurement angle parameter corresponding to the earthquake source, includes:

[0024] Determine the cosine value of the second measured angle parameter to obtain the fourth parameter value;

[0025] The third measurement intensity parameter is determined based on the fourth parameter value and the measurement intensity parameter.

[0026] In another possible implementation, determining the first correction matrix based on the first, second, and third measurement intensity parameters corresponding to the plurality of measurement intensity parameters respectively includes:

[0027] The first, second, and third measurement intensity parameters corresponding to a given measurement intensity parameter are determined to be elements of a row in the first correction matrix, resulting in multiple row elements;

[0028] The first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters are determined to be elements of a column in the first correction matrix, thereby obtaining a plurality of column elements;

[0029] The plurality of row elements and the plurality of column elements are combined to form the first correction matrix.

[0030] In another possible implementation, the step of correcting the first seismic data based on the first correction matrix to obtain the second seismic data includes:

[0031] Determine the inverse of the first correction matrix to obtain the second correction matrix;

[0032] The second seismic data is obtained by multiplying the second correction matrix with the first seismic data.

[0033] In another possible implementation, the method further includes:

[0034] Obtain a first preset angle parameter, a second preset angle parameter, and a preset intensity parameter. The first preset angle parameter and the second preset angle parameter are used to determine a preset excitation direction. The first preset angle parameter is used to represent the angle between the preset excitation direction and the first vertical component, and the second preset angle parameter is used to represent the angle between the preset excitation direction and the first horizontal component.

[0035] Based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter, and the preset intensity parameter, a third seismic data is determined. The third seismic data is the seismic data corresponding to the second seismic data under the survey line direction, the preset excitation direction, and the preset intensity parameter.

[0036] On the other hand, embodiments of this application provide a seismic data correction device, the device comprising:

[0037] The first acquisition module is used to acquire the first earthquake data obtained from multiple seismic sources.

[0038] The second acquisition module is used to acquire the first measurement angle parameter corresponding to the first seismic data, wherein the first measurement angle parameter is used to represent the angle between the survey line direction corresponding to the first seismic data and the first horizontal component.

[0039] The third acquisition module is used to acquire multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters corresponding to the multiple earthquake sources. The second measurement angle parameter is used to represent the angle between the excitation direction of the earthquake source and the first vertical component. The third measurement angle parameter is used to represent the angle between the excitation direction of the earthquake source and the first horizontal component. The measurement intensity parameter is used to represent the magnitude of the force excited by the earthquake source. One earthquake source corresponds to one second measurement angle parameter, one third measurement angle parameter, and one measurement intensity parameter.

[0040] The first determining module is used to determine, based on the first measuring angle parameter, the plurality of second measuring angle parameters, and the plurality of third measuring angle parameters, a first measuring intensity parameter on the first horizontal component and a second measuring intensity parameter on the second horizontal component, respectively, wherein the first measuring intensity parameter and the second measuring intensity parameter are used to represent the measuring intensity of the measuring intensity parameter on the first horizontal component and the measuring intensity of the measuring intensity parameter on the second horizontal component, respectively.

[0041] The second determining module is used to determine, based on the plurality of second measuring angle parameters, a third measuring intensity parameter on the first vertical component of the plurality of measuring intensity parameters, wherein the first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other, and the third measuring intensity parameter is used to represent the measuring intensity of the measuring intensity parameter on the first vertical component, and one measuring intensity parameter corresponds to the first measuring intensity parameter, the second measuring intensity parameter, and the third measuring intensity parameter respectively;

[0042] The third determining module is used to determine the first correction matrix based on the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters respectively;

[0043] The correction module is used to correct the first seismic data based on the first correction matrix to obtain the second seismic data.

[0044] In one possible implementation, the first determining module is configured to determine the sine value of the second measurement angle parameter to obtain a first parameter value; determine the difference between the first measurement angle parameter and the third measurement angle parameter to obtain a fourth angle parameter; determine the sine value of the fourth angle parameter to obtain a second parameter value; and determine the first measurement intensity parameter based on the first parameter value, the second parameter value, and the measurement intensity parameter.

[0045] In another possible implementation, the first determining module is used to determine the product of the measured intensity parameter, the first parameter value, and the second parameter value to obtain the first measured intensity parameter.

[0046] In another possible implementation, the first determining module is used to determine the cosine value of the fourth angle parameter to obtain the third parameter value; and based on the third parameter value, the first parameter value, and the measured intensity parameter, to determine the second measured intensity parameter.

[0047] In another possible implementation, the second determining module is used to determine the cosine value of the second measurement angle parameter to obtain a fourth parameter value; and based on the fourth parameter value and the measurement intensity parameter, to determine the third measurement intensity parameter.

[0048] In another possible implementation, the third determining module is used to determine that the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to a measurement intensity parameter are elements of a row in the first correction matrix, thereby obtaining multiple row elements; to determine that the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to the multiple measurement intensity parameters are elements of a column in the first correction matrix, thereby obtaining multiple column elements; and to form the first correction matrix by combining the multiple row elements and the multiple column elements.

[0049] In another possible implementation, the correction module is used to determine the inverse of the first correction matrix to obtain a second correction matrix; and to determine the product of the second correction matrix and the first seismic data to obtain the second seismic data.

[0050] In another possible implementation, the device further includes:

[0051] The fourth acquisition module is used to acquire a first preset angle parameter, a second preset angle parameter, and a preset intensity parameter. The first preset angle parameter and the second preset angle parameter are used to determine a preset excitation direction. The first preset angle parameter is used to represent the angle between the preset excitation direction and the first vertical component, and the second preset angle parameter is used to represent the angle between the preset excitation direction and the first horizontal component.

[0052] The fourth determining module is used to determine third seismic data based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter, and the preset intensity parameter. The third seismic data is the seismic data corresponding to the second seismic data under the survey line direction, the preset excitation direction, and the preset intensity parameter.

[0053] On the other hand, embodiments of this application provide a computer device, which includes a processor and a memory. The memory stores at least one piece of program code, which is loaded and executed by the processor to implement the operations performed in the seismic data correction method of this application.

[0054] On the other hand, embodiments of this application provide a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the seismic data correction method of this application.

[0055] On the other hand, embodiments of this application provide a computer program product that stores at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the seismic data correction method of this application.

[0056] The beneficial effects of the technical solutions provided in this application are:

[0057] This application provides a seismic data correction method. The method uses a first measurement angle parameter, multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters to determine a first measurement intensity parameter on a first horizontal component, a second measurement intensity parameter on a second horizontal component, and a third measurement intensity parameter on a first vertical component. Based on the first, second, and third measurement intensity parameters corresponding to the multiple measurement intensity parameters, a first correction matrix is ​​determined. This allows the actual first seismic data to be corrected using the first correction matrix to obtain theoretical seismic data, thereby improving the accuracy of the seismic data. Attached Figure Description

[0058] Figure 1 This is a flowchart of a seismic data correction method provided in an embodiment of this application;

[0059] Figure 2 This is a schematic diagram of the excitation direction of a seismic source provided in an embodiment of this application;

[0060] Figure 3 This is a schematic diagram of the structure of a seismic data correction device provided in an embodiment of this application;

[0061] Figure 4 This is a structural block diagram of a computer device provided in an embodiment of this application. Detailed Implementation

[0062] To make the technical solution and advantages of this application clearer, the embodiments of this application will be described in further detail below.

[0063] The terms "first," "second," "third," and "fourth," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish different objects, not to describe a specific order. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or apparatuses.

[0064] It should be noted that all information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, stored data, displayed data, etc.), and signals involved in this application have been authorized by the user or fully authorized by all parties, and the collection, use, and processing of related data must comply with the relevant laws, regulations, and standards of the relevant countries and regions. For example, the measurement angle parameters, measurement intensity parameters, and seismic data involved in this application were all obtained with full authorization.

[0065] This application provides a seismic data correction method, executed by a computer device, see [link to relevant documentation]. Figure 1 The method includes:

[0066] Step 101: The computer equipment acquires the first earthquake data obtained from multiple seismic sources.

[0067] Multiple sources include X-sources and Y-sources, and may also include Z-sources.

[0068] In one possible implementation, if multiple seismic sources, including an X-source and a Y-source, record the particle vibration at the receiving point using a three-component geophone, then the first seismic data is six-component seismic data.

[0069] In another possible implementation, if multiple seismic sources include an X-source, a Y-source, and a Z-source, and the vibration of the mass points is recorded at the receiving point using a three-component geophone, then the first seismic data is nine-component seismic data.

[0070] In another possible implementation, the computer device first acquires nine-component seismic data, then acquires three-component seismic data corresponding to the X source and three-component seismic data corresponding to the Y source from the nine-component seismic data, and combines the three-component seismic data corresponding to the X source and the three-component seismic data corresponding to the Y source to form the first seismic data.

[0071] Therefore, the first earthquake data can be either six-component or nine-component earthquake data, without any specific limitation.

[0072] Step 102: The computer equipment acquires the first measurement angle corresponding to the first seismic data.

[0073] The first measurement angle parameter is used to represent the angle between the survey line direction corresponding to the first seismic data and the first horizontal component.

[0074] For example, if the first horizontal component is the Y component, then the first measurement angle parameter is the angle between the survey line direction corresponding to the first seismic data and the Y component, which can be expressed as θ.

[0075] Step 103: The computer equipment acquires multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters corresponding to multiple seismic sources.

[0076] The second measurement angle parameter is used to represent the angle between the excitation direction of the source and the first vertical component. The third measurement angle parameter is used to represent the angle between the excitation direction of the source and the first horizontal component. The measurement intensity parameter is used to represent the magnitude of the force excited by the source. One source corresponds to one second measurement angle parameter, one third measurement angle parameter and one measurement intensity parameter.

[0077] Wherein, the first vertical component is the Z component, and the first horizontal component can be either the X component or the Y component. In this embodiment, only the Y component as the first horizontal component is used as an example for explanation.

[0078] In this embodiment, the coordinate axis directions can be predefined in a Cartesian three-dimensional coordinate system. For example, the first horizontal component is the Y component, the second horizontal component is the X component, the positive direction of the X component is east, the positive direction of the Y component is north, and the positive direction of the Z component is vertically upward. The excitation direction of the seismic source can be defined as a vector in the Cartesian three-dimensional spatial coordinate system, referred to here as the excitation vector. This excitation vector can be described by a second measurement angle parameter and a third measurement angle parameter.

[0079] First: The angle between the excitation vector and the Z component is represented by the second measured angle parameter, with the positive direction of the Z component being 0° and clockwise rotation being positive;

[0080] Second: The angle between the projection of the excitation vector onto the XY plane and the Y component is represented by the third measured angle parameter. The positive direction of the Y component is 0°, and when rotated clockwise towards the X component, the positive direction of the X component is 90°.

[0081] For example, see Figure 2 The second measurement angle parameter is denoted by β, and the third measurement angle parameter is denoted by α. Then the excitation direction of the earthquake source... It can be represented as (β, α). This represents the projection of the excitation vector onto the XOY plane. Indicates the direction of the survey line.

[0082] If the first earthquake data is six-component earthquake data, then the number of the second measurement angle parameter, the third measurement angle parameter, and the measurement intensity parameter are all 2. If the first earthquake data is nine-component earthquake data, then the number of the second measurement angle parameter, the third measurement angle parameter, and the measurement intensity parameter are all 3.

[0083] It should be noted that if the magnitude of the force excited by the seismic source is defined as a standard value, then the measured intensity parameter can be normalized to a unit of 1. The second and third measured angle parameters and the measured intensity parameter for each excitation can be obtained during construction. Furthermore, the computer equipment can execute steps 101-103 sequentially, or it can execute steps 102 and 103 first, then step 101, or it can execute step 103 first, then steps 101 and 102. There is no specific limitation on the order in which the computer equipment executes steps 101-103.

[0084] Step 104: The computer device determines, based on the first measurement angle parameter, multiple second measurement angle parameters, and multiple third measurement angle parameters, the first measurement intensity parameter on the first horizontal component and the second measurement intensity parameter on the second horizontal component, respectively.

[0085] The first measurement intensity parameter and the second measurement intensity parameter are used to represent the measurement intensity of the measurement intensity parameter in the first horizontal component and the measurement intensity in the second horizontal component, respectively.

[0086] If the first level component is a Y component, then the second level component is an X component. If the first level component is an X component, then the second level component is a Y component. In this embodiment, only the example of the first level component being a Y component and the second level component being an X component is used for illustration.

[0087] In this step, the process by which the computer equipment determines the first measured intensity parameter of the first horizontal component corresponding to a seismic source can be achieved through the following steps (A-1) to (A-4), including:

[0088] (A-1) The computer equipment determines the sine value of the second measurement angle parameter and obtains the value of the first parameter.

[0089] For example, if the second measured angle parameter is denoted as β, then the first parameter value is denoted as sinβ.

[0090] (A-2) The computer equipment determines the difference between the first and third measurement angle parameters to obtain the fourth angle parameter.

[0091] In this step, the computer device determines the value after subtracting the third measurement angle parameter from the first measurement angle parameter to obtain the fourth angle parameter.

[0092] For example, if the first measured angle parameter is denoted as θ and the third measured angle parameter is denoted as α, then the fourth angle parameter can be denoted as θ-α.

[0093] (A-3) The computer equipment determines the sine value of the fourth angle parameter and obtains the value of the second parameter.

[0094] If the fourth angle parameter is represented as θ-α, then the value of the second parameter can be represented as sin(θ-α).

[0095] (A-4) The computer device determines the first measurement intensity parameter based on the first parameter value, the second parameter value, and the measurement intensity parameter.

[0096] The computer equipment determines the product of the measured intensity parameter, the first parameter value, and the second parameter value to obtain the first measured intensity parameter. For example, if the measured intensity parameter is represented as g, then the first measured intensity parameter can be represented as g×sinβ×sin(θ-α).

[0097] In this step, the process by which the computer equipment determines the second measurement intensity parameter of the second horizontal component corresponding to a seismic source can be achieved through the following steps (B-1) to (B-2):

[0098] (B-1) The computer equipment determines the cosine value of the fourth angle parameter and obtains the value of the third parameter.

[0099] The value of this third parameter can be expressed as cos(θ-α).

[0100] (B-2) The computer equipment determines the second measurement intensity parameter based on the third parameter value, the first parameter value, and the measurement intensity parameter.

[0101] The computer equipment determines the product of the measured intensity parameter, the third parameter value, and the first parameter value to obtain the second measured intensity parameter, which can be expressed as g×sinβ×cos(θ-α).

[0102] In this step, the computer equipment can determine the first measurement intensity parameter first, and then the second measurement intensity parameter, or it can determine the second measurement intensity parameter first, and then the first measurement intensity parameter, with specific limitations applied. If the second measurement intensity parameter is determined first, the computer equipment first determines the fourth angle parameter based on the first and third measurement angle parameters, and then executes steps (B-1) to (B-2) to obtain the second measurement intensity parameter. In this case, when determining the first measurement intensity parameter, the sine value of the fourth angle parameter can be directly determined, without needing to determine the fourth angle parameter again.

[0103] Step 105: The computer device determines the third measurement intensity parameter on the first vertical component based on the multiple second measurement angle parameters.

[0104] The third measurement intensity parameter is used to represent the measurement intensity of the measurement intensity parameter on the first vertical component. One measurement intensity parameter corresponds to the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter.

[0105] The computer equipment determines the cosine value of the second measurement angle parameter corresponding to a seismic source, obtains the fourth parameter value, and determines the third measurement intensity parameter based on the fourth parameter value and the measurement intensity parameter corresponding to the seismic source.

[0106] In this implementation, the computer device determines the negative of the product of the fourth parameter value and the measured intensity parameter to obtain the third measured intensity parameter. For example, if the fourth parameter value is represented as cosβ, then the third intensity parameter can be represented as -g×cosβ. The computer device obtains multiple third measured intensity parameters through this method.

[0107] Step 106: The computer device determines the first correction matrix based on the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to the multiple measurement intensity parameters.

[0108] This step can be achieved through the following steps (1) to (3), including:

[0109] (1) The computer device determines that the first, second and third measurement intensity parameters corresponding to a measurement intensity parameter are elements of a row in the first correction matrix, thus obtaining multiple row elements.

[0110] If the earthquake source includes an X-source and a Y-source, then the first, second, and third measured intensity parameters corresponding to the measured intensity parameters of the X-source form the elements of one row in the first calibration matrix, and the first, second, and third measured intensity parameters corresponding to the measured intensity parameters of the Y-source form the elements of one row in the first calibration matrix, ultimately resulting in two rows of elements, with each row containing three elements.

[0111] If the earthquake source includes X source, Y source and Z source, then except for X source and Y source, the first measurement intensity parameter, second measurement intensity parameter and third measurement intensity parameter corresponding to the measurement intensity parameter of Z source form the elements of one row in the first correction matrix, and finally obtain three rows of row elements, with each row containing three elements.

[0112] The order of elements in a row can be set and changed as needed, without any specific limitation. For example, the elements in a row can be arranged from left to right in the order of X, Y, and Z components, or in the order of X, Z, and Y components, or in the order of Y, X, or Z components.

[0113] If the earthquake source includes both X-source and Y-source, then the first correction matrix consists of two rows of elements. The order of these two rows can be set and changed as needed, and there is no specific limitation on this. For example, the first row may contain elements corresponding to the X-source and the second row may contain elements corresponding to the Y-source, or the first row may contain elements corresponding to the Y-source and the second row may contain elements corresponding to the X-source.

[0114] If the earthquake sources include X-ray, Y-ray, and Z-ray sources, then the first correction consists of three rows of elements. The order of these three rows can be set and changed as needed, and there is no specific limitation on this. For example, the first row may correspond to the X-ray source, the second row to the Y-ray source, and the third row to the Z-ray source; or the first row may correspond to the Y-ray source, the second row to the X-ray source, and the third row to the Z-ray source.

[0115] (2) The computer equipment determines that the first measurement intensity parameter, the second measurement intensity parameter, and the third measurement intensity parameter corresponding to the multiple measurement intensity parameters are elements of a column in the first correction matrix, thus obtaining multiple column elements.

[0116] The computer equipment determines multiple first measurement intensity parameters as one column, multiple second measurement intensity parameters as one column, and multiple third measurement intensity parameters as one column, thus obtaining three columns of column elements.

[0117] The order of the elements in these three columns can be set and changed as needed, without any specific limitation. For example, if the earthquake source includes an X-ray source and a Y-ray source, the first column contains the second measured intensity parameter of the corresponding intensity parameters of the X-ray source and the Y-ray source in the X component; the second column contains the first measured intensity parameter of the corresponding intensity parameters of the X-ray source and the Y-ray source in the Y component; and the third column contains the third measured intensity parameter of the corresponding intensity parameters of the X-ray source and the Y-ray source in the Z component. Similarly, if the earthquake source includes an X-ray source, a Y-ray source, and a Z-ray source, then the first column contains the second measured intensity parameter of the corresponding intensity parameters of the X-ray source, the Y-ray source, and the Z-ray source in the X component; the second column contains the first measured intensity parameter of the corresponding intensity parameters of the X-ray source, the Y-ray source, and the Z-ray source in the Y component; and the third column contains the third measured intensity parameter of the corresponding intensity parameters of the X-ray source, the Y-ray source, and the Z-ray source in the Z component.

[0118] (3) The computer device forms a first correction matrix by combining multiple row elements and multiple column elements.

[0119] If the earthquake source includes both X-source and Y-source, then a row contains three elements and a column contains two elements. For the elements in a row, the computer device can arrange them from left to right in the order of X-component, Y-component, and Z-component. For the elements in a column, the computer device can arrange them from top to bottom in the order of X-source and Y-source, forming the first correction matrix.

[0120] Accordingly, the first correction matrix can be expressed as:

[0121]

[0122] Where P represents the first correction matrix, g1, β1, and α1 represent the measured intensity parameter, second measured angle parameter, and third measured angle parameter of the X-source, respectively, g1×sinβ1×cos(θ-α1), g1×sinβ1×sin(θ-α1), and -g1×cosβ1 represent the second, first, and third measured intensity parameters of the X-source, respectively, g2, β2, and α2 represent the measured intensity parameter, second measured angle parameter, and third measured angle parameter of the Y-source, respectively, g2×sinβ2×cos(θ-α2), g2×sinβ2×sin(θ-α2), and -g2×cosβ2 represent the second, first, and third measured intensity parameters of the Y-source, respectively.

[0123] It should be noted that when the seismic source includes both an X-source and a Y-source, the excitation directions of both sources are within the XOY plane. Therefore, the angle between the excitation directions of these two sources and the Z-component is 90°, meaning that the second measurement angle parameters β1 and β2 are both 90°. In this case, the first correction matrix can be expressed as:

[0124]

[0125] Computer devices may also arrange row and column elements into a first correction matrix in other orders, without specific limitations.

[0126] If the earthquake source includes X-source, Y-source, and Z-source, then a row contains three elements and a column contains three elements. For the elements in a row, the computer device can arrange them from left to right in the order of X-component, Y-component, and Z-component. For the elements in a column, the computer device can arrange them from top to bottom in the order of X-source, Y-source, and Z-source, forming the first correction matrix.

[0127] Accordingly, the first correction matrix can be expressed as:

[0128]

[0129] Where g3, β3, and α3 represent the measured intensity parameter, the second measured angle parameter, and the third measured angle parameter of the Z-source, respectively, and g3×sinβ3×cos(θ-α3), g3×sinβ3×sin(θ-α3), and -g3×cosβ3 represent the second measured intensity parameter, the first measured intensity parameter, and the third measured intensity parameter corresponding to the Z-source, respectively.

[0130] Step 107: The computer device determines the inverse matrix of the first correction matrix to obtain the second correction matrix.

[0131] If the first correction matrix is ​​represented as:

[0132]

[0133] The second correction matrix can then be expressed as:

[0134]

[0135] Among them, P -1 This represents the second correction matrix.

[0136] If the first correction matrix is ​​represented as:

[0137]

[0138] The second correction matrix can then be expressed as:

[0139]

[0140] In the embodiments of this application, the computer device can determine the inverse matrix of the first correction matrix by means of the adjoint matrix method, elementary transformation method and other methods, without making specific limitations.

[0141] It should be noted that before determining the inverse of the first correction matrix, the computer equipment first determines whether the first correction matrix is ​​invertible. If the seismic source includes an X-source and a Y-source, the process by which the computer equipment determines whether the first correction matrix is ​​invertible can be as follows: the computer equipment determines whether the excitation direction of the X-source is parallel to the excitation direction of the Y-source. If the excitation direction of the X-source is not parallel to the excitation direction of the Y-source, the computer equipment determines that the first correction matrix is ​​invertible, and then determines the inverse of the first correction matrix.

[0142] If the earthquake sources include X-source, Y-source, and Z-source, the process by which the computer equipment determines whether the first correction matrix is ​​invertible can be as follows: The computer equipment determines whether the excitation directions of X-source and Y-source are parallel, whether the excitation directions of Y-source and Z-source are parallel, and whether the excitation directions of X-source, Y-source, and Z-source are simultaneously in the same plane. If the excitation directions of X-source and Y-source are not parallel, the excitation directions of Y-source and Z-source are not parallel, and the excitation directions of X-source and Z-source are not simultaneously in the same plane, then the computer equipment determines that the first correction matrix is ​​invertible, and then determines the inverse matrix of the first correction matrix.

[0143] Step 108: The computer equipment determines the product of the first seismic data and the second correction matrix to obtain the second seismic data.

[0144] If the first earthquake data is six-component earthquake data, then the first earthquake data can be represented as:

[0145] Where U represents the first earthquake data, the row indices of the first and second rows represent the X-source and Y-source respectively, and the column indices of the first, second, and third columns represent the X-component, Y-component, and Z-component respectively. 11 U 12 and U 13 U represents the actual seismic data of the X-source in the X component, the actual seismic data of the X-source in the Y component, and the actual seismic data of the X-source in the Z component, respectively. 21 U 22 and U 23 These represent the actual earthquake data of the Y-source in the X component, the actual earthquake data of the Y-source in the Y component, and the actual earthquake data of the Y-source in the Z component, respectively.

[0146] The computer equipment determines the product of the first seismic data and the second correction matrix to obtain the theoretical second seismic data, which can be represented as: Where V represents the second earthquake data, V 11 V 12 and V 13 V represents the theoretical seismic data of the X-source in the X component, the theoretical seismic data of the X-source in the Y component, and the theoretical seismic data of the X-source in the Z component, respectively. 21 V 22 and V 23 These represent the actual earthquake data of the Y-source in the X component, the actual earthquake data of the Y-source in the Y component, and the actual earthquake data of the Y-source in the Z component, respectively.

[0147] If the first earthquake data is a nine-component earthquake data, then the first earthquake data can be represented as:

[0148] In this context, the subscript in the third row indicates the Z-source, U... 31 U 32 and U 33 These represent the actual earthquake data of the Z-source in the X component, the actual earthquake data of the Z-source in the Y component, and the actual earthquake data of the Z-source in the Z component, respectively.

[0149] The computer equipment determines the product of the first seismic data and the second correction matrix to obtain the second seismic data, which can be represented as: Among them, V 31 V 32 and V 33 These represent the theoretical seismic data of the Z-source in the X component, the theoretical seismic data of the Z-source in the Y component, and the theoretical seismic data of the Z-source in the Z component, respectively.

[0150] In this embodiment of the application, after the computer device obtains the second seismic data, it can determine seismic data with arbitrary excitation direction and arbitrary intensity based on the second seismic data. This process can be as follows: the computer device acquires a first preset angle parameter, a second preset angle parameter, and a preset intensity parameter, and determines the third seismic data based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter, and the preset intensity parameter.

[0151] The first preset angle parameter and the second preset angle parameter are used to determine the preset excitation direction. The first preset angle parameter is used to represent the angle between the preset excitation direction and the first vertical component. The second preset angle parameter is used to represent the angle between the preset excitation direction and the first horizontal component. The third seismic data is the seismic data corresponding to the second seismic data under the survey line direction, the preset excitation direction and the preset intensity parameters.

[0152] In this implementation, the computer device determines a relation matrix based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter, and the preset intensity parameter. The product of this relation matrix and the second seismic data is then determined to obtain the third seismic data. The process by which the computer device determines the relation matrix is ​​similar to the process described above for determining the first correction matrix, and will not be repeated here.

[0153] For example, if the first measurement angle parameter is still represented as θ, the first preset angle parameter is represented as β′, the second preset angle parameter is represented as a′, and the preset intensity parameter is represented as g′, then the preset excitation direction can be represented as (β′, α′).

[0154] If the second earthquake data is nine-component earthquake data, then the relation matrix can be represented as [g′×sinβ′×cos(θ-α′)g′×sinβ′×sin(θ-α′)-g′×cosβ′], and the third earthquake data can be represented as:

[0155] W = [g′×sinβ′×cos(θ-α′)g′×sinβ′×sin(θ-α′)-g′×cosβ′]V, where W represents the third earthquake data and V represents the second earthquake data.

[0156] If the second earthquake data is six-component earthquake data, then the relation matrix can be represented as [g′×cos(θ-α′)g′×sin(θ-α′)], and the third earthquake data can be represented as: W=[g′×cos(θ-α′)g′×sin(θ-α′)]V.

[0157] The second preset angle parameter, the third preset angle parameter, and the preset intensity parameter can be set and changed as needed, and no specific limitations are imposed on them.

[0158] In this embodiment of the application, after obtaining the theoretical second seismic data, the computer device can synthesize a third seismic data with three components of arbitrary excitation direction and arbitrary output intensity.

[0159] This application provides a seismic data correction method. This method uses a first measurement angle parameter, multiple second measurement angle parameters corresponding to multiple seismic sources, multiple third measurement angle parameters, and multiple measurement intensity parameters to determine a first measurement intensity parameter on a first horizontal component, a second measurement intensity parameter on a second horizontal component, and a third measurement intensity parameter on a first vertical component. Based on the first, second, and third measurement intensity parameters corresponding to the multiple measurement intensity parameters, a first correction matrix is ​​determined. This allows the actual first seismic data to be corrected using the first correction matrix to obtain theoretical seismic data, thereby improving the accuracy of the seismic data.

[0160] Figure 3 This is a schematic diagram of the structure of a seismic data correction device provided in an embodiment of this application. See also... Figure 3 The device includes:

[0161] The first acquisition module 301 is used to acquire first earthquake data based on multiple seismic sources;

[0162] The second acquisition module 302 is used to acquire the first measurement angle parameter corresponding to the first seismic data. The first measurement angle parameter is used to represent the angle between the survey line direction corresponding to the first seismic data and the first horizontal component.

[0163] The third acquisition module 303 is used to acquire multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters corresponding to multiple earthquake sources. The second measurement angle parameters are used to represent the angle between the excitation direction of the earthquake source and the first vertical component. The third measurement angle parameters are used to represent the angle between the excitation direction of the earthquake source and the first horizontal component. The measurement intensity parameters are used to represent the magnitude of the force excited by the earthquake source. One earthquake source corresponds to one second measurement angle parameter, one third measurement angle parameter, and one measurement intensity parameter.

[0164] The first determining module 304 is used to determine, based on the first measuring angle parameter, multiple second measuring angle parameters and multiple third measuring angle parameters, the first measuring intensity parameter on the first horizontal component and the second measuring intensity parameter on the second horizontal component respectively. The first measuring intensity parameter and the second measuring intensity parameter are used to represent the measuring intensity of the measuring intensity parameter on the first horizontal component and the measuring intensity on the second horizontal component respectively.

[0165] The second determining module 305 is used to determine, based on multiple second measuring angle parameters, a third measuring intensity parameter on the first vertical component of multiple measuring intensity parameters respectively. The first horizontal component, the second horizontal component, and the first vertical component are perpendicular to each other. The third measuring intensity parameter is used to represent the measuring intensity of the measuring intensity parameter on the first vertical component. One measuring intensity parameter corresponds to the first measuring intensity parameter, the second measuring intensity parameter, and the third measuring intensity parameter respectively.

[0166] The third determining module 306 is used to determine the first correction matrix based on the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to the multiple measurement intensity parameters respectively;

[0167] The correction module 307 is used to correct the first seismic data based on the first correction matrix to obtain the second seismic data.

[0168] In one possible implementation, the first determining module 304 is used to determine the sine value of the second measurement angle parameter to obtain a first parameter value; determine the difference between the first measurement angle parameter and the third measurement angle parameter to obtain a fourth angle parameter; determine the sine value of the fourth angle parameter to obtain a second parameter value; and determine the first measurement intensity parameter based on the first parameter value, the second parameter value, and the measurement intensity parameter.

[0169] In another possible implementation, the first determining module 304 is used to determine the product of the measured intensity parameter, the first parameter value, and the second parameter value to obtain the first measured intensity parameter.

[0170] In another possible implementation, the first determining module 304 is used to determine the cosine value of the fourth angle parameter to obtain the third parameter value; and based on the third parameter value, the first parameter value and the measurement intensity parameter, to determine the second measurement intensity parameter.

[0171] In another possible implementation, the second determining module 305 is used to determine the cosine value of the second measurement angle parameter to obtain the fourth parameter value; and based on the fourth parameter value and the measurement intensity parameter, to determine the third measurement intensity parameter.

[0172] In another possible implementation, the third determining module 306 is used to determine that the first, second, and third measurement intensity parameters corresponding to a measurement intensity parameter are elements of a row in the first correction matrix, thus obtaining multiple row elements; to determine that the first, second, and third measurement intensity parameters corresponding to the multiple measurement intensity parameters are elements of a column in the first correction matrix, thus obtaining multiple column elements; and to combine the multiple row elements and multiple column elements to form the first correction matrix.

[0173] In another possible implementation, the correction module 307 is used to determine the inverse of the first correction matrix to obtain the second correction matrix; and to determine the product of the second correction matrix and the first seismic data to obtain the second seismic data.

[0174] In another possible implementation, the device also includes:

[0175] The fourth acquisition module is used to acquire a first preset angle parameter, a second preset angle parameter, and a preset intensity parameter. The first preset angle parameter and the second preset angle parameter are used to determine a preset excitation direction. The first preset angle parameter is used to represent the angle between the preset excitation direction and the first vertical component, and the second preset angle parameter is used to represent the angle between the preset excitation direction and the first horizontal component.

[0176] The fourth determining module is used to determine the third seismic data based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter, and the preset intensity parameter. The third seismic data is the seismic data corresponding to the second seismic data under the survey line direction, the preset excitation direction, and the preset intensity parameter.

[0177] This application provides a seismic data correction device. The device uses a first measurement angle parameter, multiple second measurement angle parameters, multiple third measurement angle parameters, and multiple measurement intensity parameters to determine a first measurement intensity parameter on a first horizontal component, a second measurement intensity parameter on a second horizontal component, and a third measurement intensity parameter on a first vertical component. Based on the first, second, and third measurement intensity parameters corresponding to the multiple measurement intensity parameters, a first correction matrix is ​​determined. This allows the actual first seismic data to be corrected using the first correction matrix to obtain theoretical seismic data, thereby improving the accuracy of the seismic data.

[0178] It should be noted that the seismic data correction device provided in the above embodiments is only illustrated by the division of the above functional modules when correcting seismic data. In practical 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 seismic data correction device and the seismic data correction 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.

[0179] Figure 4 This illustration shows a structural block diagram of a computer device 400 provided in an exemplary embodiment of this application. The computer device 400 may be a portable mobile computer device, such as a smartphone, tablet computer, MP3 player (Moving Picture Experts Group Audio Layer III), MP4 player (Moving Picture Experts Group Audio Layer IV), laptop computer, or desktop computer. The computer device 400 may also be referred to as a user device, portable computer device, laptop computer device, desktop computer device, or other names.

[0180] Typically, computer device 400 includes a processor 401 and a memory 402.

[0181] Processor 401 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. Processor 401 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 401 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 401 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 401 may also include an AI (Artificial Intelligence) processor, which is used to handle computational operations related to machine learning.

[0182] The memory 402 may include one or more computer-readable storage media, which may be non-transitory. The memory 402 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 402 are used to store at least one line of program code, which is executed by the processor 401 to implement the seismic data correction method provided in the method embodiments of this application.

[0183] In some embodiments, the computer device 400 may optionally include a peripheral device interface 403 and at least one peripheral device. The processor 401, memory 402, and peripheral device interface 403 can be connected via a bus or signal line. Each peripheral device can be connected to the peripheral device interface 403 via a bus, signal line, or circuit board. Specifically, the peripheral device includes at least one of the following: a radio frequency circuit 404, a display screen 405, a camera assembly 406, an audio circuit 407, a positioning assembly 408, and a power supply 409.

[0184] Peripheral device interface 403 can be used to connect at least one I / O (Input / Output) related peripheral device to processor 401 and memory 402. In some embodiments, processor 401, memory 402 and peripheral device interface 403 are integrated on the same chip or circuit board; in some other embodiments, any one or two of processor 401, memory 402 and peripheral device interface 403 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0185] The radio frequency (RF) circuit 404 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 404 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 404 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals back into electrical signals. Optionally, the RF circuit 404 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, etc. The RF circuit 404 can communicate with other computer devices through at least one wireless communication protocol. This wireless communication protocol includes, but is not limited to: the World Wide Web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 404 may also include circuitry related to NFC (Near Field Communication), which is not limited in this application.

[0186] Display screen 405 is used to display a UI (User Interface). This UI may include graphics, text, icons, videos, and any combination thereof. When display screen 405 is a touch display screen, it also has the ability to collect touch signals on or above its surface. These touch signals can be input as control signals to processor 401 for processing. In this case, display screen 405 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, display screen 405 may be a single screen, disposed on the front panel of computer device 400; in other embodiments, display screen 405 may be at least two screens, disposed on different surfaces of computer device 400 or in a folded design; in other embodiments, display screen 405 may be a flexible display screen, disposed on a curved or folded surface of computer device 400. Furthermore, display screen 405 may be configured as a non-rectangular irregular shape, i.e., a non-rectangular screen. Display screen 405 may be made of materials such as LCD (Liquid Crystal Display) or OLED (Organic Light-Emitting Diode).

[0187] The camera assembly 406 is used to acquire images or videos. Optionally, the camera assembly 406 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is located on the front panel of the computer device, and the rear-facing camera is located on the back of the computer device. In some embodiments, there are at least two rear-facing cameras, which are any one of a main camera, a depth-sensing camera, a wide-angle camera, and a telephoto camera, to achieve background blurring by fusion of the main camera and the depth-sensing camera, panoramic shooting by fusion of the main camera and the wide-angle camera, VR (Virtual Reality) shooting, or other fusion shooting functions. In some embodiments, the camera assembly 406 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash is a combination of a warm-light flash and a cool-light flash, which can be used for light compensation at different color temperatures.

[0188] The audio circuit 407 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, converting the sound waves into electrical signals that are input to the processor 401 for processing, or input to the radio frequency circuit 404 for voice communication. For stereo sound acquisition or noise reduction purposes, multiple microphones may be used, each located in a different part of the computer device 400. The microphone may also be an array microphone or an omnidirectional microphone. The speaker is used to convert electrical signals from the processor 401 or the radio frequency circuit 404 into sound waves. The speaker may be a conventional diaphragm speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can convert electrical signals not only into audible sound waves but also into inaudible sound waves for purposes such as distance measurement. In some embodiments, the audio circuit 407 may also include a headphone jack.

[0189] The positioning component 408 is used to locate the current geographical location of the computer device 400 in order to enable navigation or LBS (Location Based Service). The positioning component 408 can be a positioning component based on the US GPS (Global Positioning System), China's BeiDou system, or Russia's Galileo system.

[0190] Power supply 409 is used to supply power to the various components in computer device 400. Power supply 409 can be alternating current, direct current, a disposable battery, or a rechargeable battery. When power supply 409 includes a rechargeable battery, the rechargeable battery can be a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is a battery that is charged via a wired line, while a wireless rechargeable battery is a battery that is charged via a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0191] In some embodiments, the computer device 400 further includes one or more sensors 410. The one or more sensors 410 include, but are not limited to: an accelerometer 411, a gyroscope 412, a pressure sensor 413, a fingerprint sensor 414, an optical sensor 415, and a proximity sensor 416.

[0192] Accelerometer 411 can detect the magnitude of acceleration along the three coordinate axes of a coordinate system established by computer device 400. For example, accelerometer 411 can be used to detect the components of gravitational acceleration along the three coordinate axes. Processor 401 can control display screen 405 to display the user interface in either a landscape or portrait view based on the gravitational acceleration signal acquired by accelerometer 411. Accelerometer 411 can also be used for games or for acquiring user motion data.

[0193] The gyroscope sensor 412 can detect the orientation and rotation angle of the computer device 400. The gyroscope sensor 412, in conjunction with the accelerometer sensor 411, can collect 3D motion data from the user on the computer device 400. Based on the data collected by the gyroscope sensor 412, the processor 401 can perform the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0194] The pressure sensor 413 can be disposed on the side bezel of the computer device 400 and / or on the lower layer of the display screen 405. When the pressure sensor 413 is disposed on the side bezel of the computer device 400, it can detect the user's grip signal on the computer device 400, and the processor 401 can perform left / right hand recognition or quick operation based on the grip signal collected by the pressure sensor 413. When the pressure sensor 413 is disposed on the lower layer of the display screen 405, the processor 401 can control the operable controls on the UI interface based on the user's pressure operation on the display screen 405. The operable controls include at least one of button controls, scroll bar controls, icon controls, and menu controls.

[0195] The fingerprint sensor 414 is used to collect a user's fingerprint. The processor 401 identifies the user based on the fingerprint collected by the fingerprint sensor 414, or vice versa. When the user's identity is verified as trusted, the processor 401 authorizes the user to perform relevant sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings. The fingerprint sensor 414 can be located on the front, back, or side of the computer device 400. When the computer device 400 has physical buttons or a manufacturer's logo, the fingerprint sensor 414 can be integrated with the physical buttons or the manufacturer's logo.

[0196] An optical sensor 415 is used to collect ambient light intensity. In one embodiment, the processor 401 can control the display brightness of the display screen 405 based on the ambient light intensity collected by the optical sensor 415. Specifically, when the ambient light intensity is high, the display brightness of the display screen 405 is increased; when the ambient light intensity is low, the display brightness of the display screen 405 is decreased. In another embodiment, the processor 401 can also dynamically adjust the shooting parameters of the camera assembly 406 based on the ambient light intensity collected by the optical sensor 415.

[0197] A proximity sensor 416, also known as a distance sensor, is typically mounted on the front panel of a computer device 400. The proximity sensor 416 is used to detect the distance between the user and the front of the computer device 400. In one embodiment, when the proximity sensor 416 detects that the distance between the user and the front of the computer device 400 is gradually decreasing, the processor 401 controls the display screen 405 to switch from a screen-on state to a screen-off state; when the proximity sensor 416 detects that the distance between the user and the front of the computer device 400 is gradually increasing, the processor 401 controls the display screen 405 to switch from a screen-off state to a screen-on state.

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

[0199] This application also provides a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the seismic data correction method of this application.

[0200] This application also provides a computer program product that stores at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the seismic data correction method described in this application.

[0201] In some embodiments, the computer program involved in the present application embodiments may be deployed and executed on a computer device, or executed on multiple computer devices located in one location, or executed on multiple computer devices distributed in multiple locations and interconnected through a communication network. Multiple computer devices distributed in multiple locations and interconnected through a communication network may constitute a blockchain system.

[0202] The above description is only for the purpose of enabling those skilled in the art to understand the technical solution 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 of seismic data correction, characterized by, The method comprises: obtaining first seismic data collected based on a plurality of seismic sources; obtaining a first measurement angle parameter corresponding to the first seismic data, the first measurement angle parameter being used to represent an included angle between a line direction corresponding to the first seismic data and a first horizontal component; obtaining a plurality of second measurement angle parameters, a plurality of third measurement angle parameters and a plurality of measurement intensity parameters corresponding to the plurality of seismic sources, the second measurement angle parameter being used to represent an included angle between an excitation direction of a seismic source and a first vertical component, the third measurement angle parameter being used to represent an included angle between the excitation direction of the seismic source and the first horizontal component, and the measurement intensity parameter being used to represent a magnitude of force excited by the seismic source, one seismic source corresponding to one second measurement angle parameter, one third measurement angle parameter and one measurement intensity parameter; based on the first measurement angle parameter, the plurality of second measurement angle parameters and the plurality of third measurement angle parameters, determining a first measurement intensity parameter of the plurality of measurement intensity parameters on the first horizontal component and a second measurement intensity parameter on a second horizontal component, the first measurement intensity parameter and the second measurement intensity parameter being used to represent measurement intensities of the measurement intensity parameter on the first horizontal component and on the second horizontal component respectively; based on the plurality of second measurement angle parameters, determining a third measurement intensity parameter of the plurality of measurement intensity parameters on a first vertical component, the first horizontal component, the second horizontal component and the first vertical component being perpendicular to each other, the third measurement intensity parameter being used to represent a measurement intensity of the measurement intensity parameter on the first vertical component, one measurement intensity parameter corresponding to a first measurement intensity parameter, a second measurement intensity parameter and a third measurement intensity parameter respectively; based on the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters respectively, determining a first correction matrix; based on the first correction matrix, correcting the first seismic data to obtain second seismic data.

2. The method of claim 1, wherein, The process of determining the first measurement intensity parameter of the measurement intensity parameter corresponding to a seismic source on the first horizontal component based on the first measurement angle parameter, the second measurement angle parameter and the third measurement angle parameter corresponding to the seismic source comprises: determining a sine value of the second measurement angle parameter to obtain a first parameter value; determining a difference value between the first measurement angle parameter and the third measurement angle parameter to obtain a fourth angle parameter; determining a sine value of the fourth angle parameter to obtain a second parameter value; determining the first measurement intensity parameter based on the first parameter value, the second parameter value and the measurement intensity parameter.

3. The method of claim 2, wherein, The process of determining the first measurement intensity parameter based on the first parameter value, the second parameter value and the measurement intensity parameter comprises: determining a product of the measurement intensity parameter, the first parameter value and the second parameter value to obtain the first measurement intensity parameter.

4. The method of claim 2, wherein, The process of determining a second measurement intensity parameter of the measurement intensity parameter corresponding to the seismic source on the second horizontal component based on the first measurement angle parameter, a second measurement angle parameter and a third measurement angle parameter corresponding to the seismic source comprises: determining a cosine value of the fourth angle parameter to obtain a third parameter value; determining the second measurement intensity parameter based on the third parameter value, the first parameter value and the measurement intensity parameter.

5. The method according to any one of claims 1 to 4, characterized in that, The process of determining a third measurement intensity parameter of the measurement intensity parameter corresponding to the seismic source on the first vertical component based on the second measurement angle parameter corresponding to the seismic source comprises: determining a cosine value of the second measurement angle parameter to obtain a fourth parameter value; determining the third measurement intensity parameter based on the fourth parameter value and the measurement intensity parameter.

6. The method of claim 1, wherein, The process of determining a first correction matrix based on the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters respectively comprises: determining that the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to one measurement intensity parameter are elements of one row in the first correction matrix to obtain a plurality of row elements; determining the first measurement intensity parameter corresponding to the plurality of measurement intensity parameters as one column, determining the second measurement intensity parameter corresponding to the plurality of measurement intensity parameters as one column, and determining the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters as one column to obtain a plurality of column elements; composing the plurality of row elements and the plurality of column elements into the first correction matrix.

7. The method of claim 1, wherein, The process of correcting the first seismic data based on the first correction matrix to obtain second seismic data comprises: determining an inverse matrix of the first correction matrix to obtain a second correction matrix; determining the product of the second correction matrix and the first seismic data to obtain the second seismic data.

8. The method of claim 1, wherein, The method further comprises: obtaining a first preset angle parameter, a second preset angle parameter and a preset intensity parameter, the first preset angle parameter and the second preset angle parameter being used to determine a preset excitation direction, the first preset angle parameter being used to represent the included angle between the preset excitation direction and the first vertical component, and the second preset angle parameter being used to represent the included angle between the preset excitation direction and the first horizontal component; determining third seismic data based on the first measurement angle parameter, the first preset angle parameter, the second preset angle parameter and the preset intensity parameter, the third seismic data being the seismic data corresponding to the second seismic data under the measurement line direction, the preset excitation direction and the preset intensity parameter.

9. A device for seismic data correction, characterized in that The device comprises: a first obtaining module configured to obtain first seismic data collected based on a plurality of seismic sources; a second obtaining module configured to obtain a first measurement angle parameter corresponding to the first seismic data, the first measurement angle parameter being used to represent the included angle between the measurement line direction corresponding to the first seismic data and a first horizontal component; a third obtaining module configured to obtain a second measurement angle parameter corresponding to the first seismic data, the second measurement angle parameter being used to represent the included angle between the measurement line direction corresponding to the first seismic data and a second horizontal component; The third acquisition module is configured to acquire a plurality of second measurement angle parameters, a plurality of third measurement angle parameters and a plurality of measurement intensity parameters corresponding to the plurality of seismic sources, the second measurement angle parameter is used to represent an included angle between an excitation direction of a seismic source and a first vertical component, the third measurement angle parameter is used to represent an included angle between the excitation direction of the seismic source and a first horizontal component, and the measurement intensity parameter is used to represent a magnitude of force excited by the seismic source, one seismic source corresponds to one second measurement angle parameter, one third measurement angle parameter and one measurement intensity parameter; The first determination module is configured to determine, based on the first measurement angle parameter, the plurality of second measurement angle parameters and the plurality of third measurement angle parameters, a first measurement intensity parameter of the plurality of measurement intensity parameters on the first horizontal component and a second measurement intensity parameter on a second horizontal component, the first measurement intensity parameter and the second measurement intensity parameter are respectively used to represent a measurement intensity of the measurement intensity parameter on the first horizontal component and a measurement intensity on the second horizontal component; The second determination module is configured to determine, based on the plurality of second measurement angle parameters, a third measurement intensity parameter of the plurality of measurement intensity parameters on a first vertical component, the first horizontal component, the second horizontal component and the first vertical component are perpendicular to each other, and the third measurement intensity parameter is used to represent a measurement intensity of the measurement intensity parameter on the first vertical component, one measurement intensity parameter corresponds to the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter respectively; The third determination module is configured to determine a first correction matrix based on the first measurement intensity parameter, the second measurement intensity parameter and the third measurement intensity parameter corresponding to the plurality of measurement intensity parameters respectively; The correction module is configured to correct the first seismic data based on the first correction matrix to obtain second seismic data.

10. A computer device, comprising: The computer device includes a processor and a memory, and the memory stores at least one program code, the at least one program code is loaded and executed by the processor to implement the seismic data correction method in any one of claims 1 to 8.

Images