A method, device and equipment for constructing a geological parameter model

Abstract

The present specification relates to the technical field of seismic exploration data processing, and particularly relates to a geological parameter model construction method, device and equipment. The geological parameter model construction method comprises: acquiring a plurality of first grid points of a plurality of seismic work areas in a continuous space; each seismic work area has a first geological parameter model, a line number, a point number and an initial sampling value of the first geological parameter model corresponding to each first grid point; 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; taking a grid point other than the first grid point in the local processing window as a second grid point, and determining a weight of the second grid point according to a 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 in the local processing window; and generating a second geological parameter model of the continuous 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 advances in seismic acquisition technology, key oil and gas exploration areas, particularly those with mature exploration potential, have undergone multiple rounds of repeated 3D seismic acquisition, evolving from narrow-azimuth acquisition to wide-azimuth acquisition, and finally to "two-wide, one-high" acquisition. High-quality seismic processing and interpretation, target optimization, and well site support have been conducted in each seismic block, laying a solid foundation for high-quality seismic data. With the continued advancement and intensification of oil and gas exploration and development in recent years, the need for regional geological research and large-scale oil and gas discoveries has become increasingly urgent. Domestic oil and gas field exploration and development is gradually shifting towards large-scale 3D seismic exploration. Due to the increasing cost of seismic acquisition, obtaining large-scale 3D seismic data primarily involves processing seismic blocks from different acquisition years and observation methods, which is currently the most cost-effective solution.

[0003] To efficiently complete contiguous processing, it is necessary to fully utilize the existing seismic processing results of individual seismic work areas. However, given that different seismic work areas within a contiguous area often overlap, and the overlap is even greater for secondary and higher-level 3D seismic acquisition, and that observation methods vary significantly, including complexities such as different acquisition bins, inconsistent azimuths, and irregular observation ranges, conventional methods for constructing initial parameter models for contiguous spatial seismic work areas obtain initial parameter models for velocity, density, and anisotropy by conducting 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 requirements of the oil and gas exploration and development production cycle.

[0005] Therefore, how to overcome the problems of high consumption of manpower, material and financial resources and low efficiency in the existing methods for constructing initial parameter models of contiguous spatial seismic work areas, and propose a method for constructing contiguous spatial parameter models 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 of constructing initial parameter models for contiguous spatial seismic work areas and improve the efficiency of application, thereby meeting the production cycle requirements of oil and gas exploration and development.

[0007] On the one hand, an embodiment of this specification proposes a method for constructing a geological parameter model, which includes: obtaining first grid points of multiple seismic work areas in a continuous space; each seismic work area has a first geological parameter model, and one first grid point 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 used 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 continuous space is generated.

[0008] On the other hand, a device for constructing a geological parameter model includes: an acquisition module for acquiring first grid points of multiple seismic work areas in a continuous space; each seismic work area has a first geological parameter model, and one first grid point corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; a setting module for setting 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 for 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 within the local processing window according to the distance between the second grid point and the first grid point; a calculation module for calculating 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 for generating a second geological parameter model of the continuous 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] As can be seen from the technical solutions provided in the embodiments of this specification, the method for constructing a geological parameter model provided in the embodiments of this specification can obtain first grid points of multiple seismic work areas in a contiguous space; each seismic work area has a first geological parameter model, and one first grid point corresponds to the line number, point number, and initial sampling value of multiple first geological parameter models; a local processing window of the first grid point is set according to the line number and point number corresponding to the first grid point; grid points other than the first grid point in the local processing window are used 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, a local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; and a second geological parameter model of the contiguous space is generated based on the local weighted reconstruction value of the first grid point. Compared with existing methods, the method can accurately and locally weightedly reconstruct parameters in the parameter model based on the local processing area of ​​the grid points in the contiguous space, which is faster and more efficient, and can meet the requirements of the relevant oil and gas exploration and development production cycle. BRIEF DESCRIPTION OF THE DRAWINGS

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

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

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

[0014] Figure 3 This is a partially enlarged schematic diagram of the overlapping positions of two seismic work areas within the contiguous space A 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 cross section of the initial velocity model of the reconstructed A contiguous space provided in an embodiment of this specification;

