CN114737957B

CN114737957B - Bottom hole temperature acquisition method, device, equipment and medium

Info

Publication number: CN114737957B
Application number: CN202110021670.1A
Authority: CN
Inventors: 彭庚; 郑有成; 李维; 张华礼; 周朗; 李玉飞; 唐庚; 张�林; 陆林峰; 马梓瀚
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2021-01-08
Filing date: 2021-01-08
Publication date: 2025-02-07
Anticipated expiration: 2041-01-08
Also published as: CN114737957A

Abstract

The embodiment of the present application discloses a method, device, computer equipment and medium for acquiring bottom hole temperature, which belongs to the technical field of oil and gas field well development. The method includes: acquiring natural gas data, the natural gas data including multiple attribute parameters of natural gas, the initial temperature corresponding to the natural gas well and the corresponding initial pressure; determining the target density of the natural gas according to the first relationship data; determining the first influence coefficient corresponding to the natural gas according to the second relationship data; updating the initial temperature according to the first influence coefficient and the reference pressure variable to obtain the updated temperature; in response to the difference between the target density and the initial density being less than the reference threshold, determining the bottom hole temperature of the natural gas well according to the updated temperature. While considering the influence of the natural gas's own properties on the temperature, the method also considers the change in temperature during the process of the natural gas flowing from the gas reservoir into the bottom hole, thereby improving the accuracy of the acquired bottom hole temperature.

Description

Method, device, equipment and medium for acquiring bottom hole temperature

Technical Field

The embodiment of the application relates to the technical field of oil and gas field well development, in particular to a method, a device, equipment and a medium for acquiring bottom hole temperature.

Background

As the depth of natural gas wells increases, the bottom hole temperature is also increasing, which has an important impact on wellbore pressure control, working fluid performance, downhole tools, and tubing strength. Therefore, the accurate prediction of the bottom hole temperature has important significance for safe and efficient production of ultra-deep gas wells.

In the related art, the bottom hole temperature of the natural gas well is obtained by adopting the vertical depth of the natural gas well and the ground temperature gradient data of the area where the natural gas well is located, but the ground temperature gradient is abnormal due to the fact that the natural gas well is deeper, and the influence of related parameters of the natural gas is not considered, so that the obtained bottom hole temperature is inaccurate.

Disclosure of Invention

The embodiment of the application provides a method, a device, computer equipment and a medium for acquiring bottom hole temperature, which improve the accuracy of acquired bottom hole temperature. The technical scheme is as follows:

in one aspect, a method of downhole temperature acquisition is provided, the method comprising:

acquiring natural gas data, wherein the natural gas data comprises a plurality of attribute parameters of natural gas, corresponding initial temperature and corresponding initial pressure of a natural gas well;

determining a target density of the natural gas according to first relation data, wherein the first relation data is used for representing a corresponding relation between the density of the natural gas and the natural gas data;

Determining a first influence coefficient corresponding to the natural gas according to second relation data, wherein the second relation data is used for representing the corresponding relation between the influence coefficient corresponding to the natural gas and the natural gas data, and the first influence coefficient is used for representing the influence degree of the pressure corresponding to the natural gas well on the corresponding temperature;

updating the initial temperature according to the first influence coefficient and the reference pressure variable to obtain an updated temperature;

And determining a bottom hole temperature of the natural gas well according to the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold.

In one possible implementation, the determining the bottom hole temperature of the natural gas well from the updated temperature parameter in response to the difference between the target density and the initial density being less than a reference threshold comprises:

Determining the updated temperature parameter as the bottom hole temperature in response to the difference between the target density and the initial density being less than the reference threshold, or

Updating the initial pressure parameter according to the reference pressure variable to obtain an updated pressure parameter in response to the difference between the target density and the initial density being less than the reference threshold;

Determining a second influence coefficient corresponding to the natural gas according to the second relation data, the updated temperature parameter and the updated pressure parameter;

and updating the updated temperature parameter according to the second influence coefficient and the reference pressure variable to obtain a re-updated temperature parameter, and determining the re-updated temperature parameter as the bottom hole temperature.

In another possible implementation manner, the method further includes, after updating the initial temperature according to the target influence coefficient and the reference pressure variable and obtaining an updated temperature:

updating the initial pressure according to the reference pressure variable to obtain updated pressure in response to the difference between the target density and the initial density being not less than the reference threshold;

And taking the updated pressure as initial pressure, the updated temperature as initial temperature, the target density as initial density, and continuing to execute the step of determining the target density of the natural gas according to the first relation data.

In another possible implementation, the first data relationship includes:

Wherein p is the initial pressure, ρ is the density of the natural gas, R is a universal gas constant, T is the initial temperature, α, γ, a, b, c, and d are preset parameters, and a ₀、B₀、C₀、D₀ and E ₀ are known parameters determined based on the plurality of attribute parameters.

In another possible implementation, the second data relationship includes:

wherein mu _J-T is the influence coefficient, c _p is the constant pressure molar specific heat capacity of the natural gas, T is the initial temperature, ρ is the density of the natural gas, and R is a universal gas constant.

In another possible implementation, the natural gas includes a plurality of components, and the acquiring the natural gas data includes:

for each attribute parameter, respectively acquiring attribute subparameters and corresponding mole fractions of the plurality of components;

And fusing the acquired attribute sub-parameters according to the mole fractions corresponding to the plurality of components to obtain the attribute parameters.

In another possible implementation, the acquiring natural gas data includes:

displaying a parameter setting interface;

and acquiring the natural gas data input in the parameter setting interface.

In another possible implementation, the parameter setting interface includes a static parameter setting area for setting the initial temperature, a fluid parameter setting area for setting the plurality of attribute parameters, and a production parameter setting area for setting the initial pressure.

In another aspect, there is provided a downhole temperature acquisition device, the device comprising:

the system comprises a data acquisition module, a data processing module and a data processing module, wherein the data acquisition module is used for acquiring natural gas data, and the natural gas data comprises a plurality of attribute parameters of natural gas, initial temperatures corresponding to a natural gas well and initial pressures corresponding to the natural gas well;

The density acquisition module is used for determining the target density of the natural gas according to first relation data, wherein the first relation data is used for representing the corresponding relation between the density of the natural gas and the natural gas data;

