Movatterモバイル変換


[0]ホーム

URL:


CN113640178B - Formation water type identification method, pollution rate determination method and calculation device - Google Patents

Formation water type identification method, pollution rate determination method and calculation device
Download PDF

Info

Publication number
CN113640178B
CN113640178BCN202110912286.0ACN202110912286ACN113640178BCN 113640178 BCN113640178 BCN 113640178BCN 202110912286 ACN202110912286 ACN 202110912286ACN 113640178 BCN113640178 BCN 113640178B
Authority
CN
China
Prior art keywords
density
water
formation water
pumping
fluid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110912286.0A
Other languages
Chinese (zh)
Other versions
CN113640178A (en
Inventor
周明高
冯永仁
左有祥
孔笋
孟悦新
褚晓冬
支宏旭
周艳敏
杨玉卿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Oilfield Services Ltd
Original Assignee
China Oilfield Services Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Oilfield Services LtdfiledCriticalChina Oilfield Services Ltd
Priority to CN202110912286.0ApriorityCriticalpatent/CN113640178B/en
Publication of CN113640178ApublicationCriticalpatent/CN113640178A/en
Application grantedgrantedCritical
Publication of CN113640178BpublicationCriticalpatent/CN113640178B/en
Activelegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Classifications

Landscapes

Abstract

The invention discloses a stratum water type identification method, a pollution rate determination method and computing equipment. The formation water type identification method comprises the steps of collecting pumping parameters and actual fluid density measurement values corresponding to a plurality of moments in a preset time window, generating a fluid density change curve according to the collected pumping parameters and the actual fluid density measurement values, and identifying the formation water type based on the fluid density change curve. By adopting the scheme, the stratum water type can be monitored in real time without multiple times of underground sampling to the ground for analysis. On one hand, the recognition efficiency of the formation water type is improved, a basis is provided for determining the sampling time of the formation water, and on the other hand, the sampling cost is saved and the operation risk is reduced.

Description

Formation water type identification method, pollution rate determination method and computing equipment
Technical Field
The invention relates to the field of stratum sampling of logging technology, in particular to a stratum water type identification method, a pollution rate determination method and computing equipment.
Background
The composition of the formation water has important significance for oil field exploration and production. In order to determine the composition of formation water, it is common practice to place a sampling tool downhole at a predetermined depth to obtain a sample of formation water. However, during formation sampling, the collected formation water samples are of different types, such as formation water layers, oil or gas layers, and oil or gas water layers, and the like.
In order to be able to identify the type of formation water, the prior art uses a multiple sampling method. The method comprises the steps of taking different formation water samples from underground at different moments, taking the formation water samples back to a ground system, carrying out water analysis on the obtained formation water samples, and then determining the types of the formation water samples. However, this approach of the prior art can lead to inefficient formation water type identification and increase the cost and risk of the overall exploration process.
Disclosure of Invention
The present invention has been made in view of the above problems, and has as its object to provide a formation water type identification method, a pollution rate determination method and a computing device which overcome or at least partially solve the above problems.
According to one aspect of the present invention, there is provided a formation water type identification method comprising:
collecting pumping parameters and actual fluid density measured values corresponding to a plurality of moments in a preset time window;
generating a fluid density change curve according to the collected pumping parameters and the actual fluid density measured values;
based on the fluid density profile, a formation water type is identified.
In an alternative embodiment, the identifying the formation water type based on the fluid density profile further comprises:
Judging whether a slug flow section exists in the fluid density change curve;
if yes, extracting a plurality of density maxima from a slug flow section of the fluid density change curve;
Generating a density maximum change curve according to the density maxima and pumping parameters corresponding to each density maximum;
and identifying the stratum water type according to the density maximum change curve.
In an alternative embodiment, the identifying the formation water type based on the density maxima profile further comprises:
If the smoothness of the density maximum change curve is greater than a preset threshold, determining that the stratum water type is the same layer of oil and water or the same layer of gas and water;
and if the smoothness of the density maximum change curve is smaller than or equal to a preset threshold value, determining that the stratum water type is an oil layer or a gas layer.
In an alternative embodiment, if no slugging section exists in the fluid density profile, the formation water type is determined to be a formation water layer.
In an alternative embodiment, the slug flow section is a section with regular oscillation of density values.
In an alternative embodiment, the pumping parameters and the actual fluid density measurements are acquired by a preset sampling device;
the preset sampling equipment at least comprises a pumping module and a density sensor, wherein the density sensor is positioned at the downstream of the pumping module.
According to another aspect of the present invention, there is provided a formation water pollution rate determining method, including:
identifying the formation water type by adopting the formation water type identification method;
if the formation water type is identified as the oil-water layer or the gas-water layer, extracting a plurality of density maxima of a plug flow section in the middle section of the fluid density change curve, and determining pumping parameters corresponding to each density maximum;
Calculating the density of pure stratum water according to the pumping parameters corresponding to the density maxima;
and calculating the pollution rate of the formation water according to the pure formation water density, the pure water-based drilling fluid filtrate density and the currently acquired actual fluid measurement density.
In an alternative embodiment, the calculating the water density of the pure formation according to the plurality of density maxima and the pumping parameters corresponding to each density maxima further comprises:
Taking the density maxima and pumping parameters corresponding to the density maxima as sample data;
And training a preset neural network model through the sample data to obtain the association relationship between the actual fluid density measured value and the pumping parameter of the stratum water type oil-water same layer or gas-water same layer.
And determining the density value of the pure formation water according to the preset pumping parameter and the association relation, wherein the preset pumping parameter is the pumping parameter when the pure formation water is expected to be collected.
According to yet another aspect of the present invention, there is provided a computing device comprising a processor, a memory, a communication interface and a communication bus, the processor, the memory and the communication interface completing communication with each other through the communication bus;
the storage is used for storing at least one executable instruction, and the executable instruction enables the processor to execute the operation corresponding to the formation water type identification method;
and/or the executable instructions enable the processor to execute the operations corresponding to the formation water type identification method.
According to still another aspect of the present invention, there is provided a computer storage medium having at least one executable instruction stored therein, the executable instruction causing the processor to perform operations corresponding to the above-mentioned formation water type identification method;
and/or the executable instructions enable the processor to execute the operations corresponding to the formation water type identification method.
The invention discloses a stratum water type identification method, a pollution rate determination method and computing equipment. And acquiring pumping parameters and actual fluid density measured values corresponding to a plurality of moments in a preset time window, generating a fluid density change curve according to the acquired pumping parameters and the actual fluid density measured values, and identifying the formation water type based on the fluid density change curve, so that the formation water type can be monitored in real time without repeatedly sampling from underground to the ground for analysis. On one hand, the identification efficiency of the formation water type is improved, a basis is provided for the sampling time of the formation water, and on the other hand, the sampling cost is saved and the operation risk is reduced.
And after the formation water type is identified as the oil-water layer or the gas-water layer, extracting a plurality of density maxima of a plug flow section in the middle section of the fluid density change curve, determining pumping parameters corresponding to the density maxima, calculating the pure formation water density according to the density maxima and the pumping parameters corresponding to the density maxima, and calculating the formation water pollution rate according to the pure formation water density, the pure water-based drilling fluid filtrate density and the currently acquired actual fluid measurement density, thereby accurately calculating the formation water pollution rate and providing a basis for sampling and analyzing the formation water.
The foregoing description is only an overview of the present invention, and is intended to be implemented in accordance with the teachings of the present invention in order that the same may be more clearly understood and to make the same and other objects, features and advantages of the present invention more readily apparent.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to designate like parts throughout the figures. In the drawings:
FIG. 1 shows a schematic flow chart of a method for identifying formation water types, which is provided by an embodiment of the invention;
FIG. 2a is a schematic diagram showing data acquisition of a preset sampling device for extracting fluid from top to bottom according to an embodiment of the present invention;
FIG. 2b is a schematic diagram showing data acquisition of a preset sampling device for extracting fluid from bottom to top according to an embodiment of the present invention;
FIG. 3 is a schematic flow chart for identifying formation water type based on a fluid density change curve according to an embodiment of the present invention;
FIG. 4 shows a schematic diagram of a fluid density change curve of an oil-water layer or a gas-water layer according to an embodiment of the present invention;
FIG. 5a is a schematic view showing a fluid density change curve of an oil or gas reservoir according to an embodiment of the present invention;
FIG. 5b shows a schematic representation of the fluid density profile of yet another oil or gas reservoir provided by an embodiment of the present invention;
FIG. 6 is a schematic diagram showing a fluid density change curve of a formation water layer according to an embodiment of the present invention;
FIG. 7 is a schematic flow chart of a method for determining the pollution rate of formation water according to an embodiment of the present invention;
FIG. 8 is a schematic structural diagram of a neural network model according to an embodiment of the present invention;
FIG. 9 is a schematic diagram showing the effect of density maxima predicted by a neural network model according to an embodiment of the present invention;
FIG. 10 is a schematic diagram showing the effect of a calculated formation water pollution rate provided by an embodiment of the present invention;
FIG. 11 illustrates a schematic diagram of a computing device provided by an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Fig. 1 shows a schematic flow chart of a formation water type identification method according to an embodiment of the present invention. Specifically, the method for identifying the type of the formation water provided by the embodiment is applied to a logging scene of water-based mud.
As shown in fig. 1, the method comprises the steps of:
step S110, collecting pumping parameters and actual fluid density measurement values corresponding to a plurality of moments within a preset time window.
The pumping parameters and actual fluid density measurements are acquired by a preset sampling device. The preset sampling device at least comprises a pumping module and a density sensor. In addition, the preset sampling device may also include a probe, a sample bottle, and the like.
Fig. 2a and 2b show schematic data collection diagrams of a preset sampling device for fluid extraction from top to bottom and fluid extraction from bottom to top, respectively. As can be seen from fig. 2a and 2b, the density sensor in the preset sampling device is located downstream of the pumping module. With this structure, the pumping module can separate oil-water, gas-water or gas-oil-water under gravity according to the density difference of the fluid, thereby forming a slug flow. The density sensor is arranged downstream of the pumping module, and the change can be accurately measured.
In an alternative embodiment, the predetermined sampling device may be a cable formation pressure measurement sampler (EFDT). The data acquisition process is that before logging, the cable stratum pressure measuring and sampling instrument is placed at the underground target depth, and the probe seat is sealed on the well wall. After the seat seal is successful, a pumping module is started, stratum fluid enters the pipeline through the suction port, the actual density value of the fluid is measured in real time through a density sensor, and the density value measured by the density sensor is the actual fluid density measured value. The pumping module will pump parameters and the density sensor will upload actual fluid density measurements to the surface system via remote transmission of the cable.
In another alternative embodiment, the predetermined sampling device may be a sampling while drilling Instrument (IFSA). The data acquisition process is that before logging, the sampling instrument while drilling is placed to the underground target depth, the ground system is communicated with the underground mud transmission device through the mud transmission device, and the underground mud transmission device sends a ground command to the sampling instrument while drilling. And then, the sampling instrument while drilling seats the probe on the well wall, after the seat sealing is successful, the pumping module is started, stratum fluid enters the pipeline through the suction port, the actual density value of the fluid is measured in real time through the density sensor, and the density value measured by the density sensor is the actual fluid density measured value. The pumping module uploads the pumping parameters and the density sensor to the surface system in real time via the mud delivery device.
Further, the data collected by the preset sampling device is divided at different moments, i.e. each moment has a corresponding pumping parameter and an actual fluid density measurement. The pumping parameters are in particular pumping time or pumping volume, etc.
It can be seen that, unlike the prior art, the present embodiment does not send samples collected at a plurality of moments to the surface platform for analysis, but collects pumping parameters and actual fluid density measurement values at a plurality of moments by the preset sampling device, and sends the collected data back to the surface system for analysis. Therefore, in the prior art, the sample is collected and transmitted, and the identification of the formation water type can be realized only by collecting and transmitting data in the pumping process, so that the overall identification efficiency is improved.
Step S120, generating a fluid density variation curve according to the collected pumping parameters and the actual fluid density measurement.
And taking the pumping parameter as an independent variable, and taking the actual fluid density measured value as the dependent variable to construct a fluid density and pumping parameter change curve, namely generating a fluid density change curve. For example, a fluid density-pump time curve or a fluid density-pump volume curve, etc. may be generated.
Step S130, identifying the formation water type based on the fluid density change curve.
The fluid density profile implies formation water type characteristics. The fluid density profile varies from formation water type to formation water type. The formation water type can be accurately identified through analysis of the fluid density change curve. The stratum water type comprises stratum water layer, oil layer or gas layer, oil-water layer or gas-water layer and the like.
In an alternative embodiment, to improve the accuracy of identifying the formation water type, the formation water type may be specifically identified through the steps shown in fig. 3, and in particular, as shown in fig. 3, the formation water type may be identified through the following steps S131 to S136 based on the fluid density change curve.
Step S131, judging whether a slug flow section exists in the fluid density change curve, if so, executing step S132, and if not, executing step S136.
Because of the different miscibility of water-based drilling fluid filtrate with formation water and oil or gas, slugging can occur in some formation water types. Specifically, the water-based drilling fluid filtrate is completely miscible with formation water, so that slugging is not formed when the formation water type is a formation water layer, while the water-based drilling fluid filtrate is completely immiscible with oil or gas, so that slugging is formed in the formation water type is an oil-water or gas-water layer and an oil or gas layer. Based on this feature, the present embodiment first determines whether a slug flow exists based on the fluid density change curve.
Specifically, if a slugging flow exists, the actual fluid density measurement will regularly fall and rise with pumping parameters. The fluid density change curve is embodied as a section with regular oscillation of a density value, and the section with regular oscillation of the density value in the fluid density change curve is a slug flow section. In order to avoid that the error fluctuation is considered as a slug flow section, the difference between two adjacent maximum values and minimum values in the slug flow section in the embodiment is greater than or equal to a preset density value.
Step S132, extracting a plurality of density maxima from the slug flow section of the fluid density variation curve, and generating a density maximum variation curve according to the plurality of density maxima and pumping parameters corresponding to each density maxima.
Because the formation water type is oil-water or gas-water and oil-or gas-layer, and the slug flow is formed in the oil-or gas-layer, after determining that the slug flow exists based on the fluid density change curve, the oil-water or gas-water and oil-or gas-layer is further identified through step S132 and step S133.
Specifically, although the oil-water layer or the gas-water layer and the oil layer or the gas layer both form a slug flow, the oil-water layer or the gas-water layer and the oil layer or the gas layer have different meanings of the maximum value of the plug flow section in the middle section of the fluid density change curve. The maximum value of the plug flow section in the middle section of the fluid density change curve in the oil-water layer or the gas-water layer corresponds to the density value of the water-based drilling fluid filtrate and the stratum water mixture, which is smoother along with the change of the pumping parameter, and the maximum value of the plug flow section in the middle section of the fluid density change curve in the oil layer or the gas layer corresponds to the density value of the water-based drilling fluid filtrate, which is more severe along with the change of the pumping parameter.
Based on the above, the step firstly extracts a plurality of density maxima from the plug flow section in the middle section of the fluid density change curve, and determines pumping parameters corresponding to each extracted density maxima. And finally, generating a density maximum change curve according to the plurality of density maxima and pumping parameters corresponding to each density maximum. The independent variable in the density maximum change curve is pumping parameter, and the dependent variable is density maximum. And then identifying the formation water type according to the density maximum change curve.
Step S133, judging whether the smoothness of the density maximum change curve is larger than a preset threshold, if so, executing step S134, and if not, executing step S135.
The smoothness of the density maximum change curve shows the speed of the density maximum along with the change of the pumping parameter, namely, the faster the density maximum along with the change of the pumping parameter, the lower the smoothness of the density maximum change curve, otherwise, the slower the density maximum along with the change of the pumping parameter, the higher the smoothness of the density maximum change curve.
Step S134, determining that the stratum water type is the same water-oil layer or the same gas-water layer.
If the slug flow exists based on the fluid density change curve, and the smoothness of the density maximum change curve is larger than a preset threshold value, determining that the stratum water type is the oil-water same layer or the gas-water same layer.
FIG. 4 shows a schematic diagram of a fluid density change curve of an oil-water layer or a gas-water layer. As shown in fig. 4, because the water-based drilling fluid filtrate is completely miscible with formation water and completely immiscible with oil or gas in the same layer of oil and water or gas and water, there may be a slug flow section in the fluid density profile of the same layer of oil and water. And the frequency of the slug flow is relatively stable. And in the fluid density change curve of the oil-water layer or the gas-water layer, the density maximum value corresponds to the density value of the water-based drilling fluid filtrate and the stratum water mixture. The density maxima drop smoothly and slightly with changes in pumping time or pumping volume. Specifically, the density maximum value is changed from a density value mainly comprising water-based drilling fluid filtrate at the early stage of pumping to a density value mainly comprising formation water relatively smoothly, and finally a relatively stable density value is achieved.
In step S135, the formation water type is determined to be an oil or gas reservoir.
And if the slug flow exists based on the fluid density change curve, and the smoothness of the density maximum change curve is smaller than or equal to a preset threshold value, determining that the stratum water type is an oil layer or a gas layer.
Fig. 5a shows a schematic representation of the fluid density profile of an oil or gas reservoir. As shown in fig. 5a, in an oil layer or a gas layer, since the water-based drilling fluid filtrate is completely insoluble with oil or gas, a shock section corresponding to a slug flow exists in a fluid density change curve of the same layer of oil and water. And the fluid density change curve can fluctuate dramatically up and down the oscillation section. Wherein, the minimum value of the fluid density change curve corresponds to the density of oil or gas, and the minimum value is stable. The maximum value of the fluid density profile corresponds to the density of the water-based drilling fluid filtrate, which has a large change over time, i.e. the smoothness of the density maximum profile is low. When the water-based drilling fluid filtrate invasion zone is small or the pumping time is long, the frequency of slugging becomes smaller and smaller as the pumping time or volume increases, the density maxima slowly thin and disappear, and finally only the density value of stable oil or gas is obtained. Furthermore, if the water-based drilling fluid filtrate invasion zone requires a long time to pump the water-based drilling fluid filtrate, the maximum density is a stable value for a long period of time in the early stage of pumping, as shown in fig. 5 b.
In step S136, the formation water type is determined to be the formation water layer.
And if the slug flow is not determined to exist based on the fluid density change curve, determining the formation water type as a formation water layer. As shown in fig. 6, no slugging sections exist in the formation water type fluid density profile due to the complete miscibility of the water-based drilling fluid filtrate with the formation water. The density value in the fluid density change curve of the formation water type changes smoothly with the pumping time or the volume, in particular, the density value of the water-based drilling fluid filtrate is changed smoothly to the density value of the formation water.
Therefore, in this embodiment, the pumping parameters and the actual fluid density measurement values corresponding to a plurality of moments within the preset time window are collected, the fluid density change curve is generated according to the collected pumping parameters and the actual fluid density measurement values, and finally the formation water type is identified based on the fluid density change curve. By adopting the scheme, the stratum water type can be monitored in real time without multiple times of underground sampling to the ground for analysis. On one hand, the identification efficiency of the formation water type is improved, a basis is provided for the sampling time of the formation water, and on the other hand, the sampling cost is saved and the operation risk is reduced.
Fig. 7 shows a schematic flow chart of a method for determining a pollution rate of formation water according to an embodiment of the present invention.
When the water-based drilling fluid filtrate used in the drilling process invades into the stratum deeply, the water-based drilling fluid filtrate and the original stratum water are mutually dissolved to cause stratum pollution, and the stratum water pollution rate is calculated through the follow-up steps in the embodiment, so that a foundation is provided for stratum water sampling time, stratum water component analysis and the like.
As shown in fig. 7, the method includes the steps of:
step S710, identify formation water type.
The specific implementation process of this step may refer to the description in the embodiment shown in fig. 1, and will not be repeated here.
Step S720, if the formation water type is identified as the oil-water layer or the gas-water layer, extracting a plurality of density maxima of the plug flow section in the middle section of the fluid density change curve, and determining pumping parameters corresponding to each density maximum.
If the formation water type is identified as the oil-water or gas-water layer, a fluid density change curve is generated according to the collected pumping parameters and the actual fluid density measurement value according to the method in the embodiment shown in fig. 1, and a plurality of density maxima of a plug flow section in the middle section of the fluid density change curve are extracted. And obtaining pumping parameters corresponding to any density maximum value. Wherein the density maxima correspond to the same time as the corresponding pumping parameters.
And step S730, calculating the pure formation water density according to the pumping parameters corresponding to the density maxima.
Specifically, a preset neural network model is firstly built in advance, a plurality of density maxima and pumping parameters corresponding to the density maxima are used as sample data, the preset neural network model is trained through the sample data, and the association relation between the actual fluid density measured value of the stratum water type of the oil-water same layer or the gas-water same layer and the pumping parameters is obtained. The association relation specifically refers to the association relation between the density maximum value and the pumping parameter in the slug flow section. It should be understood herein that other models such as a meditation function, a logarithmic function, an exponential function, and a hyperbolic function may be used to obtain the association relationship, which is not limited in this embodiment.
In an alternative training mode, the pumping parameters may be used as input data to the neural network model and the density maxima may be used as output data from the neural network model. The neural network model may be an ANN model. As shown in fig. 8, the pumping volume v= [ V1,V2,…,Vn ] or the pumping time t= [ t1,t2,…,tn ] is taken as an input variable of the neural network model input layer, the model parameters are assumed to be W and b in the hidden layer (the hidden layer may be a single layer or multiple layers), the output variable of the neural network is a=g (wp+b), where p is the input variable and g is a conversion function such as Sigmoid function or other functions. The output layer predicts the density ρ= [ ρ12,…,ρn ] by using the neural network model. In this prediction process, a loss function is calculatedThe parameters W and b are adjusted by gradient descent or quasi-newton to minimize the loss function globally, resulting in an optimal model ρ=f (V, W, b) or ρ=f (t, W, b). After the trained neural network model is obtained by adopting the training mode, the preset pumping parameters are used as input data, and the output data corresponding to the neural network model is the pure formation water density.
Fig. 9 shows an effect diagram of the density maxima predicted by the neural network model, and it can be seen from the figure that the density maxima predicted by the neural network model have a high matching degree with the density maxima actually measured obtained by using the EFDT, so that the neural network model in the embodiment has a higher accuracy.
In yet another alternative training mode, the data of the variation of the density maximum value along with the pumping parameter can be used as the input data of the neural network model, and the relationship between the fluid density value and the pumping parameter can be used as the output data of the neural network model. After the trained neural network model is obtained by adopting the training mode, the preset pumping parameters are used as input data, and the association relation between the fluid density value output by the neural network model and the pumping parameters is obtained. And then obtaining the pure stratum water density according to the preset pumping parameters and the association relation.
Further optionally, when the pumping parameter is pumping volume, the preset pumping parameter is pumping volume greater than a preset threshold, and when the pumping parameter is pumping time, the preset pumping parameter is pumping time greater than the preset threshold. The preset threshold may be infinite or a fixed value that is preset. For example, when the pumping parameters are approaching infinity, the training model has a limit value, the preset threshold value is infinity, and when the pumping parameters are approaching infinity, the training model has no significant limit value, the preset threshold pumping volume is 500 liters or the pumping time is the time required to calculate 500 liters of fluid to be pumped based on the pumping rate. The preset threshold may be adjusted based on formation and fluid properties and probe type, for example, the preset threshold may be smaller (e.g., 50 liters) when formation drilling fluid invasion is small and the preset threshold may be larger (e.g., 1000 liters) when formation drilling fluid invasion is large.
In addition, if the formation water type is identified as a formation water layer, pumping parameters corresponding to a plurality of moments in a preset time window and actual fluid density measurement values are used as sample data, a preset neural network model is trained through the sample data, the association relation between the actual fluid density measurement values and the pumping parameters of the formation water type as the formation water layer is obtained, and the density value of the pure formation water is determined according to the preset pumping parameters and the association relation. The model training and predicting process of the formation water layer with the formation water type can refer to the description of the corresponding part in the oil-water layer or the gas-water layer with the formation water type, and the description is omitted herein.
Step S740, calculating the pollution rate of the formation water according to the density of the pure formation water, the density of the filtrate of the pure water-based drilling fluid and the density of the currently collected actual fluid measurement.
Wherein the pure water based drilling fluid filtrate density is determined by at least one of:
the first determination mode is that when pumping is started, the fluid density is measured as the filtrate density of the pure water-based drilling fluid;
Obtaining a drilling fluid filter cake from underground to the ground, extruding filtrate, measuring the density of the fluid by using an EFDT (electronic differential pressure transducer) measuring instrument, and correcting the density to the density at the temperature and the pressure corresponding to the acquired stratum to be used as the density of the filtrate of the pure water-based drilling fluid;
and in the third determination mode, during measurement, calculating the density of the water-based drilling fluid filtrate during measurement as the density of the pure water-based drilling fluid filtrate according to the preset corresponding relation between the density of various pure drilling fluids and the temperature and pressure, the type of the drilling fluid and the corresponding temperature and pressure of the stratum during measurement.
After obtaining the pure formation water density, the pure water-based drilling fluid filtrate density, and the actual fluid measurement density currently acquired, the formation water pollution rate is obtained by the following equation 7-1:
In the formula, epsilon is the pollution rate of the filtrate of the water-based drilling fluid in the formation water, rho is the current acquired actual fluid measurement density, rhomf is the filtrate density of the pure water-based drilling fluid, and rhofw is the density of the pure formation water.
FIG. 10 shows a schematic of the effect of a calculated formation water pollution rate. As can be seen from fig. 10, the formation water pollution rate calculated in this embodiment has a high matching degree with the formation water pollution rate obtained by the conventional EFDT multiple sampling method, so that the accuracy of the formation water pollution rate calculated in this embodiment is high.
Therefore, the embodiment adopts the formation water type identification method to identify the formation water type, if the formation water type is identified as the oil-water layer or the gas-water layer, a plurality of density maxima of a plug flow section in a fluid density change curve are extracted, pumping parameters corresponding to the density maxima are determined, pure formation water density is calculated according to the density maxima and the pumping parameters corresponding to the density maxima, and the formation water pollution rate is calculated according to the pure formation water density, the pure water-based drilling fluid filtrate density and the currently acquired actual fluid measurement density. By adopting the scheme, the pollution rate of the formation water can be accurately calculated, and a foundation is provided for sampling the formation water and analyzing the composition.
FIG. 11 illustrates a schematic diagram of a computing device provided by an embodiment of the present invention. The specific embodiments of the present invention are not limited to a particular implementation of a computing device.
As shown in FIG. 11, the computing device may include a processor 1102, a communication interface (Communications Interface) 1104, memory 1106, and a communication bus 1108.
Wherein the processor 1102, the communication interface 1104, and the memory 1106 communicate with each other via a communication bus 1108. A communication interface 1104 for communicating with network elements of other devices, such as clients or other servers. The processor 1102 is configured to execute the program 1110, and may specifically perform relevant steps in the foregoing embodiment of the method for identifying a formation water type and/or the embodiment of the method for determining a contamination rate of formation water.
In particular, program 1110 may include program code including computer-operating instructions.
The processor 1102 may be a central processing unit CPU, or an Application-specific integrated Circuit ASIC (Application SPECIFIC INTEGRATED Circuit), or one or more integrated circuits configured to implement embodiments of the present invention. The computing device may include one or more processors of the same type, such as one or more CPUs, or of different types, such as one or more CPUs and one or more ASICs.
Memory 1106 for storing program 1110. The memory 1106 may include high-speed RAM memory or may also include non-volatile memory (non-volatile memory), such as at least one disk memory. The specific implementation process may refer to the description in the corresponding method embodiment, and will not be described herein.
The embodiment of the invention provides a non-volatile computer storage medium, which stores at least one executable instruction, and the computer executable instruction can execute the formation water type identification method and/or the formation water pollution rate determination method in any of the method embodiments. The specific implementation process may refer to the description in the corresponding method embodiment, and will not be described herein.
The algorithms or displays presented herein are not inherently related to any particular computer, virtual system, or other apparatus. Various general-purpose systems may also be used with the teachings herein. The required structure for a construction of such a system is apparent from the description above. In addition, embodiments of the present invention are not directed to any particular programming language. It will be appreciated that the teachings of the present invention described herein may be implemented in a variety of programming languages, and the above description of specific languages is provided for disclosure of enablement and best mode of the present invention.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the above description of exemplary embodiments of the invention, various features of the embodiments of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. This method of disclosure, however, is not to be interpreted as reflecting an intention that the claimed invention requires more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Those skilled in the art will appreciate that the modules in the apparatus of the embodiments may be adaptively changed and disposed in one or more apparatuses different from the embodiments. The modules or units or components of the embodiments may be combined into one module or unit or component and, furthermore, they may be divided into a plurality of sub-modules or sub-units or sub-components. Any combination of all features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or units of any method or apparatus so disclosed, may be used in combination, except insofar as at least some of such features and/or processes or units are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings), may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
Furthermore, those skilled in the art will appreciate that while some embodiments herein include some features but not others included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments can be used in any combination.
Various component embodiments of the invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that some or all of the functionality of some or all of the components according to embodiments of the present invention may be implemented in practice using a microprocessor or Digital Signal Processor (DSP). The present invention can also be implemented as an apparatus or device program (e.g., a computer program and a computer program product) for performing a portion or all of the methods described herein. Such a program embodying the present invention may be stored on a computer readable medium, or may have the form of one or more signals. Such signals may be downloaded from an internet website, provided on a carrier signal, or provided in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the words first, second, third, etc. do not denote any order. These words may be interpreted as names. The steps in the above embodiments should not be construed as limiting the order of execution unless specifically stated.

