Method and device for converting Q value of surface layer into Q field

Abstract

The invention provides a method and a device for converting a surface layer Q value into a Q field. The method comprises the following steps: acquiring basic data of all Q value survey points in a research area and processing basic data of a measuring network grid; performing regular gridding on the research area according to a range determined by processing the measurement grid; translating and rotating the coordinates of the Q value survey point; performing interpolation operation according to the position coordinates of all the Q value survey points after translation and rotation; and carrying out reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence. According to the method, the compensated single-shot recording frequency band is widened, the compensation effect is obvious, and the established surface Q field can be directly and efficiently applied to Q compensation processing of seismic data.

Description

Technical Field

[0001] The present invention belongs to the field of petroleum geophysical exploration, and particularly relates to the field of converting surface Q value to Q field. Specifically, it relates to a method and device for converting surface Q value to Q field. Background Art

[0002] The absorption attenuation factor Q value of the surface layer (hereinafter referred to as Q value) severely attenuates seismic waves and is one of the important factors restricting the quality of seismic data. How to establish an accurate surface Q field and effectively compensate seismic data is of great significance for improving the resolution of seismic data.

[0003] The existing methods for establishing the surface Q field include the Q value survey method, the fitting method, and the empirical formula method. The Q value survey method mainly uses the micro-logging method to survey the Q value of the surface layer in the study area to obtain the true Q value of the survey points. The fitting method for establishing the surface Q field usually fits the relationship between the true Q value of the Q value survey points and the surface velocity, establishes a formula for the true Q value and the surface velocity, and then substitutes the surface velocity of the surface Q value prediction points obtained from the first arrivals of surface seismic large guns or previous small refraction and micro-logging data in the study area into the formula for the true Q value and the surface velocity to obtain the true Q value of the Q value prediction points, forming a Q field. The empirical formula method for establishing the surface Q field usually substitutes the surface velocity of the Q value survey points and the surface Q value prediction points obtained from the first arrivals of surface seismic large guns or previous small refraction and micro-logging data in the study area into the empirical formula between the surface Q value and the surface velocity to obtain the relative Q value of the surface layer of the Q value survey points and the Q value prediction points, and then calibrates the relative Q value of the Q value survey points with the true Q value of the Q value survey points to obtain the calibration coefficient for converting the relative Q value of the surface layer to the true Q value of the surface layer. Then, the calibration coefficient is used to convert the relative Q value of the surface layer of the Q value prediction points to the true Q value of the Q value prediction points, thereby establishing the surface Q field. The surface Q field established by these methods is composed of discrete Q values and does not consider the grid information of the processing survey network in the study area. Usually, the established surface Q field cannot be directly used for the surface Q compensation processing of seismic data in the study area and needs to be further processed before it can be directly used for the compensation processing of seismic data. Summary of the Invention

[0004] The embodiments of the present application provide a method and device for converting surface Q value to Q field, which can directly and efficiently apply the established surface Q field to the Q compensation processing of seismic data.

[0005] First aspect, an embodiment of the present application provides a method for converting surface Q value to Q field, including: S1, obtaining the basic data of all Q value survey points in the study area and the basic data of processing the survey network grid; S2, regularly gridifying the study area according to the range determined by the processed survey network grid; S3, translating and rotating the coordinates of the Q value survey points; S4, performing interpolation operations according to the position coordinates of all the Q value survey points after translation and rotation; S5, performing reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

[0006] Among them, in step S1, obtaining the basic data of all Q value survey points in the study area includes: the position coordinates of the survey points and the attribute data corresponding to the position coordinates; the position coordinates (x, y) include an east coordinate sequence x and a north coordinate sequence y. The east coordinate sequence x is a vector, which is a set of east coordinates of all Q value survey points, and the north coordinate sequence y is also a vector, which is a set of north coordinates of all Q value survey points; the attribute data corresponding to the position coordinates of the Q value survey points is the true layer Q value sequence of the surface layer that is in one-to-one correspondence with the position coordinates (x, y) of the Q value survey points.

[0007] Obtaining the basic data of the processed survey network grid in the study area includes: the position coordinates of known points, the grid attribute data corresponding to the known points, and other grid data; the position coordinates (ref_x, ref_y) of the known points include a reference east coordinate ref_x and a reference north coordinate ref_y; the grid attribute data corresponding to the known points is the reference CMP line number ref_Xline and the reference CMP point number ref_inline; other grid data is the azimuth angle amz of the inline direction of the processed survey network in the study area; the CMP line number step Xline_inc and the CMP point number step inline_inc, and the CMP line number step and the CMP point number step respectively represent the actual distances dx and dy.

[0008] Among them, step S1 further includes: inputting the Q field range to be converted according to the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline.

