Method, device, electronic device and medium for obtaining sandstone thickness map

Abstract

The present invention discloses a method, device, electronic device, and medium for obtaining a sandstone thickness map. The method comprises: performing well logging sensitivity parameter analysis on corresponding wells in a study area to determine logging parameters for lithologic inversion, performing well-seismic combined lithologic inversion based on the logging parameters to obtain an inversion data volume; performing well logging lithologic interpretation on corresponding wells in the study area to determine a threshold value for the sandstone logging parameters; obtaining sandstone threshold plane data based on the threshold value using an inverse distance interpolation method, and obtaining a sandstone threshold data volume based on the vertical sampling interval of the seismic data and the threshold plane data; determining a residual value data volume based on the inversion data volume and the sandstone threshold data volume, and obtaining a sandstone thickness map by counting the number of sandstone data points in the residual value data volume. This method can improve the accuracy of sandstone thickness prediction and sandstone thickness map drawing by performing well-seismic combined lithologic inversion and determining the threshold value for the wells in the study area.

Description

Technical Field

[0001] The present invention relates to the field of petroleum geological exploration, and in particular to a method, device, electronic equipment and medium for acquiring a sandstone thickness map. Background Art

[0002] In the field of petroleum exploration, accurately determining sandstone thickness and compiling detailed and reliable sandstone thickness maps are of great significance for the study of sand body distribution and sedimentary facies in oil and gas exploration. In particular, it is crucial to achieve reservoir evaluation before drilling of the trapped target, which directly affects the evaluation results of the drilling target oil layer thickness and trapped resource volume.

[0003] At present, the field of oil exploration mainly uses the geophysical inversion method combining drilling and seismic to obtain sandstone thickness and compile sandstone thickness maps. Due to the different lithologic combinations and fluid properties revealed by drilling in different regions, there are differences in logging responses. For the convenience of calculation, the logging curves of all wells in the study area are usually standardized in the geological modeling before inversion, so that each well has a unified lithologic baseline and sandstone threshold value, and then a joint well-seismic inversion is performed; at the same time, the unified sandstone threshold value is used to "filter" the inversion results to obtain the sandstone thickness of the entire area and compile a sandstone thickness map.

[0004] Although this method is easy and simple to use a constant threshold value to calculate the sandstone thickness of the study area, it artificially processes the logging curves of each well into a uniform lithologic baseline and sandstone threshold value, changing the true rock electrical parameter values ​​at the well points. For areas with flat terrain and small lithologic changes, this method can approximately calculate the thickness of the formation sandstone. However, for areas with large structural drop and large lithologic changes, due to the large differences in rock electrical parameters between well points, the artificial standardization of the logging curves weakens the differences in rock electrical characteristics between wells and changes the rock electrical parameter values ​​at the well points. It is difficult to accurately calculate the true sandstone thickness using the standardized logging curves for geological modeling and well-seismic joint lithologic inversion, resulting in large errors in the sandstone thickness calculated based on the inversion results and unreliable sandstone thickness maps. Summary of the Invention

[0005] The present invention provides a method, device, electronic equipment and medium for obtaining a sandstone thickness map, so as to solve the problems of large errors caused by predicting sandstone thickness based on standardized logging curves and lithology inversion and unreliable compiled sandstone thickness maps.

[0006] According to one aspect of the present invention, a method for obtaining a sandstone thickness map is provided, the method comprising:

[0007] Conduct logging sensitivity parameter analysis on the corresponding wells in the study area to determine the logging parameters used for lithologic inversion. Perform well-seismic combined lithologic inversion based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizons to obtain the inversion data volume.

[0008] Performing well logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells, wherein the sandstone threshold value is a well logging parameter value used to determine the lithologic properties of the wells;

[0009] Based on the inverse distance interpolation method and according to the sandstone threshold value, threshold plane data is obtained, and a sandstone threshold data volume is obtained according to the vertical sampling interval of the seismic data and the threshold plane data, wherein the threshold plane data is plane grid data presented in the form of coordinates obtained by the inverse distance interpolation method from the sandstone threshold values ​​of the corresponding wells in the study area;

[0010] A residual value data volume is determined according to the inversion data volume and the sandstone threshold value data volume, and a sandstone thickness map is obtained by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0011] According to another aspect of the present invention, a device for obtaining a sandstone thickness map is provided, the device comprising:

[0012] An inversion data volume acquisition module is used to analyze well logging sensitivity parameters for the corresponding wells in the study area, determine the well logging parameters used for lithologic inversion, and perform well-seismic combined lithologic inversion based on the well logging parameters, the seismic data volume of the study area, and the seismic interpretation horizons to obtain the inversion data volume.

[0013] A threshold value determination module is used to perform logging lithologic interpretation on corresponding wells in the study area to determine a threshold value of the wells, wherein the sandstone threshold value is a logging parameter value used to determine the lithology;

[0014] a threshold value data volume acquisition module, configured to acquire threshold value plane data based on the sandstone threshold value using an inverse distance interpolation method, and to acquire the sandstone threshold value data volume based on the vertical sampling interval of the seismic data and the threshold value plane data, wherein the threshold value plane interpolation data is plane grid data presented in the form of coordinates obtained by using the inverse distance interpolation method from the sandstone threshold values ​​of corresponding wells in the study area;

