Method and device for dividing high-permeability layer of fractured pore type reservoir

Abstract

The invention provides a method and device for dividing a high-permeability layer of a fractured pore type reservoir, and the method comprises the steps: obtaining a full-well-section imaging logging image and logging interpretation data of a target exploration region of the fractured pore type reservoir; obtaining fracture permeability curve data according to the whole well section imaging logging image; acquiring matrix permeability curve data from the well logging interpretation data; adding the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the stratum of the whole well section; determining a permeability threshold value of a high-permeability layer according to the total permeability curve of the whole well section stratum; and according to the permeability threshold value of the high-permeability layer, performing high-permeability layer division on the single well of the target exploration area to obtain a high-permeability layer division result. According to the scheme, the technical problem that the high-permeability layer distribution of the fractured pore type reservoir cannot be effectively determined in the prior art is solved, and the technical effect of accurately and efficiently determining the high-permeability layer distribution of the fractured pore type reservoir is achieved.

Description

Technical Field

[0001] The present application belongs to the field of geological exploration technology, and in particular relates to a method and device for dividing a high permeability layer of a fracture pore type reservoir. Background Art

[0002] High permeability layer refers to the relatively high permeability layer section that is locally developed in the overall low permeability reservoir. The formation of high permeability layer is related to factors such as sedimentation and diagenesis. Clarifying the distribution of high permeability layer can give full play to the seepage capacity of reservoir and increase oil and gas production. Otherwise, improper development measures may lead to further enhancement of heterogeneity of oil and gas reservoirs, thus causing adverse consequences such as sand production and water breakthrough in oil and gas wells.

[0003] According to the different types of reservoir space, oil and gas reservoirs can be divided into: fractured reservoirs, porous reservoirs and fracture-pore reservoirs. Among them, fractured reservoirs refer to reservoirs with fractures as the main reservoir space and seepage channels. The fracture development zone is the high permeability layer. Usually, as long as the effective fracture development zone is determined by logging, seismic and other data, high-yield oil and gas flow can be obtained; porous reservoirs refer to reservoirs with pores as the main reservoir space and throats as the main seepage channels. Fractures are not developed or are developed sporadically. Through core physical property testing and logging data interpretation, combined with production dynamic data, high-permeability layers can be determined, and high oil and gas yields can be obtained; in fracture-pore reservoirs, pores are usually the main oil and gas reservoir space, and throats and fractures can be used as oil and gas seepage channels, making the formation permeability distribution more complex.

[0004] There is currently no effective solution for how to accurately determine the distribution of high permeability layers in fracture-pore reservoirs. Summary of the invention

[0005] The purpose of the present application is to provide a method and device for dividing the high permeability layer of a fractured pore type reservoir, which can accurately and efficiently determine the high permeability layer of the fractured pore type reservoir, thereby effectively carrying out oil and gas reservoir exploitation.

[0006] The present application provides a method and device for dividing a high permeability layer of a fracture pore type reservoir, which is achieved by:

[0007] A method for dividing a high permeability layer of a fracture pore type reservoir, the method comprising:

[0008] Obtain full-well imaging logging images and logging interpretation data in the target exploration area;

[0009] Obtaining fracture permeability curve data according to the full-well imaging logging image;

[0010] Acquiring matrix permeability curve data from the well logging interpretation data;

[0011] Adding the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section;

[0012] Determining a permeability threshold value of a high permeability layer according to a total permeability curve of the entire well section formation;

[0013] According to the permeability threshold value of the high permeability layer, the high permeability layer is divided for the single well in the target exploration area to obtain a high permeability layer division result.

[0014] In one embodiment, the fracture permeability curve data is obtained according to the full-well imaging logging image, including:

[0015] identifying a sinusoidal curve from the full-well imaging logging image;

[0016] Taking the dark sine curve in the well section imaging logging image as a valid fracture, and taking the bright sine curve in the well section imaging logging image as an invalid fracture, so as to obtain a fracture identification result;

[0017] According to the crack identification results, the crack line density and effective opening of each crack are determined;

[0018] The permeability of each fracture is calculated based on the fracture line density and effective aperture of each fracture;

[0019] According to the permeability of each fracture, fracture permeability curve data is generated.

[0020] In one embodiment, the permeability of each fracture is calculated based on the fracture line density and effective aperture of each fracture, including:

[0021] The permeability of each crack is calculated according to the following formula:

[0022] K=b 3 D / 12×10 15

[0023] Among them, K is the permeability of the fracture, b is the effective opening, and D is the fracture line density.

[0024] In one embodiment, determining the permeability threshold value of the high permeability layer according to the total permeability curve of the formation in the entire well section includes:

[0025] According to the total permeability curve of the formation in the whole well section, establish the intersection diagram of the total permeability of the formation section where drilling fluid loss occurs and the amount of drilling fluid loss, as well as the intersection diagram of the average total permeability of the formation of the target layer of a single well and the open flow rate;