[0009] Among them, in step S2, it includes: discretizing the CMP line number and the CMP point number according to the range determined by the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline of the processed survey network grid, with the CMP line number step Xline_inc and the CMP point number step inline_inc respectively.

[0010] The formula (1) for obtaining the CMP line number sequence Xline is:

[0011] Xline = {st_Xline, st_Xline + xline_inc, …, st_Xline + (n - 1)*xline_inc, ed_xline} (1)

[0013] The formula (2) for obtaining the CMP point number sequence inline is as follows:

[0014] inline = {st_inline, st_inline + inline_inc, …, st_inline + (m - 1)*inline_inc, ed_inline} (2)

[0016] Among them,

[0017] n ≤ (ed_xline - st_Xline) / xline_inc + 1, m ≤ (ed_inline - st_inline) / inline_inc + 1, and n, m ∈ N.

[0018] Among them, step S2 also includes: discretizing the east coordinate and the north coordinate according to the CMP line number sequence Xline corresponding to formula (1), the CMP point number sequence inline corresponding to formula (2), the reference CMP line number ref_Xline, the reference CMP point number ref_inline, and the actual distances dx and dy represented by the CMP line number step and the CMP point number step

[0019] When 45 ≤ amz ≤ 135, the formula for obtaining the east coordinate sequence x1 corresponding to the Xline sequence is:

[0020] x1 = {(st_xline - ref_Xline) / xline_inc*dx, (st_xline - ref_Xline) / xline_inc*dx + dx, …, (st_xline - ref_Xline) / xline_inc*dx + (n - 1)*dx, (ed_xline - ref_Xline) / xline_inc*dx} (3)

[0022] The formula for obtaining the north coordinate sequence y1 corresponding to the inline sequence;

[0023] y1 = (st_inline-ref_inline) / inline_inc*(-dy), (st_inline-ref_inline) / inline_inc*(-dy)+(-dy), …, (st_inline-ref_inline) / inline_inc*(-dy)+(m - 1)*(-dy), (ed_inline-ref_inline) / inline_inc*(-dy)}(4)

[0024] When -45 ≤ amz < 45, the obtained east coordinate sequence x1 is the same as formula (3), and the north coordinate sequence y1 is as follows:

[0025] y1 = (st_inline-ref_inline) / inline_inc*(dy), (st_inline-ref_inline) / inline_inc*(dy)+(dy), …, (st_inline-ref_inline) / inline_inc*(dy)+(m - 1)*(dy), (ed_inline-ref_inline) / inline_inc*(dy)}(5).

[0026] Among them, in step S2, it also includes: the east coordinate sequence x1 and the north coordinate sequence y1 form the grid node position coordinate sequence (X, Y) of a regular grid. The east coordinate sequence X of the grid node positions of the regular grid is as follows:

[0027]

[0028] The number of rows is M;

[0029] The north coordinate sequence Y of the grid node positions of the regular grid is as follows:

[0030] Y = (y1, y1,..., y1)

[0031] The number of columns is N.

[0032] Among them, in step S3, it includes: according to the position coordinate sequence (x, y) of all Q-value survey points, the translated and rotated east coordinate sequence x_out and north coordinate sequence y_out are calculated using formula (6),

[0033] Formula (6) is:

[0034]

[0035] When -45 ≤ amz < 45, a = -amz;

[0036] When 45 ≤ amz ≤ 135, a = 180 - amz

[0037] The direction of the east coordinate sequence x_out of all Q-value survey points after translation and rotation is consistent with the direction of the Xline sequence, and the direction of the north coordinate sequence y_out of all Q-value survey points after translation and rotation is consistent with the direction of the inline sequence.

[0038] Among them, in step S4, it includes: based on the fact that the Q-value sequence corresponding to the position coordinates (x_out, y_out) of all Q-value survey points after translation and rotation is the same as the Q-value sequence before translation and rotation, establishing the fitting relationship between the surface Q-value sequence Q(x_out, y_out) and the position coordinates (x_out, y_out) as follows:

[0039] Q(x_out, y_out) = F(x_out, y_out) (7)

[0040] F represents the established fitting function;

[0041] Using the relationship (7), calculate the Q values for the coordinate sequence (X, Y), and the Q-value sequence at the positions of the coordinate sequence (X, Y) is Q1(X, Y)

[0042] Q1(X, Y) = F(X, Y) (8).

[0043] Among them, in step S5, it includes: according to the coordinate sequence (X, Y), perform reverse rotation and translation according to formula (9) to calculate the new position coordinate sequence (x_amz, y_amz),

[0044]