[0015] The sandstone thickness map acquisition module is used to determine a residual value data volume according to the inversion data volume and the sandstone threshold value data volume, and obtain a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0016] According to another aspect of the present invention, an electronic device is provided, comprising:

[0017] at least one processor; and

[0018] a memory communicatively connected to the at least one processor; wherein,

[0019] The memory stores a computer program executable by the at least one processor. The computer program is executed by the at least one processor so that the at least one processor can perform the method for obtaining a sandstone thickness map according to any embodiment of the present invention.

[0020] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for obtaining a sandstone thickness map according to any embodiment of the present invention when executed.

[0021] The technical solution of the embodiment of the present invention determines logging parameters for lithologic inversion by analyzing well logging sensitivity parameters for corresponding wells in a study area. Joint well-seismic lithologic inversion is performed based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizons to obtain an inversion data volume. Lithologic interpretation is performed on corresponding wells in the study area to determine the sandstone threshold value for the wells. Threshold plane data is obtained based on the sandstone threshold value using the inverse distance interpolation method, and a sandstone threshold data volume is obtained based on the vertical sampling interval of the seismic data and the threshold plane data. A residual data volume is determined based on the inversion data volume and the sandstone threshold data volume, and a sandstone thickness map is obtained by counting the number of data points in the residual data volume that vertically reflect sandstone data. This method performs joint well-seismic lithologic inversion using logging parameters and calculates sandstone thickness using the sandstone threshold data volume, thereby improving the accuracy of sandstone thickness prediction and sandstone thickness map drawing.

[0022] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become readily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0024] Figure 1 is a flow chart of a method for obtaining a sandstone thickness map according to an embodiment of the present invention;

[0025] Figure 2is a threshold value and a threshold value plane map provided according to an embodiment of the present invention;

[0026] Figure 3 2. It is a schematic diagram of a threshold value data body provided according to an embodiment of the present invention;

[0027] Figure 4 is a flow chart of a method for obtaining a sandstone thickness map according to an embodiment of the present invention;

[0028] Figure 5 2 is a schematic diagram of a residual value data volume provided according to an embodiment of the present invention;

[0029] Figure 6 Schematic diagram of lithologic differentiation based on residual value data provided by an embodiment of the present invention;

[0030] Figure 7 is a plan view of the number of sandstone sampling points provided in an embodiment of the present invention;

[0031] Figure 8 is a sandstone thickness map in the depth domain of sandstones of corresponding wells in the study area provided by an embodiment of the present invention;

[0032] Figure 9 is a flow chart of a method for obtaining a sandstone thickness map according to an embodiment of the present invention;

[0033] Figure 10 2 is a schematic structural diagram of a device for obtaining a sandstone thickness map according to a fourth embodiment of the present invention;

[0034] Figure 11 It is a structural schematic diagram of an electronic device for implementing the method for obtaining a sandstone thickness map according to an embodiment of the present invention. DETAILED DESCRIPTION

[0035] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0036] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0037] Example 1

[0038] Figure 1 A flow chart of a method for obtaining a sandstone thickness map is provided for an embodiment of the present invention. This embodiment is applicable to the case of obtaining a sandstone thickness map. The method can be executed by a device for obtaining a sandstone thickness map. The device for obtaining a sandstone thickness map can be implemented in the form of hardware and / or software. The device for obtaining a sandstone thickness map can be configured in any electronic device with network communication function. Figure 1 As shown, the method may include:

[0039] S110, performing logging sensitive parameter analysis on corresponding wells in the study area, determining logging parameters for lithologic inversion, and performing well-seismic combined lithologic inversion based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizon to obtain an inversion data volume.

[0040] The inversion data volume can be understood as the data volume obtained by inverting the logging parameters using the well-seismic combined lithologic inversion method.

[0041] In an embodiment of the present application, the top and bottom surface seismic layer interpretations are performed on the strata owned by the corresponding wells in the study area to obtain the seismic layer information required for the subsequent well-seismic combined lithologic inversion, and the logging sensitivity parameter analysis is performed on the corresponding wells to obtain the logging parameters used for the well-seismic combined lithologic inversion. The well-seismic calibration is then performed using the acoustic wave time difference data, and geological modeling is performed in combination with the seismic layer information and the logging parameters to obtain the initial model of the lithologic inversion. Specifically, based on the seismic data of the study area and with the initial lithologic inversion model as a constraint condition, the well-seismic combined lithologic inversion is performed using the lithologic inversion commercial software algorithm to determine the inversion data body of the study area.

[0042] S120 , performing well logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells, wherein the sandstone threshold value is a well logging parameter value used to determine the lithologic properties of the wells.

[0043] Specifically, well logging lithologic interpretation is performed on the wells involved in the joint well-seismic lithologic inversion in the study area to obtain the well logging lithologic parameters. The corresponding lithologic parameters include but are not limited to natural gamma rays. The well logging lithologic parameters of the wells are used as the sandstone threshold value of the wells, and the sandstone threshold value is used to determine the lithology of the corresponding wells in the study area.

[0044] S130. Based on the inverse distance interpolation method and in accordance with the sandstone threshold value, threshold plane data is obtained, and a sandstone threshold value data volume is obtained according to the vertical sampling interval of the seismic data and the threshold plane data, wherein the threshold plane data is plane grid data presented in coordinate form obtained by the inverse distance interpolation method from the sandstone threshold value of the corresponding well in the study area.

