Reservoir characterization method and device

Abstract

A reservoir characterization method and device belongs to the field of geophysical exploration technology. The method includes: obtaining a target well logging intersection diagram and a seismic data body, the target well logging intersection diagram includes a target reservoir area and multiple well logging sample points, each well logging sample point is used to characterize at least two types of well logging data of a position point in the target area, and the seismic data body includes multiple seismic sample points, each seismic sample point is used to characterize multiple types of seismic data of a position point in the target area; setting a target index network in the target well logging intersection diagram according to the target reservoir area, the target index network includes multiple index grids and multiple grid nodes; determining the index value of the seismic sample point in the seismic data body according to the target index network to obtain an index data body, and for the seismic sample point located in the target index grid, determining the index value of the seismic sample point according to the internal and external relationship between the seismic sample point and the target reservoir area; and characterizing the target reservoir according to the index data body. This application helps to improve the efficiency of reservoir characterization.

Description

Technical Field

[0001] The present application relates to the field of geophysical exploration technology, and in particular to a reservoir characterization method and device. Background Art

[0002] A reservoir refers to a rock formation that stores oil and gas. Reservoir characterization refers to predicting the distribution range, lithologic characteristics, and spatial distribution of a reservoir by analyzing characteristic data of the geographical area where the reservoir is located.

[0003] When characterizing a reservoir, at least two types of well logging data and a seismic data volume for the geographic region where the reservoir is located are first obtained. A well logging intersection diagram is then determined based on the at least two types of well logging data. Finally, the spatial distribution of the reservoir is characterized based on the well logging intersection diagram and the seismic data volume. The well logging intersection diagram includes a plurality of well logging points, each of which is used to represent at least two types of well logging data for a location point in the geographic region. The seismic data volume includes a plurality of seismic points, each of which is used to represent multiple types of seismic data for a location point in the geographic region, wherein the attributes of at least two types of seismic data are the same as the attributes of the at least two types of well logging data. When characterizing the spatial distribution of the reservoir based on the well logging intersection diagram and the seismic data volume, the reservoir region is first determined in the well logging intersection diagram. Then, the internal and external relationships between each seismic point in the seismic data volume and the reservoir region are determined. Finally, the spatial distribution of the reservoir is characterized based on the seismic points in the seismic data volume located in the reservoir region.

[0004] However, seismic data volumes usually include a large number of seismic sample points. When characterizing the reservoir, the above scheme needs to determine the internal and external relationships between all seismic sample points in the seismic data volume and the reservoir area, resulting in a large amount of computational complexity, long time consumption, and low efficiency in reservoir characterization. Summary of the Invention

[0005] This application provides a reservoir characterization method and device that can reduce the computational complexity of reservoir characterization, shorten the time consumption of reservoir characterization, and improve the efficiency of reservoir characterization. The technical solution of this application is as follows:

[0006] In a first aspect, a reservoir characterization method is provided, the method comprising:

[0007] Acquiring a target well logging crossplot of a target area and a seismic data volume of the target area, wherein the target area is a geographical area where a target reservoir is located, the target well logging crossplot includes a target reservoir area and a plurality of well logging points, each well logging point being used to represent at least two types of well logging data for a location point in the target area, the target reservoir area being used to represent the target reservoir, the seismic data volume including a plurality of seismic points, each seismic point being used to represent multiple types of seismic data for a location point in the target area, and attributes of at least two types of seismic data being the same as attributes of the at least two types of well logging data;

[0008] A target index network is set in the target logging crossplot according to the target reservoir region, wherein a graph formed by an outer contour of the target index network covers the target reservoir region, the target index network includes a plurality of index grids and a plurality of grid nodes, the plurality of grid nodes are intersections of grid lines constituting the target index network, and each of the plurality of grid nodes has an index value;

[0009] determining an index value of each seismic sample point in the seismic data volume according to the target index grid, and obtaining an index data volume corresponding to the seismic data volume, wherein, for each seismic sample point in the seismic data volume located within a target index grid, the index value of each seismic sample point is determined according to an internal and external relationship between each seismic sample point and the target reservoir region, the multiple index grids include the target index grid, and a boundary line of the target reservoir region passes through the target index grid;

[0010] The target reservoir is characterized according to the index data volume.

[0011] Optionally, obtaining a target well logging crossplot of the target area includes:

[0012] Acquiring at least two types of well logging data of the target area;

[0013] The target logging cross-plot is obtained according to the at least two types of logging data.

[0014] Optionally, determining the index value of each seismic sample point in the seismic data volume according to the target index network includes:

[0015] For each seismic sample point in the seismic data volume that is located within a non-target index grid, an index value of each seismic sample point is determined based on an index value of a grid node corresponding to the non-target index grid, the multiple index grids including the non-target index grid, and the grid nodes corresponding to the non-target index grid are intersections of grid lines constituting the non-target index grid.

[0016] Optionally, the non-target index grid includes a first non-target index grid and a second non-target index grid, the first non-target index grid is located within the target reservoir region, and the second non-target index grid is located outside the target reservoir region, the index values ​​of the grid nodes corresponding to the first non-target index grid are all first index values, and the index values ​​of the grid nodes corresponding to the second non-target index grid are all second index values, and for each seismic sample point in the seismic data volume located within the non-target index grid, determining the index value of each seismic sample point according to the index value of the grid node corresponding to the non-target index grid includes:

[0017] For each seismic sample point in the seismic data volume and located within the first non-target index grid, determining an index value of each seismic sample point as the first index value;

[0018] For each seismic sample point in the seismic data volume that is located within the second non-target index grid, an index value of each seismic sample point is determined to be the second index value.

[0019] Optionally, the index values ​​of the grid nodes within the target reservoir region are all first index values, and the index values ​​of the grid nodes outside the target reservoir region are all second index values. For each seismic sample point in the seismic data volume located within the target index grid, determining the index value of each seismic sample point according to the internal and external relationship between each seismic sample point and the target reservoir region includes:

[0020] For each seismic sample point in the seismic data volume and located within the target index grid:

[0021] When each of the seismic sample points is located in the target reservoir area, determining the index value of each of the seismic sample points to be the first index value;

[0022] When each seismic sample point is located outside the target reservoir area, the index value of each seismic sample point is determined to be the second index value.

[0023] Optionally, the method further includes:

[0024] For each seismic sample point in the seismic data volume and located within the target index grid:

[0025] Drawing a characteristic ray with each seismic sample point as an endpoint, wherein the characteristic ray has at least one intersection point with a boundary line of the target reservoir area;

[0026] When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is an odd number, determining that each seismic sample point is located in the target reservoir region;

[0027] When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is even, it is determined that each seismic sample point is located outside the target reservoir region.

[0028] Optionally, setting a target index network in the target well logging crossplot according to the target reservoir region includes:

[0029] Setting an initial index network in the target logging crossplot according to the target reservoir region, so that a graph enclosed by an outer contour of the initial index network covers the target reservoir region, the initial index network including the plurality of index grids and the plurality of grid nodes;

[0030] According to the relationship between each grid node in the plurality of grid nodes and the target reservoir region, an index value is set for each grid node to obtain the target index network.

[0031] Optionally, setting an index value for each grid node among the plurality of grid nodes according to a relationship between the grid node and the target reservoir region includes:

[0032] For each grid node of the plurality of grid nodes:

[0033] When each grid node is located within the target reservoir region or on a boundary line of the target reservoir region, setting the index value of each grid node to a first index value;

[0034] When each grid node is located outside the target reservoir area, the index value of each grid node is set to a second index value.

[0035] Optionally, before setting a target index network in the target well logging crossmap according to the target reservoir region, the method further includes: determining the target reservoir region in the target well logging crossmap.

[0036] In a second aspect, a reservoir characterization device is provided, comprising:

[0037] an acquisition module, configured to acquire a target well logging crossplot of a target area and a seismic data volume of the target area, wherein the target area is a geographical area where a target reservoir is located, the target well logging crossplot including a target reservoir area and a plurality of well logging sample points, each well logging sample point being used to represent at least two types of well logging data at a location point in the target area, the target reservoir area being used to represent the target reservoir, the seismic data volume including a plurality of seismic sample points, each seismic sample point being used to represent multiple types of seismic data at a location point in the target area, and attributes of at least two types of seismic data being the same as attributes of the at least two types of well logging data;

[0038] a setting module, configured to set a target index network in the target logging crossplot according to the target reservoir region, wherein a graph formed by an outer contour of the target index network covers the target reservoir region, the target index network includes a plurality of index grids and a plurality of grid nodes, the plurality of grid nodes being intersection points of grid lines constituting the target index network, and each of the plurality of grid nodes having an index value;

[0039] an index determination module, configured to determine an index value of each seismic sample point in the seismic data volume based on the target index grid, and obtain an index data volume corresponding to the seismic data volume, wherein, for each seismic sample point in the seismic data volume located within a target index grid, the index value of each seismic sample point is determined based on an internal and external relationship between each seismic sample point and the target reservoir region, the multiple index grids including the target index grid, and the boundary line of the target reservoir region passing through the target index grid;

[0040] A characterization module is used to characterize the target reservoir according to the index data body.

[0041] Optionally, the index determination module is used to determine the index value of each seismic sample point in the seismic data body located in a non-target index grid according to the index value of the grid node corresponding to the non-target index grid, the multiple index grids include the non-target index grid, and the grid nodes corresponding to the non-target index grid are the intersection points of the grid lines constituting the non-target index grid.

[0042] Optionally, the non-target index grid includes a first non-target index grid and a second non-target index grid, the first non-target index grid is located within the target reservoir region, the second non-target index grid is located outside the target reservoir region, the index values ​​of the grid nodes corresponding to the first non-target index grid are all first index values, and the index values ​​of the grid nodes corresponding to the second non-target index grid are all second index values, and the index determination module is used to:

[0043] For each seismic sample point in the seismic data volume and located within the first non-target index grid, determining an index value of each seismic sample point as the first index value;

[0044] For each seismic sample point in the seismic data volume that is located within the second non-target index grid, an index value of each seismic sample point is determined to be the second index value.

[0045] Optionally, the index values ​​of the grid nodes within the target reservoir area are all first index values, and the index values ​​of the grid nodes outside the target reservoir area are all second index values, and the index determination module is used to:

[0046] For each seismic sample point in the seismic data volume and located within the target index grid:

[0047] When each of the seismic sample points is located in the target reservoir area, determining the index value of each of the seismic sample points to be the first index value;

[0048] When each seismic sample point is located outside the target reservoir area, the index value of each seismic sample point is determined to be the second index value.

[0049] Optionally, the device further includes a relationship determination module, configured to:

[0050] For each seismic sample point in the seismic data volume and located within the target index grid:

[0051] Drawing a characteristic ray with each seismic sample point as an endpoint, wherein the characteristic ray has at least one intersection point with a boundary line of the target reservoir area;

[0052] When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is an odd number, determining that each seismic sample point is located in the target reservoir region;

[0053] When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is even, it is determined that each seismic sample point is located outside the target reservoir region.

[0054] Optionally, the setting module is used to:

[0055] Setting an initial index network in the target logging crossplot according to the target reservoir region, so that a graph enclosed by an outer contour of the initial index network covers the target reservoir region, the initial index network including the plurality of index grids and the plurality of grid nodes;

[0056] According to the relationship between each grid node in the plurality of grid nodes and the target reservoir region, an index value is set for each grid node to obtain the target index network.

[0057] Optionally, the setting module is used to:

[0058] For each grid node of the plurality of grid nodes:

[0059] When each grid node is located within the target reservoir region or on a boundary line of the target reservoir region, setting the index value of each grid node to a first index value;

[0060] When each grid node is located outside the target reservoir area, the index value of each grid node is set to a second index value.

[0061] Optionally, the device further includes: a region determination module, configured to determine the target reservoir region in the target logging cross-plot before the setting module sets a target index network in the target logging cross-plot according to the target reservoir region.

[0062] In a third aspect, a computer device is provided, comprising a memory and a processor, wherein the memory stores a computer program, and the processor is configured to execute the computer program to implement the method provided in the first aspect or any optional implementation of the first aspect.

[0063] In a fourth aspect, a computer-readable storage medium is provided, in which a computer program is stored. When the computer program is executed, the method provided in the first aspect or any optional implementation of the first aspect is implemented.

[0064] In a fifth aspect, a computer program product is provided, which includes a program or code, and when the program or code is executed, it implements the method provided by the first aspect or any optional implementation of the first aspect.

[0065] The beneficial effects of the technical solution provided by this application are:

[0066] The reservoir characterization method and apparatus provided in the present application first obtain a target well logging intersection diagram and a seismic data volume of the target area during reservoir characterization. The target well logging intersection diagram includes a target reservoir area and multiple well logging sample points, and the seismic data volume includes multiple seismic sample points. Then, a target index network is set in the target well logging intersection diagram according to the target reservoir area. The target index network includes multiple index grids and multiple grid nodes, each of which has an index value. Then, the index value of each seismic sample point in the seismic data volume is determined according to the target index network to obtain an index data volume, and finally, the target reservoir is characterized according to the index data volume. The target index network includes a target index grid, and the boundary line of the target reservoir area passes through the target index grid. When determining the index data volume according to the target index network, for each seismic sample point in the seismic data volume located within the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir area. Since the internal and external relationship judgment of the target reservoir area is only performed for the seismic sample points located in the target index grid during reservoir characterization, the calculation amount and time consumption of reservoir characterization can be reduced, thereby improving the efficiency of reservoir characterization.

[0067] It should be understood that the foregoing general description and the following detailed description are exemplary only and are not restrictive of the invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0068] In order to more clearly illustrate the technical solutions in the embodiments of the present application, 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 application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0069] Figure 1 This is a flow chart of a reservoir characterization method provided in an embodiment of the present application;

[0070] Figure 2 This is an initial well logging crossplot provided by an embodiment of the present application;

[0071] Figure 3 is a schematic diagram of a target well logging crossplot provided in an embodiment of the present application;

[0072] Figure 4 This is a schematic diagram of a target reservoir region after value assignment provided in an embodiment of the present application;

[0073] Figure 5 This is a schematic diagram of an embodiment of the present application after an initial index network is set in a target well logging crossplot;

[0074] Figure 6 This is a schematic diagram of determining the index value of a grid node provided by an embodiment of the present application;

[0075] Figure 7 This is a schematic diagram of a target index network set in a target well logging crossplot provided by an embodiment of the present application;

