Method, device and equipment for constructing geological parameter model

Abstract

The invention relates to the technical field of seismic exploration data processing, in particular to a geological parameter model construction method, device and equipment. The geological parameter model construction method comprises the following steps: acquiring first grid points of a plurality of earthquake work areas in a contiguous space; each earthquake work area is provided with a first geological parameter model, and one first grid point corresponds to line numbers, point numbers and initial sampling values of a plurality of first geological parameter models; setting a local processing window of the first grid point according to the line number and the point number corresponding to the first grid point; grid points except the first grid points in the local processing window are used as second grid points, and the weight of the second grid points is determined in the local processing window according to the distance between the second grid points and the first grid points; in the local processing window, calculating a local weighted reconstruction value of the first grid point according to the weight of the second grid point; and generating a second geological parameter model of the contiguous space according to the local weighted reconstruction value of the first grid point.

Description

Technical Field

[0001] The embodiments of this specification relate to the technical field of seismic exploration data processing, and specifically to a method, device and equipment for constructing a geological parameter model. Background Art

[0002] With the advancement of seismic acquisition technology, multiple rounds of repeated 3D seismic acquisition have been carried out in key oil and gas exploration areas, especially in mature exploration areas, from narrow azimuth acquisition to wide azimuth acquisition, and then to "two wide and one high" acquisition. Each seismic block has carried out high-quality seismic processing and interpretation, target optimization, and well location support, and has a high-quality seismic data foundation. With the continuous deepening of oil and gas exploration and development and the continuous increase in exploration and development efforts in recent years, the demand for regional geological law research and large-scale oil and gas discovery has become increasingly urgent, and the exploration and development of domestic oil and gas fields is gradually developing towards large-scale 3D seismic exploration. Due to the continuous increase in seismic acquisition costs, the way to obtain large-area 3D seismic data is mainly achieved through the contiguous processing of seismic block data from different acquisition years and observation methods. This is currently the most economical and effective solution.

[0003] To efficiently complete the contiguous processing, it is necessary to make full use of the existing seismic processing results of a single seismic work area. However, considering that the scopes of different seismic work areas in the contiguous area often overlap, the overlapping range is larger for secondary and higher three-dimensional seismic acquisitions, and there are also large differences in the observation methods, including different acquisition face elements, inconsistent azimuths, irregular observation ranges and many other complex situations. The conventional method of constructing the initial parameter model of the contiguous spatial seismic work area obtains the initial parameter models of velocity, density, anisotropy, etc. by carrying out a new round of contiguous work area processing.

[0004] The conventional method of constructing the initial parameter model of a contiguous space seismic work area cannot construct the initial parameter model of the contiguous space based on the existing initial parameter models of multiple seismic work areas in the contiguous space. It requires a large investment of manpower, material and financial resources, and is time-consuming, which cannot meet the needs of the oil and gas exploration and development production cycle.

[0005] Therefore, how to overcome the problems of large consumption of manpower, material and financial resources and low efficiency in the existing methods of constructing initial parameter models of contiguous spatial seismic work areas, and propose a method for constructing a contiguous spatial parameter model with low consumption of manpower, material and financial resources, short processing cycle and high application efficiency to meet the needs of the oil and gas exploration and development production cycle is a key issue that needs to be solved urgently. Summary of the invention

[0006] The purpose of the embodiments of this specification is to provide a method, device and equipment for constructing a geological parameter model to overcome the problems of large consumption of manpower, material resources and financial resources and low efficiency in the existing methods for constructing initial parameter models for contiguous spatial seismic work areas, shorten the processing cycle for constructing initial parameter models for contiguous spatial seismic work areas and improve the efficiency of the application, thereby meeting the production cycle requirements of oil and gas exploration and development.

[0007] On the one hand, an embodiment of the present specification proposes a method for constructing a geological parameter model, which includes: obtaining first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; according to the line number and point number corresponding to the first grid point, a local processing window of the first grid point is set; the grid points other than the first grid point in the local processing window are taken as second grid points, and within the local processing window, the weight of the second grid point is determined according to the distance between the second grid point and the first grid point; within the local processing window, the local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; and according to the local weighted reconstruction value of the first grid point, a second geological parameter model of the contiguous space is generated.

[0008] On the other hand, a device for constructing a geological parameter model includes: an acquisition module, which is used to acquire first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to line numbers, point numbers and initial sampling values ​​of multiple first geological parameter models; a setting module, which is used to set a local processing window of the first grid point according to the line number and point number corresponding to the first grid point; a determination module, which is used to take the grid points other than the first grid point in the local processing window as second grid points, and determine the weight of the second grid point within the local processing window according to the distance between the second grid point and the first grid point; a calculation module, which is used to calculate the local weighted reconstruction value of the first grid point within the local processing window according to the weight of the second grid point; and a generation module, which is used to generate a second geological parameter model of the contiguous space according to the local weighted reconstruction value of the first grid point.

[0009] On the other hand, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the above-mentioned method for constructing a geological parameter model.