[0045] Specifically, the coordinates of the study area are gridded according to the geodetic coordinates. Before gridding, the line spacing and trace spacing of the seismic lines in the study area need to be determined. The line spacing is the horizontal distance between the coordinates in the grid cells that constitute the coordinate gridding of the study area, and the trace spacing is the vertical distance between the coordinates in the grid cells that constitute the coordinate gridding of the study area. The generated plane grid coordinate points correspond one-to-one with the seismic network coordinate points. On this basis, the threshold plane data is obtained based on the sandstone threshold value of the corresponding drilling well using the inverse distance interpolation method. The process and formula for obtaining the threshold plane data are as follows:

[0046]

[0047] Among them, x i and y i are the horizontal and vertical coordinates of the coordinate points of the corresponding wells in the study area; x and y are the horizontal and vertical coordinates of the unknown coordinate points outside the well coordinate points in the study area; d i is the distance between the drilling coordinate point and other unknown coordinate points in the study area, which can be determined based on the Pythagorean theorem of a right triangle; i represents the serial number of the known drilling location in the study area.

[0048]

[0049] Among them, λ i is the distance weight between an unknown coordinate point and a known well.

[0050]

[0051] Among them, M(x,y) represents the sandstone threshold value of an unknown coordinate point in the plane grid.

[0052] For example, Figure 2The coordinate positions of the three wells A, B, and C in the study area are shown. The sandstone threshold values ​​of the three wells A, B, and C are m1, m2, and m3, respectively. The sandstone threshold value of each unknown coordinate point in the plane grid of the study area is calculated according to the above formula. On this basis, based on the threshold value of the plane grid coordinate point in the study area, the vertical sampling interval of the seismic data is used as the spacing, and the top and bottom interfaces of the rock formation are used as the boundary. The plane grid coordinate point threshold value is vertically replicated at equal distances between the top and bottom interfaces to generate a sandstone threshold value data volume with spatial coordinate characteristics, as shown in the figure below. Figure 3 As shown, the vertical sampling interval of seismic data is the seismic data acquisition parameter of the study area, which is usually a fixed constant, and its meaning is how many milliseconds (ms) it takes to collect a data point in the vertical direction.

[0053] S140 , determining a residual value data volume according to the inversion data volume and the sandstone threshold value data volume, and obtaining a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0054] Specifically, the residual value is obtained by subtracting the numerical value marked at the same coordinate position of the inversion data body and the sandstone threshold value data body, and the residual value data body is obtained by replacing the data at the same coordinate position of the sandstone threshold value data body with the residual value. According to the relationship between the residual value and zero, each coordinate point on the residual value data body is reassigned to 0 or 1. If the residual value is greater than or equal to zero, the residual value of the point is reassigned to 1. If the residual value is less than zero, the residual value of the point is reassigned to 0. By determining the relationship between the residual value and the sandstone lithology of the study area, the residual value data body is statistically analyzed to vertically reflect the sandstone data. The number of data points, that is, the number of data points with a vertical value of 0 or 1 in the statistical residual value data body, if the residual value is 1, it represents sandstone, then the number of data points with a value of 1 in the vertical direction of the plane coordinate point is counted; on the contrary, if the residual value is 0, it represents sandstone, then the number of data points with a value of 0 in the vertical direction of the plane coordinate point is counted. The number of sandstone data points counted is multiplied by the vertical sampling interval (time) of the seismic data to obtain the time domain sandstone thickness map. The time domain sandstone thickness map can be converted into a depth domain sandstone thickness map through the time-depth conversion relationship of the study area, and then the sandstone thickness of the study area is determined.

[0055] The present application provides a method for obtaining a sandstone thickness map. The method comprises performing a well logging sensitivity parameter analysis on corresponding wells in a study area to determine the logging parameters used for lithologic inversion. A well-seismic combined lithologic inversion is performed based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizon to obtain an inversion data volume. A well-seismic combined lithologic inversion is performed to obtain an inversion data volume. A well logging lithologic interpretation is performed on corresponding wells in the study area to determine the sandstone threshold value of the wells. Threshold plane data is obtained based on the sandstone threshold value using an inverse distance interpolation method. A sandstone threshold data volume is obtained based on the vertical sampling interval of the seismic data and the threshold plane data. A residual value data volume is determined based on the inversion data volume and the sandstone threshold data volume. A sandstone thickness map is obtained by counting the number of data points in the residual value data volume that vertically reflect sandstone data. This method performs well-seismic combined lithologic inversion using logging parameters and obtains sandstone thickness using the sandstone threshold data volume. This method can improve the accuracy of sandstone thickness prediction and the accuracy of sandstone thickness map drawing.

[0056] Figure 4 This is a flow chart of a method for obtaining a sandstone thickness map provided by an embodiment of the present invention. This embodiment, based on the above embodiment, details the process of determining a residual value data volume based on the inversion data volume and the sandstone threshold value data volume, and obtaining a sandstone thickness map by counting the number of data points in the residual value data volume that vertically reflect sandstone data. Figure 4 As shown, the method may include:

[0057] S210: Perform seismic layer interpretation on the top and bottom surfaces of the rock strata in the study area to obtain seismic layer information.

[0058] Specifically, seismic horizon interpretation is performed on the top and bottom surfaces of the rock strata in the study area to obtain the seismic horizon information required for joint well-seismic lithologic inversion.