[0045] The new position coordinate sequence (x_amz, y_amz) corresponds one-to-one with the coordinate sequence (X, Y), and Q2(x_amz, y_amz) = Q1(X, Y).

[0046] In a second aspect, an embodiment of the present application provides a device for converting surface Q values to a Q field, including:

[0047] An acquisition module, configured to acquire the basic data of all Q-value survey points in the study area and the basic data for processing the survey network grid;

[0048] A gridification module, configured to regularly gridify the study area according to the range determined by processing the survey network grid;

[0049] A translation and rotation module, configured to translate and rotate the coordinates of the Q-value survey points;

[0050] An interpolation module, configured to perform an interpolation operation according to the position coordinates of all Q-value survey points after translation and rotation.

[0051] A reverse translation and rotation module, configured to perform reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

[0052] In a third aspect, an embodiment of the present application provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method described in any one of the above are implemented.

[0053] In a fourth aspect, an embodiment of the present application provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, the steps of the method described in any one of the above are implemented.

[0054] The method and device for converting surface Q value to Q field in the embodiment of the present application have the following beneficial effects:

[0055] The method for converting surface Q value to Q field in the present application broadens the frequency band of the compensated single-shot record, and the compensation effect is obvious. The established surface Q field can be directly and efficiently applied to the Q compensation processing of seismic data. Description of the Drawings

[0056] Figure 1 It is a schematic flowchart of the method for converting surface Q value to Q field in the embodiment of the present application;

[0057] Figure 2 It is another schematic flowchart of the method for converting surface Q value to Q field in the embodiment of the present application;

[0058] Figure 3 It is a distribution diagram of surface discrete Q values in the study area of the embodiment of the present application;

[0059] Figure 4 It is a distribution diagram of the surface discrete Q values after rotation and translation in the study area of the embodiment of the present application;

[0060] Figure 5 It is a distribution diagram of Q values at the (X, Y) positions in the study area of the embodiment of the present application;

[0061] Figure 6 It is a distribution diagram of the Q field matching the processing survey network in the study area of the embodiment of the present application;

[0062] Figure 7 It is a comparison diagram of seismic single shots before and after surface Q compensation in the embodiment of the present application;

[0063] Figure 8 It is a comparison diagram of seismic single-shot spectra before and after surface Q compensation in the embodiment of the present application;

[0064] Figure 9 Schematic diagram of the device structure for converting the surface Q value to the Q field in the embodiments of the present application. Specific embodiments

[0065] The present application will be further introduced below in conjunction with the accompanying drawings and embodiments.

[0066] In the following description, the terms "first" and "second" are for descriptive purposes only and should not be construed as indicating or implying relative importance. The following description provides multiple embodiments of the present invention, and different embodiments can be replaced or combined. Therefore, the present application can also be considered to include all possible combinations of the same and / or different embodiments described. Thus, if one embodiment includes features A, B, and C, and another embodiment includes features B and D, then the present application should also be considered to include embodiments containing all other possible combinations of features A, B, C, and D, even though such embodiments may not be explicitly described in the following content.

[0067] Embodiment 1

[0068] As Figure 1 shown, the method for converting the surface Q value to the Q field provided by the embodiments of the present application includes: S1, obtaining the basic data of all Q value survey points in the study area and the basic data for processing the survey network grid; S2, regularly gridifying the study area according to the range determined by the processed survey network grid; S3, translating and rotating the coordinates of the Q value survey points; S4, performing interpolation operations based on the position coordinates of all the Q value survey points after translation and rotation; S5, performing reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

[0069] Using the method provided by the present application, the frequency band of the compensated single-shot record is broadened, and the compensation effect is obvious. The established surface Q field can be directly and efficiently applied to the Q compensation processing of seismic data.

[0070] Embodiment 2

[0071] The method for converting the surface Q value to the Q field provided by the embodiments of the present application, as Figure 2 shown, specifically includes the following contents:

[0072] (1) Obtaining basic data

[0073] 1) Obtaining the basic data of all Q value survey points in the study area, where the basic data of the Q value survey points includes the position coordinates of the survey points and the attribute data corresponding to the position coordinates of the Q value survey points; the position coordinates (x, y) are composed of the east coordinate sequence x and the north coordinate sequence y, that is, (x, y) = {(x1, y1), (x2, y2),..., (x n , y n)}, where the east coordinate sequence x is a vector, which is the set of the east coordinates of all Q-value survey points, i.e., x = {x1, x2, …, x n}; the north coordinate sequence y is also a vector, which is the set of the north coordinates of all Q-value survey points, i.e., y = {y1, y2, …, y n}; the attribute data corresponding to the position coordinates of the Q-value survey points is the true layer Q-value sequence Q(x, y) = {Q1(x1, y1), Q2(x2, y2), …, Q n (x n , y n )} that corresponds one-to-one with the position coordinates (x, y) of the Q-value survey points;