[0010] It can be seen from the technical solutions provided in the embodiments of the present specification that the method for constructing a geological parameter model provided in the embodiments of the present specification can obtain the first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; according to the line number and point number corresponding to the first grid point, a local processing window of the first grid point is set; the grid points other than the first grid point in the local processing window are taken as second grid points, and within the local processing window, the weight of the second grid point is determined according to the distance between the second grid point and the first grid point; within the local processing window, the local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; according to the local weighted reconstruction value of the first grid point, the second geological parameter model of the contiguous space is generated. Compared with the existing method, the local weighted reconstruction of the parameters in the parameter model can be performed based on the local processing area of ​​the grid points in the contiguous space, which is faster and more efficient, and can meet the relevant oil and gas exploration and development production cycle requirements. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the technical solutions in the embodiments of this specification or the prior art, the drawings required for use in the embodiments or the prior art description are briefly introduced below.

[0012] Figure 1 It is a flow chart of a method for constructing a geological parameter model provided in an embodiment of this specification;

[0013] Figure 2 It is a schematic diagram of the spatial geographical location distribution of two earthquake work areas in the A contiguous space provided in the embodiment of this specification;

[0014] Figure 3 It is a partially enlarged schematic diagram of the overlapping positions of two seismic work areas in the A contiguous space provided in the embodiment of this specification;

[0015] Figure 4 is a partially enlarged schematic diagram of the reconstructed A contiguous space provided in the embodiment of this specification;

[0016] Figure 5 is a schematic diagram of the initial velocity model cross section of the A contiguous space after reconstruction provided in an embodiment of this specification;

[0017] Figure 6 It is a schematic diagram of the structure of a device for constructing a geological parameter model provided in an embodiment of this specification;

[0018] Figure 7 It is a schematic diagram of the structural composition of a computer device provided in an embodiment of this specification. DETAILED DESCRIPTION

[0019] The following will be combined with the drawings in the embodiments of this specification to clearly and completely describe the technical solutions in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of this specification.

[0020] Figure 1 The following is a flow chart of a method for constructing a geological parameter model provided in an embodiment of this specification. In specific implementation, the method includes the following steps:

[0021] S101: Acquire first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to line numbers, point numbers and initial sampling values ​​of multiple first geological parameter models.

[0022] In some embodiments, first geological parameter models of multiple seismic work areas in a contiguous space may be obtained; and based on the first geological parameter models of the multiple seismic work areas, the multiple seismic work areas may be divided into a first seismic work area and multiple second seismic work areas.

[0023] By obtaining the first geological parameter models of multiple seismic work areas in a contiguous space, the multiple seismic work areas can be divided into a first seismic work area and multiple second seismic work areas, and then the first geological parameter model of the second seismic work area can be adjusted based on the first geological parameter model of the first seismic work area. Compared with the second seismic work area, the seismic data related to the first seismic work area is of better quality and has a better geographical location, which can improve the speed and efficiency of the reconstruction process.

[0024] The contiguous space may refer to a whole exploration area formed by connecting multiple adjacent or close seismic work areas during the seismic exploration process. The distance between the centers of each seismic work area in the contiguous space is less than or equal to the preset distance threshold. The first geological parameter model of multiple seismic work areas in the contiguous space can be obtained. The first geological parameter model of the seismic work area may include a velocity model in the time domain or depth domain, an anisotropy model, a density model, an inclination model, and parameter updates corresponding to various models. Each seismic work area in the contiguous space is divided into grids, and each grid uniquely corresponds to the line number, point number, and initial sampling value of multiple geological parameter models. For example, for the velocity model, each grid uniquely corresponds to the line number, point number, and initial velocity value of the velocity model. A seismic work area can be selected from multiple seismic work areas in the contiguous space, and the grid of the selected seismic work area is preset as a unified grid, that is, as the grid division benchmark of the contiguous space. A first seismic work area can be selected based on the seismic data quality and geographical location of each seismic work area in the contiguous space, and the remaining seismic work areas in the contiguous space are used as the second seismic work area.

[0025] The quality of seismic data in each seismic area can be evaluated, including frequency analysis, signal-to-noise ratio analysis, coverage number analysis, and interference wave analysis. Frequency analysis is to perform frequency scanning analysis on the seismic data in each seismic area to determine the effective bandwidth and bandwidth. The clarity and continuity of reflection information in different frequency bands can be evaluated, and seismic areas with clear reflection information and good continuity can be preferred. Signal-to-noise ratio analysis is to calculate the signal-to-noise ratio of each seismic area and evaluate the continuity of the main target layer. Seismic areas with high signal-to-noise ratio and good continuity of the main target layer can be preferred. Coverage number analysis is to analyze the number of data coverage in each seismic area to ensure that the coverage number is uniform and meets the analysis requirements. Seismic areas with uniform and high coverage number can be preferred. Interference wave analysis is to identify the types of interference waves in each seismic area, such as surface waves, shallow refraction waves, multiple waves, etc. The impact of interference waves on the quality of seismic data can be evaluated, and seismic areas with less interference wave impact can be preferred.

[0026] The geographical location of each earthquake work area can be evaluated, including centrality analysis and transportation and infrastructure analysis. Centrality analysis is to determine the geographical location of each earthquake work area using geographic information system technology. The distance from each earthquake work area to the center of the contiguous range can be calculated, and earthquake work areas closer to the center can be preferred. Transportation and infrastructure analysis considers the transportation convenience and infrastructure conditions of the earthquake work area, such as roads, communications, electricity, etc. Earthquake work areas with convenient transportation and complete infrastructure can be preferred to facilitate subsequent data collection and processing.