[0076] Figure 8 This is a schematic diagram of determining the internal and external relationship between a seismic sample point and a target reservoir area provided by an embodiment of the present application;

[0077] Figure 9 This is a schematic diagram of determining the index value of a seismic sample point provided by an embodiment of the present application;

[0078] Figure 10 is a cross-sectional view of a target reservoir provided in an embodiment of the present application;

[0079] Figure 11 This is a schematic diagram of a reservoir characterization device provided in an embodiment of the present application.

[0080] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application. DETAILED DESCRIPTION

[0081] To make the objectives, technical solutions, and advantages of this application more clear, this application will be further described in detail below with reference to the accompanying drawings. Obviously, the embodiments described are only some of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making any creative efforts are within the scope of protection of this application.

[0082] In geophysical exploration, a reservoir is a rock formation with interconnected pores that allow oil and gas to store and percolate. Most of the world's discovered oil and gas reserves are found in sedimentary rock formations, with coal seams and carbonate fractured and fissured mudstone reservoirs being the most important. Industrial oil and gas reserves have also been discovered in igneous and metamorphic rock reservoirs. Reservoir characterization is a key step in geophysical exploration and development, providing effective data support for subsequent calculations of reservoir thickness and reserves.

[0083] Currently, reservoir characterization requires determining the internal and external relationships between all seismic points in a seismic data volume and the reservoir regions in well logging crossplots. However, seismic data volumes typically contain a massive number of seismic points. For example, a seismic data volume obtained in a typical work area may contain approximately 1 billion points, occupying approximately 10GB of storage space. Seismic data obtained in basin-level work areas typically occupy over 100GB of storage space. Determining the internal and external relationships between all seismic points in a seismic data volume and the reservoir regions can lead to high computational complexity, long processing time, and low efficiency in reservoir characterization. Furthermore, well logging crossplots obtained in basin-level work areas often contain multiple strata (including reservoir regions). Reservoir characterization requires determining the internal and external relationships between all seismic points in the seismic data volume and each of these strata, further increasing the computational complexity and time consumption of reservoir characterization, reducing its efficiency. In practice, high computational load often results in slow machine performance.

[0084] The embodiment of the present application provides a reservoir characterization method and device. When characterizing the reservoir, first obtain a target well logging intersection diagram and a seismic data volume of the target area. The target well logging intersection diagram includes a target reservoir area and multiple well logging sample points. The seismic data volume includes multiple seismic sample points. Then, according to the target reservoir area, a target index network is set in the target well logging intersection diagram. The graph surrounded by the outer contour of the target index network covers the target reservoir area. The target index network includes multiple index grids and multiple grid nodes. Each grid node in the multiple grid nodes has an index value. Then, according to the target index network, the index value of each seismic sample point in the seismic data volume is determined to obtain an index data volume. Finally, the target reservoir is characterized according to the index data volume. Among them, the target index network includes a target index grid. The boundary line of the target reservoir area passes through the target index grid. When determining the index data volume according to the target index network, for each seismic sample point in the seismic data volume located in the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir area. Since the internal and external relationship judgment of the target reservoir area is only performed for the seismic sample points located in the target index grid during reservoir characterization, the calculation amount and time consumption of reservoir characterization can be reduced, thereby improving the efficiency of reservoir characterization.

[0085] The technical solution of the present application is introduced below. First, an embodiment of the reservoir characterization method of the present application is introduced.

[0086] Please refer to Figure 1 , which shows a flow chart of a reservoir characterization method provided by an embodiment of the present application. The reservoir characterization method is executed by a reservoir characterization device, which is a computer device or a functional module in the computer device. Figure 1 The method includes the following steps S101 to S105.

[0087] S101. Obtain a target well logging intersection diagram and a seismic data body of the target area, wherein the target area is the geographical area where the target reservoir is located. The target well logging intersection diagram includes a target reservoir area and multiple well logging sample points, each well logging sample point is used to characterize at least two types of well logging data of a position point in the target area, the target reservoir area is used to characterize the target reservoir, and the seismic data body includes multiple seismic sample points, each seismic sample point is used to characterize multiple types of seismic data of a position point in the target area, and the attributes of at least two types of seismic data among the multiple types of seismic data are the same as the attributes of the at least two types of well logging data.

[0088] The reservoir characterization device can first acquire at least two types of well logging data for a target area, then generate an initial well logging crossplot based on the at least two types of well logging data, and finally select a target reservoir area from the initial well logging crossplot to obtain a target well logging crossplot. The target reservoir area is a closed graphic area (i.e., an area defined by a closed graphic), and the shape of the target reservoir area can be polygonal or circular, for example. The at least two types of well logging data refer to well logging data of at least two attributes, including at least two of the following: compressional wave velocity, shear wave velocity, absorption coefficient, wave impedance, density, and wellbore diameter. In an exemplary embodiment, the at least two types of well logging data are u types of well logging data, each of which has u attributes. The reservoir characterization device first constructs a u-dimensional coordinate system, wherein the u coordinate axes of the u-dimensional coordinate system correspond to the u attributes. The reservoir characterization device then plots multiple well logging points in the u-dimensional coordinate system based on the u types of well logging data to obtain an initial well logging crossplot, where u is an integer greater than 1. The reservoir characterization device then selects a target reservoir region from the initial well logging crossplot based on the distribution of the well logging points in the initial well logging crossplot. In an exemplary embodiment, the well logging points in the initial well logging crossplot are clustered in at least one region, the at least one region corresponds one-to-one to at least one rock formation, the at least one rock formation includes the target reservoir, and the at least one region includes the target reservoir region. The reservoir characterization device selects the at least one region from the initial well logging crossplot. For example, the reservoir characterization device delineates the well logging points clustered in the same area in the initial well logging crossplot to obtain at least one closed graphical region. The closed graphical region within the at least one closed graphical region corresponding to the target reservoir is the target reservoir region. Optionally, to facilitate differentiation between different rock formations, the reservoir characterization device may use different colors to depict the well logging points clustered in different areas.

[0089] In an optional embodiment, the target area includes multiple rock layers, and the multiple rock layers include a target reservoir. During the drilling process in the target area, researchers can intuitively know the types of the multiple rock layers (for example, the types of rock layers include sand layers, mud layers, reservoir layers, etc.). The researchers select multiple location points in the target reservoir and perform logging at the multiple location points to obtain a variety of logging data for each reservoir location point in the multiple location points. The researchers select at least two types of logging data that have a greater impact on the target reservoir from the multiple logging data, and the reservoir characterization device generates an initial logging intersection diagram based on the at least two types of logging data. Optionally, the researchers also select location points in other rock layers for logging to obtain logging data for multiple location points in the other rock layers, and this embodiment of the application is not limited to this. Among them, different rock formations have different properties, and data with the same property varies differently in different rock formations. For example, the range of variation of the P-wave velocity in the reservoir is 2200m / s to 3600m / s, and the range of variation of the P-wave velocity in the sand layer is 4200m / s to 6000m / s. The logging data that has a greater impact on the target reservoir refers to the data that varies greatly in the target reservoir.

