CN113914853A

CN113914853A - Method for determining water saturation of sand shale thin interbed sandstone reservoir of deep water sedimentation system

Info

Publication number: CN113914853A
Application number: CN202010655101.8A
Authority: CN
Inventors: 邓红婴; 周进高; 左国平; 邵大力; 丁梁波; 马宏霞; 庞旭; 郭沫贞; 刘艳红
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2020-07-09
Filing date: 2020-07-09
Publication date: 2022-01-11

Abstract

The invention provides a method for determining the water saturation of a sandstone reservoir of a shale thin interbed of a deepwater sedimentation system, which comprises the following steps: acquiring conventional logging data and Rt-Scanner information of a target thin interbed; processing the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters; processing the Rt-Scanner data to obtain longitudinal resistivity; and obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity. By adopting the technical scheme, the longitudinal resistivity obtained by processing the Rt-Scanner data is utilized, the influence of the layered argillaceous surrounding rock on the resistivity of the target layer is eliminated to the maximum extent, the sensitivity of thin-layer hydrocarbon reservoir identification is improved, the water saturation or the oil-gas saturation is accurately calculated, and the phenomenon that the thin-layer hydrocarbon reservoir is missed and underestimated by the conventional well logging interpretation method is avoided.

Description

Method for determining water saturation of sand shale thin interbed sandstone reservoir of deep water sedimentation system

Technical Field

The invention relates to the technical field of petroleum geological exploration and well logging, in particular to a method and a device for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of a deep water deposition system, electronic equipment and a storage medium.

Background

The water saturation is one of important parameters for evaluating oil and gas reservoirs and is a core parameter for quantitative evaluation of the reservoirs, the water saturation is required to be accurately calculated, the properties and characteristics of oil and gas reservoir fluids are truly and objectively reflected, the apparent resistivity of a stratum measured by logging is required to be close to the true resistivity of the stratum to the greatest extent possible, however, the main source of the apparent resistivity measured by the traditional resistivity logging is the contribution of the transverse resistivity of the stratum, in a sandstone thin interbed stratum, the measurement current mainly passes through the shale stratum with smaller surrounding resistance, so the apparent resistivity of the sandstone stratum obtained by measurement is mainly derived from the contribution of the surrounding shale stratum, the measurement value of the sandstone stratum is far smaller than the true resistivity of the sandstone stratum, the water saturation or oil and gas saturation of the sandstone cannot be accurately calculated, and the phenomena of missing and underestimating the thin oil and gas layer are often caused.

Disclosure of Invention

Aiming at the problems in the prior art, the invention provides a method and a device for determining the water saturation of a sandstone reservoir in a shale thin interbed of a deepwater sedimentation system, electronic equipment and a storage medium, which can at least partially solve the problems in the prior art.

In order to achieve the purpose, the invention adopts the following technical scheme:

in a first aspect, a method for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of a deepwater sedimentary system is provided, which comprises the following steps:

acquiring conventional logging data and Rt-Scanner information of a target thin interbed;

processing the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters;

processing the Rt-Scanner data to obtain longitudinal resistivity;

and obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity.

Further, the logging response characteristic parameters include: formation water resistivity, mudstone resistivity; the sandstone reservoir physical property response parameters comprise: argillaceous content and effective porosity.

Further, the conventional well log data includes: a natural gamma log, a neutron log, a density log, and a deep resistivity curve;

the conventional logging data are processed to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters, and the logging response characteristic parameters and the sandstone reservoir physical property response parameters comprise the following steps:

obtaining the mud content of the target thin interbed according to the natural gamma logging curve;

obtaining effective porosity and total porosity according to the neutron logging curve, the density logging curve and the shale content;

dividing the sandstone reservoir according to the argillaceous content and the effective porosity;

and acquiring formation water resistivity and mudstone resistivity according to the deep resistivity curve, the lithology recognition result and the sandstone reservoir division result.

Further, obtaining effective porosity and total porosity from the neutron log, the density log, and the shale content, comprising:

establishing a neutron-density logging intersection graph according to the neutron logging curve and the density logging curve;

determining a stratum skeleton parameter according to the neutron-density logging intersection map;

and calculating the effective porosity and the total porosity according to the neutron logging curve, the density logging curve, the stratum skeleton parameter and the argillaceous content.

Further, dividing the sandstone reservoir according to the argillaceous content and the effective porosity comprises:

carrying out sensitivity analysis on the shale content and the effective porosity to obtain a shale content upper limit value and a porosity lower limit value;

identifying lithology according to the shale content and the upper limit value of the shale content;

and dividing the sandstone reservoir according to the effective porosity, the lower porosity limit value and the identified lithology.

Further, before dividing the sandstone reservoir according to the shale content and the effective porosity, the method further comprises the following steps:

establishing a shale content-total porosity intersection graph according to the shale content and the total porosity;

analyzing the shale content-total porosity intersection diagram by utilizing a Thomas Stieber model to obtain a shale distribution type of the target thin interbed, wherein the shale distribution type comprises: lamellar argillaceous, dispersive argillaceous, and structural argillaceous;

and when the mud distribution type of the target thin interbed is dispersed mud or structural mud, ending the method flow.

Further, before the conventional logging data is processed to obtain the logging response characteristic parameter and the sandstone reservoir physical property response parameter, the method further comprises the following steps:

performing quality analysis on the conventional logging data;

and performing environmental correction and borehole correction on the conventional logging data according to the quality analysis result.

Further, the method for determining the water saturation of the sandstone reservoir in the shale thin interbed of the deepwater deposition system further comprises the following steps:

and acquiring the oil-gas saturation of the sandstone reservoir according to the water saturation of the sandstone reservoir.

In a second aspect, an apparatus for determining water saturation of a sandstone reservoir in a thin interbed of sandstone in a deepwater deposition system is provided, which includes:

the data acquisition module is used for acquiring conventional logging data and Rt-Scanner data of the target thin interbed;

the conventional logging data processing module is used for processing the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters;

the longitudinal resistivity acquisition module is used for processing the Rt-Scanner data to obtain longitudinal resistivity;

and the water saturation calculation module is used for obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity.

the conventional well logging data processing module comprises:

the mud content acquisition unit is used for acquiring the mud content of the target thin interbed according to the natural gamma logging curve;

the porosity obtaining unit is used for obtaining effective porosity and total porosity according to the neutron logging curve, the density logging curve and the shale content;

the sandstone reservoir dividing unit is used for dividing the sandstone reservoir according to the argillaceous content and the effective porosity;

and the resistivity acquisition unit is used for acquiring formation water resistivity and mudstone resistivity according to the deep resistivity curve, the lithology recognition result and the sandstone reservoir division result.

Further, the porosity obtaining unit includes:

a neutron-density logging intersection graph establishing subunit, which is used for establishing a neutron-density logging intersection graph according to the neutron logging curve and the density logging curve;

a stratum skeleton parameter determining subunit, which determines a stratum skeleton parameter according to the neutron-density logging intersection map;

and the porosity calculation operator unit is used for calculating the effective porosity and the total porosity according to the neutron logging curve, the density logging curve, the stratum skeleton parameter and the shale content.

Further, the sandstone reservoir partition unit comprises:

a limit value obtaining subunit, which performs sensitivity analysis on the shale content and the effective porosity to obtain a shale content upper limit value and a porosity lower limit value;

a lithology identifying subunit, for identifying lithology according to the shale content and the upper limit value of the shale content;

and the reservoir dividing subunit is used for dividing the sandstone reservoir according to the effective porosity, the lower limit value of the porosity and the identified lithology.

Further, the conventional well logging data processing module further comprises:

a shale content-total porosity intersection map establishing unit for establishing a shale content-total porosity intersection map according to the shale content and the total porosity;

a argillaceous distribution type analysis unit, which analyzes the argillaceous content-total porosity intersection diagram by using a Thomas Stieber model to obtain the argillaceous distribution type of the target thin interbed, wherein the argillaceous distribution type comprises: lamellar argillaceous, dispersive argillaceous, and structural argillaceous;

and the flow control unit is used for finishing the flow of the method when the mud distribution type of the target thin interbed is dispersed mud or structural mud.

Further, the device for determining the water saturation of the sandstone reservoir in the thin interbed of the sandstone of the deepwater deposition system further comprises:

the quality analysis module is used for carrying out quality analysis on the conventional logging data;

and the data correction module is used for carrying out environmental correction and borehole correction on the conventional logging data according to the quality analysis result.

and the oil-gas saturation calculation module is used for acquiring the oil-gas saturation of the sandstone reservoir according to the water saturation of the sandstone reservoir.

In a third aspect, an electronic device is provided, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, and the processor executes the computer program to implement the steps of the method for determining the water saturation of a sandstone reservoir in a sandstone thin interbed of a deepwater sedimentation system.

In a fourth aspect, a computer readable storage medium is provided, on which a computer program is stored, which computer program, when being executed by a processor, is adapted to carry out the steps of the method for determining the water saturation of a sandstone reservoir in a sandstone thin interbed of a deepwater sedimentary system as described above.

The invention provides a method and a device for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of a deepwater sedimentation system, electronic equipment and a storage medium, wherein the method comprises the following steps: acquiring conventional logging data and Rt-Scanner information of a target thin interbed; processing the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters; processing the Rt-Scanner data to obtain longitudinal resistivity; and obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity. By adopting the technical scheme, the longitudinal resistivity obtained by processing the Rt-Scanner data is utilized, the influence of the layered argillaceous surrounding rock on the resistivity of the target sandstone reservoir is eliminated to the maximum extent, the sensitivity of thin-layer hydrocarbon reservoir identification is improved, the water saturation or the oil-gas saturation is accurately calculated, the phenomenon that the thin-layer hydrocarbon reservoir is omitted and underestimated by the conventional well logging interpretation method is avoided, and the interpretation precision and accuracy of the thin-layer hydrocarbon reservoir are greatly improved.

In order to make the aforementioned and other objects, features and advantages of the invention comprehensible, preferred embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts. In the drawings:

FIG. 1 is a schematic diagram of an architecture between a server S1 and a client device B1 according to an embodiment of the present invention;

FIG. 2 is a first schematic flow chart of a method for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of the deepwater sedimentary system in the embodiment of the invention;

FIG. 3 shows a plot of shale content versus measured resistivity;

FIG. 4 is a schematic diagram of the invention, the left diagram showing the path of the longitudinal resistivity measurement current through the path; the right graph is a conventional measuring current passing path graph;

FIG. 5 is a schematic flow chart of a sandstone reservoir water saturation determination method in a sandstone-shale thin interbed of the deepwater deposition system in the embodiment of the invention;

FIG. 6 is a third schematic flow chart of a method for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of the deepwater sedimentary system in the embodiment of the invention;

fig. 7 shows the specific steps of step S200 in fig. 2 or fig. 5 or fig. 6;

fig. 8 shows the detailed steps of step S220 in fig. 7;

fig. 9 shows the detailed steps of step S230 in fig. 7;

FIG. 10 is a technical flow diagram of an implementation of the present invention;