[0027] In some embodiments, a comprehensive evaluation index system for seismic work areas can be established based on the results of seismic data quality and geographical location evaluation. The index system may include effective bandwidth, signal-to-noise ratio, coverage times, interference wave impact, geographical location centrality, traffic convenience and infrastructure conditions, etc. The comprehensive evaluation index system can be used to comprehensively evaluate each seismic work area. According to the evaluation results, the seismic work areas are sorted, and the seismic work area with the best comprehensive evaluation results is selected as the first seismic work area. According to the seismic data grid of the first seismic work area selected, the direction, vertices, coordinate data, grid vertical and horizontal step lengths, point numbers and line numbers of the unified grid of the contiguous space are determined. The unified grid is a grid range composed of four-point coordinates, and the four-point coordinates of the work area specifically include 1) minimum line number, minimum point number, X and Y; 2) minimum line number, maximum point number, X and Y; 3) maximum line number, minimum point number, X and Y; 4) maximum line number, maximum point number, X and Y. The maximum and minimum line numbers and point numbers can be adjusted according to the numerical range, and their values ​​are usually adjusted to positive integers greater than 0.

[0028] In some embodiments, the fourth grid point of the second seismic working area can be adjusted according to the third grid point of the first seismic working area; the first grid point of multiple seismic working areas in the contiguous space can be determined according to the third grid point of the first seismic working area and the fourth grid points of the multiple second seismic working areas.

[0029] By dividing the three-dimensional seismic work area in a grid manner, a plurality of grid points of the three-dimensional seismic work area can be obtained. The third grid point is a grid point obtained by dividing the first seismic work area in a grid manner. The fourth grid point is a grid point obtained by dividing the second seismic work area in a grid manner.

[0030] By adjusting the grid points of the second seismic work area and determining the first grid points of multiple seismic work areas in the contiguous space, the initial unification of the grid data in the contiguous space was achieved, laying a data foundation for the further reconstruction of the grid data in the contiguous space.

[0031] According to the first seismic work area selected, the line number, point number and initial sampling value of the fourth grid point of the second seismic work area can be adjusted according to the line number, point number and initial sampling value of the third grid point of the first seismic work area. The third grid point of the first seismic work area and the adjusted fourth grid points of the plurality of second seismic work areas can be considered to be based on the same grid reference. Therefore, according to the third grid point of the first seismic work area and the adjusted fourth grid points of the plurality of second seismic work areas, the first grid points of the plurality of seismic work areas in the contiguous space can be determined.

[0032] In some embodiments, the line number and point number of the fourth grid point of the second seismic work area may be adjusted according to the line number and point number of the third grid point of the first seismic work area.

[0033] By adjusting the line number and point number of the fourth grid point in the second seismic work area, a data foundation is laid for the unification and reconstruction of grid data in the contiguous space.

[0034] According to the line number and point number of the third grid point of the first seismic work area, the four-point coordinates of the grid area of ​​the first seismic work area can be determined. According to the four-point coordinates of the grid area of ​​the first seismic work area, the line number and point number of the fourth grid point of the second seismic work area can be regenerated, and then the four-point coordinates of the unified grid of the entire contiguous space can be determined. The line number and point number of the four-point coordinates of the unified grid of the entire contiguous space are integers. If the maximum and minimum line numbers and point numbers become negative numbers, the line number and point number need to be adjusted to make their values ​​become positive numbers greater than 0. The specific method is to add a positive integer greater than the absolute value of the negative number, which will not be repeated here.

[0035] In some embodiments, the initial sampling value of the fourth grid point of the second seismic working area may be adjusted according to the initial sampling value of the third grid point of the first seismic working area.

[0036] By adjusting the initial sampling value of the fourth grid point in the second seismic work area, a data foundation is laid for the unification and reconstruction of grid data in the contiguous space.

[0037] Taking the value range of the initial sampling value of the preferred first seismic work area as the standard, the value range of the initial sampling value in the initial parameter model of the multiple second seismic work areas can be adjusted, and finally the value range of all seismic work areas in the contiguous space is made consistent. The adjustment method is specifically to obtain the maximum and minimum values ​​of the value range of the initial sampling value of the preferred first seismic work area, compare the initial sampling value in the initial parameter model of the multiple second seismic work area models with the aforementioned maximum and minimum values, if it is less than the aforementioned minimum value, then the initial sampling value in the initial parameter model of the second seismic work area is set to the aforementioned minimum value, if it is greater than the aforementioned maximum value, then the initial sampling value in the initial parameter model of the second seismic work area is set to the aforementioned maximum value.

[0038] In some embodiments, line numbers and point numbers of the plurality of first grid points may be generated according to line numbers and point numbers of the third grid points and the fourth grid points.

[0039] By determining the first grid point, the initial unification of grid data in the contiguous space is achieved, laying a data foundation for further reconstruction of grid data in the contiguous space.

[0040] According to the third grid point of the first seismic work area and the adjusted fourth grid points of the multiple second seismic work areas, the first grid points of the multiple seismic work areas in the contiguous space can be determined. The third grid point and the adjusted fourth grid point are determined based on the same standard. Therefore, according to the line number, point number and initial sampling value of the third grid point of the first seismic work area and the line number, point number and initial sampling value of the fourth grid point of the multiple second seismic work areas, the line number, point number and initial sampling value of all the first grid points in the contiguous space can be determined, that is, the third grid point of the first seismic work area and the fourth grid points of the multiple second seismic work areas can be directly used as the corresponding first grid points. Thus, the first grid points of the multiple seismic work areas in the contiguous space include the third grid point of the first seismic work area and the adjusted fourth grid points of each second seismic work area.