the influence coefficient acquisition module is used for determining a first influence coefficient corresponding to the natural gas according to second relation data, the second relation data are used for representing the corresponding relation between the influence coefficient corresponding to the natural gas and the natural gas data, and the first influence coefficient is used for representing the influence degree of the pressure corresponding to the natural gas well on the corresponding temperature;

the temperature updating module is used for updating the initial temperature according to the first influence coefficient and the reference pressure variable to obtain an updated temperature;

And the bottom hole temperature acquisition module is used for responding to the difference value between the target density and the initial density to be smaller than a reference threshold value and determining the bottom hole temperature of the natural gas well according to the updated temperature.

In one possible implementation, the bottom hole temperature acquisition module is configured to:

In another possible implementation, the apparatus further includes:

a pressure updating module, configured to update the initial pressure according to the reference pressure variable in response to the difference between the target density and the initial density being not less than the reference threshold, and obtain an updated pressure;

and the data acquisition module is further used for taking the updated pressure as initial pressure, taking the updated temperature as initial temperature, taking the target density as initial density, and continuing to execute the step of determining the target density of the natural gas according to the first relation data.

In another possible implementation, the first data relationship includes:

Wherein p is the initial pressure, ρ is the density of the natural gas, R is a universal gas constant, T is the initial temperature, a, γ, a, b, c, and d are preset parameters, and a ₀、B₀、C₀、D₀ and E ₀ are known parameters determined based on the plurality of attribute parameters.

In another possible implementation, the second data relationship includes:

In another possible implementation, the natural gas includes a plurality of components, and the data acquisition module is configured to:

In another possible implementation manner, the data acquisition module is configured to:

displaying a parameter setting interface;

and acquiring the natural gas data input in the parameter setting interface.

In another aspect, a computer device is provided that includes a processor and a memory having at least one computer program stored therein, the at least one computer program loaded and executed by the processor to perform the operations performed in the bottom hole temperature acquisition method as described in the above aspect.

In another aspect, a computer readable storage medium having at least one computer program stored therein is provided, the at least one computer program loaded and executed by a processor to perform the operations performed in the downhole temperature acquisition method as described in the above aspect.

In another aspect, a computer program product or a computer program is provided, the computer program product or the computer program comprising computer program code stored in a computer readable storage medium, the computer program code being read from the computer readable storage medium by a processor of a computer device, the computer program code being executed by the processor such that the computer device implements the operations performed in the downhole temperature acquisition method as described in the above aspect.

The technical scheme provided by the embodiment of the application has the beneficial effects that at least:

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the description of the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

FIG. 1 is a flow chart of a method for acquiring bottom hole temperature according to an embodiment of the present application;

FIG. 2 is a flow chart of another method of acquiring bottom hole temperature provided by an embodiment of the present application;

FIG. 3 is a schematic illustration of a bottom hole construction of a natural gas well provided in accordance with an embodiment of the present application;

FIG. 4 is a flow chart of another method of acquiring bottom hole temperature provided by an embodiment of the present application;

FIG. 5 is a schematic diagram of simulation software provided by an embodiment of the present application;

FIG. 6 is a schematic illustration of a change in bottom hole temperature provided by an embodiment of the present application;

FIG. 7 is a schematic diagram of a gas reservoir temperature distribution provided by an embodiment of the present application;

FIG. 8 is a schematic diagram of a bottom hole temperature acquisition device according to an embodiment of the present application;

FIG. 9 is a schematic diagram of another embodiment of a bottom hole temperature acquisition device;

Fig. 10 is a schematic structural diagram of a terminal according to an embodiment of the present application;

Fig. 11 is a schematic structural diagram of a server according to an embodiment of the present application.

Detailed Description

For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the following detailed description of the embodiments of the present application will be given with reference to the accompanying drawings.

It is to be understood that the terms "first," "second," and the like, as used herein, may be used to describe various concepts, but are not limited by these terms unless otherwise specified. These terms are only used to distinguish one concept from another. For example, the first relationship data may be referred to as second relationship data and the second relationship data may be referred to as first relationship data without departing from the scope of the application.

The terms "at least one", "a plurality", "each", "any" and the like as used herein, at least one includes one, two or more, a plurality includes two or more, each means each of the corresponding plurality, and any one means any of the plurality. For example, the plurality of attribute parameters includes 3 attribute parameters, and each attribute parameter refers to each of the 3 attribute parameters, and any one refers to any one of the 3 attribute parameters, which may be the first, the second, or the third.

For a better understanding of the method provided by the present application, the following explains the keywords related to the present application:

Joule-Thomson effect (Joule-Thomson effect), which refers to the phenomenon of a change in temperature of a gas after it has flowed from a high pressure zone to a low pressure zone, is characterized by the Joule-Thomson coefficient. A gas Joule-Thomson coefficient greater than 0 indicates a decrease in temperature of the gas after the pressure decrease, a gas Joule-Thomson coefficient less than 0 indicates an increase in temperature of the gas after the pressure decrease, and a gas Joule-Thomson coefficient equal to 0 indicates a constant temperature of the gas after the pressure decrease. The curve formed by the points where the Joule-Thomson coefficient is equal to 0 is a Joule-Thomson inversion curve, and both sides of the curve represent conditions for increasing and decreasing the gas temperature, respectively.

The drilling rate determining method provided by the embodiment of the application can be applied to computer equipment, wherein the computer equipment is a terminal or a server. The terminal is a portable, pocket-sized, hand-held terminal of various types, such as a mobile phone, a computer, a tablet personal computer and the like. The server is a server, or a server cluster formed by a plurality of servers, or a cloud computing service center.

FIG. 1 is a flow chart of a method for acquiring bottom hole temperature according to an embodiment of the present application. The execution subject of the embodiment of the application is a computer device, referring to fig. 1, the method includes the following steps:

101. Natural gas data is acquired, the natural gas data including a plurality of property parameters of the natural gas, a corresponding initial temperature of the natural gas well, and a corresponding initial pressure.

102. And determining the target density of the natural gas according to first relation data, wherein the first relation data is used for representing the corresponding relation between the density of the natural gas and the natural gas data.

103. And determining a first influence coefficient corresponding to the natural gas according to second relation data, wherein the second relation data is used for representing the corresponding relation between the influence coefficient corresponding to the natural gas and the natural gas data, and the first influence coefficient is used for representing the influence degree of the pressure corresponding to the natural gas well on the corresponding temperature.