[0074] Among them, the subscript in the formula represents the Q-value survey points at different positions. The coordinates in this article are all plane projection coordinates, with the east coordinate being x and the north coordinate being y.

[0075] 2) Obtain the basic data of the processing survey network grid in the study area. The basic data of the processing survey network grid includes the position coordinates of known points, the grid attribute data corresponding to the known points, and other grid data; the position coordinates (ref_x, ref_y) of the known points are composed of the reference east coordinate ref_x and the reference north coordinate ref_y; the grid attribute data corresponding to the known points is the reference CMP line number ref_Xline and the reference CMP point number ref_inline; the other grid data is the azimuth angle amz of the inline direction of the processing survey network in the study area (the angle rotated clockwise from the grid north direction to the inline direction), the CMP line number step Xline_inc and the CMP point number step inline_inc, and the actual distances dx and dy represented by the CMP line number step and the CMP point number step respectively.

[0076] 3) Input the Q-field range to be transformed, i.e., the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline.

[0077] (2) Regular grid division of the study area

[0078] 1) Discretize the CMP line number and the CMP point number according to the range determined by the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline of the processing survey network grid, with the CMP line number step Xline_inc and the CMP point number step inline_inc respectively, to obtain the CMP line number sequence Xline represented by formula (1) and the CMP point number sequence inline represented by formula (2)

[0079] Xline = {st_Xline, st_Xline + xline_inc, …, st_Xline + (n - 1)*xline_inc, ed_xline} (1)

[0081] inline = {st_inline, st_inline + inline_inc, …, st_inline + (m - 1)*inline_inc, ed_inline} (2)

[0083] wherein,

[0084] n ≤ (ed_xline - st_Xline) / xline_inc + 1, m ≤ (ed_inline - st_inline) / inline_inc + 1, and n, m ∈ N

[0085] 2) Discretize the east coordinate and north coordinate according to the CMP line number sequence Xline corresponding to formula (1), the CMP point number sequence inline corresponding to formula (2), the reference CMP line number ref_Xline, the reference CMP point number ref_inline, and the actual distances dx and dy represented by the CMP line number step and CMP point number step.

[0086] When 45 ≤ amz ≤ 135, obtain the east coordinate sequence x1 corresponding one-to-one to the Xline sequence represented by formula (3) and the north coordinate sequence y1 corresponding one-to-one to the inline sequence represented by formula (4) respectively.

[0087] x1 = {(st_xline - ref_Xline) / xline_inc*dx, (st_xline - ref_Xline) / xline_inc*dx + dx, …, (st_xline - ref_Xline) / xline_inc*dx + (n - 1)*dx, (ed_xline - ref_Xline) / xline_inc*dx} (3)

[0089] y1 = (st_inline-ref_inline) / inline_inc*(-dy), (st_inline-ref_inline) / inline_inc*(-dy)+(-dy), …, (st_inline-ref_inline) / inline_inc*(-dy)+(m - 1)*(-dy), (ed_inline-ref_inline) / inline_inc*(-dy)} (4)

[0091] When -45 ≤ amz < 45, the obtained east coordinate sequence x1 is the same as formula (3), and the north coordinate sequence y1 is formula (5).

[0092] y1 = (st_inline-ref_inline) / inline_inc*(dy), (st_inline-ref_inline) / inline_inc*(dy)+(dy), …, (st_inline-ref_inline) / inline_inc*(dy)+(m - 1)*(dy), (ed_inline-ref_inline) / inline_inc*(dy)} (5)

[0094] 3) The grid node position coordinate sequence (X, Y) of the regular grid can be composed of the east coordinate sequence x1 and the north coordinate sequence y1. Among them, the east coordinate sequence X of the grid node position of the regular grid is as follows:

[0095] The number of rows is M,

[0096] The north coordinate sequence Y of the grid node position of the regular grid is as follows:

[0097] Y = (y1, y1,..., y1), and the number of columns is N.

[0098] (3) Coordinate translation and rotation of the Q-value survey points

[0099] For the position coordinate sequences (x, y) of all Q-value survey points, use formula (6) to subtract the reference east coordinate ref_x from all the east coordinate sequences x of the Q-value survey points, subtract the reference north coordinate ref_y from all the north coordinate sequences y of the Q-value survey points, and rotate clockwise by an angle a to obtain the east coordinate sequence x_out and north coordinate sequence y_out of all the Q-value survey points after translation and rotation; the direction of the east coordinate sequence x_out of all the Q-value survey points after translation and rotation is consistent with the direction of the Xline sequence, and the direction of the north coordinate sequence y_out of all the Q-value survey points after translation and rotation is consistent with the direction of the inline sequence.