[0059] S220 , performing well-seismic calibration based on the acoustic wave time difference data and performing geological modeling according to the seismic layer information and logging parameters to obtain an initial lithologic inversion model.

[0060] Among them, the logging parameters are parameters that can reflect the sandstone lithology in the study area, including but not limited to natural gamma logging. The initial lithology inversion model is a logging parameter data model obtained through real logging data for predicting the sandstone lithology in the study area.

[0061] Specifically, well-seismic calibration is performed based on acoustic wave time difference data, and geological modeling is performed in combination with seismic layer information and logging parameters to form an initial lithologic inversion model.

[0062] S230. Performing well-seismic combined lithologic inversion based on the seismic data of the study area and the initial lithologic inversion model to generate a lithologic inversion data volume for the study area. Specifically, based on the seismic data of the study area and with the initial lithologic inversion model as a constraint, performing well-seismic combined lithologic inversion using a commercial lithologic inversion software algorithm to generate a lithologic inversion data volume for the study area.

[0063] S240 , performing well logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells, wherein the sandstone threshold value is a well logging parameter value used to determine the lithologic properties of the wells.

[0064] S250. Based on the inverse distance interpolation method and in accordance with the sandstone threshold value, threshold plane data is obtained, and a sandstone threshold value data body is obtained according to the vertical sampling interval of the seismic data and the threshold plane data, wherein the threshold plane data is plane grid data presented in coordinate form obtained by the inverse distance interpolation method from the sandstone threshold value of the corresponding well in the study area.

[0065] S260 , determining a residual value data volume according to the inversion data volume and the sandstone threshold value data volume, and obtaining a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0066] Specifically, the residual value data body is determined based on the inversion data body and the sandstone threshold value data body obtained in S210-S250, and the number of data points vertically reflecting the sandstone data in the residual value data body on the spatial coordinate system is counted to generate a plane diagram of the number of sandstone data points and multiply the data in the plane diagram of the number of sandstone data points with the vertical sampling interval of the seismic data to generate a sandstone thickness map in the time domain. The sandstone thickness map in the time domain is converted into a sandstone thickness map in the depth domain through the time-depth conversion relationship of the study area, thereby determining the sandstone thickness of the study area.

[0067] Optionally, the residual value data volume is volume data that presents the residual values ​​of the inversion data volume and the threshold sandstone data volume in the form of coordinates; each coordinate point in the inversion data volume, threshold value data volume and residual value data volume has a one-to-one correspondence.

[0068] As an optional but not limiting implementation, determining the residual value data volume according to the inversion data volume and the threshold value data volume may include steps A1-A7:

[0069] Step A1: sequentially subtract the lithologic inversion parameter values ​​of the inversion data volume from the threshold values ​​of the corresponding coordinate positions of the sandstone threshold value data volume to obtain the residual value of each data point, and form the residual value data volume based on the residual value of each data point.

[0070] Specifically, the values ​​of each coordinate point of the inversion data volume and the sandstone threshold value data volume in the spatial coordinate system are traversed respectively to determine the values ​​of each coordinate point of the inversion data volume and the sandstone threshold value data volume in the spatial coordinate system, the lithologic inversion parameter value obtained on the inversion data volume is subtracted from the threshold value obtained on the sandstone threshold value data volume to obtain the residual value, and the residual value is marked at the corresponding coordinate position on the residual data volume spatial coordinate system to obtain the residual value data volume.

[0071] For example, the lithologic inversion parameter values ​​marked at the coordinate points of the inversion data body are traversed respectively, and the threshold values ​​marked at the corresponding coordinate points of the sandstone threshold value data body are subtracted to obtain the residual values, and the residual values ​​are re-marked at the corresponding coordinate points in the spatial coordinate system of the residual value data body, so as to obtain the following: Figure 5 The residual value data body is shown.

[0072] Step A2: Determine the relationship between the residual value displayed at each coordinate position in the residual value data body and zero; if the residual value displayed at the coordinate position in the residual value data body is greater than or equal to zero, assign the residual value to 1.

[0073] Specifically, the relationship between the residual values ​​displayed at all coordinate positions on the residual value data body and 0 is determined, the residual values ​​displayed at each coordinate position on the residual value data body are traversed in turn, and it is determined whether the residual value is greater than or equal to 0. If so, the residual value is assigned a value of 1.

[0074] Step A3: If the residual value displayed by the coordinate position in the residual value data body is less than zero, the residual value is assigned to 0.

[0075] Specifically, if it is determined that the residual value is less than 0, the residual value is assigned to 0.

[0076] For example, Figure 6 The residual value data body represented by 0 and 1 obtained after assignment is shown, which serves as data support for counting the number of 0 and 1 in the residual value data body.

[0077] As an optional but not limiting implementation, obtaining a sandstone thickness map by counting the number of vertically reflecting sandstone data points in the residual value data volume may include steps B1-B4:

[0078] Step B1: Determine the relationship between the lithologic inversion parameter value and the sandstone threshold value in the study area.

[0079] Specifically, the relationship between the lithologic inversion parameter value and the sandstone threshold value at each coordinate point of the inversion data volume in the spatial coordinate system is determined.

[0080] Step B2: If the lithologic inversion parameter value is greater than or equal to the sandstone threshold value, which indicates sandstone, then the number of data points with a vertical residual value of 1 in the residual value data volume is counted.