104. And updating the initial temperature according to the first influence coefficient and the reference pressure variable to obtain the updated temperature.

105. And determining a bottom hole temperature of the natural gas well from the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold.

According to the method provided by the embodiment of the application, when the bottom hole temperature is acquired, not only the attribute parameters of the natural gas are adopted, but also the initial temperature and the initial pressure are combined, namely, the influence of the attribute of the natural gas on the temperature is considered, meanwhile, the temperature change in the process of flowing the natural gas from the gas reservoir into the bottom hole is considered, the temperature is updated according to the influence coefficient, and the accuracy of the acquired bottom hole temperature is improved.

In one possible implementation, determining a bottom hole temperature of the natural gas well from the updated temperature parameter in response to the difference between the target density and the initial density being less than a reference threshold comprises:

determining an updated temperature parameter as a bottom hole temperature in response to the difference between the target density and the initial density being less than a reference threshold, or

In response to the difference between the target density and the initial density being less than the reference threshold, updating the initial pressure parameter according to the reference pressure variable to obtain an updated pressure parameter;

In another possible implementation manner, the method further includes, after updating the initial temperature according to the target influence coefficient and the reference pressure variable and obtaining the updated temperature:

in response to the difference between the target density and the initial density being not less than a reference threshold, updating the initial pressure according to a reference pressure variable to obtain an updated pressure;

And taking the updated pressure as initial pressure, taking the updated temperature as initial temperature, taking the target density as initial density, and continuously executing the step of determining the target density of the natural gas according to the first relation data.

In another possible implementation, the first data relationship includes:

Where p is the initial pressure, ρ is the density of natural gas, R is the universal gas constant, T is the initial temperature, α, γ, a, b, c, and d are preset parameters, and a ₀、B₀、C₀、D₀ and E ₀ are known parameters determined based on a plurality of attribute parameters.

In another possible implementation, the second data relationship includes:

wherein mu _J-T is an influence coefficient, c _p is a constant pressure molar specific heat capacity of natural gas, T is an initial temperature, ρ is the density of the natural gas, and R is a universal gas constant.

In another possible implementation, the natural gas includes a plurality of constituents, and acquiring the natural gas data includes:

for each attribute parameter, respectively acquiring attribute subparameters and corresponding mole fractions of various components;

And fusing the acquired attribute sub-parameters according to the mole fractions corresponding to the components to obtain the attribute parameters.

In another possible implementation, acquiring natural gas data includes:

displaying a parameter setting interface;

and acquiring natural gas data input in a parameter setting interface.

In another possible implementation, the parameter setting interface includes a static parameter setting area for setting an initial temperature, a fluid parameter setting area for setting a plurality of attribute parameters, and a production parameter setting area for setting an initial pressure.

FIG. 2 is a flow chart of another method for acquiring bottom hole temperature provided by an embodiment of the present application. The execution subject of the embodiment of the application is a computer device, referring to fig. 2, the method includes the following steps:

201. Natural gas data is acquired.

In the embodiment of the application, in order to meet the actual situation of the natural gas well, the natural gas distribution in the gas reservoir is considered to be uniform, the natural gas is uniformly and stably flowed from the gas reservoir to the bottom of the well, and the expansion of the natural gas is adiabatic and irreversible, so that the bottom-hole temperature of the natural gas well is predicted on the basis.

The natural gas data comprises a plurality of attribute parameters of the natural gas, wherein the attribute parameters of the natural gas comprise mole fraction, equivalent molar mass, critical temperature, critical density, critical pressure of each component of the natural gas, binary interaction coefficient of the natural gas and the like.

In one possible implementation, the natural gas includes a plurality of components, and the attribute parameters of the natural gas are determined according to attribute sub-parameters of the plurality of components. For each attribute parameter of the natural gas, respectively obtaining attribute subparameters and corresponding mole fractions of various components of the natural gas; and fusing the acquired attribute sub-parameters according to the mole fractions corresponding to the plurality of components to obtain the attribute parameters.

For example, natural gas contains n components, and the equivalent molar mass of natural gas is determined using the following formula:

Wherein M represents an equivalent molar mass of natural gas, n represents the number of constituent components contained in natural gas, x _i represents a molar fraction of the ith constituent component, and M _i represents an equivalent molar mass of the ith constituent component.

The critical temperature of natural gas is determined using the following formula:

Wherein T _c represents the critical temperature of natural gas, n represents the number of constituent components contained in natural gas, x _i represents the mole fraction of the ith constituent component, and T _ci represents the critical temperature of the ith constituent component.

The critical pressure of natural gas is determined using the following formula:

Where P _c represents the critical pressure of natural gas, n represents the number of constituent components contained in natural gas, x _i represents the mole fraction of the ith constituent component, and P _ci represents the critical pressure of the ith constituent component.

The natural gas data also comprises an initial temperature corresponding to the natural gas well and an initial pressure corresponding to the natural gas well, wherein the initial temperature refers to a gas reservoir temperature (formation temperature), and the initial pressure refers to a gas reservoir pressure (formation pressure).

For example, referring to the schematic bottom hole structure of FIG. 3, there is shown the relationship between bottom hole temperature and bottom hole pressure and reservoir temperature and reservoir pressure, where the temperature and pressure may change during the flow of natural gas from the reservoir into the well bore of the natural gas well. Wherein P _res refers to the reservoir pressure, T _res refers to the reservoir temperature, P _wf refers to the bottom hole temperature, and T _wf refers to the bottom hole temperature.

202. And determining the target density of the natural gas according to the first relation data.

The first relation data are used for representing the corresponding relation between the density of the natural gas and the natural gas data. Optionally, the first relationship data is a natural gas state equation, where the natural gas state equation is shown in the following formula (1):

In one possible implementation, for the above formula (1), under the condition of keeping the initial pressure unchanged, deriving T, and obtaining the following deformed first relationship data:

Wherein C ₁、C₂、C₃ and C ₄ are the following formulas (3), (4), (5) and (6), respectively:

In one possible implementation manner, since the expression about density in the above formula (1) or formula (2) is a hidden function, an iteration method is required to solve the density of the natural gas, and the above formula (1) or formula (2) is deformed to obtain the following formula (7):