FIG. 11 is a comparison graph of the interpretation and evaluation of a gas production zone of a conventional method applied to a deep water sedimentation system A well P sand group and the method of the invention;

FIG. 12 is a first block diagram illustrating the configuration of a sandstone reservoir water saturation determination device in a sandstone-shale thin interbed of a deepwater sedimentary system according to an embodiment of the present invention;

FIG. 13 is a block diagram II of the structure of a sandstone reservoir water saturation determination device in a sandstone-shale thin interbed of the deepwater sedimentary system in the embodiment of the invention;

FIG. 14 is a block diagram of a third structural diagram of a sandstone reservoir water saturation determination device in a sandstone-shale thin interbed of the deepwater deposition system in the embodiment of the invention;

FIG. 15 shows a specific structure of the conventional well log data processing module 20 in any of FIGS. 12 to 14;

fig. 16 shows a specific structure of the porosity obtaining unit 22 in fig. 15;

fig. 17 shows a specific structure of the sandstone reservoir partitioning unit 23 in fig. 15;

fig. 18 is a block diagram of an electronic device according to an embodiment of the present invention.

Detailed Description

In order to make the technical solutions better understood by those skilled in the art, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only partial embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

It should be noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of this application and the above-described drawings, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.

It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.

The evaluation of the thin interbed oil and gas reservoir is a difficult problem of oil and gas exploration and development all the time, along with the deepening of the oil and gas exploration and development at home and abroad, the identification of the lithology, pore structure and seepage characteristics of the oil and gas reservoirs such as the thin interbe, shale, compact sandstone and the like and the evaluation of the oil and gas reservoir are increasingly difficult, and meanwhile, because the thin-layer logging response is greatly influenced by surrounding rocks, the traditional logging method and the traditional evaluation technology cannot accurately and effectively identify the oil and gas reservoir and calculate the oil and gas saturation of the oil and gas reservoir, and the phenomena of omission and underestimation of the thin oil and gas reservoir are often caused.

In order to at least partially solve the technical problems in the prior art, the embodiment of the invention provides a method for determining the water saturation of a sandstone reservoir in a thin interbed of sandstone and shale of a deep water sedimentary system, which utilizes the longitudinal resistivity obtained by Rt-Scanner data processing to eliminate the influence of the layered shale surrounding rock on the resistivity of a target layer to the maximum extent, improves the sensitivity of thin hydrocarbon layer identification, accurately calculates the water saturation or the oil gas saturation, avoids the phenomenon that the thin hydrocarbon layer is omitted and underestimated by the conventional well logging interpretation method, greatly improves the interpretation precision of the thin hydrocarbon layer, and provides a feasible new way for the well logging interpretation and reservoir evaluation of the complex hydrocarbon reservoir of the thin interbed of sandstone and shale of the deep water sedimentary system.

In view of the above, the present application provides an apparatus for determining water saturation of sandstone reservoir in thin interbed of sandstone in deepwater sedimentary system, which may be a server S1, see fig. 1, where the server S1 may be communicatively connected to at least one client device B1, the client device B1 may transmit conventional logging data and Rt-Scanner data of a target thin interbed to the server S1, and the server S1 may receive the conventional logging data and Rt-Scanner data of the target thin interbed online. The server S1 can carry out online or offline preprocessing on the acquired conventional logging data and Rt-Scanner data of the target thin interbed, and the conventional logging data are processed to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters; processing the Rt-Scanner data to obtain longitudinal resistivity; and obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity. The server S1 may then send the sandstone reservoir water saturation online to the client device B1. The client device B1 may receive the sandstone reservoir water saturation online.

Based on the above, the client device B1 may have a display interface so that a user can view the sandstone reservoir water saturation sent by the server S1 according to the interface.

It is understood that the client device B1 may include a smart phone, a tablet electronic device, a network set-top box, a portable computer, a desktop computer, a Personal Digital Assistant (PDA), a vehicle-mounted device, a smart wearable device, etc. Wherein, intelligence wearing equipment can include intelligent glasses, intelligent wrist-watch, intelligent bracelet etc..

In practical applications, the part of the determination of the water saturation of the sandstone reservoir in the sandstone thin interbed of the deepwater sedimentary system may be performed on the side of the server S1 as described above, that is, the architecture shown in fig. 1, or all the operations may be completed in the client device B1. Specifically, the selection may be performed according to the processing capability of the client device B1, the limitation of the user usage scenario, and the like. This is not a limitation of the present application. If all operations are completed in the client device B1, the client device B1 may further include a processor for performing specific processing for sandstone reservoir water saturation determination in a sandstone thin interbed of the deepwater sedimentary system.

The server and the client device may communicate using any suitable network protocol, including network protocols not yet developed at the filing date of this application. The network protocol may include, for example, a TCP/IP protocol, a UDP/IP protocol, an HTTP protocol, an HTTPS protocol, or the like. Of course, the network Protocol may also include, for example, an RPC Protocol (Remote Procedure Call Protocol), a REST Protocol (Representational State Transfer Protocol), and the like used above the above Protocol.

FIG. 2 is a first schematic flow chart of a method for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of the deepwater sedimentary system in the embodiment of the invention; as shown in fig. 2, the method for determining the water saturation of the sandstone reservoir in the sandstone thin interbed of the deepwater sedimentary system can comprise the following steps:

step S100: acquiring conventional logging data and Rt-Scanner information of a target thin interbed;

specifically, the target thin interbed is used as a target interval of the scheme and comprises a sandstone reservoir and a non-sandstone reservoir; the purpose of the scheme is to determine the water saturation of the sandstone reservoir in the thin interbed.