[0081] Specifically, if the lithologic inversion parameter value of a coordinate point in the spatial coordinate system is greater than or equal to the sandstone threshold value, the lithologic type represented by the coordinate point is determined to be sandstone, and the number of data points with a residual value of 1 in the vertical direction of the coordinate point is counted.

[0082] Step B3: If the lithologic inversion parameter value is less than the sandstone threshold value, which indicates sandstone, then the number of data points with a vertical residual value of 0 in the residual value data volume is counted.

[0083] Specifically, if the lithologic inversion parameter value of a certain coordinate point in the spatial coordinate system is less than the sandstone threshold value, the lithologic type represented by the coordinate point is determined to be sandstone, and the number of data points with a residual value of 0 in the vertical direction of the coordinate point is counted.

[0084] Step B4: Draw a plane diagram of the number of sandstone data points based on the number of data points with vertical residual values ​​of 1 or the number of data points with residual values ​​of 0 in the residual value data volume, wherein the plane coordinate points of the plane diagram of the number of sandstone data points correspond one-to-one to the plane coordinate points of the sandstone threshold value.

[0085] In the embodiment of the present application, by counting the number of data points with vertical residual values ​​of 1 or the number of data points with residual values ​​of 0 in the residual value data volume, a plane diagram of the number of sandstone data points is drawn, such as Figure 7 As shown, the sandstone thickness map in the time domain can be determined through the plane map of the number of sandstone data points.

[0086] Optionally, a contour map is drawn based on the value of each coordinate point in the plane diagram of the number of sandstone data points multiplied by the vertical sampling interval of the seismic data to obtain a sandstone thickness map in the time domain; and the sandstone thickness map in the time domain is converted into a sandstone thickness map in the depth domain by using the time and depth conversion relationship of the study area.

[0087] In the embodiment of the present application, based on the number of data points representing sandstone on the plane graph of the number of sandstone data points, that is, by counting the number of data points whose vertical coordinates of the residual value data volume are 1 or 0, the number of data points representing sandstone is multiplied by the vertical sampling interval of the seismic data, and a contour mapping operation is performed to obtain a sandstone thickness map in the time domain. At the same time, based on the time and depth conversion relationship of the study area, the time domain sandstone thickness map is converted into a sandstone thickness map in the depth domain. For example, Figure 8 Sandstone thickness map of the study area at depth is shown.

[0088] An embodiment of the present application proposes a method for obtaining a sandstone thickness map, which comprises the following steps: performing logging sensitivity parameter analysis on corresponding wells in a study area to determine logging parameters for lithologic inversion; performing well-seismic combined lithologic inversion based on the logging parameters, a seismic data volume in the study area, and a seismic interpretation layer to obtain an inversion data volume; performing logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value for the wells; obtaining threshold plane data based on an inverse distance interpolation method and in accordance with the sandstone threshold value; obtaining a sandstone threshold data volume based on a vertical sampling interval of seismic data and the threshold plane data; determining a residual data volume based on the inversion data volume and the sandstone threshold data volume; obtaining a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual data volume; counting the number of data points representing sandstone on the residual data volume and drawing a plane map of the number of sandstone data points; and determining a sandstone time thickness map by multiplying the number of data points by the vertical sampling interval of seismic data. The sandstone thickness map can be determined by the time-depth conversion relationship of the study area. This method performs joint well-seismic lithologic inversion using logging parameters reflecting sandstone lithology, and obtains sandstone thickness using sandstone threshold data volumes, which can improve the accuracy of sandstone thickness prediction and sandstone thickness mapping.

[0089] Figure 9 This is a flowchart of a method for obtaining a sandstone thickness map provided by an embodiment of the present invention. This embodiment, based on the above embodiment, details the process of obtaining threshold plane data based on the inverse distance interpolation method and the sandstone threshold value, and obtaining a sandstone thickness map based on the vertical sampling interval of the seismic data and the threshold plane data. Figure 9 As shown, the method may include:

[0090] S310, performing logging sensitive parameter analysis on corresponding wells in the study area to determine logging parameters for lithologic inversion, and performing well-seismic combined lithologic inversion based on the logging parameters, the seismic data volume of the study area and the seismic interpretation horizon to obtain an inversion data volume.

[0091] S320: Perform logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells, wherein the sandstone threshold value is a logging parameter value used to determine the lithologic properties of the wells.

[0092] S330, based on the line spacing and trace spacing of the seismic survey lines in the study area as grid units, coordinate gridding is performed on the study area to generate a plane coordinate grid that corresponds one-to-one to the coordinate points of the seismic survey network.

[0093] S340 , assigning values ​​to data points on the plane coordinate grid based on the sandstone threshold value of the study area and the inverse distance interpolation method to generate plane data of the sandstone threshold value of the study area.

[0094] Specifically, based on the sandstone threshold values ​​of known wells in the study area, the inverse distance interpolation method is used to assign values ​​to the data points on the plane coordinate grid to generate the plane data of the sandstone threshold values ​​in the study area.

[0095] S350: Using the vertical sampling interval of the seismic data as the spacing and the top and bottom interfaces of the rock formation as the boundary, the sandstone threshold plane data is vertically replicated at equal distances between the top and bottom interfaces of the rock formation to determine the sandstone threshold data volume.