Assuming that the iteration number is k times, iteration initial values ρ ⁽⁰⁾ and ρ ⁽¹⁾ are preset, and the following formula (8) is adopted as an iteration equation to perform iterative calculation:

in the above formula (8), ρ ^(k+1) can be obtained by solving for ρ ^(k-1) and ρ ^(k).

203. And determining a first influence coefficient corresponding to the natural gas according to the second relation data.

The second relation data are used for representing the corresponding relation between the influence coefficient corresponding to the natural gas and the natural gas data, and the first influence coefficient is used for representing the influence degree of the pressure corresponding to the natural gas well on the corresponding temperature. Alternatively, the influence coefficient is a Joule-Thomson coefficient.

In one possible implementation, the second relationship data is the following formula (9):

Wherein mu _J-T is an influence coefficient (Joule-Thomson coefficient), c _p is a constant pressure molar specific heat capacity of natural gas, T is an initial temperature, ρ is a density of the natural gas, and R is a general gas constant.

Optionally, the second relationship data is deformed according to the first relationship data shown in the above formula (2), so as to obtain second relationship data shown in the following formula (10):

optionally, the natural gas has a constant pressure molar specific heat capacity of:

c_p＝1243+3.14T+7.931×10^-4T²-6.881×10^-7T³;

wherein T is the initial temperature.

204. And updating the initial temperature according to the first influence coefficient and the reference pressure variable to obtain the updated temperature.

In the embodiment of the application, the influence of the Joule-Thomson effect on the bottom hole temperature is considered, so that in the process of determining the temperature, the initial temperature is updated according to the first influence coefficient and the reference pressure variable to obtain the updated temperature, and whether the updated temperature can be used as the bottom hole temperature or not is determined according to whether the acquired density meets the condition or not. If the acquired target density satisfies the condition, step 205 is performed, and if the target density does not satisfy the condition, step 206 is performed.

In one possible implementation, the initial temperature is updated using the following equation (11):

T′=T-ΔP·μ_J-T (11)

wherein T' is the updated temperature, T is the initial temperature, deltaP is the reference pressure variable, and mu _J-T is the first influence coefficient. Where Δp is any number, for example, Δp is 1, 2, or other number.

205. And determining a bottom hole temperature of the natural gas well from the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold.

In the embodiment of the application, when the difference between the target density and the initial density is smaller than the reference threshold value, the target density is considered to meet the condition. The target density is the density obtained in the current iteration process, and the initial density is the density obtained in the last iteration process of the current iteration process, or if the target density is the density obtained by the first calculation, the initial density is the preset density.

In one possible implementation, the difference between the target density and the initial density needs to satisfy the following condition:

|ρ^(k+1)-ρ^(k)|＜10^-5

Where k+1 represents the current iteration process, k represents the last iteration process, and the reference threshold is 10 ^-5. The reference threshold may also be other values, a smaller reference threshold may be set if a higher accuracy is required, or a larger reference threshold may be set if an accuracy is not required to be too high.

For example, the initial density of iterations is preset to ρ ⁽⁰⁾ =0,

In one possible implementation, the updated temperature is determined to be the bottom hole temperature of the natural gas well if the obtained target density meets the condition.

In another possible implementation, in response to the difference between the target density and the initial density being less than the reference threshold, updating the initial pressure parameter according to the reference pressure variable to obtain an updated pressure parameter, determining a second influence coefficient corresponding to the natural gas according to the second relationship data, the updated temperature parameter and the updated pressure parameter, updating the updated temperature parameter according to the second influence coefficient and the reference pressure variable to obtain a re-updated temperature parameter, and determining the re-updated temperature parameter as the bottom hole temperature. Wherein the second influence coefficient is similar to the first influence coefficient.

That is, after the updated temperature and the updated pressure are obtained, the updated temperature can be calculated again according to the updated temperature and the updated pressure, so that the corresponding temperatures under different pressures are obtained, and the corresponding temperatures under different pressures are obtained, namely, the temperature distribution of the natural gas well is obtained.

Optionally, the initial pressure is updated using the following equation (12):

P′=P-ΔP (12)

Wherein P' is the updated pressure, P is the initial pressure, and ΔP is the reference pressure variable.

206. And in response to the difference between the target density and the initial density not being less than the reference threshold, updating the initial pressure according to the reference pressure variable to obtain an updated pressure.

207. The process continues to step 202 with the updated pressure as the initial pressure, the updated temperature as the initial temperature, and the target density as the initial density.

In the embodiment of the application, when the difference between the target density and the initial density is not smaller than the reference threshold, the target density is considered to be unsatisfied with the condition, and iteration is further required to be continued until the obtained target density meets the condition, and the obtained temperature can be used as the bottom hole temperature. At this time, the updated temperature and the updated pressure are set as the initial temperature and the initial pressure, respectively, and the obtained target density is set as the initial density, and step 202 is executed again until the obtained target density satisfies the condition.

It should be noted that, in the embodiment of the present application, only one update process is taken as an example, and in another embodiment, after the updated temperature and the updated pressure are obtained, the process of step 202 and step 207 can be further adopted, and the temperature and the pressure are continuously updated, so as to obtain the temperatures under different pressures, i.e. obtain the bottom hole temperature distribution.

It should be noted that, in another embodiment, the above formula is merely described by way of example, and in another embodiment, the above formula representing the first relationship data and the second relationship data may also be in other forms, which is not limited in this embodiment of the present application.

And moreover, the bottom hole temperature is accurately controlled, scientific basis is provided for selection of ultra-deep well downhole tools, working fluid and pipe columns, design level of ultra-deep well exploitation equipment is improved, and therefore development of ultra-deep gas reservoirs is facilitated. Wherein, ultra-deep well refers to natural gas well with depth greater than 7000 meters.

Additionally, in one possible implementation, referring to the flow chart of the bottom hole temperature acquisition method shown in fig. 4, the bottom hole temperature acquisition process is referred to as follows:

401. an initial temperature and an initial pressure are obtained.

402. A plurality of attribute parameters of the natural gas are determined based on the attribute sub-parameters of the plurality of constituent components of the natural gas.

403. An initial density is obtained.

404. And calculating the target density of the natural gas according to the first relation data.