[0017] Figure 6 This 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 the computer device provided in the 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 embodiments described 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 making creative efforts should fall within the scope of protection of this specification.

[0020] Figure 1 This 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 first grid point 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 for multiple seismic areas within a contiguous space, the areas can be divided into a first seismic area and multiple second seismic areas. The first geological parameter models for the second seismic areas can then be adjusted based on the first geological parameter models for the first seismic area. Compared to the second seismic areas, the first seismic area has higher-quality seismic data and a better geographical location, which can improve the speed and efficiency of the reconstruction process.

[0024] A contiguous space can refer to a single exploration area formed by connecting multiple adjacent or nearby seismic work areas during seismic exploration. The distance between the centers of each seismic work area in the contiguous space is less than or equal to a preset distance threshold. First geological parameter models for multiple seismic work areas within the contiguous space can be obtained. These first geological parameter models for the seismic work areas can include velocity models, anisotropy models, density models, dip models, and corresponding parameter updates for each model in the time or depth domain. Each seismic work area within the contiguous space is divided into grids, with each grid uniquely corresponding to the line number, point number, and initial sampling values ​​of multiple geological parameter models. For example, for a velocity model, each grid uniquely corresponds to the line number, point number, and initial velocity value of the velocity model. One seismic work area can be selected from the multiple seismic work areas within the contiguous space, and the grid of the selected seismic work area can be preset as a unified grid, serving as the grid division basis for the contiguous space. A first seismic work area can be selected based on the seismic data quality and geographic location of each seismic work area within the contiguous space, with the remaining seismic work areas within the contiguous space serving as second seismic work areas.

[0025] Seismic data quality can be assessed for each seismic site, including frequency analysis, signal-to-noise ratio analysis, coverage analysis, and interference wave analysis. Frequency analysis involves performing a frequency sweep analysis of seismic data from each seismic site to determine the effective bandwidth and bandwidth. The clarity and continuity of reflection information within different frequency bands can be assessed, prioritizing seismic sites with clear and continuous reflection information. Signal-to-noise ratio analysis calculates the signal-to-noise ratio for each seismic site and assesses the continuity of the primary target layer. Seismic sites with high signal-to-noise ratios and good continuity of the primary target layer can be prioritized. Coverage analysis analyzes the coverage of data from each seismic site to ensure uniform coverage and meet analysis requirements. Seismic sites with uniform and high coverage can be prioritized. Interference wave analysis identifies the types of interference waves in each seismic site, such as surface waves, shallow refraction waves, and multiple waves. The impact of interference waves on seismic data quality can be assessed, prioritizing seismic sites with minimal interference waves.

[0026] The geographic location of each seismic work area can be evaluated, including centrality analysis and transportation and infrastructure analysis. Centrality analysis uses GIS technology to determine the geographic location of each seismic work area. The distance from each seismic work area to the center of the contiguous area can be calculated, with seismic work areas closer to the center being prioritized. Transportation and infrastructure analysis considers accessibility and infrastructure conditions within the seismic work area, such as roads, communications, and electricity. Seismic work areas with convenient transportation and comprehensive infrastructure can be prioritized 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 geographic location assessments. This index system may include effective bandwidth, signal-to-noise ratio, coverage times, interference wave impact, geographic location centrality, transportation convenience, and infrastructure conditions. This comprehensive evaluation index system can be used to conduct a comprehensive assessment of each seismic work area. Based on the assessment results, the seismic work areas are ranked, and the seismic work area with the best comprehensive evaluation results is selected as the first seismic work area. Based on the seismic data grid of the selected first seismic work area, parameters such as the direction, vertices, coordinate data, vertical and horizontal grid step sizes, point numbers, and line numbers of the contiguous space unified grid are determined. The unified grid consists of a grid range composed of four-point coordinates. 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; and 4) maximum line number, maximum point number, X and Y. The maximum and minimum line and point numbers are adjustable according to a numerical range, and their values ​​are typically adjusted to positive integers greater than 0.

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