[0026] Draw a trend line based on the cross-plot of the total formation permeability and the amount of drilling fluid loss in the layer where drilling fluid loss occurs, and the cross-plot of the average total formation permeability and the open flow rate of the target layer of a single well;

[0027] The total formation permeability corresponding to the point where the slope of the trend line undergoes an obvious mutation is used as the permeability threshold value of the high permeability layer.

[0028] In one embodiment, dividing the high permeability layer of a single well in the target exploration area according to the permeability threshold value of the high permeability layer includes:

[0029] The layer section whose total permeability in each single well of the target exploration area is higher than the permeability threshold value of the high permeability layer is regarded as the high permeability layer;

[0030] According to the distribution of high permeability layers in each single well, the thickness of high permeability layers in the whole well section and each layer section of each single well and the proportion of high permeability layers in the total formation thickness are obtained statistically.

[0031] In one embodiment, after dividing the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the method further includes:

[0032] Obtaining a high permeability layer thickness threshold value of water wells in the target exploration area;

[0033] According to the high permeability layer thickness threshold value of the water-breaking wells and the high permeability layer thickness corresponding to each well in the target exploration area, the wells with water-breaking risk are determined and early warning information is generated.

[0034] In one embodiment, after dividing the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the method further includes:

[0035] According to the single well core phase marks, lithology and lithofacies combinations, imaging logging images, and logging curve characteristics of the target exploration area, each single well in the target exploration area is divided into sedimentary phases, and a planar distribution map of the sedimentary phases is formed;

[0036] Taking the thickness of the high permeability layer of each single well as hard data and the planar distribution of the sedimentary phase as the qualitative constraint condition, the contour map of the high permeability layer thickness in the target exploration period was drawn to obtain the distribution law of the high permeability layer thickness on the plane.

[0037] A device for dividing a high permeability layer of a fracture pore type reservoir, comprising:

[0038] The first acquisition module is used to acquire full-well imaging logging images and logging interpretation data of the target exploration area;

[0039] A generating module, used for obtaining fracture permeability curve data according to the full-well imaging logging image;

[0040] A second acquisition module is used to acquire matrix permeability curve data from the well logging interpretation data;

[0041] A processing module, used for adding the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section formation;

[0042] A determination module, used for determining a permeability threshold value of a high permeability layer according to a total permeability curve of the formation in the entire well section;

[0043] The division module is used to divide the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer to obtain the high permeability layer division result.

[0044] An electronic device comprises a processor and a memory for storing instructions executable by the processor, wherein the steps of the above method are implemented when the processor executes the instructions.

[0045] A computer-readable storage medium stores a computer program / instruction, which implements the steps of the above method when executed by a processor.

[0046] The method for dividing the high permeability layer of the fractured pore type reservoir provided in the present application obtains the fracture permeability curve data through the full-well imaging logging image of the target exploration area, obtains the matrix permeability curve data through the logging interpretation data, and then adds the fracture permeability curve data and the matrix permeability curve data to obtain the total permeability curve of the formation of the whole well section; based on the total permeability curve, the permeability threshold value of the high permeability layer is determined, so that the high permeability layer can be accurately determined. The above scheme solves the existing technical problem of being unable to effectively determine the distribution of the high permeability layer of the fractured pore type reservoir, and achieves the technical effect of accurately and efficiently determining the distribution of the high permeability layer of the fractured pore type reservoir, thereby increasing the production of oil and gas reservoirs. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0048] Figure 1 This is a method flow chart of an embodiment of a method for dividing a high permeability layer of a fracture pore type reservoir provided by the present application;

[0049] Figure 2 This is a schematic diagram of the natural fracture identification results of the Paleogene 5279.5-5283.0m of the DN2-20 well provided in this application;

[0050] Figure 3This is a schematic diagram of the Paleogene permeability curve, sedimentary phase and high permeability layer identification and division results of the DN2-20 well provided in this application;

[0051] Figure 4 It is a cross-plot of total formation permeability and drilling fluid loss provided by the present application;

[0052] Figure 5 It is a cross-plot of total formation permeability and open flow provided by this application;

[0053] Figure 6 This is the planar distribution diagram of the thickness of the high permeability layer of the Paleogene in the Dina 2 gas field provided by this application;

[0054] Figure 7 This is a hardware structure block diagram of an electronic device for a method for dividing a high permeability layer of a fracture pore type reservoir provided in the present application;

[0055] Figure 8 It is a schematic diagram of the module structure of an embodiment of a device for dividing a high permeability layer of a fracture pore type reservoir provided in the present application. DETAILED DESCRIPTION

[0056] In order to enable those skilled in the art to better understand the technical solutions in this application, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part 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 creative work should fall within the scope of protection of this application.