[0090] Take the above-mentioned at least two types of logging data as an example, wherein the attributes of the two types of logging data are longitudinal wave velocity and shear wave velocity. The initial logging cross-plot generated by the reservoir characterization device based on the two types of logging data can be as follows: Figure 2 As shown, the well logging cross-plot (i.e., the target well logging cross-plot) after the reservoir characterization device picks up the target reservoir area in the initial well logging cross-plot can be as follows: Figure 3 As shown. Figure 2 and Figure 3 In the well logging cross-plot shown in the figure, the horizontal axis represents the P-wave velocity (in Figure 2 and Figure 3 The vertical axis represents the shear wave velocity (in Figure 2 and Figure 3 S_VEL is used in the diagram), the units of the longitudinal wave velocity and the shear wave velocity are both meters per second (m / s), and the small black dots represent logging sample points. Both the initial logging intersection diagram and the target logging intersection diagram include multiple logging sample points, and the multiple logging sample points are clustered in the same area. Each of the multiple logging sample points is used to characterize the longitudinal wave velocity and shear wave velocity of a position point in the target reservoir of the target area. The area where the multiple logging sample points are clustered is the target reservoir area. The longitudinal wave velocity range corresponding to the target reservoir area is a~b, and the shear wave velocity range corresponding to the target reservoir area is c~d. Figure 3 As shown, the shape of the target reservoir area picked up in the initial well logging intersection map is a polygon. Figure 3 The shape of the target reservoir region shown is only an example. In other examples, the shape of the target reservoir region may be other shapes, such as a circle.

[0091] After the reservoir characterization device obtains the target well logging intersection map, the reservoir characterization device can also assign a value to the target reservoir area in the target well logging intersection map. For example, the reservoir characterization device picks up multiple closed graphic areas in the initial well logging intersection map to obtain the target well logging intersection map. The target reservoir area is one of the multiple closed graphic areas. The reservoir characterization device assigns values ​​to the multiple closed graphic areas and the areas other than the multiple closed graphics in the target well logging intersection map. The values ​​assigned by the reservoir characterization device to the multiple closed graphic areas are different, and the values ​​assigned by the reservoir characterization device to the multiple closed graphic areas are different from the values ​​assigned by the reservoir characterization device to the areas other than the multiple closed graphics in the target well logging intersection map. For example, Figure 4 This embodiment of the present application provides a method for Figure 3 The target reservoir area in the target well logging crossplot and the area other than the target reservoir area in the target well logging crossplot are respectively assigned values. The reservoir characterization device assigns a first index value to the target reservoir area and a second index value to the area other than the target reservoir area. Figure 4 Take the first index value as "1" and the second index value as "0" as an example. Figure 4 For the sake of simplicity, the well logging sample points in the target well logging cross-plot are not shown.

[0092] In an optional embodiment, a seismic test is conducted in the target area to obtain seismic data for the target area, and the reservoir characterization device obtains a seismic data volume for the target area based on the seismic data of the target area. Each seismic sample point in the seismic data volume is used to represent multiple types of seismic data at a location in the target area, and the attributes of at least two types of seismic data are the same as the attributes of the at least two types of well logging data represented by the well logging sample points in the target well logging crossplot. For example, if the two types of well logging data include compressional wave velocity and shear wave velocity, then the multiple types of seismic data include compressional wave velocity and shear wave velocity.

[0093] S102. A target index network is set in the target logging intersection diagram according to the target reservoir area. The figure enclosed by the outer contour of the target index network covers the target reservoir area. The target index network includes multiple index grids and multiple grid nodes. The multiple grid nodes are intersection points of grid lines constituting the target index network. Each of the multiple grid nodes has an index value.

[0094] The reservoir characterization device can set an initial index network in the target logging intersection diagram according to the target reservoir area, so that the figure surrounded by the outer contour of the initial index network covers the target reservoir area. The initial index network includes multiple index grids and multiple grid nodes. The multiple grid nodes are the intersection points of the grid lines constituting the initial index network; the reservoir characterization device sets an index value for each grid node in the multiple grid nodes according to the relationship between each grid node and the target reservoir area, thereby obtaining the target index network.

[0095] In an optional embodiment, the index grids in the target index network are square in shape. The reservoir characterization device first determines the grid width of the index grid; the reservoir characterization device determines the number of rows and columns of the target index network based on the grid width of the index grid and the target reservoir area; and the reservoir characterization device sets an initial index network in the target well logging crossplot based on the grid width of the index grid, the number of rows and columns of the target index network. For example, the reservoir characterization device determines the grid width of the index grid based on the area of ​​the target reservoir area and the target number, such that the number of logging points included in each index grid is not less than the target number. The reservoir characterization device determines the number of rows of the target index network as the ratio of the projected length of the target reservoir area on the vertical axis of the target well logging crossplot to the grid width, and determines the number of columns of the target index network as the ratio of the projected length of the target reservoir area on the horizontal axis of the target well logging crossplot to the grid width. For example, the projection length of the target reservoir area on the vertical axis of the target logging crossplot is the difference between the maximum and minimum values ​​of the vertical coordinate corresponding to the target reservoir area; the projection length of the target reservoir area on the horizontal axis of the target logging crossplot is the difference between the maximum and minimum values ​​of the horizontal coordinate corresponding to the target reservoir area.

[0096] As an example, see Figure 5 , which shows a method provided by an embodiment of the present application Figure 4 The schematic diagram of the target well logging intersection diagram after setting the initial index network is shown. Figure 5 As shown in FIG, the horizontal axis of the target well logging crossplot represents the longitudinal wave velocity, the vertical axis represents the shear wave velocity, the projection length of the target reservoir area on the vertical axis of the target well logging crossplot is dc, and the projection length of the target reservoir area on the horizontal axis of the target well logging crossplot is ba. Assuming that the reservoir characterization device determines that the grid width of the index grid is h, the reservoir characterization device determines that the number of rows of the initial index network is (dc) / h and the number of columns is (ba) / h. The reservoir characterization device determines the initial index network according to the grid width h, the number of rows of the initial index network (dc) / h and the number of columns of the target index grid (ba) / h. Figure 4 The target well logging intersection diagram shown in the figure is used to set the initial index network. Figure 5 The well logging crossplot is shown. Figure 5Taking an initial index network having 10 rows and 10 columns as an example, a graph formed by the outer contour of the initial index network covers the target reservoir area. The initial index network includes a plurality of index grids 041 and a plurality of grid nodes 043. The plurality of grid nodes 043 are intersection points of the grid lines 042 constituting the initial index network.

[0097] As previously described, the initial index network includes multiple index grids and multiple grid nodes. The reservoir characterization device can set an index value for each grid node in the multiple grid nodes based on its relationship with the target reservoir region. For example, for each grid node in the multiple grid nodes, when the grid node is located within the target reservoir region or on the boundary of the target reservoir region, the reservoir characterization device sets the index value of the grid node to a first index value; when the grid node is located outside the target reservoir region, the reservoir characterization device sets the index value of the grid node to a second index value. In actual practice, due to the large number of index grids in the initial index network, the reservoir characterization device can first determine the index values ​​of the grid nodes corresponding to the index grids that the boundary of the target reservoir region passes through (for ease of description, the index grids that the boundary of the target reservoir region passes through are referred to as target index grids). Then, the reservoir characterization device can determine the index values ​​of the grid nodes corresponding to non-target index grids (i.e., index grids in the initial index network other than the target index grid).

[0098] When determining the index value of the grid node corresponding to the target index grid, the reservoir characterization device can process each boundary line of the target reservoir area in the form of linear parameter equation stepping. For example, the first boundary line is a boundary line of the target reservoir area. For the first boundary line: first, the reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the first endpoint of the first boundary line (any endpoint of the first boundary line can be used as the first endpoint) is located, based on the internal and external relationship between the grid node corresponding to the target index grid and the target reservoir area; then, the reservoir characterization device takes the first endpoint as the starting point and uses the first linear equation (the first linear equation is the linear equation of the first boundary line) to step the first distance along the first boundary line to the first step point. The reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the first step point is located, based on the internal and external relationship between the grid node corresponding to the target index grid where the first step point is located and the target reservoir area; then, the reservoir characterization device uses the first endpoint as the starting point and the first linear equation (the first linear equation is the linear equation of the first boundary line) to step the first distance along the first boundary line to the first step point. The first step point is used as the starting point, and the first linear equation is used to continue stepping along the first boundary line for a second distance to the second step point. The reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the second step point is located based on the internal and external relationship between the grid node corresponding to the target index grid where the second step point is located and the target reservoir area. Then, the reservoir characterization device uses the second step point as the starting point and continues stepping along the first boundary line for a third distance to the third step point using the first linear equation. The reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the third step point is located based on the internal and external relationship between the grid node corresponding to the target index grid where the third step point is located and the target reservoir area. This process is repeated until the reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the second endpoint of the first boundary line is located. The first distance, the second distance, and the third distance may be equal or unequal. The first linear equation may include a horizontal axis linear equation and a vertical axis linear equation. The horizontal axis linear equation is used to express the horizontal axis data of the first boundary line, and the vertical axis linear equation is used to express the vertical axis data of the first boundary line. The reservoir characterization device performs the above operation process on each boundary line of the target reservoir area to determine the index value of the grid node corresponding to each target index grid.For example, the first endpoint, the second endpoint and the above-mentioned step points are all called operating points. For each operating point among the first endpoint, the second endpoint and the above-mentioned step points, the reservoir characterization device determines the index value of the grid node corresponding to the target index grid where each operating point is located, including: the reservoir characterization device draws out a characteristic ray with each grid node of the target index grid as an endpoint, and the characteristic ray has at least one intersection with the boundary line of the target reservoir area; when the reservoir characterization device determines that the number of intersections of the characteristic ray and the boundary line of the target reservoir area is an odd number, the reservoir characterization device determines that each grid node is located within the target reservoir area, and then the reservoir characterization device sets the index value of each grid node to the first index value; when the reservoir characterization device determines that the number of intersections of the characteristic ray and the boundary line of the target reservoir area is an even number, the reservoir characterization device determines that each grid node is located outside the target reservoir area, and then sets the index value of each grid node to the second index value.

[0099] For example, the embodiment of the present application takes the first boundary line of the target reservoir area as Figure 5 Take the line segment PQ in the figure as an example. Figure 6 As shown in FIG, it shows a schematic diagram of determining the index value of the grid node corresponding to the target index grid that the line segment PQ passes through according to an embodiment of the present application. The two endpoints of the line segment PQ are endpoint P and endpoint Q, and the coordinates of endpoint P are (X P ,Y P ), the coordinates of the endpoint Q are (X Q ,Y Q ), the thread equation of line segment PQ (i.e., the first linear equation) includes the horizontal axis linear equation and the vertical axis linear equation. The horizontal axis linear equation is X=X P +et, the linear equation of the vertical axis is Y=Y P +ft,e=X Q -X P , f=Y Q -Y P , X represents the horizontal coordinate of the step point, Y represents the vertical coordinate of the step point, X P Indicates the horizontal coordinate of the endpoint P, Y P represents the vertical coordinate of the endpoint P, X Q Indicates the horizontal coordinate of endpoint Q, Y Q represents the vertical coordinate of the endpoint Q, and t represents the step parameter. Figure 6As shown, the endpoint P is taken as the first endpoint for explanation. First, the reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the endpoint P is located according to the internal and external relationship between the grid node corresponding to the target index grid where the endpoint P is located and the target reservoir area; then, the reservoir characterization device takes the endpoint P as the starting point, and uses the above-mentioned first linear equation to step a first distance along the first boundary line PQ to the first step point A1, the coordinates of the first step point are (X1, Y1), and the reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the first step point A1 is located according to the internal and external relationship between the grid node corresponding to the target index grid where the first step point A1 is located and the target reservoir area; then, the reservoir characterization device takes the first step point A1 as the starting point, and uses the above-mentioned first linear equation to continue stepping a second distance along the first boundary line PQ to the second step point A2, the second step point The coordinates of step point A2 are (X2, Y2). The reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the second step point A2 is located based on the internal and external relationship between the grid node corresponding to the target index grid where the second step point A2 is located and the target reservoir area; thereafter, the reservoir characterization device takes the second step point A2 as the starting point and uses the above-mentioned first linear equation to continue stepping a third distance along the first boundary line PQ to the third step point A3. The coordinates of the third step point A3 are (X3, Y3). The reservoir characterization device determines the index value of the grid node corresponding to the target index grid where the third step point A3 is located based on the internal and external relationship between the grid node corresponding to the target index grid where the third step point A3 is located and the target reservoir area; and so on, until the reservoir characterization device determines the index value of the grid node of the index grid where the endpoint Q is located.

[0100] When determining the index value of the grid node corresponding to the non-target index grid, the reservoir characterization device scans the grid lines of the initial index network row by row or column by column to determine the index value of the grid node on each grid line. For each grid line in the initial index network: when there is a grid node with an index value set on each grid line, the reservoir characterization device sets an index value for other grid nodes on each grid line according to the index value of the grid node with the first index value set; when there is no grid node with the first index value set on each grid line, the reservoir characterization device sets the index value of the grid node on each grid line to the second index value. Among them, when there is a grid node with a first index value set on each grid line, the reservoir characterization device sets index values ​​for other grid nodes on each grid line according to the index value of the grid node with an index value set, including: when the number of grid nodes with the first index value set on each grid line is 1, the reservoir characterization device sets the index values ​​of other grid nodes on each grid line to the second index value; when the number of grid nodes with the first index value set on each grid line is 2, the reservoir characterization device sets the index values ​​of the grid nodes on each grid line located between the two grid nodes (i.e., the two grid nodes with the first index value set) between the two grid nodes (i.e., the two grid nodes with the first index value set) to the first index value, and sets the index values ​​of other grid nodes on each grid line to the second index value.

[0101] As an example, see Figure 7 , which shows a method provided by an embodiment of the present application Figure 5 The schematic diagram of the initial index network after the network nodes are indexed (i.e., the schematic diagram of the target index network) is shown. This embodiment of the application uses the first index value of "1" and the second index value of "0" as an example. When a grid node is located exactly on the boundary line of the target reservoir area, the index value of the grid node is set to "1."

[0102] S103. Determine the index value of each seismic sample point in the seismic data volume according to the target index network, and obtain an index data volume corresponding to the seismic data volume, wherein, for each seismic sample point in the seismic data volume located within the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir area, the multiple index grids include the target index grid, and the boundary line of the target reservoir area passes through the target index grid.

[0103] Optionally, for each seismic sample point in the seismic data volume that is within a target index grid, the reservoir characterization device determines an index value for each seismic sample point based on the internal and external relationship between each seismic sample point and the target reservoir region; and for each seismic sample point in the seismic data volume that is within a non-target index grid, the reservoir characterization device determines an index value for each seismic sample point based on the index value of the grid node corresponding to the non-target index grid. The multiple index grids of the target index network include target index grids and non-target index grids, the boundary line of the target reservoir region passes through the target index grid, and the non-target index grids are index grids in the target index network other than the target index grid.

[0104] In an optional embodiment, the non-target index grid includes a first non-target index grid and a second non-target index grid, wherein the first non-target index grid is located within the target reservoir region, and the second non-target index grid is located outside the target reservoir region. The index values ​​of the grid nodes corresponding to the first non-target index grid are all first index values, and the index values ​​of the grid nodes corresponding to the second non-target index grid are all second index values. For each seismic sample point in the seismic data volume located within the first non-target index grid, the reservoir characterization device determines the index value of each seismic sample point to be the first index value. For each seismic sample point in the seismic data volume located within the second non-target index grid, the reservoir characterization device determines the index value of each seismic sample point to be the second index value. For example, the orthographic projection of the seismic data body in the target well logging interaction diagram coincides with the figure enclosed by the outer contour of the target index network. The reservoir characterization device maintains a first corresponding relationship, which is the correspondence between the coordinates of the seismic sample points and the seismic data. The coordinates of each seismic sample point are the coordinates of the seismic sample point in the seismic data body. Each seismic sample point corresponds to multiple seismic data. The multiple seismic data corresponding to each seismic sample point are multiple seismic data of the position point represented by each seismic sample point in the seismic data body. The reservoir characterization device determines the position of each seismic sample point in the target index network based on the first corresponding relationship, that is, determines which index grid of the target index network each seismic sample point is in. For example, the multiple seismic data corresponding to each seismic sample point include P-wave velocity and P-wave velocity, the reservoir characterization device determines the difference between the P-wave velocity corresponding to each seismic sample point and the minimum P-wave velocities corresponding to all seismic sample points in the seismic data body (referred to as the P-wave velocity difference), and the reservoir characterization device determines the ratio of the P-wave velocity difference to the grid width of the index grid of the target index network as the horizontal coordinate of each seismic sample point in the target index network; and the reservoir characterization device determines the difference between the S-wave velocity corresponding to each seismic sample point and the minimum S-wave velocities corresponding to all seismic sample points in the seismic data body (referred to as the S-wave velocity difference), and the reservoir characterization device determines the ratio of the S-wave velocity difference to the grid width of the index grid of the target index network as the vertical coordinate of each seismic sample point in the target index network.

[0105] For example, the first correspondence is shown in Table 1. Each seismic sample point corresponds to a longitudinal wave velocity and a transverse wave velocity. It is assumed that among the m×n×r seismic sample points, the minimum longitudinal wave velocity corresponding to all seismic sample points is S0, and the minimum transverse wave velocity corresponding to the seismic sample point is W0. For a seismic sample point with coordinates (x1, y1, z1), the horizontal coordinate of the seismic sample point in the target index network is (S1-S0) / h, and the vertical coordinate of the seismic sample point in the target index network is (W1-W0) / h; for a seismic sample point with coordinates (x1, y1, z2), the horizontal coordinate of the seismic sample point in the target index network is (S2-S0) / h, and the vertical coordinate of the seismic sample point in the target index network is (W2-W0) / h; and so on. m ,y n ,z r ) of the earthquake sample point, the horizontal coordinate of the earthquake sample point in the target index network is (S3-S0) / h, and the vertical coordinate of the earthquake sample point in the target index network is (W3-W0) / h.