[0100]

[0101] 1) When -45 ≤ amz < 45, a = -amz

[0102] 2) When 45 ≤ amz ≤ 135, a = 180 - amz.

[0103] (4) Regular grid Q-value interpolation

[0104] According to the Q-value sequence corresponding to the position coordinates (x_out, y_out) of all the Q-value survey points after translation and rotation described in (3), it is the same as the Q-value sequence before translation and rotation, that is, Q(x_out, y_out) = Q(x, y), where (x, y) and (x_out, y_out) are in one-to-one correspondence.

[0105] Establish a fitting relationship formula (7) between the surface Q-value sequence Q(x_out, y_out) and the position coordinates (x_out, y_out)

[0106] Q(x_out, y_out) = F(x_out, y_out) (7)

[0107] In formula (7), F represents the established fitting function.

[0108] Using the relationship formula (7), calculate the Q-values for the coordinate sequence (X, Y). The Q-value sequence at the positions of the coordinate sequence (X, Y) is Q1(X, Y)

[0109] Q1(X, Y) = F(X, Y) (8)

[0110] (5) Reverse translation and reverse rotation of regular grid coordinates

[0111] Perform reverse rotation and translation on the coordinate sequence (X, Y) according to formula (9) to obtain a new position coordinate sequence (x_amz, y_amz);

[0112]

[0113] The new position coordinate sequence (x_amz, y_amz) corresponds one-to-one with the coordinate sequence (X, Y). Therefore, the CMP line number sequence uniquely corresponding to the new position coordinate sequence (x_amz, y_amz) is the Xline sequence, the CMP point number sequence uniquely corresponding to it is the inline sequence; the Q value sequence uniquely corresponding to it is Q2(x_amz, y_amz), where Q2(x_amz, y_amz) = Q1(X, Y).

[0114] Therefore, the newly obtained Q value sequence Q2(x_amz, y_amz) is the Q field to be established in the study area. This Q field corresponds one-to-one with the position coordinates (x_amz, y_amz) of the study area, as well as the CMP line number, CMP point number, and azimuth angle of the processing survey network, that is, the establishment of a Q field matching the processing survey network is achieved.

[0115] In this application, through coordinate translation, rotation, and interpolation operations, the established surface Q field is matched with the processing survey network, and without other related processing, the established surface Q field can be directly and efficiently applied to the Q compensation processing of seismic data.

[0116] Example Three

[0117] As Figure 2 shown, the method for converting the surface Q value to a Q field in this application form includes:

[0118] (1) Obtain basic data

[0119] 1) Obtain the basic data of all Q value survey points in a certain study area. The basic data of the Q value survey points includes the position coordinates of the survey points and the attribute data corresponding to the position coordinates of the Q value survey points; the position coordinates (x, y) are composed of an east coordinate sequence x and a north coordinate sequence y, that is, (x, y) = {(x1, y1), (x2, y2)…, (x n , y n )}, where the east coordinate sequence x is a vector, which is the set of the east coordinates of all Q value survey points, that is, x = {x1, x2, …, x n}, and the north coordinate sequence y is also a vector, which is the set of the north coordinates of all Q value survey points, that is, y = {y1, y2, …, y n}; the attribute data corresponding to the position coordinates of the Q value survey points is the true layer Q value sequence Q(x, y) = {Q1(x1, y1), Q2(x2, y2), …, Q n (x n , y n )} that corresponds one-to-one with the position coordinates (x, y) of the Q value survey points; the subscripts in the above formulas represent Q value survey points at different positions; asFigure 3 As shown, for the surface discrete Q - value distribution map of the study area, the abscissa is the east coordinate and the ordinate is the north coordinate, and different colors represent the magnitude of the Q - value.

[0120] 2) Obtain the basic data of the processing survey network grid for a certain study area. The basic data of the processing survey network grid includes the coordinates of known positions, the grid attribute data corresponding to the coordinates of known positions, and other grid data. The coordinates (ref_x, ref_y) of the known position are composed of the reference east coordinate ref_x and the reference north coordinate ref_y. The reference CMP line number ref_Xline in the grid attribute data corresponding to the coordinates of the known position is 1, and the reference CMP point number ref_inline is 1. In the other grid data, the azimuth angle amz of the processing survey network of the study area is 131.9°, the CMP line number step Xline_inc is 1, the CMP point number step inline_inc is 1, and the actual distances represented by the CMP line number step and the CMP point number step are both 25m.

[0121] 3) Input the range of the Q - field to be transformed. The starting CMP line number st_Xline is 3500, the ending CMP line number ed_Xline is 4500, the starting CMP point number st_inline is 3000, and the ending CMP point number ed_inline is 4000.