[0057] Figure 1 It is a method flow chart of an embodiment of the method for dividing the high permeability layer of the fracture pore type reservoir provided by the present application. Although the present application provides the method operation steps or device structure as shown in the following embodiments or drawings, more or fewer operation steps or module units may be included in the method or device based on routine or no creative labor. In the steps or structures where there is no necessary causal relationship logically, the execution order of these steps or the module structure of the device is not limited to the execution order or module structure described in the embodiments of the present application and shown in the drawings. When the method or module structure is applied in an actual device or terminal product, it can be connected according to the method or module structure shown in the embodiments or drawings for sequential execution or parallel execution (for example, a parallel processor or a multi-threaded processing environment, or even a distributed processing environment).

[0058] Specifically, Figure 1 As shown, the method for dividing the high permeability layer of the above-mentioned fracture pore type reservoir may include the following steps:

[0059] Step 101: Acquire full-well imaging logging images and logging interpretation data of the target exploration area;

[0060] Step 102: obtaining fracture permeability curve data according to the full-well imaging logging image;

[0061] Considering that natural fractures appear as sinusoidal curves on imaging logging images, dark sinusoidal curves are high-conductivity fractures, i.e., effective fractures, and bright sinusoidal curves are high-resistance fractures, i.e., invalid fractures. Therefore, when obtaining fracture permeability curve data based on the full-well imaging logging image, the sinusoidal curves can be identified from the full-well imaging logging image; the dark sinusoidal curves in the well imaging logging image are used as effective fractures, and the bright sinusoidal curves in the well imaging logging image are used as invalid fractures to obtain fracture identification results; based on the fracture identification results, the fracture line density and effective aperture of each fracture are determined; based on the fracture line density and effective aperture of each fracture, the permeability of each fracture is calculated; based on the permeability of each fracture, the fracture permeability curve data is formed.

[0062] The above crack line density is the number of cracks per unit length in the normal direction of the crack surface, and the unit is m. -1 ; For high-conductivity cracks (i.e., effective cracks), the effective opening is the width between the two walls of the crack, and for high-resistance cracks (i.e., invalid cracks), the effective opening of the crack is 0.

[0063] For example, the permeability of each fracture can be calculated according to the following formula:

[0064] K=b 3 D / 12×10 15

[0065] Among them, K is the permeability of the fracture, b is the effective opening, and D is the fracture line density.

[0066] Step 103: obtaining matrix permeability curve data from the well logging interpretation data;

[0067] Step 104: Add the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section formation;

[0068] Step 105: determining a permeability threshold value of a high permeability layer according to the total permeability curve of the formation in the entire well section;

[0069] Specifically, based on the total permeability curve of the formation in the entire well section, an intersection diagram of the total formation permeability of the layer section where drilling fluid loss occurs and the drilling fluid loss amount, as well as an intersection diagram of the average total formation permeability of the single well target layer and the open flow rate can be established; based on the intersection diagram of the total formation permeability of the layer section where drilling fluid loss occurs and the drilling fluid loss amount, and the intersection diagram of the average total formation permeability of the single well target layer and the open flow rate, a trend line is drawn; the total formation permeability corresponding to the point where the slope of the trend line undergoes an obvious mutation is used as the permeability threshold value of the high permeability layer.

[0070] Step 106: According to the permeability threshold value of the high permeability layer, the high permeability layer is divided for the single well in the target exploration area to obtain a high permeability layer division result.

[0071] In the actual division of high permeability layers, the layer sections in each single well in the target exploration area whose total permeability is higher than the permeability threshold value of the high permeability layer can be regarded as high permeability layers; based on the determined distribution of high permeability layers in each single well, the thickness of the high permeability layers in the entire well section and each layer section of each single well and the proportion of the total formation thickness are statistically obtained.

[0072] That is, according to the above-mentioned high permeability layer division standard and based on the total formation permeability curve, the high permeability layer of a single well can be identified and divided, and according to the formation stratification, the total thickness of the high permeability layer of the entire well section and each layer section and the proportion of the total formation thickness can be counted. For wells that have not drilled through the target layer, the high permeability layer division situation of the adjacent wells can be used for fitting and inference.

[0073] Furthermore, on the basis of regional geological background and sedimentary phase research, the planar distribution of high permeability layers in fracture-pore reservoirs can be characterized to predict wells that may experience abnormal water breakthrough, and then formulate reasonable development measures to ensure the efficient development of oil and gas reservoirs.

[0074] For example, after dividing the single wells in the target exploration area into high permeability layers according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the high permeability layer thickness threshold value of the water-seeking wells in the target exploration area can be obtained; based on the high permeability layer thickness threshold value of the water-seeking wells, the high permeability layer thickness corresponding to each well of each single well in the target exploration area is determined to determine the wells with water-seeking risks, and generate early warning information.

[0075] For example, if the thickness of the high permeability layer of all wells in the study area can be counted, the confirmed water-seeping wells can be compared with normal wells to determine the threshold value of the high permeability layer thickness of water-seeping wells under current conditions. For example, the threshold value of the high permeability layer thickness of the Paleogene water-seeping wells in a gas field is 65m. The high permeability layer thicknesses of wells DN2-2 and DN2-23 in the gas field are 105m and 80m respectively, both exceeding 65m. Therefore, it can be determined that wells DN2-2 and DN2-23 are wells that may see water, and waterproof protection measures need to be taken.