[0041] S102: Setting a local processing window of the first grid point according to the line number and the point number corresponding to the first grid point.

[0042] By setting a local processing window for the first grid point, the first grid point can be better characterized based on the set local processing window, thus laying a foundation for accurate reconstruction of the first grid point.

[0043] For any first grid point in multiple seismic work areas in a contiguous space, the position of the first grid point can be determined based on the line number, point number and initial sampling value of the first grid point. The line number and point number are dimensions in two spatial directions, and the initial sampling value is the dimension in the vertical direction. A grid point and its position can be uniquely determined based on a set of line numbers, point numbers and initial sampling values. A local processing window containing A parameter model data channels and B grid points can be set with the position of the first grid point as the center. A and B are both preset positive integers. For example, for a given first grid point, A can be set to 10, B can be set to 4, and a local processing window containing 10 geological parameter model data channels and 4 grid points can be set with the position of the given first grid point as the center.

[0044] S103: Taking grid points other than the first grid point in the local processing window as second grid points, and determining the weight of the second grid point in the local processing window according to the distance between the second grid point and the first grid point.

[0045] In some embodiments, the weight of the second grid point can be calculated using a preset weight model according to the distance between the second grid point and the first grid point; the weights of all second grid points in the local processing window satisfy the formula Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; ω i,j The spatial coordinates are (x i ,y i ) is the weight of the second grid point.

[0046] The weight of the second grid point is calculated by using the distance between the second grid point and the first grid point, which lays a data foundation for the accurate reconstruction of the first grid point.

[0047] The azimuth of all seismic work areas in the contiguous space can be rotated to 0 degrees, and then the line number and point number of the first grid point can be subtracted from the line number and point number of any grid point in the contiguous space after the rotation. The obtained difference can be multiplied by the face element of the line number space direction and the point number space direction of the first grid point, and the product of the line number space direction and the point number space direction of the first grid point can be added to the spatial coordinates corresponding to the line number and point number of any grid point in the contiguous space, to obtain the coordinates of all seismic work areas in the contiguous space when the azimuth is 0 degrees. The coordinates of the first grid point obtained when the azimuth of all seismic work areas in the contiguous space is 0 degrees are rotated to the azimuth of the seismic work area where the first grid point is located, that is, the spatial coordinates corresponding to the line number and point number of the first grid point are obtained.

[0048] Activation functions such as Sigmod function, ReLU function and LeakyReLU function can be used as preset weight models. A weight model can also be trained based on machine learning models such as KNN, BP, random forest and decision tree, which will not be repeated here. The distance between each second grid point and the first grid point in the local processing window can be calculated, and the calculated distance between each second grid point and the first grid point can be input into the preset weight model, and the preset weight model can output the initial weight coefficient of each second grid point. The initial weight coefficients of each second grid point in the local processing window can be added to obtain the initial weight sum, and then the initial weight coefficients of each second grid point in the local processing window can be divided by the initial weight sum to obtain the weight of each second grid point. The weight of each second grid point in the local processing window satisfies the following formula:

[0049]