[0122] (2) Regular grid - making of the study area

[0123] 1) Discretize the grid of the processing survey network according to the range determined by the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline, with the CMP line number step Xline_inc and the CMP point number step inline_inc to obtain the CMP line number sequence Xline represented by formula (1) and the CMP point number sequence inline represented by formula (2).

[0124] Xline = {3500, 3501, …, 4499, 4500}

[0125] inline = {3000, 3001, …, 3999, 4000}

[0126] 2) According to the CMP line number sequence Xline, CMP point number sequence inline described in 1), with the reference CMP line number ref_Xline being 1, the reference CMP point number ref_inline being 1, and the actual distances represented by the CMP line number step and CMP point number step both being 25m, discretize the east coordinate and north coordinate to obtain the east coordinate sequence x1 corresponding one-to-one with the Xline sequence and the north coordinate sequence y1 corresponding one-to-one with the inline sequence respectively;

[0127] x1 = {3499 * 25, 3500 * 25, …, 4498 * 25, 4499 * 25}

[0128] y1 = {-2999 * 25, -3000 * 25, … -3998 * 25, -3999 * 25}

[0129] 3) The east coordinate sequence x1 and north coordinate sequence y1 can form a regular grid, and the coordinate sequence of the grid nodes of the regular grid is (X, Y)

[0130]

[0131] Among them, the east coordinate sequence X of the regular grid is

[0132] The number of rows is 1001,

[0133] The north coordinate sequence Y of the regular grid is

[0134] The number of columns is 1001

[0135] (3) Coordinate translation and rotation of Q-value survey points

[0136] Since the azimuth of the processing survey network in this study area is 131.9 degrees, then a = 180 - 131.9 = 48.1; use formula (10) for the position coordinate sequence (x, y) of all Q-value survey points to subtract the reference east coordinate ref_x from the east coordinate sequence x of all Q-value survey points, subtract the reference north coordinate ref_y from the north coordinate sequence y of all Q-value survey points, and perform a clockwise rotation of 48.1 degrees to obtain the east coordinate sequence x_out and north coordinate sequence y_out of all Q-value survey points after translation and rotation; the direction of the east coordinate sequence x_out of all Q-value survey points after translation and rotation is consistent with the direction of the Xline sequence, and the direction of the north coordinate sequence y_out of all Q-value survey points after translation and rotation is consistent with the direction of the inline sequence; the Q-value distribution map after translation and rotation is as Figure 4 shown, with the abscissa being the east coordinate, the ordinate being the north coordinate, and different colors representing the magnitude of the Q value.

[0137]

[0138] (4) Regular grid Q-value interpolation

[0139] The Q-value sequence Q(x_out, y_out) at the position determined by the east coordinate x_out and north coordinate y_out described in (3) corresponds one-to-one with the Q-value sequence Q(x, y) before rotation and translation, that is, Q(x_out, y_out) = Q(x, y).

[0140] As Figure 5 shown, in this embodiment, the MATLAB's griddata function is used to obtain the Q value at the position of coordinates (X, Y).

[0141] (5) Reverse translation and reverse rotation of regular grid coordinates

[0142] The position coordinate sequence (X, Y) is reversely rotated and translated according to formula (11) to obtain a new position coordinate sequence (x_amz, y_amz).

[0143]

[0144] The new position coordinate sequence (x_amz, y_amz) corresponds one-to-one with the coordinate sequence (X, Y). Therefore, the Xline sequence is the unique CMP line number sequence corresponding to the new position coordinate sequence (x_amz, y_amz), the inline sequence is the unique CMP point number sequence corresponding to it; the Q-value sequence corresponding to it is Q2(x_amz, y_amz), where Q2(x_amz, y_amz) = Q1(X, Y).

[0145] Therefore, the newly obtained Q-value sequence Q2(x_amz, y_amz) is the Q field to be established in the study area. This Q field corresponds one-to-one with the position coordinates (x_amz, y_amz) of the study area, as well as the CMP line number, CMP point number, and azimuth angle of the processed survey grid, that is, the establishment of a Q field matching the processed survey grid is achieved. The Q field matching the processed survey grid is as Figure 6 shown.