[0076] After dividing the high permeability layers of the single wells in the target exploration area according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the sedimentary phases of each single well in the target exploration area can be divided according to the core phase marks, lithology and lithofacies combinations, imaging logging images, and logging curve characteristics of the single wells in the target exploration area, and a planar distribution map of the sedimentary phases is formed; the thickness of the high permeability layer of each single well is used as hard data and the planar distribution of the sedimentary phase is used as a qualitative constraint condition to draw a contour map of the high permeability layer thickness in the target exploration period to obtain the distribution law of the high permeability layer thickness on the plane.

[0077] The above method is described below in conjunction with a specific embodiment. However, it should be noted that this specific embodiment is only for better illustrating the present application and does not constitute an improper limitation on the present application.

[0078] In this example, a method for dividing high permeability layers of fracture-pore type reservoirs and early warning of water-seeping wells is provided, which may include the following steps:

[0079] S1: Use imaging logging data to identify fractures and calculate fracture permeability;

[0080] Among them, crack identification and permeability calculation can include:

[0081] 1) On the full-well imaging logging image, natural fractures are identified using the human-computer interaction mode in the Borview module of GeoFrame software. Natural fractures appear as sinusoidal curves on the imaging logging image. The dark sinusoidal curves are high-conductivity fractures, i.e., effective fractures, and the bright sinusoidal curves are high-resistance fractures, i.e., invalid fractures.

[0082] 2) Based on the crack identification results, quantitative parameters such as crack line density and effective opening are derived from the software. The crack line density is defined as the number of cracks per unit length in the normal direction of the crack surface, with the unit of m. -1 ; For high-conductivity cracks, the effective opening is the width between the two walls of the crack, and for high-resistance cracks, the effective opening of the crack is zero.

[0083] 3) Calculate the fracture permeability curve of the entire well section according to the following fracture permeability calculation formula:

[0084] K=b 3 D / 12×10 15

[0085] Where K is the permeability of the fracture, in mD, b is the effective opening, in m, and D is the fracture line density, in m -1 .

[0086] S2: The total permeability of the formation is obtained by adding the fracture permeability and the matrix permeability interpreted by logging. The permeability threshold of the high permeability layer is determined by using the relationship between the total permeability and the drilling fluid loss and open flow data to divide the high permeability layer of a single well.

[0087] Among them, the calculation of total formation permeability and the division of high permeability layers can include:

[0088] 1) Add the fracture permeability curve data and the matrix permeability curve data to obtain the total permeability curve of the entire well section.

[0089] 2) Establish an intersection diagram of the total formation permeability and the amount of drilling fluid loss in the layer section where drilling fluid loss occurs, as well as an intersection diagram of the average total formation permeability and the open flow rate of the target layer of a single well, and draw a trend line. The total formation permeability corresponding to the point where the slope of the trend line undergoes an obvious mutation is the permeability threshold value of the high permeability layer.

[0090] 3) According to the high permeability layer classification standard (i.e., the layer section with a total permeability higher than this permeability threshold value is a high permeability layer), based on the total formation permeability curve, the high permeability layer of a single well is identified and classified, and according to the formation stratification, the total thickness of the high permeability layer of the entire well section and each layer section and the proportion of the total formation thickness are counted. For wells that have not penetrated the target layer, inferences can be made based on the high permeability layer classification of adjacent wells.

[0091] S3: Characterize the thickness planar distribution law of the high permeability layer in combination with the sedimentary phase characteristics;

[0092] Among them, the plane distribution characterization of the high permeability layer may include:

[0093] 1) Combined with the regional sedimentary background, the single well sedimentary facies is identified and divided according to the single well core phase marks, lithology and lithofacies combinations, imaging logging images, and conventional logging curve characteristics, and then the planar distribution map of the sedimentary facies is drawn.

[0094] 2) Taking the total thickness of the high-permeability layer as hard data and the planar distribution of the sedimentary phase as the qualitative constraint condition, a contour map of the high-permeability layer thickness is drawn to obtain the distribution law of the high-permeability layer thickness on the plane.

[0095] S4: Provide early warning for wells that may experience water breakthrough.

[0096] Among them, water well warning may include:

[0097] 1) Count the thickness of the high permeability layer of all wells in the target study area, compare the water-seeping wells with normal wells, and determine the threshold value of the high permeability layer thickness of the water-seeping wells.

[0098] 2) Based on the thickness threshold of the high permeability layer in water-bearing wells, predict the wells that may experience water-bearing conditions among the currently normal wells so that corresponding protective measures can be taken.

[0099] Through the above method, the distribution of high permeability layers in fracture-pore reservoirs can be quantitatively evaluated, and wells that may experience water breakthrough can be predicted. This method has the advantages of being systematic and quantitative, with sufficient geological basis, strong operability and high accuracy.