405. And according to the second relation data, obtaining an influence coefficient corresponding to the natural gas.

406. And updating the initial temperature according to the influence coefficient and the reference pressure variable to obtain the updated temperature.

407. And updating the initial pressure according to the influence coefficient to obtain the updated pressure.

408. It is determined whether the difference between the target density and the initial density is less than a reference threshold, if the difference is less than the reference threshold, step 406 is performed, and if the difference is not less than the reference threshold, step 404 is performed again.

409. The updated temperature is determined to be the bottom hole temperature.

In addition, in one possible implementation, simulation software is installed in the computer device, where the simulation software includes a data storage module, an input-output module, an algorithm module, and a graphics module.

The system comprises a parameter setting interface, a data storage module, an input/output module and an algorithm module, wherein the parameter setting interface is used for storing natural gas data, namely a user inputs the natural gas data in the parameter setting interface, the computer equipment acquires the natural gas data input in the parameter setting interface, the input/output module is used for reading and outputting the natural gas data from a memory, the algorithm module is used for processing the natural gas data by adopting the bottom hole temperature acquisition process shown in the steps 202-207, the obtained bottom hole temperature graphic module is used for processing the obtained data, and the obtained data is displayed to the user in the form of icons.

In one possible implementation, the natural gas data is acquired through a parameter setting interface that includes a static parameter setting area, a fluid parameter setting area, and a production parameter setting area. Wherein the static parameter setting area is used for setting an initial temperature, the fluid parameter setting area is used for setting a plurality of attribute parameters, and the production parameter setting area is used for setting an initial pressure.

For example, referring to the schematic diagram of simulation software shown in fig. 5, it includes a data storage module, an input-output module, a graphics module, and an algorithm module. The system comprises a data storage module, an input/output module, a graphic module, an algorithm module and a control module, wherein the data storage module stores static data, fluid data and production data, the input/output module comprises a natural gas database and a gas reservoir database, the graphic module comprises a temperature prediction graphic and a gas reservoir distribution graphic, and the algorithm module comprises an established model and a set algorithm. The data storage module and the input/output module interact with each other, the input/output module interacts with the algorithm module, and the input/output module interacts with the graphic module, i.e. the input/output module can send the obtained output to the graphic module, and the graphic module generates a corresponding graphic.

For example, the bottom hole temperature is obtained using the natural gas data shown in tables 1 and 2 below, the property parameters of the natural gas are shown in table 1 below, and the gas reservoir data of the natural gas are shown in table 2 below:

TABLE 1

Name of the name	Numerical value	Unit (B)	Type(s)
				CH ₄ mole fraction	0.91	/	F
C ₂H₆ mole fraction	0.09	/	F
				CH ₄ molar mass	16.042	g/mol	F
Molar mass of C ₂H₆	30.069	g/mol	F
				Critical temperature of CH ₄	-82.595	°C	F
C ₂H₆ critical temperature	32.172	°C	F
				Critical density of CH ₄	162.658	kg/mol	F
Critical density of C ₂H₆	206.179	kg/mol	F
				Critical pressure of CH ₄	4.5992	MPa	F
Critical pressure of C ₂H₆	4.8722	MPa	F
				Binary interaction coefficient	0.01	/	F
Bottom hole pressure	48	MPa	P

TABLE 2

Name of the name	Numerical value	Unit (B)	Type(s)
				Gas reservoir temperature	155	°C	S
Pressure of gas reservoir	98	MPa	S

The natural gas data in tables 1 and 2 above were input into simulation software to obtain the bottom hole temperature. For example, referring to the schematic diagram of the bottom hole temperature shift shown in fig. 6, the corresponding Joule-Thomson coefficient on the left side of the Joule-Thomson effect inversion curve in fig. 5 is less than 0, and the corresponding Joule-Thomson coefficient on the left side of the Joule-Thomson effect inversion curve is greater than 0, i.e. the natural gas releases heat and absorbs heat as the bottom hole pressure decreases, and as can be seen in fig. 6, for three different natural gas wells, the corresponding temperatures of the bottom hole pressure increase and decrease as the bottom hole pressure decreases. The three different natural gas wells have different input parameters, namely, the specific conditions of the natural gas wells are different, so that the change degree of the bottom hole temperature is also different.

See, for example, the schematic diagram of the temperature profile shown in fig. 7. As can be seen from fig. 6, the closer the temperature and pressure of the natural gas in the same natural gas well is to the Joule-Thomson effect inversion curve, the higher the corresponding bottom hole temperature.

In addition, in the related art, a geothermal gradient method is adopted to obtain the bottom hole temperature, and the geothermal gradient method is used for determining the bottom hole temperature according to the vertical depth of a gas well and the geothermal gradient data. In an actual application scene, the bottom hole temperature is calculated by selecting an average value of the ground temperature gradients in a certain area, the natural gas distribution area is mostly in a ground temperature gradient abnormal zone, and the method does not consider the Joule-Thomson effect in the process that natural gas flows into the bottom hole from a reservoir, so that the ground temperature gradient method is only suitable for calculating the bottom hole temperature of a shallow well, and the difference between the bottom hole temperature calculated by the method and the actual bottom hole temperature is larger and the accuracy is lower for an ultra-deep well. In the embodiment of the application, the influence of the Joule-Thomson effect on the bottom hole temperature is considered, so that the obtained bottom hole temperature is more accurate.

In another related art, a ground inverse algorithm is used to obtain the bottom hole temperature, and the ground inverse algorithm is based on the data of the wellhead temperature, pressure, flow rate and the like of the gas producing well, and a thermodynamic equation and a material balance equation are used to inversely calculate the bottom hole temperature. In the practical application process, the solving process of the method is complex, a large number of known parameters are required to be input, when the yield of the gas well changes, the prediction calculation of the bottom hole temperature is more complex, the method provided by the application embodiment requires the attribute parameters of the natural gas, the gas reservoir temperature and the gas reservoir pressure, the parameters required to be input are fewer, the calculation process is simpler, and the influence of the Joule-Thomson effect is considered, so that the obtained bottom hole temperature is more accurate.