[0106] Table 1

[0107] Earthquake sample point coordinates P-wave velocity data Shear wave velocity data <![CDATA[(x1,y1,z1)]]> S1 W1 <![CDATA[(x1,y1,z2)]]> S2 W2 ... ... ... <![CDATA[(x m ,y n ,z r )]]> S3 W3

[0108] Optionally, the reservoir characterization device maintains a second correspondence between seismic sample point coordinates and seismic sample point index values. The reservoir characterization device determines the index value of each seismic sample point based on its position in the target index network, thereby obtaining the second correspondence. The reservoir characterization device generates a feature data volume based on the second correspondence.

[0109] In an optional embodiment, the reservoir characterization device calculates the position of each seismic sample point in the target index network. For each seismic sample point in the first non-target index grid in the target index network in the seismic data volume, the index value of the index grid where each seismic sample point is located is the first index value, and the reservoir characterization device determines the index value of each seismic sample point to be the first index value; for each seismic sample point in the second non-target index grid in the target index network in the seismic data volume, the index value of the index grid where each seismic sample point is located is the second index value, and the reservoir characterization device determines the index value of each seismic sample point to be the second index value.

[0110] Optionally, the reservoir characterization device calculates the position of each seismic sample point in the target index network. For each seismic sample point in the seismic data volume that is located within the first target index grid, when the seismic sample point is located within the boundary line of the target reservoir area, the reservoir characterization device determines the index value of the seismic sample point to be the first index value; when the seismic sample point is located outside the boundary line of the target reservoir area, the reservoir characterization device determines the index value of the seismic sample point to be the first index value. In an optional embodiment, for each seismic sample point in the seismic data volume that is located within the first target index grid, a characteristic ray is drawn with the seismic sample point as an endpoint. The characteristic ray has at least one intersection with the boundary line of the target reservoir area. When the number of the intersections is an odd number, the reservoir characterization device determines that the seismic sample point is located within the target reservoir area. When the number of the intersections is an even number, the reservoir characterization device determines that the seismic sample point is located outside the target reservoir area. For an example, please refer to Figure 8 , which shows a schematic diagram of determining the internal and external relationship between a seismic sample point and a target reservoir area provided in the implementation of this application. Figure 8 As shown, point J and point G are the position points of two seismic sample points among the multiple seismic sample points in the target index network. A characteristic ray is drawn with point J as the endpoint. The characteristic ray passes through the target reservoir area and has one intersection with the target reservoir area. The number of the intersections is an odd number. The reservoir characterization device determines that the seismic sample point J is located in the target reservoir area; a characteristic ray is drawn with point G as the endpoint. The characteristic ray passes through the target reservoir area and has two intersections with the target reservoir area. The number of the intersections is an even number. The reservoir characterization device determines that the seismic sample point G is located outside the target reservoir area.