[0100] Taking the Paleogene of the Dina 2 gas reservoir in the Kuche Depression of the Tarim Basin as an example, the high permeability layer division of the fracture-pore type reservoir and the early warning method of water breakthrough wells provided in this example can include:

[0101] S1: Crack identification and permeability calculation:

[0102] 1) On the full-well imaging logging image, the human-computer interaction mode is used in the Borview module of GeoFrame software to identify the following Figure 2 The natural fractures shown in the figure are shown as sinusoidal curves on the imaging logging image, wherein the dark sinusoidal curves are high-conductivity fractures, i.e., effective fractures, and the bright sinusoidal curves are high-resistance fractures, i.e., ineffective fractures. For example, Figure 2 The sinusoidal curve in the figure is the crack position.

[0103] 2) Based on the crack identification results, quantitative parameters such as crack line density and effective opening are derived from GeoFrame software. The crack line density is defined as the number of cracks per unit length in the normal direction of the crack surface, with the unit of m. -1 ; For high-conductivity cracks, the effective opening is the width between the two walls of the crack, and for high-resistance cracks, the effective opening of the crack is zero.

[0104] 3) According to the calculation formula of fracture permeability, calculate the fracture permeability curve of the entire well section and import the result data into ResForm software:

[0105] K=b 3 D / 12×10 15

[0106] Where K is the permeability of the fracture, in mD, b is the effective opening, in m, and D is the fracture line density, in m -1 .

[0107] S2: Calculation of total formation permeability and division of high permeability layers:

[0108] 1) In ResForm software, the fracture permeability curve data and the matrix permeability curve data obtained by logging interpretation were resampled at intervals of 0.125 m and exported to Excel. The two were added to obtain the total formation permeability data. The total permeability data was imported into ResForm software to obtain the following: Figure 3 The total permeability curve of the formation in the entire well section is shown.

[0109] 2) Establish the intersection diagram of the total permeability of the formation and the amount of drilling fluid loss in the layer section where drilling fluid loss occurs, as well as the intersection diagram of the average total permeability of the formation of the target layer of a single well and the open flow rate, and draw a trend line. The total permeability of the formation corresponding to the point where the slope of the trend line changes significantly is the permeability threshold of the high permeability layer. The layer section with a total permeability higher than this permeability threshold is the high permeability layer. For example, the permeability threshold of the high permeability layer of the Paleogene in the Dina 2 gas field is as follows: Figure 4 and Figure 5 As shown, it is 7mD. Correspondingly, the layer with a total permeability higher than 7mD is the high permeability layer. Figure 4 is the cross-plot of total permeability and drilling fluid loss. The lower left corner curve corresponds to: y = 23.1x + 12.407, R 2 =0.8479, R 2 Represents the correlation coefficient of the formula; the upper right corner curve corresponds to: y = 137.74x-688.08, R 2 =0.5883. Figure 5 It is the intersection diagram of total permeability and open flow rate. The lower left corner curve corresponds to: y = 4.0109x + 96.759, R 2 =0.9305; the upper right corner curve corresponds to: y = 83.812x-379.32, R 2 =0.6513.

[0110] 3) According to the high permeability layer classification standard, the high permeability layer of a single well is classified in the ResForm software based on the total permeability curve of the formation. Figure 3 The identification divisions are shown, and the total thickness of the hyperpermeable layer is calculated.

[0111] S3: Planar distribution characterization of the hyperosmotic layer:

[0112] 1) Combined with the regional sedimentary background, the single well sedimentary facies is identified and divided according to the single well core phase marks, lithology and lithofacies combinations, imaging logging images, and conventional logging curve characteristics, and then the planar distribution map of the sedimentary facies is drawn, in which the sedimentary phase is divided into subfacies and below.

[0113] 2) Taking the total thickness of the high permeability layer of a single well as hard data and the planar distribution of the sedimentary phase as a qualitative constraint, a contour map of the high permeability layer thickness is drawn to obtain the distribution law of the high permeability layer thickness on the plane. The sedimentary phase has a significant correlation with the thickness of the high permeability layer. The Suweiyi Formation of the Paleogene in the Dina 2 gas field develops fan delta plain phase, fan delta front phase and coastal shallow lake phase. The thickness of the high permeability layer in the fan delta plain phase is the largest, followed by the fan delta front phase, and the coastal shallow lake phase is the smallest. With this as a constraint, the following is drawn: Figure 6 The plane contour map of the high permeability layer thickness is shown.

[0114] S4: See water well warning:

[0115] 1) Count the thickness of the high permeability layer of all wells in the study area, compare the confirmed water-seeping wells with normal wells, and determine the threshold value of the high permeability layer thickness of water-seeping wells under current conditions. Currently, the threshold value of the high permeability layer thickness of the water-seeping wells in the Paleogene of the Dina 2 gas field is 65m.