Claims (6)

Translated fromChinese
1.一种地层水污染率确定方法,其特征在于,包括:1. A method for determining formation water contamination rate, characterized by comprising:采集预设时间窗口内多个时刻对应的泵抽参数以及实际流体密度测量值;Collecting pumping parameters and actual fluid density measurements corresponding to multiple moments within a preset time window;根据采集到的所述泵抽参数以及所述实际流体密度测量值,生成流体密度变化曲线;Generating a fluid density variation curve according to the collected pumping parameters and the actual fluid density measurement value;从所述流体密度变化曲线的段塞流区段中提取出多个密度极大值,根据多个密度极大值以及每个密度极大值对应的泵抽参数生成密度极大值变化曲线,若密度极大值变化曲线的平滑度大于预设阈值,则确定地层水类型为油水同层或气水同层;若密度极大值变化曲线的平滑度小于或等于预设阈值,则确定地层水类型为油层或气层;其中,在油水同层或气水同层的流体密度变化曲线中,密度极大值对应于水基钻井液滤液和地层水混合物的密度值,密度极大值由泵抽早期时以水基钻井液滤液为主的密度值平滑地变化到以地层水为主的密度值;而油层或气层中流体密度变化曲线中密度极大值对应的是水基钻井液滤液的密度值,其随泵抽参数的变化剧烈,平滑度低,随着泵抽时间或体积的增加,段塞流的频率变得愈来愈小,密度极大值变稀并消失,最终只有稳定的油或气的密度值;Extract multiple density maxima from the slug flow section of the fluid density change curve, generate a density maximum change curve according to the multiple density maxima and the pumping parameters corresponding to each density maximum, if the smoothness of the density maximum change curve is greater than a preset threshold, determine the type of formation water to be oil-water in the same layer or gas-water in the same layer; if the smoothness of the density maximum change curve is less than or equal to the preset threshold, determine the type of formation water to be an oil layer or a gas layer; wherein, in the fluid density change curve of the oil-water in the same layer or the gas-water in the same layer, the density maximum corresponds to the density value of the mixture of water-based drilling fluid filtrate and formation water, and the density maximum changes smoothly from the density value dominated by water-based drilling fluid filtrate in the early stage of pumping to the density value dominated by formation water; and the density maximum in the fluid density change curve in the oil layer or gas layer corresponds to the density value of the water-based drilling fluid filtrate, which changes dramatically with the pumping parameters and has low smoothness. As the pumping time or volume increases, the frequency of the slug flow becomes smaller and smaller, and the density maximum becomes thinner and disappears, and finally only a stable density value of oil or gas remains;若地层水类型为油水同层或气水同层,确定每个密度极大值对应的泵抽参数,将所述多个密度极大值及每个密度极大值对应的泵抽参数,作为样本数据;通过所述样本数据训练预设的神经网络模型,得到地层水类型为油水同层或气水同层的实际流体密度测量值和泵抽参数之间的关联关系;If the formation water type is oil-water in the same layer or gas-water in the same layer, determine the pumping parameters corresponding to each density maximum, and use the multiple density maximums and the pumping parameters corresponding to each density maximum as sample data; train a preset neural network model with the sample data to obtain the correlation between the actual fluid density measurement value and the pumping parameters when the formation water type is oil-water in the same layer or gas-water in the same layer;根据预设泵抽参数和所述关联关系,确定纯地层水密度;其中,所述预设泵抽参数是预计采集到纯地层水时的泵抽参数;Determine the density of pure formation water according to the preset pumping parameters and the correlation relationship; wherein the preset pumping parameters are pumping parameters expected when pure formation water is collected;根据所述纯地层水密度、纯水基钻井液滤液密度以及当前采集的实际流体测量密度,计算地层水污染率。The formation water contamination rate is calculated based on the pure formation water density, the pure water-based drilling fluid filtrate density and the actual fluid measurement density currently collected.2.根据权利要求1所述的方法,其特征在于,若所述流体密度变化曲线中不存在段塞流区段,则确定地层水类型为地层水层。2 . The method according to claim 1 , wherein if there is no slug flow section in the fluid density variation curve, the formation water type is determined to be a formation water layer.3.根据权利要求1所述的方法,其特征在于,所述段塞流区段为密度值规律震荡的区段。3. The method according to claim 1, characterized in that the slug flow section is a section where the density value oscillates regularly.4.根据权利要求1-3中任一项所述的方法,其特征在于,所述泵抽参数以及所述实际流体密度测量值由预设取样设备采集获得;4. The method according to any one of claims 1 to 3, characterized in that the pumping parameters and the actual fluid density measurement values are acquired by a preset sampling device;其中,所述预设取样设备至少包括:泵抽模块以及密度传感器;所述密度传感器位于泵抽模块的下游。Wherein, the preset sampling device at least includes: a pumping module and a density sensor; the density sensor is located downstream of the pumping module.5.一种计算设备,其特征在于,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;5. A computing device, comprising: a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other through the communication bus;所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如权利要求1-4中任一项所述的地层水污染率确定方法对应的操作。The memory is used to store at least one executable instruction, and the executable instruction enables the processor to perform operations corresponding to the formation water contamination rate determination method according to any one of claims 1-4.6.一种计算机存储介质,其特征在于,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如权利要求1-4中任一项所述的地层水污染率确定方法对应的操作。6. A computer storage medium, characterized in that at least one executable instruction is stored in the storage medium, and the executable instruction enables a processor to execute operations corresponding to the formation water contamination rate determination method as described in any one of claims 1-4.
CN202110912286.0A2021-08-102021-08-10 Formation water type identification method, pollution rate determination method and calculation deviceActiveCN113640178B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN202110912286.0ACN113640178B (en)2021-08-102021-08-10 Formation water type identification method, pollution rate determination method and calculation device

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN202110912286.0ACN113640178B (en)2021-08-102021-08-10 Formation water type identification method, pollution rate determination method and calculation device