[0050] Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; ω i,j The spatial coordinates are (x i ,y i The sum of the weights of the second grid points in the local processing window is 1, which helps to ensure the rationality and accuracy of the initial sampling values ​​corresponding to the first grid points of the subsequent reconstruction.

[0051] For example, a Sigmod function can be used as a preset weight model. The distance between each second grid point and the first grid point in the local processing window can be calculated, and the calculated distance between each second grid point and the first grid point is input into the Sigmod function, and the Sigmod function can output the initial weight coefficient of each second grid point. The initial weight coefficients of each second grid point in the local processing window can be added to obtain an initial weight sum, and then the initial weight coefficients of each second grid point in the local processing window can be divided by the initial weight sum to obtain the weight of each second grid point.

[0052] S104: Calculate local weighted reconstruction values ​​of the first grid points according to the weights of the second grid points within the local processing window.

[0053] In some embodiments, the formula may be used based on the local weighted sampling value of the second grid point and the weight of the second grid point. Calculate the local weighted reconstruction value of the first grid point. Wherein, N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X coordinate and Y coordinate corresponding to the first grid point, respectively; R(x, y) is the local weighted reconstruction value of the first grid point; x i and j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,j The spatial coordinates are (x i ,y i )’s second grid point weight; R i,j The spatial coordinates are (x i ,y i ) is the local weighted sampling value of the second grid point.

[0054] By calculating the local weighted reconstruction value of the first grid point using the local weighted sampling value and weight of the second grid point in the local processing window, accurate sparsification of data encrypted locations in the contiguous space and accurate encryption of data sparse locations in the contiguous space can be achieved.

[0055] According to the local weighted sampling values ​​and weights of all the second grid points in the local processing window, the local weighted reconstruction value of the first grid point can be calculated using the following formula:

[0056]

[0057] Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X-coordinate and Y-coordinate corresponding to the first grid point, respectively; R(x, y) is the local weighted reconstruction value of the first grid point; x i and j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,jThe spatial coordinates are (x i ,y i )’s second grid point weight; R i,j The spatial coordinates are (x i ,y i ). The local weighted reconstruction value of the first grid point is calculated using the local weighted sampling value of the second grid point in the local processing window and the corresponding weight. The local weighted reconstruction value calculated in this way can not only make full use of the first geological parameter model data, but also generate accurate and reasonable local weighted reconstruction values ​​based on the local processing window.

[0058] In some embodiments, the L second grid points closest to the second grid point in the local processing window can be obtained as the fifth grid point; L is a preset positive integer; based on the initial sampling values ​​of the L fifth grid points, the local weighted sampling value of the second grid point can be determined.

[0059] By determining the local weighted sampling values ​​of the second grid points, the influence of abnormal initial sampling values ​​and data noise on the subsequent reconstruction process is further alleviated.

[0060] For each second grid point in the local processing window, the distance between each second grid point and the remaining second grid points in the local processing window can be calculated. Take the nearest L second grid points as the fifth grid point corresponding to the current grid point. Activation functions such as Sigmod function, ReLU function and LeakyReLU function can be used as the preset weight model. A weight model can also be trained based on machine learning models such as KNN, BP, random forest and decision tree. The weight of each fifth grid point can be calculated based on the preset weight model. Specifically, the calculated distance between each fifth grid point and the second grid point can be input into the preset weight model, and the preset weight model can output the initial weight coefficient of each fifth grid point. The initial weight coefficients of each fifth grid point in the local processing window can be added to obtain the initial weight sum, and then the initial weight coefficients of each fifth grid point in the local processing window can be divided by the initial weight sum to obtain the weight of each fifth grid point. According to the initial sampling value and weight corresponding to each fifth grid point, the weighted sum of all fifth grid points can be calculated, and then the calculated weighted sum can be used as the local weighted sampling value of each second grid point.

[0061] S105: Generate a second geological parameter model of the contiguous space according to the local weighted reconstruction values ​​of the first grid points.

[0062] In some embodiments, the first geological parameter models of the multiple seismic work areas can be updated according to the local weighted reconstruction values ​​of the first grid points; and the second geological parameter model of the contiguous space can be generated according to the updated first geological parameter models of the multiple seismic work areas.

[0063] The second geological parameter model of the contiguous space generated by the local weighted reconstruction values ​​of the first grid points is faster and more efficient, and can meet the needs of the relevant oil and gas exploration, development and production cycle.

[0064] According to the local weighted reconstruction value of the first grid point, the first geological parameter model of the multiple seismic work areas can be updated, and then the second geological parameter model of the contiguous space can be generated based on the updated first geological parameter model of the multiple seismic work areas. Specifically, the local weighted reconstruction value of the first grid point can be used as the parameter value (such as the initial sampling value) of the first geological parameter model of the contiguous space at the corresponding first grid point, so as to update the initial sampling value of the first geological parameter model; and then the updated multiple first geological parameter models can be merged to obtain the second geological parameter model of the contiguous space. The above multiple first geological parameter models are the first geological parameter models of multiple seismic work areas in the contiguous space. The geological parameter model of the contiguous space generated by the local weighted reconstruction can be smoothed to eliminate the abnormal points of the local weighted reconstruction. The smoothing method can be a mean filtering method, a median filtering method, and a diffusion filtering method. For example, the geological parameter model of the contiguous space generated by the local weighted reconstruction can be median filtered to eliminate the abnormal points of the local weighted reconstruction. Specifically, a smoothing window containing an odd number of points can be set in the grid area of ​​the contiguous space, and this window can be scanned on the grid area of ​​the entire contiguous space. The parameter values ​​of the grid points contained in the window are arranged in ascending or descending order, and the parameter value in the middle is taken to replace the parameter value of the point, thereby achieving the elimination of abnormal local weighted reconstruction points.

[0065] A specific embodiment of this description is provided below:

[0066] 1. Obtain the initial velocity parameter model of the two seismic work areas in the contiguous space A, and then obtain the range, bin size, work area azimuth, value range of the parameter model and other information of the two seismic work areas. The range of the two seismic work areas is rectangular, and their spatial geographical location distribution is as follows Figure 2 As shown, they are in an oblique relationship, where the horizontal axis is the line number direction and the vertical axis is the point number direction. Figure 3This is a partial enlarged view of the overlapping position of the two seismic work areas. Analysis shows that the initial velocity model traces are unevenly distributed in space and the trace spacing is irregular. Among them, the bin size of seismic work area 1 is 20m×20m, the azimuth of the work area is 90 degrees (the north direction is 0 degrees, the angle between the main survey line and the horizontal axis), the range of the velocity model is [2000m / s, 5650m / s], and the total area of ​​the work area is 280km 2 , the line number range of the work area is 100-800, and the point number range is 100-1100. The area of ​​the seismic work area 2 is 25m×25m, the azimuth of the work area is 45 degrees, the range of the velocity model is [1680m / s, 7100m / s], and the total area of ​​the work area is 200km 2 , the work area line number range is 200-1000, and the point number range is 300-700.

[0067] 2. Select the grid points of seismic work area 1 from the A contiguous space and define it as the unified grid of the A contiguous space, because the seismic data quality of this seismic work area is the best, and its geographical location is relatively central among all seismic work areas in the contiguous range, and its range occupies the largest area. The grid of seismic work area 1 is a grid range composed of four-point coordinates, including 1) minimum line number, minimum point number, X and Y; 2) minimum line number, maximum point number, X and Y; 3) maximum line number, minimum point number, X and Y; 4) maximum line number, maximum point number, X and Y;