[0116] 2) According to the threshold value of the thickness of the high permeability layer of the water-seeping wells, the wells that may see water among the normal wells are predicted so that corresponding protective measures can be taken. The thickness of the high permeability layer of the DN2-2 and DN2-23 wells in the Indina 2 gas field are 105m and 80m respectively, both exceeding 65m. Therefore, it can be determined that the DN2-2 and DN2-23 wells are wells that may see water and need to take waterproof protection measures.

[0117] The method embodiments provided in the above embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking running on an electronic device as an example, Figure 7 This is a hardware structure block diagram of an electronic device for a method of dividing a high permeability layer of a fracture pore type reservoir provided in this application. Figure 7 As shown, the electronic device 10 may include one or more (only one is shown in the figure) processors 02 (the processor 02 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA), a memory 04 for storing data, and a transmission module 06 for communication functions. It can be understood by those skilled in the art that Figure 7 The structure shown is only for illustration and does not limit the structure of the above electronic device. Figure 7 More or fewer components as shown, or with Figure 7 Different configurations are shown.

[0118] The memory 04 can be used to store software programs and modules of application software, such as program instructions / modules corresponding to the method for dividing the high permeability layer of the fracture pore type reservoir in the embodiment of the present application. The processor 02 executes various functional applications and data processing by running the software programs and modules stored in the memory 04, that is, the method for dividing the high permeability layer of the fracture pore type reservoir of the above-mentioned application is realized. The memory 04 may include a high-speed random access memory, and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 04 may further include a memory remotely arranged relative to the processor 02, and these remote memories may be connected to the electronic device 10 via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0119] The transmission module 06 is used to receive or send data via a network. The specific example of the above network may include a wireless network provided by a communication provider of the electronic device 10. In one example, the transmission module 06 includes a network adapter (Network Interface Controller, NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission module 06 can be a radio frequency (Radio Frequency, RF) module, which is used to communicate with the Internet wirelessly.

[0120] At the software level, the high permeability layer division device of the above-mentioned fracture pore type reservoir can be as follows Figure 8 As shown, including:

[0121] The first acquisition module 801 is used to acquire full-well imaging logging images and logging interpretation data of the target exploration area;

[0122] A generating module 802 is used to obtain fracture permeability curve data according to the full-well imaging logging image;

[0123] The second acquisition module 803 is used to acquire matrix permeability curve data from the well logging interpretation data;

[0124] A processing module 804 is used to add the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section formation;

[0125] A determination module 805 is used to determine the permeability threshold value of the high permeability layer according to the total permeability curve of the formation in the entire well section;

[0126] The division module 806 is used to divide the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer to obtain a high permeability layer division result.

[0127] In one embodiment, the above-mentioned generation module 802 is specifically used to identify sinusoidal curves from the full-section imaging logging image; use the dark sinusoidal curves in the section imaging logging image as effective fractures, and use the bright sinusoidal curves in the section imaging logging image as invalid fractures to obtain fracture identification results; determine the fracture line density and effective aperture of each fracture based on the fracture line density and effective aperture of each fracture; calculate the permeability of each fracture based on the fracture line density and effective aperture of each fracture; and form fracture permeability curve data based on the permeability of each fracture.

[0128] In one embodiment, the permeability of each fracture may be calculated according to the following formula:

[0129] K=b 3 D / 12×1015

[0130] Among them, K is the permeability of the fracture, b is the effective opening, and D is the fracture line density.

[0131] In one embodiment, the determination module 805 can be specifically used to establish an intersection diagram of the total formation permeability of the layer section where drilling fluid loss occurs and the drilling fluid loss amount, as well as an intersection diagram of the average total formation permeability of the single well target layer and the open flow rate based on the total permeability curve of the entire well section; draw a trend line based on the intersection diagram of the total formation permeability of the layer section where drilling fluid loss occurs and the drilling fluid loss amount, and the intersection diagram of the average total formation permeability of the single well target layer and the open flow rate; and use the total formation permeability corresponding to the point where the slope of the trend line undergoes an obvious mutation as the permeability threshold value of the high permeability layer.

[0132] In one embodiment, the division module 806 can be specifically used to regard the layer segments in each single well in the target exploration area whose total permeability is higher than the permeability threshold value of the high permeability layer as the high permeability layer; based on the determined distribution of the high permeability layer of each single well, the thickness of the high permeability layer of the entire well segment and each layer segment of each single well and the proportion of the total formation thickness are statistically obtained.

[0133] In one embodiment, the device for dividing the high permeability layer of the above-mentioned fractured pore reservoir can also be used to divide the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer, and after obtaining the high permeability layer division result, obtain the high permeability layer thickness threshold value of the water-seeking wells in the target exploration area; according to the high permeability layer thickness threshold value of the water-seeking wells, the high permeability layer thickness corresponding to each well of each single well in the target exploration area is determined to determine the wells with water risk, and generate early warning information.

[0134] In one embodiment, the high permeability layer division device of the above-mentioned fracture pore type reservoir can also be used to divide the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer. After obtaining the high permeability layer division result, the sedimentary phases of each single well in the target exploration area are divided according to the core phase marks, lithology and lithofacies combinations, imaging logging images, and logging curve characteristics of the single well in the target exploration area, and a planar distribution map of the sedimentary phases is formed; the thickness of the high permeability layer of each single well is used as hard data and the planar distribution of the sedimentary phase is used as a qualitative constraint condition to draw a contour map of the high permeability layer thickness in the target exploration period to obtain the distribution law of the high permeability layer thickness on the plane.

[0135] The embodiment of the present application also provides a specific implementation of an electronic device capable of implementing all the steps in the method for dividing a high permeability layer of a fracture pore type reservoir in the above embodiment, and the electronic device specifically includes the following contents: a processor, a memory, a communication interface, and a bus; wherein the processor, the memory, and the communication interface communicate with each other through the bus; the processor is used to call a computer program in the memory, and when the processor executes the computer program, all the steps in the method for dividing a high permeability layer of a fracture pore type reservoir in the above embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0136] Step 1: Obtain full-well imaging logging images and logging interpretation data of the target exploration area;