The Rt-Scanner data are obtained by measuring by using an Rt-Scanner (resistivity scanning logging) method, the three-component resistivity can be obtained by processing the measurement result of the Rt-Scanner method, and the Rt-Scanner method has the greatest advantage that the measurement current is forced to pass through the stratums of all layers, so that the influence of the shale surrounding rock on the apparent resistivity value of the sandstone stratum is reduced to the greatest extent by the measurement result. The relationship between the mud content of the thin interbed and the measured resistivity is shown in fig. 3. FIG. 4 is a schematic diagram of the invention, the left diagram showing the path of the longitudinal resistivity measurement current through the path; the right graph is a conventional measurement current through path graph.

Additionally, conventional well log data includes: the natural gamma-ray logging curve GR, the borehole diameter logging curve CAL, the neutron logging curve NHPI, the density logging curve RHOB, the sound wave time difference logging curve DTCO, the mud filtrate resistivity logging curve RMF, the shallow resistivity logging curve Rs, the deep resistivity logging curve Rd, the longitudinal resistivity logging curve Rv, the transverse resistivity logging curve Rh and the like.

Step S200: processing the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters;

specifically, the logging response characteristic parameters comprise: formation water resistivity, mudstone resistivity; the sandstone reservoir physical property response parameters comprise: argillaceous content, total porosity, and effective porosity.

Step S300: processing the Rt-Scanner data to obtain longitudinal resistivity;

specifically, the longitudinal resistivity can be obtained by processing the three-component resistivity.

Step S400: and obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the physical property response parameter of the sandstone reservoir and the longitudinal resistivity.

The longitudinal resistivity obtained by Rt-Scanner data processing eliminates the influence of the layered argillaceous surrounding rock on the resistivity of a target layer to the maximum extent, improves the sensitivity of thin-layer hydrocarbon reservoir identification, accurately calculates the water saturation or oil-gas saturation of a sandstone reservoir, avoids the phenomenon that a previous well logging interpretation method omits and underestimates the thin hydrocarbon reservoir, greatly improves the interpretation precision of the thin hydrocarbon reservoir, and is suitable for the sand shale thin interbed hydrocarbon reservoir of a deep water sedimentary system.

In an alternative embodiment, the calculation of water saturation may be performed using the following equation:

wherein Sw is the water saturation of the sandstone reservoir in the target thin interbed; rv is the longitudinal resistivity, ohmm; rw is formation water resistivity, ohmm; rsh is mudstone resistivity, ohmm; vsh is the sludge content, decimal fraction; phi is effective porosity).

By adopting the calculation mode, the water saturation of the sandstone reservoir in the target thin interbed can be accurately calculated, and the subsequent exploration and development are facilitated.

In an alternative embodiment, referring to fig. 5, the method for determining the water saturation of the sandstone reservoir in the sandstone thin interbed of the deepwater deposition system may further comprise:

step S500: and acquiring the oil-gas saturation of the sandstone reservoir according to the water saturation of the sandstone reservoir.

Specifically, the fluid saturations include a water saturation and a hydrocarbon saturation; the oil and gas saturation is equal to 1 minus the water saturation in order to accurately and efficiently evaluate sandstone reservoir and fluid properties.

In an alternative embodiment, referring to fig. 6, the method for determining the water saturation of the sandstone reservoir in the sandstone thin interbed of the deepwater deposition system may further comprise:

step S600: performing quality analysis on the conventional well logging data;

step S700: and performing environmental correction and borehole correction on the conventional logging data according to the quality analysis result.

Wherein, through carrying out environmental correction and well correction to conventional logging data to guarantee the authenticity and the accuracy of logging information.

In an alternative embodiment, referring to fig. 7, this step S200 may include the following:

step S210: obtaining the mud content of a target thin interbed according to the natural gamma logging curve;

specifically, the argillaceous content is calculated by the following formula:

wherein Vsh is the argillaceous content; GR is the natural gamma measurement, API; GR_minAPI for dealing with the natural gamma minimum of the well section; GR_maxAPI for processing the natural gamma maximum of the well section; the new formation (third family and renewed formation) GCUR is 3.7 and the old formation GCUR is 2.0.

Step S220: obtaining effective porosity and total porosity according to the neutron logging curve, the density logging curve and the shale content;

step S230: dividing the sandstone reservoir according to the argillaceous content and the effective porosity;

step S240: and acquiring formation water resistivity and mudstone resistivity according to the deep resistivity curve, the lithology recognition result and the sandstone reservoir division result.

In an alternative embodiment, referring to fig. 8, this step S220 may include the following:

step S221: establishing a neutron-density logging intersection map according to the neutron logging curve and the density logging curve;

step S222: determining a stratum skeleton parameter according to the neutron-density logging intersection map;

wherein, the formation skeleton parameters may include: rock skeleton density, rock skeleton neutron value, mudstone neutron value and mudstone density value.

Step S223: and calculating the effective porosity and the total porosity according to the neutron logging curve, the density logging curve, the stratum skeleton parameters and the shale content.

Specifically, the porosity calculation can be implemented using the following formula:

wherein PHIE represents effective porosity; RHOB refers to the value of a sampling point on the density log; RHOB_maIs the density of the rock skeleton in g/cm³；RHOB_fIs the formation fluid density, g/cm³；RHOB_shIs mudstone density in g/cm³(ii) a NPHI refers to the value of a sampling point on the neutron log; NPHI_maThe neutron value and decimal number of the rock skeleton; NPHI_fIn the formation fluidSub-value, decimal; NPHI_shIs the mudstone neutron value, decimal; vsh is the sludge content, decimal fraction.