[0111] As an example, see Figure 9 , which shows a schematic diagram of determining the index value of a seismic sample point provided by an embodiment of the present application, where point K, point J, point G and point I are the positions of four seismic sample points in the target index network among the multiple seismic sample points, and the positions of the four points in the target index grid encompass the above-mentioned multiple situations. Figure 9As shown, point J is located in the first non-target index grid, and the index values ​​of the grid nodes of the index grid where point J is located are all "1", so the reservoir characterization device determines the index value of point J to be "1"; point G is located in the second non-target index grid, and the index values ​​of the grid nodes of the index grid where point G is located are all "0", so the reservoir characterization device determines the index value of point G to be "0"; point K is located in the first target index grid, and the index values ​​of the grid nodes of the index grid where point K is located are "1", "0", "0" and "0", so according to the internal and external relationship between point K and the polygon, it is judged that point K is located outside the polygon, and the reservoir characterization device determines the index value of point K to be "0"; point I is located in the first target index grid, and the index values ​​of the grid nodes of the index grid where point I is located are "1", "1", "1" and "0", so according to the internal and external relationship between point I and the polygon, it is judged that point I is located inside the polygon, and the reservoir characterization device determines the index value of point K to be "1".

[0112] Optionally, the reservoir characterization device performs the above steps on each seismic sample point in the seismic data body, determines the index value corresponding to each seismic sample point in the seismic data body, and obtains a second corresponding relationship. For example, the second corresponding relationship is shown in Table 2. The reservoir characterization device represents each seismic sample point in the seismic data body with the index value based on the second corresponding relationship, thereby obtaining the seismic data body.