In another related technology, a test temperature-pressure distribution coupling prediction model is adopted to obtain the bottom hole temperature, and the prediction model comprehensively considers factors influencing temperature and pressure according to the basic principles of momentum, energy conservation law and mass and heat transfer theory. In this predictive model, the heat transfer in the wellbore is considered to be steady-state heat transfer, while the heat transfer in the formation is non-steady-state heat transfer, taking into account the effect of wellbore heat loss on temperature. The four-order Longgar-Kutta method is adopted to carry out coupling iterative solution on the prediction model, the well bottom is used as a reference to respectively carry out sensitivity analysis on the well bore temperature and pressure under the conditions of different yields, different relative gas densities and different production time, and the prediction results of the coupling model and the linear model are compared and analyzed. The prediction result of the prediction model shows that the well bore pressure is reduced and the temperature is increased along with the increase of the yield, the well bore pressure is reduced and the temperature is increased along with the increase of the gas density, the well bore pressure is kept unchanged and the temperature is increased along with the increase of the production time, and the temperature in the well bore changes along with the well depth in a nonlinear rule. However, the prediction model has excessive input parameters and complex calculation procedures, and the influence of the reverse throttling effect at the bottom of the well is not considered, but the method provided by the embodiment of the application only needs the attribute parameters of the natural gas, the gas reservoir temperature and the gas reservoir pressure, the parameters needing to be input are fewer, the calculation process is simpler, and the influence of the Joule-Thomson effect is considered, so that the obtained bottom of the well temperature is more accurate.

In another related technology, aiming at the problems that the conventional fluid physical property parameter calculation method is difficult to accurately apply in a high-temperature high-pressure gas reservoir, the conventional steady-state temperature pressure model is difficult to accurately predict the pressure and temperature distribution of a shaft, and the multiphase pipe flow pressure calculation model under different water-gas ratios is difficult to select, a calculation model of deviation factors and viscosity of a high-temperature high-pressure gas reservoir is adopted, and an unsteady-state temperature model in the well opening and closing process is established by using a shaft-stratum unsteady-state heat transfer model and a fluid thermodynamic model. The prediction result of the model shows that the unsteady state temperature model can accurately predict the temperature distribution higher by Wen Qijing than the steady state temperature model, and the unsteady state heat transfer model considering the well closing shaft temperature can correct the pressure recovery abnormal curve in the well closing process, so that the well test interpretation result can be closer to the measured value. However, the method needs to carry out additional correction, a large amount of field measured data needs to be acquired during correction, and the method provided by the embodiment of the application is difficult to widely apply, does not need a large amount of measured data, needs less workload and has a wider application range.

Fig. 8 is a schematic structural diagram of a bottom hole temperature acquiring device according to an embodiment of the present application. Referring to fig. 8, the apparatus includes:

the data acquisition module 801 is configured to acquire natural gas data, where the natural gas data includes a plurality of attribute parameters of natural gas, an initial temperature corresponding to a natural gas well, and a corresponding initial pressure;

A density acquisition module 802, configured to determine a target density of the natural gas according to first relationship data, where the first relationship data is used to represent a correspondence between the density of the natural gas and the natural gas data;

the influence coefficient obtaining module 803 is configured to determine a first influence coefficient corresponding to natural gas according to second relationship data, where the second relationship data is used to represent a correspondence between the influence coefficient corresponding to natural gas and natural gas data, and the first influence coefficient is used to represent a degree of influence of pressure corresponding to a natural gas well on a corresponding temperature;

A temperature updating module 804, configured to update an initial temperature according to the first influence coefficient and the reference pressure variable, to obtain an updated temperature;

a bottom hole temperature acquisition module 805 for determining a bottom hole temperature of the natural gas well from the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold.

According to the device provided by the embodiment of the application, when the bottom hole temperature is acquired, not only the attribute parameters of the natural gas are adopted, but also the initial temperature and the initial pressure are combined, namely, the influence of the attribute of the natural gas on the temperature is considered, meanwhile, the temperature change in the process of flowing the natural gas from the gas reservoir into the bottom hole is considered, the temperature is updated according to the influence coefficient, and the accuracy of the acquired bottom hole temperature is improved.

In one possible implementation, the bottom hole temperature acquisition module 805 is configured to:

In another possible implementation, referring to fig. 9, the apparatus further includes:

a pressure updating module 806, configured to update the initial pressure according to the reference pressure variable to obtain an updated pressure in response to the difference between the target density and the initial density being not less than the reference threshold;

The data obtaining module 801 is further configured to continuously perform the step of determining the target density of the natural gas according to the first relationship data, with the updated pressure as the initial pressure, the updated temperature as the initial temperature, and the target density as the initial density.

In another possible implementation, the first data relationship includes:

In another possible implementation, the second data relationship includes:

In another possible implementation, the natural gas includes a plurality of components, and the data acquisition module 801 is configured to:

In another possible implementation, the data acquisition module 801 is configured to:

displaying a parameter setting interface;

and acquiring natural gas data input in a parameter setting interface.

Any combination of the above optional solutions may be adopted to form an optional embodiment of the present application, which is not described herein.

It should be noted that, when the bottom hole temperature acquiring device provided in the above embodiment acquires the bottom hole temperature, only the division of the above functional modules is used for illustration, in practical application, the above functional allocation may be completed by different functional modules according to needs, that is, the internal structure of the computer device is divided into different functional modules to complete all or part of the functions described above. In addition, the bottom hole temperature acquiring device and the bottom hole temperature acquiring method provided in the foregoing embodiments belong to the same concept, and specific implementation processes thereof are detailed in the method embodiments and are not described herein again.

The embodiment of the application also provides a computer device, which comprises a processor and a memory, wherein at least one computer program is stored in the memory, and the at least one computer program is loaded and executed by the processor to realize the operation executed in the bottom hole temperature acquisition method of the embodiment.

Optionally, the computer device is provided as a terminal. Fig. 10 is a schematic structural diagram of a terminal 1000 according to an embodiment of the present application. The terminal 1000 can be a portable mobile terminal such as a smart phone, tablet, MP3 player (Moving Picture Experts Group Audio Layer III, MPEG audio layer 3), MP4 (Moving Picture Experts Group Audio Layer IV, MPEG audio layer 4) player, notebook, or desktop. Terminal 1000 can also be referred to by other names of user equipment, portable terminal, laptop terminal, desktop terminal, etc.

Terminal 1000 can include a processor 1001 and a memory 1002.