The density of the formation fluid and the neutron value of the formation fluid are preset parameters, and the engineering application is determined according to the actual engineering environment, for example, 0.9-1.1, such as 0.98, 1, 1.02, and the like, which is not limited in the embodiment of the present invention.

PHIT＝PHIE+Vsh×PHIT_sh

Wherein PHIT represents total porosity; PHIT_ShIs the total porosity of the mudstone; RHOB_DShIs the density of dry mudstone in g/cm³；RHOB_ShIs mudstone density in g/cm³；RHOB_WIs the density of the water of the stratum in g/cm³(ii) a Vsh is the sludge content, decimal fraction.

In an alternative embodiment, referring to fig. 9, this step S230 may include the following:

step S231: carrying out sensitivity analysis on the shale content and the effective porosity to obtain a shale content upper limit value and a porosity lower limit value;

step S232: identifying lithology according to the shale content and the upper limit value of the shale content;

wherein lithology can be preliminarily identified as mudstone of sandstone according to the two limitations.

Step S233: and dividing the sandstone reservoir according to the effective porosity, the lower porosity limit value and the identified lithology.

In an optional embodiment, before step S230, the following steps may be further included:

step I: establishing a shale content-total porosity intersection graph according to the shale content and the total porosity;

step II: analyzing the shale content-total porosity intersection diagram by utilizing a Thomas Stieber model to obtain a shale distribution type of the target thin interbed, wherein the shale distribution type comprises: lamellar argillaceous, dispersive argillaceous, and structural argillaceous;

step III: and when the mud distribution type of the target thin interbed is dispersed mud or structural mud, ending the method flow.

It is worth to be noted that the method for determining the water saturation of the sandstone reservoir in the thin interbed of the sandstone of the deep water deposition system, which is provided by the embodiment of the invention, is suitable for the thin interbed with the shale distribution type of lamellar shale.

In order to make the present application better understood by those skilled in the art, the method steps provided by the embodiment of the present invention are illustrated with reference to fig. 10:

firstly, curve quality analysis and curve correction are carried out on conventional logging data, the mud content and the porosity are calculated after pretreatment, the water saturation is calculated by combining Rt-Scanner data, the CutOff value is determined by combining adjacent region parameters and the like, and in addition, quantitative fine evaluation of a hydrocarbon reservoir is realized by combining lithology based on conventional logging data identification.

For example, the method for determining the water saturation of the sandstone reservoir in the shale thin interbed of the deepwater deposition system can comprise the following steps:

(1) and loading conventional logging data of the interval to be evaluated of the research well, wherein the conventional logging data comprise a natural gamma logging curve GR, a caliper logging curve CAL, a neutron logging curve NHPI, a density logging curve RHOB, a sound wave time difference logging curve DTCO, a mud filtrate resistivity logging curve RMF, a shallow resistivity logging curve Rs, a deep resistivity logging curve Rd, a longitudinal resistivity logging curve Rv and a transverse resistivity logging curve Rh.

(2) And (3) performing quality analysis on the conventional logging data loaded in the step (1), and performing environment correction and borehole correction according to the quality analysis result to ensure the authenticity and accuracy of logging data.

(3) And (3) calculating the mud content Vsh of the stratum (thin interbed) to be evaluated by using the natural gamma logging curve corrected in the step (2).

(4) And (3) establishing a neutron-density logging intersection map by using the neutron logging curve and the density logging curve corrected in the step (2), determining the parameters of the stratum skeleton, and quantitatively calculating the stratum effective porosity PHIE and the total porosity by using the neutron logging curve, the density logging curve and the stratum shale content Vsh calculated in the step (3).

(5) And (4) carrying out sensitivity analysis on the effective porosity PHIE and the argillaceous content Vsh of the stratum to be evaluated obtained in the steps (3) and (4), determining a porosity lower limit value (Cutoff-1) and an argillaceous content upper limit value (Cutoff-2) for dividing the sandstone reservoir, and identifying the lithology according to the porosity lower limit value and the argillaceous content upper limit value.

(6) And (4) establishing a shale content-total porosity intersection diagram by using the shale content and the total porosity obtained in the steps (3) and (4), and analyzing the shale distribution type in the stratum to be evaluated by using a Thomas Stieber model.

(7) And (3) applying the effective porosity PHIE obtained in the step (4), the porosity cutoff-1 and the argillaceous content cutoff-1 value determined in the step (5) and the identified lithology to divide the effective reservoir (namely the sandstone reservoir).

(8) And (6) reading longitudinal resistivity Rv data obtained by Rt-Scanner data processing, reading the pure water layer formation resistivity Ro of the formation to be evaluated according to the depth resistivity curve and the sandstone reservoir division result, and reading the actual measurement lithoelectric parameters of the formation to be evaluated or the actual measurement lithoelectric parameters of an adjacent well if the result of the analysis is the layered argillaceous quality.

(9) And determining the water saturation of the sandstone reservoir of the sandstone-shale thin interbed based on the constraint of the formation test data and according to the parameters, and accurately and effectively evaluating the sandstone reservoir and the fluid property of the sandstone reservoir.

(10) And collecting and analyzing the upper limit value of water saturation Cutoff-3 of the adjacent region drilling distinguishing oil-gas reservoir.

(11) And establishing a logging evaluation chart according to the water saturation Sw value, carrying out comprehensive interpretation evaluation on the effective sandstone reservoir and quantitatively judging the oil-gas content of the sandstone reservoir.

It is worth noting that steps 10 and 11 serve as subsequent procedures for calculating the water saturation.