[0068] 3. Taking the value range of seismic area 1 as the standard, adjust the value range of the initial parameter model of seismic area 2 from the value range of [1680m / s, 7100m / s] to [2000m / s, 5650m / s], and finally make the value ranges of each seismic area consistent. The adjustment method is to obtain the maximum and minimum values ​​of the value range of the preferred seismic area 1, compare the value in the initial parameter model of each seismic area model with the aforementioned maximum and minimum values, if it is less than the aforementioned minimum value, then set the value in the initial parameter model to the aforementioned minimum value, if it is greater than the aforementioned maximum value, then set the value in the initial parameter model to the aforementioned maximum value.

[0069] 4. Determine that the spatial range of the initial parameter model of the contiguous earthquake work area is greater than or equal to the spatial range that includes all earthquake work areas, and regenerate the unified grid of the contiguous range to obtain the new maximum and minimum line numbers and point numbers of the contiguous range. After processing, the line number becomes 100-1230, and the point number becomes 100-1193. The maximum and minimum line numbers and point numbers do not become negative integers, so there is no need to adjust according to the specific numerical range of the line number and point number.

[0070] 5. According to the regenerated contiguous spatial range, that is, the grid range determined by the new maximum and minimum line numbers and point numbers, the grid position is determined one by one for each seismic work area, line number, point number, and sampling point, and the data range for local weighted reconstruction processing is preset, where the line number and point number are the dimensions in two spatial directions, and the sampling point is the dimension in the vertical direction. The local in local weighted reconstruction refers to the line number, point number, and sampling point of the current position to be processed as the center, and contains at least 5 initial parameter model data channels and 1 sampling point.

[0071] 6. Determine the X and Y coordinates of the position of the grid to be processed. Rotate the azimuth of the unified grid of the contiguous range to 0 degrees; make a difference between the line number and point number of the position to be determined and any line number and point number in the unified grid of the contiguous space range, and multiply the difference with the face element in the spatial direction of the line number and point number; add the product of the spatial direction of the line number and point number to the spatial coordinates corresponding to any line number and point number in the unified grid of the contiguous space range, and obtain the coordinates when the azimuth of the unified grid of the contiguous range is 0 degrees; rotate the coordinates in the opposite direction to obtain the spatial coordinates corresponding to the line number and point number of the grid position to be processed.

[0072] 7. Complete the local weighted reconstruction process of the grid position to be processed. The specific formula is as follows:

[0073]

[0074] Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X-coordinate and Y-coordinate of the grid position to be processed, respectively; R(x,y) is the local weighted reconstruction value of the grid position to be processed; x i and j are the X-coordinates and Y-coordinates of all positions in the local processing window except the grid position to be processed; ω i,j The spatial coordinates are (x i ,y i ) is the weight of all positions in the local processing window except the grid position to be processed; R i,j The spatial coordinates are (x i ,y i ) is the local weighted sampling value of all positions except the grid position to be processed in the local processing window.

[0075] 8. The initial parameter model after local weighted reconstruction is smoothed to eliminate interpolation anomalies. Specific smoothing methods include mean filtering method, median filtering method, diffusion filtering method, etc., thereby completing the establishment of a unified grid initial parameter model for the contiguous spatial seismic work area. Figure 4 This is a local magnified image of the spatial position after the unified grid regularization reconstruction of the contiguous spatial range. Figure 5The initial velocity model profile of the overlapping position of the two work areas after unified grid regularization reconstruction in the contiguous spatial range. It can be seen from the analysis that the technical solution of the embodiment of this specification solves the problem of reconstructing the initial velocity model of different seismic work areas in the contiguous spatial range, and the processed velocity model profile data is regularly distributed without abnormal interference, thus meeting the application requirements of establishing the initial parameter model for contiguous seismic processing.

[0076] Based on the above-mentioned method for constructing a geological parameter model, this specification also proposes an embodiment of a device for constructing a geological parameter model. Figure 6 As shown, the geological parameter model construction device 600 may specifically include the following modules:

[0077] The acquisition module 601 can be used to acquire the first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to the line numbers, point numbers and initial sampling values ​​of multiple first geological parameter models.

[0078] The setting module 602 may be configured to set a local processing window of the first grid point according to a line number and a point number corresponding to the first grid point.

[0079] The determination module 603 may be configured to use the grid points other than the first grid points in the local processing window as second grid points, and determine the weight of the second grid points according to the distance between the second grid points and the first grid points in the local processing window.

[0080] The calculation module 604 may be configured to calculate the local weighted reconstruction value of the first grid point according to the weight of the second grid point within the local processing window.

[0081] The generating module 605 may be used to generate a second geological parameter model of the contiguous space according to the local weighted reconstruction values ​​of the first grid points.