[0146] In this application Figure 6 is the Q field distribution map within the range where the starting CMP line number is 3500, the ending CMP line number is 4500, the starting CMP point number is 3000, and the ending CMP point number is 4000 for the processed survey grid of a certain study area. This Q field matches the processed survey grid and can be directly used for Q compensation processing of seismic data. Figure 7 is obtained by using Figure 6The comparison diagram of the surface Q field before and after the seismic single-shot Q compensation processing shows that qualitatively, the reflection event axis of the single-shot record becomes thinner after the surface Q compensation, indicating that the frequency band of the single-shot record has been broadened after the compensation. Figure 8 For Figure 7 the frequency spectrum comparison diagram of the seismic single-shot before and after the surface Q compensation in [area], it can be seen that the main frequency of the seismic single-shot expands from 17 Hz to 29 Hz after the surface Q field compensation, indicating that the compensation effect is obvious and also proving the effectiveness of the Q value to Q field method proposed in this patent.

[0147] Example 4

[0148] As Figure 9 shown, the embodiment of the device for converting surface Q value to Q field provided by this application includes:

[0149] An acquisition module 201, configured to acquire the basic data of all Q value survey points in the study area and the basic data of the processing survey network grid;

[0150] A gridification module 202, configured to regularly gridify the study area according to the range determined by the processing survey network grid;

[0151] A translation and rotation module 203, configured to translate and rotate the coordinates of the Q value survey points;

[0152] An interpolation module 204, configured to perform interpolation operations according to the position coordinates of all the Q value survey points after translation and rotation;

[0153] A reverse translation and rotation module 205, configured to perform reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

[0154] In this application, the device embodiment for converting surface Q value to Q field is basically similar to the method embodiment for converting surface Q value to Q field. For the relevant parts, please refer to the introduction of the method embodiment for converting surface Q value to Q field.

[0155] This application also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method steps for converting the above-mentioned surface Q value to Q field are implemented.

[0156] An embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored. When the program is executed by a processor, the method steps of the above-mentioned conversion of the surface Q value to the Q field are implemented. Among them, the computer-readable storage medium may include, but is not limited to, any type of disk, including floppy disks, optical disks, DVDs, CD-ROMs, microdrives, and magneto-optical disks, ROMs, RAMs, EPROMs, EEPROMs, DRAMs, VRAMs, flash memory devices, magnetic cards or optical cards, nanosystems (including molecular memory ICs), or any type of medium or device suitable for storing instructions and / or data.

[0157] The above introduction is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for converting surface Q value to Q field, characterized in that, Including: S1. Obtain the basic data of all Q-value survey points in the study area and the basic data for processing the survey network grid; S2. Regularly grid the study area according to the range determined by the processed survey network grid; S3. Translate and rotate the coordinates of the Q-value survey points; S4. Perform interpolation operations based on the position coordinates of all Q-value survey points after translation and rotation; S5. Reverse-rotate and translate the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

2. The method for converting the surface Q value to a Q field according to claim 1, characterized in that In step S1, obtaining the basic data of all Q-value survey points in the study area includes: the position coordinates of the survey points and the attribute data corresponding to the position coordinates; the position coordinates (x, y) include an east coordinate sequence x and a north coordinate sequence y. The east coordinate sequence x is a vector, which is the set of east coordinates of all Q-value survey points, and the north coordinate sequence y is also a vector, which is the set of north coordinates of all Q-value survey points; the attribute data corresponding to the position coordinates of the Q-value survey points is the true layer Q-value sequence of the surface layer that corresponds one-to-one with the position coordinates (x, y) of the Q-value survey points; Obtaining the basic data of the processed survey network grid in the study area includes: the position coordinates of known points, the grid attribute data corresponding to the known points, and other grid data; the position coordinates (ref_x, ref_y) of the known points include a reference east coordinate ref_x and a reference north coordinate ref_y; the grid attribute data corresponding to the known points is the reference CMP line number ref_Xline and the reference CMP point number ref_inline; other grid data is the azimuth angle amz of the inline direction of the processed survey network in the study area; the CMP line number step Xline_inc and the CMP point number step inline_inc, and the actual distances dx and dy represented by the CMP line number step and the CMP point number step respectively.

3. The method for converting the surface Q value to a Q field according to claim 2, characterized in that, Step S1 also includes: inputting the Q-field range to be converted according to the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline.

4. The method for converting the surface Q value to the Q field according to any one of claims 1-3, characterized in that In step S2, it includes: discretizing the CMP line number and the CMP point number respectively with the CMP line number step Xline_inc and the CMP point number step inline_inc according to the range determined by the starting CMP line number st_Xline, the ending CMP line number ed_Xline, the starting CMP point number st_inline, and the ending CMP point number ed_inline for the grid of the processed survey network, The formula (1) for obtaining the CMP line number sequence Xline is: Xline = {st_Xline, st_Xline + xline_inc, …, st_Xline + (n - 1)*xline_inc, ed_xline} (1) The formula (2) for obtaining the CMP point number sequence inline is: inline = {st_inline, st_inline + inline_inc, …, st_inline + (m - 1)*inline_inc, ed_inline} (2) where, n ≤ (ed_xline - st_Xline) / xline_inc + 1, m ≤ (ed_inline - st_inline) / inline_inc + 1, and n, m ∈ N.