[0137] Step 2: obtaining fracture permeability curve data according to the full-well imaging logging image;

[0138] Step 3: Obtaining matrix permeability curve data from the well logging interpretation data;

[0139] Step 4: Add the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section;

[0140] Step 5: Determine the permeability threshold value of the high permeability layer according to the total permeability curve of the formation in the entire well section;

[0141] Step 6: According to the permeability threshold value of the high permeability layer, the high permeability layer is divided for the single well in the target exploration area to obtain a high permeability layer division result.

[0142] The embodiment of the present application also provides a computer-readable storage medium capable of implementing all the steps in the method for dividing a high permeability layer of a fractured pore type reservoir in the above embodiment. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, all the steps in the method for dividing a high permeability layer of a fractured pore type reservoir in the above embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0143] Step 1: Obtain full-well imaging logging images and logging interpretation data of the target exploration area;

[0144] Step 2: obtaining fracture permeability curve data according to the full-well imaging logging image;

[0145] Step 3: Obtaining matrix permeability curve data from the well logging interpretation data;

[0146] Step 4: Add the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section;

[0147] Step 5: Determine the permeability threshold value of the high permeability layer according to the total permeability curve of the formation in the entire well section;

[0148] Step 6: According to the permeability threshold value of the high permeability layer, the high permeability layer is divided for the single well in the target exploration area to obtain a high permeability layer division result.

[0149] From the above description, it can be seen that the embodiment of the present application obtains fracture permeability curve data through the full-well imaging logging image of the target exploration area, obtains matrix permeability curve data through logging interpretation data, and then adds the fracture permeability curve data and the matrix permeability curve data to obtain the total permeability curve of the formation in the whole well section; based on the total permeability curve, the permeability threshold value of the high permeability layer is determined, so that the high permeability layer can be accurately determined. The above scheme solves the existing technical problem of being unable to effectively determine the distribution of high permeability layers in fracture pore type reservoirs, and achieves the technical effect of accurately and efficiently determining the distribution of high permeability layers in fracture pore type reservoirs, thereby increasing the production of oil and gas reservoirs.

[0150] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the hardware + program embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0151] The above is a description of a specific embodiment of the specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0152] Although the present application provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-creative labor. The order of steps listed in the embodiments is only one way of executing the order of many steps and does not represent the only execution order. When the actual device or client product is executed, it can be executed in the order of the method shown in the embodiments or the drawings or in parallel (for example, in a parallel processor or multi-threaded processing environment).

[0153] Although the present specification embodiment provides the method operation steps as described in the embodiment or flow chart, more or less operation steps may be included based on conventional or non-creative means. The order of steps listed in the embodiment is only one way in the order of execution of many steps, and does not represent a unique execution order. When the device or terminal product in practice is executed, it can be executed in sequence or in parallel (such as a parallel processor or a multi-threaded processing environment, or even a distributed data processing environment) according to the method shown in the embodiment or the accompanying drawings. The term "include", "comprise" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, product or equipment including a series of elements not only includes those elements, but also includes other elements not clearly listed, or also includes elements inherent to such process, method, product or equipment. In the absence of more restrictions, it is not excluded that there are other identical or equivalent elements in the process, method, product or equipment including the elements.

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

[0155] Those skilled in the art also know that, in addition to implementing the controller in a purely computer-readable program code, the controller can be made to implement the same function in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded microcontrollers by logically programming the method steps. Therefore, such a controller can be considered as a hardware component, and the devices for implementing various functions included therein can also be considered as structures within the hardware component. Or even, the devices for implementing various functions can be considered as both software modules for implementing the method and structures within the hardware component.

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

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

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