[0113] Table 2

[0114] Earthquake sample point coordinates Index value <![CDATA[(x1,y1,z1)]]> 0 <![CDATA[(x1,y1,z2)]]> 1 ... ... <![CDATA[(x m ,y n ,z r )]]> 0

[0115] S104. Characterize the target reservoir according to the index data volume.

[0116] The reservoir characterization device displays all seismic sample points with the first index value in the index data body as one color, and displays all seismic sample points with the second index value in the index data body as another color. The area formed by all seismic sample points with the first index value is the target reservoir. For an example, please refer to Figure 10 , which shows a cross-sectional view of a target reservoir provided in an embodiment of the present application. Taking the target reservoir as a coal seam as an example, the horizontal axis is the plane position of the target reservoir cross-section within the target area, and the vertical axis is the time axis. The reservoir characterization device sets the seismic sample points with an index value of "1" to black, and sets the seismic sample points with an index value of "0" to white. The black part is the target reservoir obtained by characterization. After the index data body is color-set, the cross-sectional structure of the target reservoir is obtained at any position. Figure 10 A cross-sectional view of the target reservoir is shown.

[0117] In summary, the reservoir characterization method provided by the embodiment of the present application. When characterizing the reservoir, first obtain the target well logging intersection diagram of the target area and the seismic data body of the target area, the target well logging intersection diagram includes the target reservoir area and multiple well logging sample points, and the seismic data body includes multiple seismic sample points. Then, according to the target reservoir area, a target index network is set in the target well logging intersection diagram. The graph surrounded by the outer contour of the target index network covers the target reservoir area. The target index network includes multiple index grids and multiple grid nodes. Each grid node in the multiple grid nodes has an index value. Then, according to the target index network, the index value of each seismic sample point in the seismic data body is determined to obtain an index data body, and finally, the target reservoir is characterized according to the index data body. Among them, the target index network includes a target index grid, and the boundary line of the target reservoir area passes through the target index grid. When determining the index data body according to the target index network, for each seismic sample point in the seismic data body located in the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir area. Since the internal and external relationship judgment of the target reservoir area is only performed for the seismic sample points located in the target index grid during reservoir characterization, the calculation amount and time consumption of reservoir characterization can be reduced, thereby improving the efficiency of reservoir characterization.

[0118] The above is an introduction to the method embodiments of the present application. The following describes the device embodiments of the present application, which can be used to perform the method of the present application. For details not disclosed in the device embodiments, please refer to the method embodiments.

[0119] Please refer to Figure 11 , which shows a schematic diagram of a reservoir characterization device 100 provided in an embodiment of the present application, the reservoir characterization device 100 can be used to perform Figure 1 The embodiment shown provides a reservoir characterization method. Figure 11 The reservoir characterization device 100 includes: an acquisition module 101 , a setting module 102 , an index determination module 103 and a characterization module 104 .

[0120] An acquisition module 101 is configured to acquire a target well logging crossplot of a target region and a seismic data volume of the target region, wherein the target region is a geographical region where a target reservoir is located, the target well logging crossplot including a target reservoir region and a plurality of well logging points, each well logging point being used to represent at least two types of well logging data at a location point in the target region, the target reservoir region being used to represent the target reservoir, the seismic data volume including a plurality of seismic points, each seismic point being used to represent multiple types of seismic data at a location point in the target region, and attributes of at least two types of seismic data being the same as attributes of the at least two types of well logging data;

[0121] A setting module 102 is configured to set a target index network in the target logging crossplot according to the target reservoir region, wherein the target index network is surrounded by an outer contour that covers the target reservoir region, and the target index network includes a plurality of index grids and a plurality of grid nodes, wherein the plurality of grid nodes are intersections of grid lines constituting the target index network, and each of the plurality of grid nodes has an index value;

[0122] An index determination module 103 is configured to determine an index value of each seismic sample point in the seismic data volume based on the target index grid, thereby obtaining an index data volume corresponding to the seismic data volume. For a seismic sample point in the seismic data volume that is within a target index grid, the index value of the seismic sample point is determined based on an internal / external relationship between the seismic sample point and the target reservoir region, wherein the multiple index grids include the target index grid, and the boundary line of the target reservoir region passes through the target index grid.

[0123] The characterization module 104 is configured to characterize the target reservoir according to the index data volume.

[0124] Optionally, the index determination module 103 is used to determine the index value of each seismic sample point in the seismic data body located in a non-target index grid according to the index value of the grid node corresponding to the non-target index grid, the multiple index grids include the non-target index grid, and the grid nodes corresponding to the non-target index grid are the intersection points of the grid lines constituting the non-target index grid.

[0125] Optionally, the non-target index grid includes a first non-target index grid and a second non-target index grid, the first non-target index grid is located in the target reservoir area, and the second non-target index grid is located outside the target reservoir area, the index values ​​of the grid nodes corresponding to the first non-target index grid are all first index values, and the index values ​​of the grid nodes corresponding to the second non-target index grid are all second index values, and the index determination module 103 is used to: for each seismic sample point in the seismic data body located in the first non-target index grid, determine that the index value of each seismic sample point is the first index value; for each seismic sample point in the seismic data body located in the second non-target index grid, determine that the index value of each seismic sample point is the second index value.

[0126] Optionally, the index values ​​of the grid nodes within the target reservoir area are all first index values, and the index values ​​of the grid nodes outside the target reservoir area are all second index values. The index determination module 103 is used to: for each seismic sample point in the seismic data volume that is located within the target index grid: when each seismic sample point is located within the target reservoir area, determine the index value of each seismic sample point to be the first index value; when each seismic sample point is located outside the target reservoir area, determine the index value of each seismic sample point to be the second index value.

[0127] Optionally, the reservoir characterization device 100 also includes: a relationship determination module 105, which is used to: for each seismic sample point in the seismic data body located in the target index grid: draw a characteristic ray with each seismic sample point as an endpoint, and the characteristic ray has at least one intersection with the boundary line of the target reservoir area; when the number of intersections between the characteristic ray and the boundary line of the target reservoir area is an odd number, determine that each seismic sample point is located within the target reservoir area; when the number of intersections between the characteristic ray and the boundary line of the target reservoir area is an even number, determine that each seismic sample point is located outside the target reservoir area.

[0128] Optionally, a setting module 102 is used to: set an initial index network in the target logging intersection diagram according to the target reservoir area, so that the figure surrounded by the outer contour of the initial index network covers the target reservoir area, and the initial index network includes multiple index grids and multiple grid nodes; according to the relationship between each grid node in the multiple grid nodes and the target reservoir area, set an index value for each grid node to obtain the target index network.

[0129] Optionally, the setting module 102 is used to: for each grid node among the multiple grid nodes: when each grid node is located within the target reservoir area or on the boundary line of the target reservoir area, set the index value of each grid node to a first index value; when each grid node is located outside the target reservoir area, set the index value of each grid node to a second index value.