[0082] In some embodiments, the above-mentioned acquisition module 601 can be specifically used to obtain the first geological parameter models of multiple seismic work areas in a contiguous space; divide the multiple seismic work areas into a first seismic work area and multiple second seismic work areas according to the first geological parameter models of the multiple seismic work areas; adjust the fourth grid point of the second seismic work area according to the third grid point of the first seismic work area; determine the first grid point of the multiple seismic work areas in the contiguous space according to the third grid point of the first seismic work area and the fourth grid point of the multiple second seismic work areas.

[0083] In some embodiments, the above-mentioned acquisition module 601 can also be specifically used to adjust the line number and point number of the fourth grid point of the second seismic working area according to the line number and point number of the third grid point of the first seismic working area; according to the initial sampling value of the third grid point of the first seismic working area, adjust the initial sampling value of the fourth grid point of the second seismic working area.

[0084] In some embodiments, the acquisition module 601 may be further configured to generate line numbers and point numbers of a plurality of first grid points according to line numbers and point numbers of the third grid points and the fourth grid points.

[0085] In some embodiments, the determination module 603 can be specifically used to calculate the weight of the second grid point using a preset weight model according to the distance between the second grid point and the first grid point; the weights of all second grid points in the local processing window satisfy the formula Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; ω i,j The spatial coordinates are (x i ,y i ) is the weight of the second grid point.

[0086] In some embodiments, the calculation module 604 may be specifically configured to calculate the weight of the second grid point using the formula: Calculate the local weighted reconstruction value of the first grid point. Wherein, N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X coordinate and Y coordinate corresponding to the first grid point, respectively; R(x, y) is the local weighted reconstruction value of the first grid point; x i and j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,j The spatial coordinates are (x i ,y i )’s second grid point weight; R i,j The spatial coordinates are (x i ,y i ) is the local weighted sampling value of the second grid point.

[0087] In some embodiments, the above-mentioned calculation module 604 can also be specifically used to obtain the L second grid points closest to the second grid point in the local processing window as the fifth grid point; L is a preset positive integer; and the local weighted sampling value of the second grid point is determined based on the initial sampling values ​​of the L fifth grid points.

[0088] In some embodiments, the above-mentioned generation module 605 can be specifically used to update the first geological parameter models of the multiple seismic work areas according to the local weighted reconstruction values ​​of the first grid points; and generate the second geological parameter model of the contiguous space according to the updated first geological parameter models of the multiple seismic work areas.

[0089] It should be noted that the units, devices or modules described in the above embodiments can be implemented by computer chips or entities, or by products with certain functions. For the convenience of description, the above devices are described separately by functions divided into various modules. Of course, when implementing this specification, the functions of each module can be implemented in the same or more software and / or hardware, or the modules that implement the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0090] As can be seen from the above, based on the construction device of the geological parameter model provided in the embodiment of this specification, the first grid points of multiple seismic work areas in the contiguous space can be obtained; each seismic work area has a first geological parameter model, and one of the first grid points corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; according to the line number and point number corresponding to the first grid point, the local processing window of the first grid point is set; the grid points other than the first grid point in the local processing window are taken as the second grid points, and the weight of the second grid point is determined according to the distance between the second grid point and the first grid point within the local processing window; within the local processing window, the local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; according to the local weighted reconstruction value of the first grid point, the second geological parameter model of the contiguous space is generated. Compared with the existing method, the local weighted reconstruction of the parameters in the parameter model can be accurately performed based on the local processing area of ​​the grid points in the contiguous space, which is faster and more efficient, and can meet the relevant oil and gas exploration and development production cycle requirements.

[0091] The embodiment of the present specification also provides a computer device for a method of constructing a geological parameter model, including a processor and a memory for storing instructions executable by the processor. When the processor is implemented, the following steps can be performed according to the instructions: the first grid points of multiple seismic work areas in a contiguous space can be obtained; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; according to the line number and point number corresponding to the first grid point, a local processing window of the first grid point is set; the grid points other than the first grid point in the local processing window are taken as second grid points, and within the local processing window, the weight of the second grid point is determined according to the distance between the second grid point and the first grid point; within the local processing window, the local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; according to the local weighted reconstruction value of the first grid point, a second geological parameter model of the contiguous space is generated.

[0092] In order to complete the above instructions more accurately, refer to Figure 7 As shown, the embodiment of this specification also provides another specific computer device 700, wherein the computer device 700 includes a network communication port 701, a processor 702 and a memory 703, and the above structures are connected through internal cables so that each structure can perform specific data interaction.

[0093] The processor 702 can be specifically used to obtain the first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; according to the line number and point number corresponding to the first grid point, a local processing window of the first grid point is set; the grid points other than the first grid point in the local processing window are taken as second grid points, and within the local processing window, the weight of the second grid point is determined according to the distance between the second grid point and the first grid point; within the local processing window, the local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; and according to the local weighted reconstruction value of the first grid point, a second geological parameter model of the contiguous space is generated.

[0094] The memory 703 may be specifically used to store corresponding instruction programs.

[0095] In this embodiment, the network communication port 701 can be a virtual port that is bound to different communication protocols so that different data can be sent or received. For example, the network communication port can be a port responsible for web data communication, a port responsible for FTP data communication, or a port responsible for email data communication. In addition, the network communication port can also be a physical communication interface or communication chip. For example, it can be a wireless mobile network communication chip, such as GSM, CDMA, etc.; it can also be a Wifi chip; it can also be a Bluetooth chip.