[0096] Based on the plane grid coordinate point values ​​of the sandstone threshold value in the study area, with the vertical sampling interval of seismic data as the spacing, and the top and bottom interfaces of the rock layer as the boundary, the plane grid coordinate point values ​​are vertically replicated at equal distances between the top and bottom interfaces to determine the sandstone threshold value data body with spatial coordinate characteristics.

[0097] S360: Determine a residual value data volume according to the inversion data volume and the sandstone threshold value data volume, and obtain a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0098] The present application provides a method for obtaining a sandstone thickness map. The method comprises performing well logging sensitivity parameter analysis on corresponding wells in a study area to determine logging parameters for lithologic inversion, performing well-seismic combined lithologic inversion based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizon to obtain an inversion data volume. A well logging lithologic interpretation is performed on corresponding wells in the study area to determine the sandstone threshold value of the wells. Threshold plane data is obtained based on the inverse distance interpolation method and the sandstone threshold value, and a sandstone threshold value data volume is obtained based on the vertical sampling interval of the seismic data and the threshold plane data. A residual value data volume is determined based on the inversion data volume and the sandstone threshold value data volume, and a sandstone thickness map is obtained by counting the number of data points in the residual value data volume that vertically reflect sandstone data. This method performs well-seismic combined lithologic inversion using logging parameters that reflect sandstone lithology and obtains sandstone thickness using the sandstone threshold value data volume. This method can improve the accuracy of sandstone thickness prediction and the accuracy of sandstone thickness map drawing.

[0099] Example 4

[0100] Figure 10 This is a schematic diagram of the structure of a device for obtaining a sandstone thickness map provided in the fourth embodiment of the present invention. The device can execute the method for obtaining a sandstone thickness map provided in any embodiment of the present application, and has the corresponding functional modules and beneficial effects of the execution method. Figure 10 As shown, the device may include:

[0101] The inversion data volume acquisition module 410 is used to perform well logging sensitivity parameter analysis on the corresponding wells in the study area to determine the well logging parameters used for lithologic inversion, and perform well-seismic combined lithologic inversion based on the well logging parameters, the seismic data volume of the study area, and the seismic interpretation horizons to obtain the inversion data volume;

[0102] A threshold value determination module 420 is configured to perform well logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells, wherein the sandstone threshold value is a well logging parameter value used to determine the lithologic properties of the wells;

[0103] A threshold data volume acquisition module 430 is configured to acquire threshold plane data based on the sandstone threshold value using an inverse distance interpolation method, and to acquire a sandstone threshold data volume based on the vertical sampling interval of the seismic data and the threshold plane data, wherein the threshold plane data is plane grid data presented in coordinate form obtained by using the inverse distance interpolation method based on the sandstone threshold values ​​of corresponding wells in the study area;

[0104] The sandstone thickness map acquisition module 440 is configured to determine a residual value data volume based on the inversion data volume and the sandstone threshold value data volume, and acquire a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume.

[0105] Furthermore, the residual value data volume is volume data that presents the residual values ​​of the inversion data volume and the sandstone threshold value data volume in the form of coordinates; each coordinate point in the inversion data volume, threshold value data volume and residual value data volume has a one-to-one correspondence.

[0106] Furthermore, the sandstone thickness map acquisition module 440 includes:

[0107] a residual value data volume determining unit, configured to sequentially subtract the lithologic inversion parameter values ​​of the inversion data volume from the threshold values ​​of the corresponding coordinate positions of the sandstone threshold value data volume to obtain a residual value for each data point, and form the residual value data volume based on the residual value of each data point;

[0108] a first residual determination unit, configured to determine a relationship between a residual value displayed at each coordinate position in the residual value data body and zero, and assigning a value of 1 to the residual value if the residual value displayed at the coordinate position in the residual value data body is greater than or equal to zero;

[0109] The second residual determination unit is configured to assign the residual value to 0 if the residual value displayed by the coordinate position in the residual value data body is less than zero.

[0110] Furthermore, the sandstone thickness map acquisition module 440 further includes:

[0111] A parameter and threshold value relationship determination unit, used to determine the relationship between the lithologic inversion parameter value and the sandstone threshold value in the study area;

[0112] a first coordinate point number counting unit, configured to count the number of data points with a residual value of 1 in a vertical direction in the residual value data volume if the lithologic inversion parameter value is greater than or equal to the sandstone threshold value, representing sandstone;

[0113] a second coordinate point number counting unit, configured to count the number of data points with a residual value of 0 in a vertical direction in the residual value data volume if the lithologic inversion parameter value is less than the sandstone threshold value, indicating sandstone;

[0114] A plane map acquisition unit is used to draw a plane map of the number of sandstone data points based on the number of data points with vertical residual values ​​of 1 or the number of data points with residual values ​​of 0 in the residual value data volume, wherein the plane coordinate points of the plane map of the number of sandstone data points correspond one-to-one to the plane coordinate points of the sandstone threshold value.

[0115] Furthermore, the sandstone thickness map acquisition module 440 further includes:

[0116] A time thickness map acquisition unit is used to obtain a time domain sandstone thickness map by performing contour mapping based on the value of each coordinate point in the plane map of the number of sandstone data points multiplied by the vertical sampling interval of the seismic data;

[0117] The sandstone thickness map acquisition unit is used to convert the sandstone thickness map in the time domain into a sandstone thickness map in the depth domain by using the time-depth conversion relationship of the study area.