Publications (2)

Publication NumberPublication Date
CN113640178A CN113640178A (en)2021-11-12
CN113640178Btrue CN113640178B (en)2025-03-25

Family

ID=78420427

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN202110912286.0AActiveCN113640178B (en)2021-08-102021-08-10 Formation water type identification method, pollution rate determination method and calculation device

Country Status (1)

CountryLink
CN (1)CN113640178B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN115855872A (en)*2022-12-192023-03-28中国海洋石油集团有限公司 Fluid type identification method, device and sampler based on downhole spectral measurement

Citations (2)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN111411949A (en)*2020-03-312020-07-14中国海洋石油集团有限公司Pumping device, sampling system and sampling method
CN112761631A (en)*2021-01-122021-05-07中国海洋石油集团有限公司Method for determining density, sampling method and pollution rate of pure formation water

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN101498215B (en)*2008-02-012014-12-10普拉德研究及开发股份有限公司Enhanced downhole fluid analysis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN111411949A (en)*2020-03-312020-07-14中国海洋石油集团有限公司Pumping device, sampling system and sampling method
CN112761631A (en)*2021-01-122021-05-07中国海洋石油集团有限公司Method for determining density, sampling method and pollution rate of pure formation water

Also Published As

Publication numberPublication date
CN113640178A (en)2021-11-12