[0159] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0160] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0161] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0162] Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems or computer program products. Therefore, the embodiments of this specification may take the form of complete hardware embodiments, complete software embodiments or embodiments combining software and hardware. Moreover, the embodiments of this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0163] The present specification embodiments may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present specification embodiments may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0164] Each embodiment in this specification is described in a progressive manner, and the same and similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment. In the description of this specification, the description of the reference terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the embodiment of this specification. In this specification, the schematic representation of the above terms does not necessarily target the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner. In addition, in the absence of contradiction, a person skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0165] The above is only an example of the embodiment of the present specification and is not intended to limit the embodiment of the present specification. For those skilled in the art, the embodiment of the present specification may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiment of the present specification shall be included in the scope of the claims of the embodiment of the present specification.

Claims

1. A method for dividing a high permeability layer of a fractured pore reservoir, characterized in that: The method comprises: Obtain full-well imaging logging images and logging interpretation data for the target exploration area of ​​fracture-pore reservoirs; Obtaining fracture permeability curve data according to the full-well imaging logging image; Acquiring matrix permeability curve data from the well logging interpretation data; Adding the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section; Determining a permeability threshold value of a high permeability layer according to a total permeability curve of the entire well section formation; According to the permeability threshold value of the high permeability layer, the high permeability layer is divided for the single well in the target exploration area to obtain a high permeability layer division result.

2. The method according to claim 1, characterized in that According to the full-well imaging logging image, fracture permeability curve data is obtained, including: identifying a sinusoidal curve from the full-well imaging logging image; Taking the dark sine curve in the well section imaging logging image as a valid fracture, and taking the bright sine curve in the well section imaging logging image as an invalid fracture, so as to obtain a fracture identification result; According to the crack identification results, the crack line density and effective opening of each crack are determined; The permeability of each fracture is calculated based on the fracture line density and effective aperture of each fracture; According to the permeability of each fracture, fracture permeability curve data is generated.

3. The method according to claim 2, characterized in that According to the fracture line density and effective aperture of each fracture, the permeability of each fracture is calculated, including: The permeability of each crack is calculated according to the following formula: K=b 3 D / 12×10 15 Among them, K is the permeability of the fracture, b is the effective opening, and D is the fracture line density.

4. The method according to claim 1, characterized in that: Determining the permeability threshold value of the high permeability layer according to the total permeability curve of the formation in the entire well section includes: According to the total permeability curve of the formation in the whole well section, establish the intersection diagram of the total permeability of the formation in the layer section where drilling fluid loss occurs and the amount of drilling fluid loss, as well as the intersection diagram of the average total permeability of the formation of the target layer in a single well and the open flow rate; Draw a trend line based on the cross-plot of the total formation permeability and the amount of drilling fluid loss in the layer where drilling fluid loss occurs, and the cross-plot of the average total formation permeability and the open flow rate of the target layer of a single well; The total formation permeability corresponding to the point where the slope of the trend line undergoes an obvious mutation is used as the permeability threshold value of the high permeability layer.

5. The method according to claim 1, characterized in that According to the permeability threshold value of the high permeability layer, the high permeability layer of the single well in the target exploration area is divided, including: The layer section whose total permeability in each single well of the target exploration area is higher than the permeability threshold value of the high permeability layer is regarded as the high permeability layer; According to the distribution of high permeability layers in each single well, the thickness of high permeability layers in the whole well section and each layer section of each single well and the proportion of high permeability layers in the total formation thickness are obtained statistically.

6. The method according to claim 5, characterized in that After dividing the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the method further includes: Obtaining a high permeability layer thickness threshold value of water wells in the target exploration area; According to the high permeability layer thickness threshold value of the water-breaking wells and the high permeability layer thickness corresponding to each well in the target exploration area, the wells with water-breaking risk are determined and early warning information is generated.

7. The method according to claim 5, characterized in that After dividing the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer and obtaining the high permeability layer division result, the method further includes: According to the single well core phase marks, lithology and lithofacies combinations, imaging logging images, and logging curve characteristics of the target exploration area, each single well in the target exploration area is divided into sedimentary phases, and a planar distribution map of the sedimentary phases is formed; Taking the thickness of the high permeability layer of each single well as hard data and the planar distribution of the sedimentary phase as the qualitative constraint condition, the contour map of the high permeability layer thickness in the target exploration period was drawn to obtain the distribution law of the high permeability layer thickness on the plane.

8. A device for dividing a high permeability layer of a fractured pore reservoir, characterized in that: include: The first acquisition module is used to acquire full-well imaging logging images and logging interpretation data of the target exploration area of ​​the fracture-pore type reservoir; A generating module, used for obtaining fracture permeability curve data according to the full-well imaging logging image; A second acquisition module is used to acquire matrix permeability curve data from the well logging interpretation data; A processing module, used for adding the fracture permeability curve data and the matrix permeability curve data to obtain a total permeability curve of the entire well section formation; A determination module, used for determining a permeability threshold value of a high permeability layer according to a total permeability curve of the formation in the entire well section; The division module is used to divide the high permeability layer of the single well in the target exploration area according to the permeability threshold value of the high permeability layer to obtain the high permeability layer division result.

9. An electronic device comprising a processor and a memory for storing instructions executable by the processor, characterized in that: When the processor executes the instructions, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instructions are executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images