[0118] Furthermore, the threshold data acquisition module 430 includes:

[0119] The coordinate grid forming unit is used to coordinate grid the study area based on the line spacing and trace spacing of the seismic survey lines in the study area as grid units, and generate a plane coordinate grid that corresponds one-to-one with the coordinate points of the seismic survey network;

[0120] A grid intersection assignment unit is used to assign values ​​to data points on the plane coordinate grid based on the sandstone threshold value of the study area and the inverse distance interpolation method to generate plane data of the sandstone threshold value of the study area;

[0121] The threshold value data volume determination unit is used to use the vertical sampling interval of the seismic data as the spacing, the top and bottom interfaces of the rock layer as the boundary, and vertically replicate the sandstone threshold value plane data at equal distances between the top and bottom interfaces of the rock layer to determine the sandstone threshold value data volume.

[0122] Furthermore, the inversion data volume acquisition module 410 includes:

[0123] A seismic layer information acquisition unit, configured to perform seismic layer interpretation on the top and bottom surfaces of the rock layers in the study area to obtain seismic layer information;

[0124] An inversion initial model construction unit is used to perform well-seismic calibration based on acoustic wave time difference data and perform geological modeling based on the seismic layer information and logging parameters to obtain a lithologic inversion initial model;

[0125] The inversion data volume determination unit performs well-seismic combined lithologic inversion based on the seismic data of the study area and the lithologic inversion initial model to generate the lithologic inversion data volume of the study area.

[0126] The device for obtaining a sandstone thickness map provided in an embodiment of the present invention can execute the method for obtaining a sandstone thickness map provided in any embodiment of the present invention, and has the corresponding functions and beneficial effects of executing the method for obtaining a sandstone thickness map. For detailed processes, please refer to the relevant operations of the method for obtaining a sandstone thickness map in the aforementioned embodiment.

[0127] Example 5

[0128] Figure 11 A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0129] like Figure 11 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. Various programs and data required for the operation of the electronic device 10 can also be stored in the RAM 13. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0130] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0131] The processor 11 can be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any other suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the method for obtaining a sandstone thickness map.

[0132] In some embodiments, the sandstone thickness map acquisition method can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the sandstone thickness map acquisition method described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to execute the sandstone thickness map acquisition method in any other suitable manner (e.g., via firmware).

[0133] Various embodiments of the systems and techniques described herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system comprising at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0134] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0135] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0136] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0137] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0138] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0139] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in the present invention can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved. This is not limited herein.

[0140] The above specific embodiments do not limit the scope of protection of the present invention. Those skilled in the art will appreciate that various modifications, combinations, sub-combinations, and substitutions may be made based on design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention are intended to be included within the scope of protection of the present invention.

Claims

1. A method for obtaining a sandstone thickness map, characterized in that: include: Conduct logging sensitivity parameter analysis on the corresponding wells in the study area to determine the logging parameters used for lithologic inversion, and perform well-seismic combined lithologic inversion based on the logging parameters, the seismic data volume of the study area, and the seismic interpretation horizons to obtain the inversion data volume; Performing well logging lithologic interpretation on corresponding wells in the study area to determine a sandstone threshold value of the wells; wherein the sandstone threshold value is a well logging parameter value used to determine the lithologic properties of the wells; Based on the inverse distance interpolation method and the sandstone threshold value, threshold plane data is obtained, and a sandstone threshold data volume is obtained according to the vertical sampling interval of the seismic data and the threshold plane data; wherein the threshold plane data is plane grid data presented in the form of coordinates obtained by the inverse distance interpolation method from the sandstone threshold values ​​of the corresponding wells in the study area; determining a residual value data volume according to the inversion data volume and the sandstone threshold value data volume, and obtaining a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume; The residual value data volume is volume data in which the residual values ​​of the inversion data volume and the threshold value data volume are presented in the form of coordinates; each coordinate point in the inversion data volume, the threshold value data volume and the residual value data volume has a one-to-one correspondence. The method of determining a residual value data body according to the inversion data body and the sandstone threshold value data body comprises: sequentially subtracting the lithologic inversion parameter value of the inversion data body from the threshold value of the corresponding coordinate position of the sandstone threshold value data body to obtain a residual value of each data point, and forming the residual value data body according to the residual value of each data point; determining a relationship between a residual value displayed at each coordinate position in the residual value data body and zero, and if the residual value displayed at the coordinate position in the residual value data body is greater than or equal to zero, assigning the residual value to 1; and if the residual value displayed at the coordinate position in the residual value data body is less than zero, assigning the residual value to 0; The sandstone thickness map is obtained by counting the number of data points vertically reflecting sandstone data in the residual value data volume, including: Determine the relationship between the lithologic inversion parameter values ​​and the sandstone threshold values ​​in the study area; If the lithologic inversion parameter value is greater than or equal to the sandstone threshold value, representing sandstone, then counting the number of data points with a vertical residual value of 1 in the residual value data volume; If the lithologic inversion parameter value is less than the sandstone threshold value, which indicates sandstone, then counting the number of data points with a vertical residual value of 0 in the residual value data volume; Draw a sandstone data point quantity plane diagram based on the number of data points with vertical residual values ​​of 1 or the number of data points with vertical residual values ​​of 0 in the residual value data volume; wherein the plane coordinate points of the sandstone data point quantity plane diagram correspond one-to-one to the plane coordinate points of the sandstone threshold value; The sandstone thickness map in the time domain is obtained by multiplying the value of each coordinate point in the plane diagram of the number of sandstone data points by the vertical sampling interval of the seismic data to perform contour mapping; The sandstone thickness map in the time domain is converted into a sandstone thickness map in the depth domain by using the time-depth conversion relationship of the study area.