FIG. 11 is a comparison graph of the interpretation and evaluation of a gas production zone of a conventional method applied to a deep water sedimentation system A well P sand group and the method of the invention; therefore, the invention has the following characteristics that:

(1) the method is suitable for sand shale thin interbed oil and gas reservoirs of deep water sedimentation systems, and is not suitable for dispersed or structural mudstone profiles.

(2) The longitudinal resistivity is used as a main parameter for calculating the water saturation, the influence of the thin-layer argillaceous surrounding rock on the formation resistivity is eliminated to a large extent, the calculated sandstone reservoir water saturation reflects the oil-gas-containing property and state of the formation to the maximum extent, and the method has high accuracy and practicability.

By applying the water saturation calculation method established by the patent, on the basis of determining the shale in the stratum as the layer distribution by utilizing a Thomas Stieber model for the logging information of a deep water sedimentation body system of a well A in a certain work area, longitudinal resistivity data obtained by logging is adopted to obtain a graph 11, and the graph shows: the sensitivity of longitudinal resistivity to oil-gas containing of a stratum is high, the logging response of the longitudinal resistivity of an upper gas layer is 5-13.3ohmm and is 9.1ohmm on average, and the response of the longitudinal resistivity of a lower gas layer is 7.5-22.8ohmm and is 14.9ohmm on average; and the response value of conventional deep resistivity logging is that the upper gas layer is only 4.4-5.1ohmm, the average is 4.45ohmm, and the response value of conventional deep resistivity logging of the lower gas layer is 4.8-6.3ohmm, the average is 5.46 ohmm. The accuracy of the calculated water saturation is greatly improved, the water saturation upper gas layer obtained by the method is 36.3 percent, the water saturation lower gas layer is 27.3 percent, and the water saturation upper gas layer obtained by the conventional deep resistivity method is 67.1 percent and the water saturation lower gas layer is 54.5 percent. Thirdly, identifying that the gas layer is more practical, explaining that the well sections with the well depths of X503.6-X536.2 m and X536.2-X566.0 m are both gas layers, and conforming to the practical test result, wherein the practical test result is that the well sections with the well depths of X503.6-X536.2 m are used for obtaining natural gas 34MMscf produced in daily life; the lower gas layer X536.2-X566 well section is obtained by daily natural gas production 50MMscf, while the conventional well logging explains that the X503.6-X543.0 meter well section is a gas layer, and the X543.0-X566.0 meter well section is a gas layer, and the two methods are compared: the gas layer thickness is improved from conventionally explained 23 meters (effective thickness 19.6 meters) to 62.6 meters (effective thickness 59.2 meters), with an upper gas layer thickness of 32.8 meters (effective thickness 31.55 meters) and a lower gas layer thickness of 29.8 meters (effective thickness 28.65 meters), and the effective gas layer thickness is increased by nearly a factor of two. The above examples demonstrate the correctness and applicability of the present invention.

In summary, according to the method for determining the water saturation of the sandstone reservoir in the thin interbed of the deepwater sedimentary system, provided by the embodiment of the invention, the water saturation of the sandstone reservoir in the thin interbed section of the sandstone is accurately calculated by using the longitudinal resistivity data obtained by processing the Rt-Scanner data and the calculation formula of the saturation of the logging fluid of the thin interbed of the deepwater sedimentary system based on the Rt-Scanner technology innovatively established by the invention, and finally, the thin sandstone reservoir and the oil-gas layer in the deepwater sedimentary system are accurately identified and divided.

The method provided by the invention is applied to sandstone reservoir and gas reservoir evaluation of deep water sand shale thin interbed sedimentary system drilling in a certain sea area, and DST trial production data verifies that the method has higher accuracy and better applicability, thereby providing a new idea for comprehensive interpretation of thin interbed oil and gas reservoir logging.

Based on the same inventive concept, the embodiment of the application also provides a device for determining the water saturation of the sandstone reservoir in the sandstone and shale thin interbed of the deepwater sedimentary system, which can be used for realizing the method described in the embodiment, and is described in the following embodiment. The principle of solving the problems by the sandstone reservoir water saturation determining device in the deepwater sedimentary system sandstone-shale thin interbed is similar to that of the method, so the implementation of the sandstone reservoir water saturation determining device in the deepwater sedimentary system sandstone-shale thin interbed can refer to the implementation of the method, and repeated parts are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.

Fig. 12 is a first structural block diagram of an apparatus for determining water saturation of a sandstone reservoir in a sandstone thin interbed of a deepwater deposition system in an embodiment of the present invention. As shown in fig. 12, the device for determining the water saturation of the sandstone reservoir in the sandstone-shale thin interbed of the deep water sedimentary system specifically comprises: the system comprises a data acquisition module 10, a conventional logging data processing module 20, a longitudinal resistivity acquisition module 30 and a water saturation calculation module 40.

The data acquisition module 10 acquires conventional logging data and Rt-Scanner data of a target thin interbed;

specifically, the target thin interbed is used as a target interval of the scheme and comprises a sandstone reservoir, a sandstone non-reservoir and a mudstone layer; the purpose of the scheme is to determine the water saturation of a sandstone reservoir in a thin interbed.

The Rt-Scanner data are obtained by measuring by using an Rt-Scanner method, the three-component resistivity can be obtained by processing the measurement result of the Rt-Scanner method, and the Rt-Scanner method has the greatest advantage that the measurement current is forced to pass through the stratums of all layers, so that the influence of mudstone on the apparent resistivity value of the sandstone stratum is reduced to the greatest extent by the measurement result. The relationship between the mud content of the thin interbed and the measured resistivity is shown in fig. 3. FIG. 4 is a schematic diagram of the invention, the left diagram showing the path of the longitudinal resistivity measurement current through the path; the right graph is a conventional measurement current through path graph.