[0029] By dividing the three-dimensional seismic work area using a grid method, multiple grid points in 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 using a grid method. The fourth grid point is a grid point obtained by dividing the second seismic work area using a grid method.

[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 preliminary 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] Based on the selected first seismic work area, the line number, point number, and initial sampling value of the fourth grid point in the second seismic work area can be adjusted based on the line number, point number, and initial sampling value of the third grid point in the first seismic work area. The third grid point in the first seismic work area and the adjusted fourth grid points in the multiple second seismic work areas can be considered to be based on the same grid reference. Therefore, based on the third grid point in the first seismic work area and the adjusted fourth grid points in the multiple second seismic work areas, the first grid point of multiple 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 working area may be adjusted according to the line number and point number of the third grid point of the first seismic working 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] Based on the line number and point number of the third grid point in the first seismic work area, the four-point coordinates of the first seismic work area grid area can be determined. Based on the four-point coordinates of the first seismic work area grid area, the line number and point number of the fourth grid point in 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 number and point number become negative, the line number and point number need to be adjusted to a positive number 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] Based on the value range of the initial sampling values ​​of the preferred first seismic working area, the value range of the initial sampling values ​​in the initial parameter models of the multiple second seismic working areas can be adjusted, and finally the value range of all seismic working areas in the contiguous space is made consistent. The adjustment method specifically comprises obtaining the maximum and minimum values ​​of the value range of the initial sampling values ​​of the preferred first seismic working area, comparing the initial sampling values ​​in the initial parameter models of the multiple second seismic working area models with the aforementioned maximum and minimum values, and if the initial sampling values ​​are less than the aforementioned minimum value, setting the initial sampling values ​​in the initial parameter models of the second seismic working area to the aforementioned minimum value; if the initial sampling values ​​are greater than the aforementioned maximum value, setting the initial sampling values ​​in the initial parameter models of the second seismic working area 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 the 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 the further reconstruction of grid data in the contiguous space.

[0040] Based on the third grid point of the first seismic work area and the adjusted fourth grid points of multiple second seismic work areas, the first grid points of multiple seismic work areas in the contiguous space can be determined. The third grid point and the adjusted fourth grid point are both determined based on the same standard. Therefore, based on 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 points of multiple second seismic work areas, the line number, point number, and initial sampling value of all 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 multiple second seismic work areas can be directly used as the corresponding first grid points. Therefore, the first grid points of 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 for the first grid point according to the line number and 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, laying the 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 the dimensions of the two spatial directions, and the initial sampling value is the dimension of the vertical direction. Based on a set of line numbers, point numbers and initial sampling values, a grid point and its location can be uniquely determined. 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 and B can be set to 4. 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 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.

[0045] In some embodiments, the weight of the second grid point can be calculated using a preset weight model based on 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 coordinate is (x i ,y i ) is the weight of the second grid point.

[0046] By using the distance between the second grid point and the first grid point to calculate the weight of the second grid point, a data foundation is laid for achieving accurate reconstruction of the first grid point.