2. The method according to claim 1, characterized in that The method includes: obtaining threshold plane data based on the inverse distance interpolation method and the sandstone threshold value, and obtaining a sandstone threshold data volume based on the vertical sampling interval of the seismic data and the threshold plane data, including: Based on the line spacing and trace spacing of the seismic survey lines in the study area as grid units, the study area is gridded to generate a plane coordinate grid that corresponds one-to-one to the coordinate points of the seismic survey network. Assigning values ​​to data points on the plane coordinate grid based on the sandstone threshold value of the study area and the inverse distance interpolation method to generate plane data of the sandstone threshold value of the study area; The vertical sampling interval of seismic data is used as the spacing, the top and bottom interfaces of the rock formation are used as the boundary, and the sandstone threshold plane data is vertically replicated at equal distances between the top and bottom interfaces of the rock formation to determine the sandstone threshold data volume.

3. The method according to claim 1, characterized in that Perform logging sensitivity parameter analysis on the corresponding wells in the study area to determine the logging parameters used for lithologic inversion. Based on the logging parameters, the seismic data volume of the study area and the seismic interpretation horizons, perform well-seismic combined lithologic inversion to obtain the inversion data volume, including: Conduct seismic layer interpretation on the top and bottom surfaces of the rock formations in the study area to obtain seismic layer information; Perform well-seismic calibration based on acoustic time difference data and perform geological modeling based on the seismic layer information and logging parameters to obtain an initial lithologic inversion model; Based on the seismic data of the study area and combined with the initial lithologic inversion model, well-seismic joint lithologic inversion is performed to generate the lithologic inversion data volume of the study area.

4. A device for obtaining a sandstone thickness map, characterized in that: include: The inversion data volume acquisition module is used to replicate the sandstone threshold value plane data vertically and equidistantly between the top and bottom interfaces of the rock formation with the top and bottom interfaces of the rock formation as the boundary, thereby determining the sandstone threshold value data volume; A threshold value determination module is used to perform well logging lithology interpretation on corresponding wells in the study area to determine the sandstone threshold value of the wells, where the sandstone threshold value is a well logging parameter value used to determine the lithology; a threshold value data volume acquisition module, configured to acquire threshold value plane data based on the sandstone threshold value using an inverse distance interpolation method, and to acquire a sandstone threshold value data volume based on the vertical sampling interval of the seismic data and the threshold value plane data; wherein the threshold value plane data is plane grid data presented in coordinate form, acquired by using the inverse distance interpolation method from the sandstone threshold values ​​of corresponding wells in the study area; a sandstone thickness map acquisition module, configured to determine a residual value data volume based on the inversion data volume and the sandstone threshold value data volume, and acquire a sandstone thickness map by counting the number of data points vertically reflecting sandstone data in the residual value data volume; Wherein, the threshold value determination module includes: a residual value data volume determining unit, configured to sequentially subtract the lithologic inversion parameter values ​​of the inversion data volume from the threshold values ​​of the corresponding coordinate positions of the sandstone threshold value data volume to obtain a residual value for each data point, and form the residual value data volume based on the residual value of each data point; a first residual determination unit, configured to determine a relationship between a residual value displayed at each coordinate position in the residual value data body and zero, and assigning a value of 1 to the residual value if the residual value displayed at the coordinate position in the residual value data body is greater than or equal to zero; a second residual determination unit, configured to assign the residual value to 0 if the residual value displayed by the coordinate position in the residual value data body is less than zero; A parameter and threshold value relationship determination unit is used to determine the relationship between the lithologic inversion parameter value and the sandstone threshold value in the study area; a first coordinate point number counting unit, configured to count the number of data points with a residual value of 1 in a vertical direction in the residual value data volume if the lithologic inversion parameter value is greater than or equal to the sandstone threshold value, representing sandstone; a second coordinate point number counting unit, configured to count the number of data points with a residual value of 0 in a vertical direction in the residual value data volume if the lithologic inversion parameter value is less than the sandstone threshold value, indicating sandstone; a plane map acquiring unit, configured to draw a plane map of the number of sandstone data points based on the number of data points with a vertical residual value of 1 or the number of data points with a residual value of 0 in the residual value data volume; wherein the plane coordinate points of the plane map of the number of sandstone data points correspond one-to-one to the plane coordinate points of the sandstone threshold value; A time thickness map acquisition unit is used to obtain a time domain sandstone thickness map by performing contour mapping based on the value of each coordinate point in the plane map of the number of sandstone data points multiplied by the vertical sampling interval of the seismic data; The sandstone thickness map acquisition unit is used to convert the sandstone thickness map in the time domain into a sandstone thickness map in the depth domain by using the time-depth conversion relationship of the study area.

5. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor. The computer program is executed by the at least one processor to enable the at least one processor to perform the method for obtaining a sandstone thickness map according to any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for obtaining a sandstone thickness map according to any one of claims 1 to 3 when executed.

Images