[0130] Optionally, the reservoir characterization device 100 further includes: a region determination module 106, configured to determine the target reservoir region in the target logging crossmap before the setting module 102 sets the target index network in the target logging crossmap according to the target reservoir region.

[0131] In summary, the reservoir characterization device provided in the embodiment of the present application first obtains a target well logging intersection diagram and a seismic data volume of the target area, wherein the target well logging intersection diagram includes a target reservoir area and multiple well logging sample points, and the seismic data volume includes multiple seismic sample points. Then, the reservoir characterization device sets a target index network in the target well logging intersection diagram according to the target reservoir area, wherein the outer contour of the target index network forms a graph covering the target reservoir area, and the target index network includes multiple index grids and multiple grid nodes, wherein each grid node in the multiple grid nodes has an index value. Then, the reservoir characterization device determines the index value of each seismic sample point in the seismic data volume according to the target index network to obtain an index data volume, and finally, the reservoir characterization device characterizes the target reservoir according to the index data volume. In this case, the target index network includes a target index grid, and the boundary line of the target reservoir area passes through the target index grid. When determining the index data volume according to the target index network, for each seismic sample point in the seismic data volume located within the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir area. Since the internal and external relationship judgment of the target reservoir area is only performed for the seismic sample points located in the target index grid during reservoir characterization, the calculation amount and time consumption of reservoir characterization can be reduced, thereby improving the efficiency of reservoir characterization.

[0132] The embodiment of the present application provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor is configured to execute the computer program to implement the following Figure 1 The method shown.

[0133] The present invention provides a computer-readable storage medium in which a computer program is stored. When the computer program is executed, the following is achieved: Figure 1 The method shown.

[0134] The embodiment of the present application provides a computer program product, which includes a program or code, and when the program or code is executed, realizes the following Figure 1 The method shown.

[0135] It should be understood that, for the sake of clarity, in this application, terms such as "first," "second," and "third" are used to distinguish identical or similar items with substantially the same functions and effects. Those skilled in the art will understand that terms such as "first," "second," and "third" do not limit the quantity or order of execution.

[0136] Those skilled in the art will understand that all or part of the steps to implement the above embodiments may be accomplished by hardware, or by a program to instruct the relevant hardware, and the program may be stored in a computer-readable storage medium, which may be a read-only memory, a disk, or an optical disk, etc.

[0137] The above description is merely an exemplary embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A reservoir characterization method, characterized in that: The method comprises: Acquiring a target well logging crossplot of a target area and a seismic data volume of the target area, wherein the target area is a geographical area where a target reservoir is located, the target well logging crossplot includes a target reservoir area and a plurality of well logging points, each well logging point being used to represent at least two types of well logging data for a location point in the target area, the target reservoir area being used to represent the target reservoir, the seismic data volume including a plurality of seismic points, each seismic point being used to represent multiple types of seismic data for a location point in the target area, and attributes of at least two types of seismic data being the same as attributes of the at least two types of well logging data; A target index network is set in the target logging crossplot according to the target reservoir region, wherein a graph formed by an outer contour of the target index network covers the target reservoir region, the target index network includes a plurality of index grids and a plurality of grid nodes, the plurality of grid nodes are intersections of grid lines constituting the target index network, and each of the plurality of grid nodes has an index value; determining an index value of each seismic sample point in the seismic data volume according to the target index grid, and obtaining an index data volume corresponding to the seismic data volume, wherein, for each seismic sample point in the seismic data volume located within a target index grid, the index value of each seismic sample point is determined according to an internal and external relationship between each seismic sample point and the target reservoir region, the multiple index grids include the target index grid, and a boundary line of the target reservoir region passes through the target index grid; The target reservoir is characterized according to the index data volume.

2. The method according to claim 1, characterized in that Determining the index value of each seismic sample point in the seismic data volume according to the target index network includes: For each seismic sample point in the seismic data volume that is located within a non-target index grid, an index value of each seismic sample point is determined based on an index value of a grid node corresponding to the non-target index grid, the multiple index grids including the non-target index grid, and the grid nodes corresponding to the non-target index grid are intersections of grid lines constituting the non-target index grid.

3. The method according to claim 2, characterized in that The non-target index grid includes a first non-target index grid and a second non-target index grid, wherein the first non-target index grid is located within the target reservoir region, and the second non-target index grid is located outside the target reservoir region, and the index values ​​of the grid nodes corresponding to the first non-target index grid are all first index values, and the index values ​​of the grid nodes corresponding to the second non-target index grid are all second index values. The method of determining, for each seismic sample point in the seismic data volume and located in a non-target index grid, each seismic sample point according to an index value of a grid node corresponding to the non-target index grid, includes: For each seismic sample point in the seismic data volume and located within the first non-target index grid, determining an index value of each seismic sample point as the first index value; For each seismic sample point in the seismic data volume that is located within the second non-target index grid, an index value of each seismic sample point is determined to be the second index value.

4. The method according to claim 1, wherein The index values ​​of the grid nodes within the target reservoir region are all first index values, and the index values ​​of the grid nodes outside the target reservoir region are all second index values. For each seismic sample point in the seismic data volume located within the target index grid, the index value of each seismic sample point is determined according to the internal and external relationship between each seismic sample point and the target reservoir region, including: For each seismic sample point in the seismic data volume and located within the target index grid: When each of the seismic sample points is located in the target reservoir area, determining the index value of each of the seismic sample points to be the first index value; When each seismic sample point is located outside the target reservoir area, the index value of each seismic sample point is determined to be the second index value.

5. The method according to claim 4, characterized in that The method further comprises: For each seismic sample point in the seismic data volume and located within the target index grid: Drawing a characteristic ray with each seismic sample point as an endpoint, wherein the characteristic ray has at least one intersection point with a boundary line of the target reservoir area; When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is an odd number, determining that each seismic sample point is located in the target reservoir region; When the number of intersections between the characteristic ray and the boundary line of the target reservoir region is even, it is determined that each seismic sample point is located outside the target reservoir region.

6. The method according to claim 1, characterized in that Setting a target index network in the target well logging crossplot according to the target reservoir area includes: Setting an initial index network in the target logging crossplot according to the target reservoir region, so that a graph enclosed by an outer contour of the initial index network covers the target reservoir region, the initial index network including the plurality of index grids and the plurality of grid nodes; According to the relationship between each grid node in the plurality of grid nodes and the target reservoir region, an index value is set for each grid node to obtain the target index network.

7. The method according to claim 6, characterized in that The step of setting an index value for each grid node among the plurality of grid nodes according to a relationship between the grid node and the target reservoir region includes: For each grid node of the plurality of grid nodes: When each grid node is located within the target reservoir region or on a boundary line of the target reservoir region, setting the index value of each grid node to a first index value; When each grid node is located outside the target reservoir area, the index value of each grid node is set to a second index value.

8. The method according to any one of claims 1 to 7, characterized in that Before setting a target index network in the target well logging crossmap according to the target reservoir region, the method further includes: The target reservoir area in the target well logging cross-plot is determined.

9. A reservoir characterization device, characterized in that: The method comprises modules for executing the method according to any one of claims 1 to 8.

10. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and the processor is configured to execute the computer program to implement the method according to any one of claims 1 to 8.

Images