5. The method for converting the surface Q value to a Q field according to claim 4, characterized in that, In step S2, it further includes: discretizing the east coordinate and north coordinate according to the CMP line number sequence Xline corresponding to formula (1), the CMP point number sequence inline corresponding to formula (2), the reference CMP line number ref_Xline, the reference CMP point number ref_inline, and the actual distances dx and dy represented by the CMP line number step and CMP point number step. When 45 ≤ amz ≤ 135, the formula for obtaining the east coordinate sequence x1 corresponding one-to-one to the Xline sequence is: x1 = {(st_xline - ref_Xline) / xline_inc*dx, (st_xline - ref_Xline) / xline_inc*dx + dx, …, (st_xline - ref_Xline) / xline_inc*dx + (n - 1)*dx, (ed_xline - ref_Xline) / xline_inc*dx} (3) The formula for obtaining the north coordinate sequence y1 corresponding one-to-one to the inline sequence; y1 = (st_inline - ref_inline) / inline_inc*(-dy), (st_inline - ref_inline) / inline_inc*(-dy) + (-dy), …, (st_inline - ref_inline) / inline_inc*(-dy) + (m - 1)*(-dy), (ed_inline - ref_inline) / inline_inc*(-dy)} (4) When -45 ≤ amz < 45, the obtained east coordinate sequence x1 is the same as formula (3), and the north coordinate sequence y1 is: y1 = (st_inline - ref_inline) / inline_inc*(dy), (st_inline - ref_inline) / inline_inc*(dy) + (dy), …, (st_inline - ref_inline) / inline_inc*(dy) + (m - 1)*(dy), (ed_inline - ref_inline) / inline_inc*(dy)} (5).

6. The method for converting the surface Q value to a Q field according to claim 4, wherein In step S2, it further includes: the east coordinate sequence x1 and the north coordinate sequence y1 form the grid node position coordinate sequence (X, Y) of a regular grid. The east coordinate sequence X of the grid node positions of the regular grid is: The number of rows is M; The north coordinate sequence Y of the grid node positions of the regular grid is: Y = (y1, y1,..., y1) The number of columns is N.

7. The method for converting the surface Q value to the Q field according to any one of claims 1-3, characterized in that, Step S3 includes: According to the position coordinate sequences (x, y) of all Q-value survey points, use formula (6) to calculate the translated and rotated east coordinate sequence x_out and north coordinate sequence y_out. Formula (6) is: When -45 ≤ amz < 45, a = -amz; When 45 ≤ amz ≤ 135, a = 180 - amz The direction of the east coordinate sequence x_out of all Q-value survey points after translation and rotation is the same as the direction of the Xline sequence, and the direction of the north coordinate sequence y_out of all Q-value survey points after translation and rotation is the same as the direction of the inline sequence.

8. The method for converting the surface Q value to the Q field according to any one of claims 1-3, characterized in that, Step S4 includes: Based on the fact that the Q-value sequence corresponding to the position coordinates (x_out, y_out) of all Q-value survey points after translation and rotation is the same as the Q-value sequence before translation and rotation, establish the fitting relationship between the surface Q-value sequence Q(x_out, y_out) and the position coordinates (x_out, y_out) as: Q(x_out, y_out) = F(x_out, y_out) (7) F represents the established fitting function; Using relationship (7), obtain the Q-values for the coordinate sequence (X, Y). The Q-value sequence at the positions of the coordinate sequence (X, Y) is Q1(X, Y) Q1(X, Y) = F(X, Y) (8).

9. The method for converting the surface Q value to the Q field according to any one of claims 1-3, characterized in that, Step S5 includes: According to the coordinate sequence (X, Y), perform reverse rotation and translation according to formula (9) to calculate the new position coordinate sequence (x_amz, y_amz). The new position coordinate sequence (x_amz, y_amz) corresponds one-to-one with the coordinate sequence (X, Y), and Q2(x_amz, y_amz) = Q1(X, Y).

10. A device for converting surface Q value to Q field, characterized in that, Includes: An acquisition module for acquiring the basic data of all Q-value survey points in the study area and the basic data for processing the survey network grid; A gridification module for regularly gridifying the study area according to the range determined by processing the survey network grid; A translation and rotation module for translating and rotating the coordinates of the Q-value survey points; An interpolation module for performing interpolation operations according to the position coordinates of all Q-value survey points after translation and rotation; A reverse translation and rotation module for performing reverse rotation and translation on the coordinate sequence (X, Y) to obtain a new position coordinate sequence.

Images