The processor 1001 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so on. The processor 1001 may be implemented in at least one hardware form of DSP (DIGITAL SIGNAL Processing), FPGA (Field-Programmable gate array), PLA (Programmable Logic Array ). The processor 1001 may also include a main processor for processing data in the awake state, which is also called a CPU (Central Processing Unit ), and a coprocessor for processing data in the standby state, which is a low-power-consumption processor. In some embodiments, the processor 1001 may be integrated with a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 1001 may also include an AI (ARTIFICIAL INTELLIGENCE ) processor for processing computing operations related to machine learning.

Memory 1002 may include one or more computer-readable storage media, which may be non-transitory. Memory 1002 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 1002 is used to store at least one computer program for execution by processor 1001 to implement the downhole temperature acquisition method provided by the method embodiments of the present application.

Terminal 1000 can also optionally include a peripheral interface 1003 and at least one peripheral in some embodiments. The processor 1001, the memory 1002, and the peripheral interface 1003 may be connected by a bus or signal line. The various peripheral devices may be connected to the peripheral device interface 1003 via a bus, signal wire, or circuit board. Specifically, the peripheral devices include at least one of radio frequency circuitry 1004, a display 1005, a camera assembly 1006, audio circuitry 1007, a positioning assembly 1008, and a power supply 1009.

Peripheral interface 1003 may be used to connect I/O (Input/Output) related at least one peripheral to processor 1001 and memory 1002. In some embodiments, processor 1001, memory 1002, and peripheral interface 1003 are integrated on the same chip or circuit board, and in some other embodiments, either or both of processor 1001, memory 1002, and peripheral interface 1003 may be implemented on separate chips or circuit boards, as this embodiment is not limiting.

Radio Frequency circuit 1004 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. Radio frequency circuitry 1004 communicates with a communication network and other communication devices via electromagnetic signals. The radio frequency circuit 1004 converts an electrical signal into an electromagnetic signal for transmission, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuitry 1004 includes an antenna system, an RF transceiver, one or more amplifiers, tuners, oscillators, digital signal processors, codec chipsets, subscriber identity module cards, and so forth. Radio frequency circuitry 1004 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to, the world wide web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and/or WiFi (WIRELESS FIDELITY ) networks. In some embodiments, the radio frequency circuit 1004 may further include NFC (NEAR FIELD Communication) related circuits, which is not limited by the present application.

The display screen 1005 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display 1005 is a touch screen, the display 1005 also has the ability to capture touch signals at or above the surface of the display 1005. The touch signal may be input to the processor 1001 as a control signal for processing. At this time, the display 1005 may also be used to provide virtual buttons and/or virtual keyboards, also referred to as soft buttons and/or soft keyboards. In some embodiments, display 1005 may be one, disposed on the front panel of terminal 1000, in other embodiments, display 1005 may be at least two, disposed on different surfaces of terminal 1000 or in a folded configuration, respectively, and in other embodiments, display 1005 may be a flexible display disposed on a curved surface or a folded surface of terminal 1000. Even more, the display 1005 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The display 1005 may be made of LCD (Liquid CRYSTAL DISPLAY), OLED (Organic Light-Emitting Diode) or other materials.

The camera assembly 1006 is used to capture images or video. Optionally, camera assembly 1006 includes a front camera and a rear camera. The front camera is arranged on the front panel of the terminal, and the rear camera is arranged on the back of the terminal. In some embodiments, the at least two rear cameras are any one of a main camera, a depth camera, a wide-angle camera and a tele camera, so as to realize that the main camera and the depth camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize a panoramic shooting and Virtual Reality (VR) shooting function or other fusion shooting functions. In some embodiments, camera assembly 1006 may also include a flash. The flash lamp can be a single-color temperature flash lamp or a double-color temperature flash lamp. The dual-color temperature flash lamp refers to a combination of a warm light flash lamp and a cold light flash lamp, and can be used for light compensation under different color temperatures.

The audio circuit 1007 may include a microphone and a speaker. The microphone is used for collecting sound waves of users and environments, converting the sound waves into electric signals, and inputting the electric signals to the processor 1001 for processing, or inputting the electric signals to the radio frequency circuit 1004 for voice communication. For purposes of stereo acquisition or noise reduction, the microphone may be multiple, each located at a different portion of terminal 1000. The microphone may also be an array microphone or an omni-directional pickup microphone. The speaker is used to convert electrical signals from the processor 1001 or the radio frequency circuit 1004 into sound waves. The speaker may be a conventional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only the electric signal can be converted into a sound wave audible to humans, but also the electric signal can be converted into a sound wave inaudible to humans for ranging and other purposes. In some embodiments, audio circuit 1007 may also include a headphone jack.

The location component 1008 is used to locate the current geographic location of terminal 1000 to enable navigation or LBS (Location Based Service, location-based services). The positioning component 1008 may be a positioning component based on the united states GPS (Global Positioning System ), the beidou system of china, the russian graver positioning system, or the galileo positioning system of the european union.

Power supply 1009 is used to power the various components in terminal 1000. The power source 1009 may be alternating current, direct current, disposable battery or rechargeable battery. When the power source 1009 includes a rechargeable battery, the rechargeable battery may be a wired rechargeable battery or a wireless rechargeable battery. The wired rechargeable battery is a battery charged through a wired line, and the wireless rechargeable battery is a battery charged through a wireless coil. The rechargeable battery may also be used to support fast charge technology.

In some embodiments, terminal 1000 can further include one or more sensors 1010. The one or more sensors 1010 include, but are not limited to, an acceleration sensor 1011, a gyroscope sensor 1012, a pressure sensor 1013, a fingerprint sensor 1014, an optical sensor 1015, and a proximity sensor 1016.

The acceleration sensor 1011 can detect the magnitudes of accelerations on three coordinate axes of the coordinate system established with the terminal 1000. For example, the acceleration sensor 1011 may be used to detect components of gravitational acceleration in three coordinate axes. The processor 1001 may control the display screen 1005 to display a user interface in a landscape view or a portrait view according to the gravitational acceleration signal acquired by the acceleration sensor 1011. The acceleration sensor 1011 may also be used for the acquisition of motion data of a game or a user.