Similar Documents

PublicationPublication DateTitle
CN110056348B (en) A method and system for determining formation fluid composition and properties
US20230053236A1 (en)Multivariate statistical contamination prediction using multiple sensors or data streams
CN110114552A (en)Method and system for determining drilling cuttings depth
EP3397833A1 (en)Machine learning for production prediction
EP1984860A2 (en)Methods, systems, and computer-readable media for fast updating of oil and gas field production models with physical and proxy simulators
CN111783847B (en)Low-contrast hydrocarbon reservoir identification method, device, equipment and system
CN116644284A (en) Method, system, electronic equipment and medium for determining characteristic factors of stratigraphic classification
CN113640178B (en) Formation water type identification method, pollution rate determination method and calculation device
CN112761631B (en)Density determination method, sampling method and pollution rate determination method for pure formation water
WO2023066549A1 (en)Virtual flow meter for well using a sucker rod pump
EP4208623A1 (en)Multi-sensor data assimilation and predictive analytics for optimizing well operations
CN113627640A (en)Productivity well testing prediction method and system for fracture-cavity type oil reservoir oil and gas well
US12400425B2 (en)Automated thermal maturation estimation from palynological sample images
CN119333123A (en) Downhole parameter prediction method and device
CN113392538B (en)Method and device for treating stratum water pollution rate under water-based mud condition
CN112377175B (en)Method and system for optimizing drilling mud and rapidly identifying low-resistance oil-gas layer
US20230306164A1 (en)Method for predicting sand production in a formation
Fan et al.Integration of image recognition and expert system for real-time wellbore stability analysis
CN119437809B (en)Method and device for predicting pollution degree of oil-based drilling fluid filtrate in real time
CN117496090B (en) A virtual borehole modeling method and system based on fine outcrop characterization
US20250116188A1 (en)Method and system for locating a downhole tool in steel-cased holes
CN119538469A (en) Formation interface picking method, device, electronic equipment and storage medium
WO2025080890A1 (en)Method and system for locating an untethered downhole tool in steel-cased wellbores
CN120608676A (en) Predrilling formation pressure prediction method, system, equipment and medium
CN118426055A (en)Porosity inversion method, device, equipment and medium based on support vector machine

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant

[8]ページ先頭

©2009-2025 Movatter.jp