[0047] The azimuths 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 differences can be multiplied by the bins of the line number spatial direction and the point number spatial direction of the first grid point, and the products of the line number spatial direction and the point number spatial 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 obtained coordinates of the first grid point 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 described 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. 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 coordinate is (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 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 within 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. The Sigmod function can output the initial weight coefficient for each second grid point. The initial weight coefficients of each second grid point within the local processing window can be added together to obtain an initial weight sum, and then the initial weight coefficients of each second grid point within the local processing window can be divided by the initial weight sum to obtain the weight of each second grid point.

[0052] S104: 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.

[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. Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X and Y coordinates corresponding to the first grid point, respectively; R(x,y) is the local weighted reconstruction value of the first grid point; x i and y j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,j The spatial coordinate is (x i ,y i )’s second grid point weight; R i,j The spatial coordinate is (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 areas in the contiguous space and accurate encryption of data sparse areas in the contiguous space can be achieved.

[0055] According to the local weighted sampling values ​​and weights of all 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 of the first grid point, respectively; R(x,y) is the local weighted reconstruction value of the first grid point; x i and y j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,jThe spatial coordinate is (x i ,y i )’s second grid point weight; R i,j The spatial coordinate is (x i ,y i ). The local weighted reconstructed value of the first grid point is calculated using the local weighted reconstructed value of the second grid point within the local processing window and the corresponding weight. The local weighted reconstructed value calculated in this way can fully utilize the first geological parameter model data and generate an accurate and reasonable local weighted reconstructed value 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. The L closest second grid points are selected as the fifth grid points corresponding to the current grid point. An activation function such as a Sigmoid function, a ReLU function, or a LeakyReLU function can be used as a 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, which can output the initial weight coefficient for each fifth grid point. The initial weight coefficients of each fifth grid point in the local processing window can be added together to obtain an 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. Based on the initial sampling values ​​and weights corresponding to each fifth grid point, a weighted sum of all fifth grid points can be calculated, and the calculated weighted sum can be used as the local weighted sampling value for 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 based on the local weighted reconstruction values ​​of the first grid points; and the second geological parameter model of the contiguous space can be generated based on 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 requirements of the relevant oil and gas exploration and development production cycle.

[0064] Based on the locally weighted reconstructed values ​​of the first grid points, the first geological parameter models of the multiple seismic work areas can be updated, and then a second geological parameter model of the contiguous space can be generated based on the updated first geological parameter models of the multiple seismic work areas. Specifically, the locally weighted reconstructed values ​​of the first grid points can be used as parameter values ​​(e.g., initial sampling values) of the first geological parameter model of the contiguous space at the corresponding first grid points, thereby updating the initial sampling values ​​of the first geological parameter model. The multiple updated first geological parameter models can then be merged to obtain the second geological parameter model of the contiguous space. The multiple first geological parameter models described above are the first geological parameter models of the multiple seismic work areas within the contiguous space. The contiguous space geological parameter model generated by the local weighted reconstruction can be smoothed to eliminate outliers in the local weighted reconstruction. Smoothing methods can include mean filtering, median filtering, and diffusion filtering. For example, the contiguous space geological parameter model generated by the local weighted reconstruction can be median filtered to eliminate outliers in 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 over 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, and value range of the parameter model 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 coordinate is the line number direction and the vertical coordinate is the point number direction. Figure 3This is a partial zoomed-in image of the overlapping location of the two seismic work areas. Analysis shows that the initial velocity model traces are unevenly distributed and the trace spacing is irregular. The bin size of seismic work area 1 is 20m×20m, the work area azimuth is 90 degrees (with due north as 0 degrees, the angle between the main survey line and the horizontal axis), the velocity model range is [2000m / s, 5650m / s], and the total work area is 280km. 2 The work area line number range is 100-800, and the point number range is 100-1100. The cell size of seismic work area 2 is 25m×25m, the work area azimuth is 45 degrees, the velocity model range is [1680m / s, 7100m / s], and the total 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 area 1 from contiguous space A and define it as the unified grid of contiguous space A. This is because the seismic data quality of this seismic area is the best, its geographical location is relatively central among all seismic areas in the contiguous area, and its range occupies the largest area. The grid of seismic area 1 consists of a grid range composed of four coordinates: 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. Using the range of earthquake zone 1 as a standard, adjust the range of the initial parameter model of earthquake zone 2 from [1680 m / s, 7100 m / s] to [2000 m / s, 5650 m / s], ultimately making the ranges of all earthquake zones consistent. The adjustment method is to obtain the maximum and minimum values ​​of the optimal range of earthquake zone 1, compare the value in the initial parameter model of each earthquake zone model with the maximum and minimum values, and if the value is less than the minimum value, set the value in the initial parameter model to the minimum value; if the value is greater than the maximum value, set the value in the initial parameter model to the maximum value.

[0069] 4. Determine the spatial extent of the initial parameter model for the contiguous earthquake zone to be greater than or equal to the spatial extent encompassing all earthquake zones. Regenerate a unified grid for the contiguous zone, and then obtain new maximum and minimum line and point numbers for the contiguous zone. This processing changes the line numbers to 100-1230, and the point numbers to 100-1193. The maximum and minimum line and point numbers do not become negative integers, so no adjustments are required based on their specific numerical ranges.

[0070] 5. Based on the regenerated contiguous spatial range (i.e., the grid range defined by the new maximum and minimum line and point numbers), the grid locations are determined for each seismic work area, line number, point number, and sampling point. The data range for local weighted reconstruction processing is also preset, where the line and point numbers represent the dimensions in two spatial directions, and the sampling point represents the vertical dimension. The term "local" in local weighted reconstruction refers to a region centered around the line number, point number, and sampling point of the current location to be processed, and includes at least five initial parameter model data channels and one sampling point.

[0071] 6. Determine the X and Y coordinates of the grid position to be processed. Rotate the azimuth of the unified grid within the contiguous range to 0 degrees; subtract the line number and point number of the position to be processed from any line number and point number in the unified grid within the contiguous spatial range, and multiply the difference by the bin in the spatial direction of the line number and point number; add the product of the spatial directions of the line number and point number to the spatial coordinates corresponding to the line number and point number in the unified grid within the contiguous spatial range to obtain the coordinates when the azimuth of the unified grid within 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 y j are the X-coordinates and Y-coordinates of all positions in the local processing window except the position of the grid to be processed; ω i,j The spatial coordinate is (x i ,y i ) in the local processing window except the grid position to be processed; R i,j The spatial coordinate is (x i ,y i ) is the local weighted sampling value of all positions except the grid position to be processed within the local processing window of .

[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 processing of the continuous spatial range. Figure 5The initial velocity model profile at the overlapping location of two seismic work areas after unified grid regularization and reconstruction within a contiguous spatial range is shown. Analysis shows that the technical solutions of the embodiments of this specification solve the problem of reconstructing the initial velocity model for different seismic work areas within a contiguous spatial range. The processed velocity model profile data traces are regularly distributed and free of abnormal interference, thus meeting the application requirements for establishing initial parameter models 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 obtain the first grid points of multiple seismic work areas in a continuous space; each seismic work area has a first geological parameter model, and one first grid point corresponds to the line number, point number and initial sampling value of multiple first geological parameter models.

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

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

[0080] The calculation module 604 may be 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.

[0081] The generating module 605 may be configured 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 continuous 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 continuous 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 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; and adjust the initial sampling value 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.

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

[0085] In some embodiments, the determination module 603 may be specifically configured 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 within 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 coordinate is (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. Where N and M are half of the horizontal and vertical sizes of the local processing window, respectively; x and y are the X and Y coordinates corresponding to the first grid point, respectively; R(x,y) is the local weighted reconstruction value of the first grid point; x i and y j are the X-coordinate and Y-coordinate of the second grid point in the local processing window respectively; ω i,j The spatial coordinate is (x i ,y i )’s second grid point weight; R i,j The spatial coordinate is (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; the 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 based on the local weighted reconstruction values ​​of the first grid points; and generate the second geological parameter model of the contiguous space based on 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 in terms of functions and are divided into various modules and described separately. Of course, when implementing this specification, the functions of each module can be implemented in the same or multiple software and / or hardware, or the module that implements the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods, 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, the geological parameter model construction device provided in the embodiments of this specification can obtain first grid points of multiple seismic work areas within a contiguous space; each seismic work area has a first geological parameter model, and each first grid point corresponds to the line number, point number, and initial sampling value of multiple first geological parameter models; a local processing window for the first grid point is set based on the line number and point number corresponding to the first grid point; grid points other than the first grid point within the local processing window are used as second grid points, and within the local processing window, the weight of the second grid point is determined based on the distance between the second grid point and the first grid point; within the local processing window, a local weighted reconstruction value of the first grid point is calculated based on the weight of the second grid point; and a second geological parameter model of the contiguous space is generated based on the local weighted reconstruction value of the first grid point. Compared with existing methods, this method can accurately and locally weightedly reconstruct parameters in the parameter model based on the local processing area of ​​the grid points within the contiguous space, which is faster and more efficient, and can meet the requirements of the relevant oil and gas exploration and development production cycle.

[0091] The embodiments of this specification also provide a computer device for a method of constructing a geological parameter model, including a processor and a memory for storing processor-executable instructions. When the processor is specifically implemented, the following steps can be performed according to the instructions: first grid points of multiple seismic work areas in a contiguous space can be obtained; each seismic work area has a first geological parameter model, and one first grid point corresponds to the line number, point number and initial sampling value of multiple first geological parameter models; a local processing window of the first grid point is set according to the line number and point number corresponding to the first grid point; grid points other than the first grid point in the local processing window are used 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, a local weighted reconstruction value of the first grid point is calculated according to the weight of the second grid point; and a second geological parameter model of the contiguous space is generated according to the local weighted reconstruction value of the first grid point.

[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 continuous 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 used 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 continuous 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, thereby being capable of sending or receiving different data. 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. Furthermore, 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 or CDMA; it can also be a Wi-Fi chip; or it can be a Bluetooth chip.

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

[0097] In this embodiment, the memory 703 includes volatile memory and non-volatile memory. The memory 703 can include multiple levels. In digital systems, anything that can store binary data can be considered a memory. In integrated circuits, a circuit with a storage function that does not have a physical form is also called a memory, such as RAM and FIFO. In systems, a physical storage device is also called a memory, such as a memory stick or TF card.

[0098] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. The present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. 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 magnetic 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 flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, 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 flowcharts and / or block diagrams. 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 that can direct a computer or other programmable data processing device to work 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 The 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 operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function 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: Obtaining first geological parameter models of multiple seismic work areas in a contiguous space; dividing the plurality of seismic working areas into a first seismic working area and a plurality of second seismic working areas according to the first geological parameter models of the plurality of seismic working areas; Adjusting the line number and point number of the fourth grid point in the second seismic working area according to the line number and point number of the third grid point in the first seismic working area; Adjusting the initial sampling value 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; Generate line numbers and point numbers of a plurality of first grid points according to the line numbers and point numbers of the third grid point and the fourth grid point; Setting a local processing window for the first grid point according to the line number and point number corresponding to the first grid point; Taking the grid points other than the first grid point within the local processing window as second grid points, and determining 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; Calculating a local weighted reconstruction value of the first grid point according to the weight of the second grid point within the local processing window; updating first geological parameter models of the plurality of seismic working 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.

2. The method according to claim 1, characterized in that The determining the weight of the second grid point according to the distance between the second grid point and the first grid point includes: Calculating 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, and are half of the horizontal and vertical sizes of the local processing window respectively; The spatial coordinates are The weight of the second grid point.

3. 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, and are half of the horizontal and vertical sizes of the local processing window respectively; and are the X coordinate and Y coordinate corresponding to the first grid point respectively; is the local weighted reconstruction value of the first grid point; and are the X-coordinate and Y-coordinate corresponding to the second grid point in the local processing window respectively; The spatial coordinates are The weight of the second grid point; The spatial coordinates are The local weighted sampling value of the second grid point.

4. The method according to claim 3, characterized in that The method further comprises: Get the nearest grid point in the local processing window to the second grid point The second grid point is used as the fifth grid point; is a preset positive integer; According to the The initial sampling value of the fifth grid point is used to determine the local weighted sampling value of the second grid point.

5. A device for constructing a geological parameter model, characterized in that: The device comprises: An acquisition module is configured to acquire first geological parameter models of multiple seismic working areas within a contiguous space; divide the multiple seismic working areas into a first seismic working area and multiple second seismic working areas based on the first geological parameter models of the multiple seismic working areas; adjust the line number and point number of the fourth grid point of the second seismic working area based on the line number and point number 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 based on the initial sampling value of the third grid point of the first seismic working area; and generate the line number and point number of the multiple first grid points based on the line number and point number of the third grid point and the fourth grid point; a setting module, configured to set a local processing window for 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 a grid point other than the first grid point within the local processing window as a second grid point, and determine a weight of the second grid point within the local processing window according to a distance between the second grid point and the first grid point; 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 generation module is 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 a second geological parameter model of the contiguous space according to the updated first geological parameter models of the multiple seismic work areas.

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

Images