The gyro sensor 1012 may detect the body direction and the rotation angle of the terminal 1000, and the gyro sensor 1012 may collect the 3D motion of the user to the terminal 1000 in cooperation with the acceleration sensor 1011. The processor 1001 can realize functions such as motion sensing (e.g., changing a UI according to a tilting operation of a user), image stabilization at photographing, game control, and inertial navigation, based on data acquired by the gyro sensor 1012.

Pressure sensor 1013 may be disposed on a side frame of terminal 1000 and/or on an underlying layer of display 1005. When the pressure sensor 1013 is provided at a side frame of the terminal 1000, a grip signal of the terminal 1000 by a user can be detected, and the processor 1001 performs right-and-left hand recognition or quick operation according to the grip signal collected by the pressure sensor 1013. When the pressure sensor 1013 is provided at the lower layer of the display screen 1005, the processor 1001 controls the operability control on the UI interface according to the pressure operation of the user on the display screen 1005. The operability controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.

The fingerprint sensor 1014 is used to collect a fingerprint of the user, and the processor 1001 identifies the identity of the user based on the fingerprint collected by the fingerprint sensor 1014, or the fingerprint sensor 1014 identifies the identity of the user based on the collected fingerprint. Upon recognizing that the user's identity is a trusted identity, the processor 1001 authorizes the user to perform relevant sensitive operations including unlocking the screen, viewing encrypted information, downloading software, paying for and changing settings, etc. Fingerprint sensor 1014 may be disposed on the front, back, or side of terminal 1000. When a physical key or vendor Logo is provided on terminal 1000, fingerprint sensor 1014 may be integrated with the physical key or vendor Logo.

The optical sensor 1015 is used to collect ambient light intensity. In one embodiment, the processor 1001 may control the display brightness of the display screen 1005 based on the ambient light intensity collected by the optical sensor 1015. Specifically, the display luminance of the display screen 1005 is turned up when the ambient light intensity is high, and the display luminance of the display screen 1005 is turned down when the ambient light intensity is low. In another embodiment, the processor 1001 may dynamically adjust the shooting parameters of the camera module 1006 according to the ambient light intensity collected by the optical sensor 1015.

Proximity sensor 1016, also known as a distance sensor, is disposed on the front panel of terminal 1000. Proximity sensor 1016 is used to collect the distance between the user and the front of terminal 1000. In one embodiment, processor 1001 controls display 1005 to switch from the on-screen state to the off-screen state when proximity sensor 1016 detects a gradual decrease in the distance between the user and the front of terminal 1000, and processor 1001 controls display 1005 to switch from the off-screen state to the on-screen state when proximity sensor 1016 detects a gradual increase in the distance between the user and the front of terminal 1000.

Those skilled in the art will appreciate that the structure shown in fig. 10 is not limiting and that terminal 1000 can include more or fewer components than shown, or certain components can be combined, or a different arrangement of components can be employed.

Optionally, the computer device is provided as a server. Fig. 11 is a schematic structural diagram of a server according to an embodiment of the present application, where the server 1100 may have a relatively large difference due to configuration or performance, and may include one or more processors (Central Processing Units, CPU) 1101 and one or more memories 1102, where the memories 1102 store at least one computer program, and the at least one computer program is loaded and executed by the processors 1101 to implement the methods according to the above-mentioned method embodiments. Of course, the server may also have a wired or wireless network interface, a keyboard, an input/output interface, and other components for implementing the functions of the device, which are not described herein.

The embodiment of the application also provides a computer readable storage medium, wherein at least one computer program is stored in the computer readable storage medium, and the at least one computer program is loaded and executed by a processor to realize the operations performed in the bottom hole temperature acquisition method of the embodiment.

Embodiments of the present application also provide a computer program product or computer program comprising computer program code stored in a computer readable storage medium. The processor of the computer device reads the computer program code from the computer readable storage medium, and the processor executes the computer program code so that the computer device realizes the operations performed in the bottom hole temperature acquisition method of the above embodiment.

It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program for instructing relevant hardware, where the program may be stored in a computer readable storage medium, and the above storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.

The foregoing is merely an alternative embodiment of the present application and is not intended to limit the embodiment of the present application, and any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the embodiment of the present application should be included in the protection scope of the present application.

Claims

1. A method of acquiring a bottom hole temperature, the method comprising:

determining a bottom hole temperature of the natural gas well from the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold;

wherein the first relationship data comprises:

Wherein p is the initial pressure, ρ is the density of the natural gas, R is a universal gas constant, T is the initial temperature, α, γ, a, b, c, and d are preset parameters, and a ₀、B₀、C₀、D₀ and E ₀ are known parameters determined based on the plurality of attribute parameters;

The second relationship data includes:

2. The method of claim 1, wherein the determining the bottom hole temperature of the natural gas well from the updated temperature parameter in response to the difference between the target density and the initial density being less than a reference threshold comprises:

Updating an initial pressure parameter according to the reference pressure variable to obtain an updated pressure parameter in response to the difference between the target density and the initial density being less than the reference threshold;

3. The method of claim 1, wherein the updating the initial temperature based on the first influence coefficient and the reference pressure variable, after which the updated temperature is obtained, further comprises:

4. The method of claim 1, wherein the natural gas comprises a plurality of constituents, the acquiring the natural gas data comprising:

5. The method of claim 1, wherein the acquiring natural gas data comprises:

displaying a parameter setting interface;

and acquiring the natural gas data input in the parameter setting interface.

6. The method of claim 5, wherein the parameter setting interface comprises a static parameter setting area for setting the initial temperature, a fluid parameter setting area for setting the plurality of attribute parameters, and a production parameter setting area for setting the initial pressure.

7. A downhole temperature acquisition device, the device comprising:

a bottom hole temperature acquisition module for determining a bottom hole temperature of the natural gas well from the updated temperature in response to the difference between the target density and the initial density being less than a reference threshold;

wherein the first relationship data comprises:

The second relationship data includes:

8. A computer device comprising a processor and a memory, wherein the memory has stored therein at least one computer program that is loaded and executed by the processor to perform the operations performed in the bottom hole temperature acquisition method of any one of claims 1 to 6.

9. A computer readable storage medium having stored therein at least one computer program loaded and executed by a processor to perform the operations performed in the bottom hole temperature acquisition method of any one of claims 1 to 6.