[0096] In this embodiment, the processor 702 may be implemented in any appropriate manner. For example, the processor may take the form of a microprocessor or processor and a computer-readable medium storing a computer-readable program code (such as software or firmware) executable by the (micro)processor, a logic gate, a switch, an application specific integrated circuit (ASIC), a programmable logic controller, and an embedded microcontroller, etc. This specification does not limit this.

[0097] In this embodiment, the memory 703 includes a volatile memory and a non-volatile memory. The memory 703 may include multiple levels. In a digital system, anything that can store binary data can be a memory; in an integrated circuit, a circuit with a storage function without a physical form is also called a memory, such as RAM, FIFO, etc.; in a system, a storage device with a physical form is also called a memory, such as a memory stick, TF card, etc.

[0098] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. The present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0099] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0100] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0101] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0102] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for constructing a geological parameter model, characterized in that: include: Obtain the first grid point of multiple seismic work areas in a contiguous space; Each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to a plurality of line numbers, point numbers and initial sampling values ​​of the first geological parameter model; According to the line number and point number corresponding to the first grid point, a local processing window of the first grid point is set; Taking the grid points other than the first grid point in the local processing window as second grid points, and determining the weight of the second grid point in the local processing window according to the distance between the second grid point and the first grid point; In the local processing window, calculating a local weighted reconstruction value of the first grid point according to the weight of the second grid point; A second geological parameter model of the contiguous space is generated according to the local weighted reconstruction values ​​of the first grid points.

2. The method according to claim 1, characterized in that: The step of obtaining the first grid points of a plurality of seismic work areas in the contiguous space includes: Obtaining the first geological parameter model of multiple seismic work areas in a contiguous space; According to the first geological parameter models of the multiple seismic working areas, the multiple seismic working areas are divided into a first seismic working area and multiple second seismic working areas; According to the third grid point of the first seismic working area, adjusting the fourth grid point of the second seismic working area; According to the third grid point of the first seismic working area and the fourth grid points of the plurality of second seismic working areas, the first grid point of the plurality of seismic working areas in the contiguous space is determined.

3. The method according to claim 2, characterized in that: The adjusting the fourth grid point of the second seismic working area according to the third grid point of the first seismic working area comprises: According to the line number and point number of the third grid point of the first seismic working area, adjust the line number and point number of the fourth grid point of the second seismic working area; According to the initial sampling value of the third grid point of the first seismic working area, the initial sampling value of the fourth grid point of the second seismic working area is adjusted.

4. The method according to claim 2, characterized in that: The step of determining the first grid point of the plurality of seismic working areas in the contiguous space according to the third grid point of the first seismic working area and the fourth grid points of the plurality of second seismic working areas comprises: The line numbers and point numbers of the plurality of first grid points are generated according to the line numbers and point numbers of the third grid points and the fourth grid points.

5. The method according to claim 1, characterized in that: The step of determining the weight of the second grid point according to the distance between the second grid point and the first grid point includes: Calculate the weight of the second grid point using a preset weight model according to the distance between the second grid point and the first grid point; The weights of all second grid points in the local processing window satisfy the following formula: Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; ω i,j The spatial coordinates are (x i ,y i ) is the weight of the second grid point.

6. The method according to claim 1, characterized in that: The calculating the local weighted reconstruction value of the first grid point according to the weight of the second grid point includes: According to the local weighted sampling value of the second grid point and the weight of the second grid point, the local weighted reconstruction value of the first grid point is calculated using the following formula: Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X-coordinate and Y-coordinate corresponding to the first grid point, respectively; R(x, y) is the local weighted reconstruction value of the first grid point; x i and j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,j The spatial coordinates are (x i ,y i )’s second grid point weight; R i,j The spatial coordinates are (x i ,y i ) is the local weighted sampling value of the second grid point.

7. The method according to claim 6, characterized in that: The method further comprises: Acquire L second grid points closest to the second grid point in the local processing window as the fifth grid point; L is a preset positive integer; The local weighted sampling value of the second grid point is determined according to the initial sampling values ​​of the L fifth grid points.

8. The method according to claim 1, characterized in that: The step of generating a second geological parameter model of the contiguous space according to the local weighted reconstruction values ​​of the first grid points comprises: updating the first geological parameter models of the plurality of seismic work areas according to the local weighted reconstruction values ​​of the first grid points; A second geological parameter model of the contiguous space is generated based on the updated first geological parameter models of the plurality of seismic work areas.

9. A device for constructing a geological parameter model, characterized in that: The device comprises: An acquisition module is used to acquire first grid points of multiple seismic work areas in a contiguous space; each of the seismic work areas has a first geological parameter model, and one of the first grid points corresponds to line numbers, point numbers and initial sampling values ​​of multiple first geological parameter models; A setting module, configured to set a local processing window of the first grid point according to a line number and a point number corresponding to the first grid point; A determination module, configured to use the grid points other than the first grid points in the local processing window as second grid points, and determine the weight of the second grid points according to the distance between the second grid points and the first grid points in the local processing window; A calculation module, configured to calculate, within the local processing window, a local weighted reconstruction value of the first grid point according to the weight of the second grid point; A generating module is used to generate a second geological parameter model of the contiguous space according to the local weighted reconstruction values ​​of the first grid points.

10. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 8 is implemented.

Images