The conventional logging data processing module 20 processes the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters;

The longitudinal resistivity acquisition module 30 processes the Rt-Scanner data to obtain longitudinal resistivity;

And the water saturation calculation module 40 obtains the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the sandstone reservoir physical property response parameter and the longitudinal resistivity.

The longitudinal resistivity obtained by Rt-Scanner data processing eliminates the influence of the layered argillaceous surrounding rock on the resistivity of a target layer to the maximum extent, improves the sensitivity of thin-layer hydrocarbon reservoir identification, accurately calculates the water saturation or the oil gas saturation, avoids the phenomenon that the conventional well logging interpretation method omits and underestimates the thin hydrocarbon reservoir, greatly improves the interpretation precision of the thin hydrocarbon reservoir, and is suitable for the sand-shale thin interbed hydrocarbon reservoir of a deep water deposition system.

In an alternative embodiment, referring to fig. 13, the sandstone reservoir water saturation determination device in the sandstone thin interbed of the deepwater deposition system may further comprise: a quality analysis module 50 and a data correction module 60.

The quality analysis module 50 performs quality analysis on the conventional well logging data;

the data correction module 60 performs environmental correction and borehole correction on the conventional well log data according to the quality analysis result.

In an alternative embodiment, referring to fig. 14, the sandstone reservoir water saturation determination device in the sandstone thin interbed of the deepwater deposition system may further comprise: and a hydrocarbon saturation calculation module 70.

And the oil and gas saturation calculation module 70 acquires the oil and gas saturation of the sandstone reservoir according to the water saturation of the sandstone reservoir.

In an alternative embodiment, referring to FIG. 15, the conventional well log data processing module comprises: the device comprises a argillaceous content acquisition unit 21, a porosity acquisition unit 22, a sandstone reservoir division unit 23 and a resistivity acquisition unit 24.

The shale content obtaining unit 21 obtains the shale content of the target thin interbed according to the natural gamma logging curve;

the porosity obtaining unit 22 obtains effective porosity and total porosity according to the neutron logging curve, the density logging curve and the shale content;

the sandstone reservoir dividing unit 23 divides the sandstone reservoir according to the argillaceous content and the effective porosity;

the resistivity obtaining unit 24 obtains formation water resistivity and mudstone resistivity according to the deep resistivity curve, the lithology recognition result and the sandstone reservoir division result.

In an alternative embodiment, referring to fig. 16, the porosity obtaining unit includes: the neutron-density well logging intersection map building subunit 22a, the stratum skeleton parameter determining subunit 22b and the porosity operator subunit 22 c.

The neutron-density well logging crosspoint map establishing subunit 22a establishes a neutron-density well logging crosspoint map according to the neutron well logging curve and the density well logging curve;

the formation skeleton parameter determining subunit 22b determines a formation skeleton parameter according to the neutron-density logging intersection map;

the porosity operator unit 22c calculates the effective porosity and the total porosity according to the neutron log, the density log, the formation framework parameters and the shale content.

In an alternative embodiment, referring to fig. 17, the sandstone reservoir partitioning unit includes: a limit acquisition subunit 23a, a lithology identification subunit 23b, and a reservoir zoning subunit 23 c.

The limit value obtaining subunit 23a performs sensitivity analysis on the shale content and the effective porosity to obtain a shale content upper limit value and a porosity lower limit value;

the lithology identifying subunit 23b identifies lithology according to the shale content and the upper limit value of the shale content;

and the reservoir dividing unit 23c divides the sandstone reservoir according to the effective porosity, the lower porosity limit value and the identified lithology.

In an alternative embodiment, the conventional well log data processing module may further include: the device comprises a shale content-total porosity intersection map establishing unit, a shale distribution type analyzing unit and a process control unit.

a argillaceous distribution type analysis unit, which analyzes the argillaceous content-total porosity intersection map by using a Thomas Stieber model to obtain the argillaceous distribution type of the target thin interbed, wherein the argillaceous distribution type comprises: lamellar argillaceous, dispersive argillaceous, and structural argillaceous;

The apparatuses, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or implemented by a product with certain functions. A typical implementation device is an electronic device, which may be, for example, a personal computer, a laptop computer, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

In a typical example, the electronic device specifically includes a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor when executing the program performing the above-described steps of determining water saturation of a sandstone reservoir in a sandstone thin interbed of a water deposition system.

Referring now to FIG. 18, shown is a schematic diagram of an electronic device 600 suitable for use in implementing embodiments of the present application.

As shown in fig. 18, the electronic apparatus 600 includes a Central Processing Unit (CPU)601 that can perform various appropriate works and processes according to a program stored in a Read Only Memory (ROM)602 or a program loaded from a storage section 608 into a Random Access Memory (RAM)) 603. In the RAM603, various programs and data necessary for the operation of the system 600 are also stored. The CPU601, ROM602, and RAM603 are connected to each other via a bus 604. An input/output (I/O) interface 605 is also connected to bus 604.

The following components are connected to the I/O interface 605: an input portion 606 including a keyboard, a mouse, and the like; an output portion 607 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage section 608 including a hard disk and the like; and a communication section 609 including a network interface card such as a LAN card, a modem, or the like. The communication section 609 performs communication processing via a network such as the internet. The driver 610 is also connected to the I/O interface 605 as needed. A removable medium 611 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 610 as necessary, so that a computer program read out therefrom is mounted as necessary on the storage section 608.

In particular, according to an embodiment of the present invention, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, an embodiment of the invention includes a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the above-described method for determining water saturation of a sandstone reservoir in a sandstone thin interbed of a water deposition system.

In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 609, and/or installed from the removable medium 611.

Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, 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 disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.

For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.

As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.

The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims

1. A method for determining the water saturation of a sandstone reservoir in a sandstone-shale thin interbed of a deepwater deposition system is characterized by comprising the following steps:

processing the Rt-Scanner data to obtain longitudinal resistivity;

2. The method for determining the water saturation of a sandstone reservoir in a shale thin interbed of a deepwater sedimentary system according to claim 1, wherein the logging response characteristic parameters comprise: formation water resistivity, mudstone resistivity; the sandstone reservoir physical property response parameters comprise: argillaceous content and effective porosity.

3. The method for determining the water saturation of a sandstone reservoir in a shale thin interbed of a deepwater sedimentary system according to claim 2, wherein the conventional well logging data comprises: a natural gamma log, a neutron log, a density log, and a deep resistivity curve;

the processing of the conventional logging data to obtain logging response characteristic parameters and sandstone reservoir physical property response parameters comprises the following steps:

obtaining the mud content of a target thin interbed according to the natural gamma logging curve;

4. The method for determining the water saturation of a sandstone reservoir in a shale thin interbed of a deepwater sedimentary system according to claim 3, wherein the obtaining of the effective porosity and the total porosity according to the neutron log, the density log and the shale content comprises:

establishing a neutron-density logging intersection map according to the neutron logging curve and the density logging curve;

and calculating the effective porosity and the total porosity according to the neutron logging curve, the density logging curve, the stratum skeleton parameters and the shale content.

5. The method for determining the water saturation of the sandstone reservoir in the thin interbed of the sandstone in the deepwater deposition system according to claim 3, wherein the dividing the sandstone reservoir according to the shale content and the effective porosity comprises:

6. The method for determining the water saturation of the sandstone reservoir in the thin interbed of the sandstone in the deepwater deposition system according to claim 3, wherein before the step of dividing the sandstone reservoir according to the shale content and the effective porosity, the method further comprises the following steps of:

7. The method for determining the water saturation of the sandstone reservoir in the sandstone-shale thin interbed of the deepwater sedimentary system according to claim 1, wherein before the conventional logging data are processed to obtain the logging response characteristic parameter and the sandstone reservoir physical property response parameter, the method further comprises the following steps:

performing quality analysis on the conventional well logging data;

8. The method for determining the water saturation of the sandstone reservoir in the sandstone thin interbed of the deepwater deposition system according to any of claims 1 to 7, further comprising the following steps:

9. A device for determining the water saturation of a sandstone reservoir in a thin interbed of sandstone in a deepwater deposition system is characterized by comprising:

and the water saturation calculation module is used for obtaining the water saturation of the sandstone reservoir according to the logging response characteristic parameter, the sandstone reservoir physical property response parameter and the longitudinal resistivity.

10. The device for determining the water saturation of a sandstone reservoir in a shale thin interbed of a deepwater sedimentary system according to claim 9, wherein the logging response characteristic parameters comprise: formation water resistivity, mudstone resistivity; the sandstone reservoir physical property response parameters comprise: argillaceous content and effective porosity.

11. The apparatus for determining water saturation of a sandstone reservoir in a thin interbed of sandstone in a deepwater deposition system according to claim 9, wherein the conventional well log data comprises: a natural gamma log, a neutron log, a density log, and a deep resistivity curve;

the conventional well logging data processing module comprises:

and the resistivity obtaining unit is used for obtaining the formation water resistivity and the mudstone resistivity according to the deep resistivity curve, the lithology recognition result and the sandstone reservoir division result.

12. The device for determining the water saturation of the sandstone reservoir in the thin interbed of sandstone in a deepwater sedimentary system according to claim 11, wherein the porosity acquiring unit comprises:

the neutron-density logging intersection graph establishing subunit establishes a neutron-density logging intersection graph according to the neutron logging curve and the density logging curve;

the stratum skeleton parameter determining subunit determines stratum skeleton parameters according to the neutron-density logging intersection map;

and the porosity operator unit is used for calculating the effective porosity and the total porosity according to the neutron logging curve, the density logging curve, the stratum skeleton parameters and the shale content.

13. The device for determining the water saturation of the sandstone reservoir in the thin interbed of sandstone and mudstone of a deepwater sedimentary system according to claim 11, wherein the sandstone reservoir dividing unit comprises:

the limit value acquisition subunit is used for carrying out sensitivity analysis on the shale content and the effective porosity to obtain a shale content upper limit value and a porosity lower limit value;

a lithology identification subunit for identifying lithology according to the shale content and the upper limit value of the shale content;

and the reservoir dividing subunit is used for dividing the sandstone reservoir according to the effective porosity, the lower porosity limit value and the identified lithology.

14. The device for determining the water saturation of a sandstone reservoir in a thin interbed of sandstone in a deepwater sedimentary system according to claim 11, wherein the conventional logging data processing module further comprises:

15. The device for determining the water saturation of a sandstone reservoir in a thin interbed of sandstone in a deepwater deposition system according to claim 9, further comprising:

16. The sandstone reservoir water saturation determination device in a sandstone thin interbed of a deepwater deposition system according to any of claims 9 to 15, further comprising:

17. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the method for determining water saturation of a sandstone reservoir in a sandstone thin interbed of a deepwater deposition system according to any of claims 1 to 8.

18. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method for determining the water saturation of a sandstone reservoir in a sandstone thin interbed of a deepwater sedimentary system according to any of claims 1 to 8.