A single-phase flow-multiphase flow coupling calculation method and device for gas field gathering and transportation pipeline network

Abstract

This article relates to the field of oil and gas field development engineering, and in particular to a single-phase flow-multiphase flow coupling calculation method and device for a gas field gathering and transportation network. This method includes performing hydraulic and thermodynamic coupling iterations based on the basic parameters of the gas field gathering and transportation network, the set initial value of the pipeline volume flow of the gas gathering network, the set initial value of the pipeline mass flow of the gas production network, the set initial value of the node temperature, and the physical properties of natural gas, until the pipeline volume flow, pipeline mass flow, and node temperature converge, and then using the obtained node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature as the coupling calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network. Through the embodiments of this article, the hydraulic and thermodynamic coupling calculation of the multiphase flow gas production network and the single-phase flow gas gathering network is realized, thereby guiding the exploitation of the gas field based on the obtained node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature.

Description

Technical Field

[0001] This article relates to the field of oil and gas field development engineering, and in particular to a single-phase flow-multiphase flow coupling calculation method and device for a gas field gathering and transportation pipeline network. Background Art

[0002] Gas field gathering and transportation systems are crucial infrastructure for gas field development. Their pipeline networks are characterized by large scale, numerous energy-consuming nodes, and a complex network structure. Their production and operating parameters are numerous, rapidly changing, complex, and poorly controllable. Furthermore, as gas field production progresses, the landscape evolves, and multiple development approaches coexist. As most gas fields enter the middle and late stages of development, boosting projects increase, leading to a year-on-year increase in energy consumption. Compressors are the primary energy consumers. Therefore, to optimize and reduce energy consumption in gas field gathering and transportation systems, steady-state simulation analysis combining the gathering and transportation pipeline network and compressors is essential.

[0003] Reciprocating compressors are the preferred choice for gas field boosting projects due to their high exhaust pressure, stable exhaust, and ability to achieve high-pressure conditions with small gas volumes, making them suitable for actual gas field production. However, existing simulation models are mostly targeted at centrifugal compressors, lacking models and calculation methods suitable for gas field gathering and transmission pipeline boosting projects. Furthermore, existing steady-state simulation technology suffers from weaknesses in pipeline network hydraulic calculations, such as weak computational stability and sensitivity to iterative initial values. It also simplifies the treatment of equipment control equations, limiting calculations to fixed equipment inlet and outlet pressures or pressure ratios. This makes it unsuitable for adjusting compressor speed, exhaust volume, and other various operating conditions. Furthermore, research on the coupled calculations of reciprocating compressors and gathering and transmission pipelines is lacking, making it difficult to provide precise guidance for gas field production.

[0004] There is an urgent need for a single-phase flow-multiphase flow coupling calculation method for gas field gathering and transportation pipelines, so as to solve the problem that the existing technology lacks research on the coupling calculation of reciprocating compressors and gathering and transportation pipelines, and cannot provide accurate guidance for gas field production. Summary of the Invention

[0005] In order to solve the problems in the existing technology, the embodiments of this article provide a single-phase flow-multiphase flow coupling calculation method and device for a gas field gathering and transportation pipeline network, and establish a set of gas field gathering and transportation pipeline network calculation methods that include single-phase flow and multiphase flow in the pipeline network. Secondly, a unified model solution algorithm is used for the multiphase flow gas production pipeline network and the single-phase flow gas gathering pipeline network, and a coupled simulation method for the gas-liquid two-phase flow of the gas production pipeline network and the single-phase flow of the gas gathering pipeline network is proposed, which improves the correlation of the pipeline network simulation parameters and makes the simulation results closer to production reality.

[0006] In order to solve the above technical problems, the specific technical solutions of this article are as follows:

[0007] On the one hand, the embodiment of this invention provides a single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation network, including:

[0008] Calculating the pipeline friction and admittance matrix of the gas field gathering and transportation network according to the basic parameters of the gas field gathering and transportation network, the set initial value of the pipeline volume flow of the gas gathering and transportation network, the set initial value of the pipeline mass flow of the gas production and transportation network, the set initial value of the node temperature, and the physical properties of natural gas, and generating a coefficient matrix;

[0009] Obtaining the node pressure and pipeline volume flow of the gas gathering network, and the node pressure and pipeline mass flow of the gas production network according to the pipeline friction resistance, admittance matrix and coefficient matrix;

[0010] Performing thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network;

[0011] Recalculating the natural gas physical properties according to the BWRS state equation and the node pressures of the gas gathering network, the node pressures of the gas production network, and the node temperatures;

[0012] Determining whether a difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether a difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether a difference between the node temperature and the initial value of the node temperature is less than a third threshold;

[0013] If not, taking the pipeline volume flow as the pipeline volume flow initial value, taking the pipeline mass flow as the pipeline mass flow initial value, taking the node temperature as the node temperature initial value, and calculating the pipeline friction resistance and admittance matrix of the gas field gathering and transportation pipeline network according to the pipeline network basic parameters, the pipeline volume flow initial value, the pipeline mass flow initial value, the set node temperature initial value, and the recalculated natural gas physical property parameters, and generating a coefficient matrix;

[0014] If so, the obtained node pressure, node pressure and node temperature of the gas gathering network are used as the coupling calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network, so as to guide the production of the gas field according to the coupling calculation results.

[0015] Furthermore, the basic parameters of the pipe network include the diameter of the gas collection pipe D 集 , Gas collecting pipe length L 集 , gas gathering network friction coefficient λ 集 , Gas pipe diameter D 采 , Gas pipe length L 采 , friction coefficient of gas pipeline network λ 采, a pipe section node association matrix, the natural gas physical property parameters include the natural gas compressibility factor Z, the natural gas relative density Δ, and the natural gas average temperature T, wherein the natural gas average temperature T is calculated based on the set node temperature initial value and the pipe section node association matrix, the pipeline friction resistance of the gas field gathering and transportation network includes the gas gathering network friction resistance and the gas production network friction resistance, the admittance matrix includes the gas gathering network admittance matrix and the gas production network admittance matrix, and the coefficient matrix includes the gas gathering network coefficient matrix and the gas production network coefficient matrix.

[0016] Furthermore, the steps of calculating the pipeline friction and admittance matrix of the gas field gathering and transportation network according to the basic parameters of the gas field gathering and transportation network, the set initial value of the pipeline volume flow of the gas gathering and transportation network, the set initial value of the pipeline mass flow of the gas production and transportation network, the set initial value of the node temperature, and the natural gas physical properties, and generating a coefficient matrix include:

[0017] According to the formula Calculate the friction resistance of the gas gathering network, where S 集 is the friction resistance of the gas collecting pipe network, and C is a constant;

[0018] According to the formula Calculate the gas collection network admittance matrix, where G 集 is the admittance matrix of the gas gathering network, Q is the initial value of the pipeline volume flow rate, and α is the flow index;

[0019] According to formula Y 集 =A×G 集 ×A T Calculate the gas collection network coefficient matrix, where Y 集 is the gas gathering network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose;

[0020] According to the formula Calculate the friction resistance of the gas production network, where S 采 is the friction resistance of the gas production network, C1 is a constant, x is the gas mass content, ρ tp is the density of the mixture in the gas pipeline network, ρ g is the gas phase density, M represents the initial value of the pipeline mass flow rate, and p is the average absolute pressure of the pipeline;

[0021] According to the formula Calculate the gas production network admittance matrix, where G 采 is the admittance matrix of the gas production network, M is the initial value of the pipeline mass flow rate, and α is the flow index;

[0022] According to formula Y 采 =A×G 采 ×A TCalculate the gas production network coefficient matrix, where Y 采 is the gas production network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose.

[0023] Furthermore, the natural gas physical property parameters also include the compression factor of the natural gas at the compressor inlet, the natural gas density, the specific heat capacity of the gas at the inflow node, and the specific heat capacity of the gas at the compressor outlet.

[0024] Furthermore, the step of obtaining the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction, admittance matrix and coefficient matrix includes:

[0025] According to the formula:

[0026] Y 集 P 集 =-q,

[0027] Calculate the node pressure and pipeline volume flow of the gas gathering network, where P 集 is the node pressure of the gas gathering network, q is the pipeline volume flow rate;

[0028] According to the formula:

[0029] Y 采 P 采 =-m,

[0030] Calculate the node pressure and pipeline mass flow of the gas production network, where P 采 is the node pressure of the gas pipeline network, and m is the pipeline mass flow rate.

[0031] Furthermore, performing thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network further includes:

[0032] Constructing a correlation matrix between pipeline ends and nodes according to the pipeline segment node correlation matrix and the pipeline volume flow;

[0033] According to the correlation matrix between the pipeline end and the node and the pipeline mass flow, the formula is used:

[0034]

[0035] Calculate the node temperature of the gas pipeline network, where T i,采 is the node temperature of the gas pipeline network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,采 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, Tin,j,采 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,采 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,采 It is calculated based on the mass flow of the pipeline, and the outflow node flow m out,i,采 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,采 The inflow node temperature at the end of the pipeline includes a set known node temperature and an unknown node temperature, wherein the unknown node temperature is the set initial value of the node temperature;

[0036] According to the correlation matrix between the pipe end and the node and the pipe volume flow, the formula is used:

[0037]

[0038] Calculate the node temperature of the gas pipeline network, where T i,集 is the node temperature of the gas gathering network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,集 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, T in,j,集 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,集 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,集 It is calculated based on the pipeline volume flow and the natural gas density in the pipeline, and the outflow node flow m out,i,集 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,集 The inflow node temperature at the end of the pipeline includes the set known node temperature and unknown node temperature, wherein the unknown node temperature is the gas network node temperature T i,采 .

[0039] Furthermore, the step of constructing a pipeline end-node association matrix based on the pipeline segment node association matrix and the pipeline volume flow includes:

[0040] Determining that the pipeline volume flow rate of the middle pipeline of the gas field gathering and transportation pipeline network is less than 0;

[0041] If so, the element of the pipeline in the pipeline segment node association matrix is ​​used as the element of the pipeline in the pipeline end and node association matrix;

[0042] If not, the inverse numbers of the elements of the pipeline in the pipeline segment node association matrix are used as the elements of the pipeline in the pipeline end and node association matrix;

[0043] According to the structure of the gas field gathering and transportation pipeline network, the value of the element corresponding to the inlet pipeline in the pipeline-node association matrix is ​​set to 1.

[0044] On the other hand, the embodiment of this invention also provides a single-phase flow-multiphase flow coupling calculation device for a gas field gathering and transportation network, comprising:

[0045] It includes: iterative parameter calculation unit, pipeline flow calculation unit, node temperature calculation unit, natural gas physical parameter calculation unit and result verification unit;

[0046] The iterative parameter calculation unit is used to calculate the pipeline friction resistance and admittance matrix of the gas field gathering and transportation pipeline network according to the basic parameters of the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical property parameters, and generate a coefficient matrix;

[0047] The pipeline flow calculation unit is used to obtain the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction resistance, admittance matrix and coefficient matrix;

[0048] The node temperature calculation unit is used to perform thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network;

[0049] The natural gas physical property parameter calculation unit is used to recalculate the natural gas physical property parameters according to the BWRS state equation and the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature;

[0050] The result verification unit is used to determine whether the difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether the difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether the difference between the node temperature and the initial value of the node temperature is less than a second threshold; if not, the pipeline volume flow is used as the initial value of the pipeline volume flow, the pipeline mass flow is used as the initial value of the pipeline mass flow, and the node temperature is used as the initial value of the node temperature, and the pipeline friction and admittance matrix of the gas field gathering and transportation network is calculated according to the basic parameters of the pipeline network, the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, the set initial value of the node temperature, and the recalculated natural gas physical properties, and a coefficient matrix is ​​generated. If so, the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature are used as the coupled calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network, so as to guide the exploitation of the gas field according to the coupled calculation results.

[0051] On the other hand, an embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above method when executing the computer program.

[0052] Finally, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

[0053] Using the embodiment of this article, first, according to the basic parameters of the pipeline network in the gas field gathering and transportation network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical properties, the pipeline friction resistance and admittance matrix of the gas field gathering and transportation pipeline network are calculated, and a coefficient matrix is ​​generated. Then, hydraulic calculations are performed on the gas gathering pipeline network and the gas production pipeline network in the gas field gathering and transportation pipeline network respectively to obtain the node pressure and pipeline volume flow of the gas gathering pipeline network, and the node pressure and pipeline mass flow of the gas production pipeline network. Then, the pipeline volume flow of the gas gathering pipeline network and the pipeline mass flow of the gas production pipeline network are used to perform thermal calculations to obtain the node temperature of the gas field gathering and transportation pipeline network, and the calculated node pressure of the gas gathering pipeline network, the node pressure of the gas production pipeline network and the node The natural gas physical properties are recalculated based on the temperature. Then, it is determined whether the calculated pipeline volume flow, pipeline mass flow, and node pressure are less than their respective preset threshold values. If not, the calculated pipeline volume flow, pipeline mass flow, and node temperature are used as the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, and the initial value of the node temperature. The hydraulic calculation of the gas gathering network and the gas production network and the thermal calculation of the gas field gathering network are re-performed until the obtained pipeline volume flow, pipeline mass flow, and node pressure are greater than or equal to their respective preset threshold values. The hydraulic and thermal coupling iteration of the gas field gathering network is terminated. Finally, the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature are used as the single-phase flow-multiphase flow coupling calculation results of the gas field gathering network. A coupled solution method for the multiphase flow gas production network and the unidirectional flow gas gathering network is realized, which improves the correlation of the network calculation parameters and makes the calculation results closer to the actual production. Finally, the node temperature of the gas field gathering network is calculated according to the pipeline mass flow of the multiphase flow gas production network and the pipeline volume flow of the single-phase flow gas gathering network, and the thermal calculation of the gas field gathering network is realized based on the results of the hydraulic calculation of the gas field gathering network, so as to guide the gas field exploitation according to the obtained node pressure of the gas gathering network, node pressure of the gas production network and node temperature. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the embodiments of this article or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of this article. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0055] Figure 1 Schematic diagram of an implementation system of a single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation network in an embodiment of this invention;

[0056] Figure 2Schematic diagram of a flow chart of a single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation network according to an embodiment of this invention;

[0057] Figure 3 The figure shows the process of constructing the correlation matrix between the pipeline end and the node according to the pipeline segment node correlation matrix and the pipeline volume flow in the embodiment of this article;

[0058] Figure 4 FIG2 is a schematic diagram of the structure of a single-phase flow-multiphase flow coupling calculation device for a gas field gathering and transportation network according to an embodiment of the present invention;

[0059] Figure 5 Shown is a schematic structural diagram of a computer device according to an embodiment of this invention.

[0060]

Description of the accompanying drawings

[0061] 101. Terminal;

[0062] 102. Server;

[0063] 401, iterative parameter calculation unit;

[0064] 402. Pipeline flow calculation unit;

[0065] 403. Node temperature calculation unit;

[0066] 404. Natural gas physical property parameter calculation unit;

[0067] 405. Result verification unit;

[0068] 502. Computer equipment;

[0069] 504, processor;

[0070] 506. Memory;

[0071] 508, driving mechanism;

[0072] 510, input / output module;

[0073] 512. Input devices;

[0074] 514. Output device;

[0075] 516. Presentation equipment;

[0076] 518. Graphical User Interface;

[0077] 520, network interface;

[0078] 522, communication link;

[0079] 524. Communication bus. DETAILED DESCRIPTION

[0080] The following will be combined with the accompanying drawings to clearly and completely describe the technical solutions in the embodiments of this document. Obviously, the embodiments described are only part of the embodiments of this document, not all of the embodiments. Based on the embodiments of this document, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this document.

[0081] It should be noted that the terms "first," "second," and the like in the specification and claims herein and in the accompanying drawings are used to distinguish similar objects and are not necessarily used to describe a particular order or precedence. It should be understood that the terms used in this manner are interchangeable where appropriate so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having," as well as any variations thereof, are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product, or device comprising a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units that are not explicitly listed or that are inherent to these processes, methods, products, or devices.

[0082] like Figure 1 The figure shows a schematic diagram of an implementation system of a single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation network according to an embodiment of the present invention, which may include a terminal 101 and a server 102. A communication connection is established between the terminal 101 and the server 102 to enable data interaction. The terminal 101 can input the structure of the gas field gathering and transportation pipeline network, the basic parameters of the pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, and the set initial value of the node temperature to the server 102, and also input the threshold value required for judgment. The server 102 iteratively calculates the pipeline flow and node temperature of the gas field gathering and transportation pipeline network based on the basic parameters of the pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, and the set initial value of the node temperature, and judges whether the calculated pipeline volume flow, pipeline mass flow and node temperature meet the requirements based on the corresponding threshold value. If so, the obtained pipeline volume flow, pipeline mass flow and node temperature are used as the coupled calculation results of the unidirectional flow-multiphase flow of the field gathering and transportation pipeline network, thereby guiding the exploitation of the gas field.

[0083] In the embodiments of this specification, the server 102 can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers. It can also be a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, content delivery networks (CDN, Content Delivery Network), and big data and artificial intelligence platforms.

[0084] In an optional embodiment, the terminal 101 may include but is not limited to electronic devices such as desktop computers, tablet computers, and laptop computers. Optionally, the operating system running on the electronic device may include but is not limited to Android, iOS, Linux, and Windows.

[0085] In addition, it should be noted that Figure 1 What is shown is only one application environment provided by the present disclosure. In actual applications, other application environments may also be included, and this specification does not limit them.

[0086] In order to solve the problems existing in the prior art, the embodiment of this article provides a single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation pipeline network, and establishes a set of gas field gathering and transportation pipeline network calculation methods that include single-phase flow and multiphase flow in the pipeline network. Secondly, a unified model solution algorithm is used for the multiphase flow gas production pipeline network and the single-phase flow gas gathering pipeline network, and a coupled simulation method for the gas-liquid two-phase flow of the gas production pipeline network and the single-phase flow of the gas gathering pipeline network is proposed, which improves the correlation of the pipeline network simulation parameters and makes the simulation results closer to production reality. Figure 2 Shown is a flow chart of a method for calculating the coupling of single-phase flow and multi-phase flow in a gas field gathering and transportation network according to an embodiment of this invention. This figure describes the process of calculating the coupling of single-phase flow and multi-phase flow in a gas field gathering and transportation network, but it may include more or fewer operation steps based on conventional or non-creative work. The order of steps listed in the embodiment is only one way of executing the steps among many, and does not represent the only execution order. When the system or device product is actually executed, it may be executed in sequence or in parallel according to the method shown in the embodiment or the accompanying drawings. Specifically, as Figure 2 As shown, the method may be executed by a processing unit on a server, and may include:

[0087] Step 201: Calculating the pipeline friction and admittance matrices of the gas field gathering and transportation network based on the basic parameters of the gas field gathering and transportation network, the set initial values ​​of the pipeline volume flow of the gas gathering and transportation network, the set initial values ​​of the pipeline mass flow of the gas production and transportation network, the set initial values ​​of the node temperature, and the natural gas physical properties, and generating a coefficient matrix;

[0088] In this step, the basic parameters of the pipe network include the gas collection pipe diameter D集 , Gas collecting pipe length L 集 , gas gathering network friction coefficient λ 集 , Gas pipe diameter D 采 , Gas pipe length L 采 , friction coefficient of gas pipeline network λ 采 , a pipe section node association matrix, the natural gas physical property parameters include the natural gas compressibility factor Z, the natural gas relative density Δ, and the natural gas average temperature T, wherein the natural gas average temperature T is calculated based on the set node temperature initial value and the pipe section node association matrix, the pipeline friction resistance of the gas field gathering and transportation network includes the gas gathering network friction resistance and the gas production network friction resistance, the admittance matrix includes the gas gathering network admittance matrix and the gas production network admittance matrix, and the coefficient matrix includes the gas gathering network coefficient matrix and the gas production network coefficient matrix.

[0089] In this step, the natural gas physical property parameters also include the compression factor of the natural gas at the compressor inlet, the natural gas density, the specific heat capacity of the gas at the inflow node, and the specific heat capacity of the gas at the compressor outlet.

[0090] The steps of calculating the pipeline friction and admittance matrix of the gas field gathering and transportation pipeline network according to the basic parameters of the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical properties, and generating a coefficient matrix include:

[0091] According to formula (1):

[0092]

[0093] Calculate the friction resistance of the gas gathering network, where S 集 is the friction resistance of the gas collecting pipe network, and C is a constant;

[0094] According to formula (2):

[0095]

[0096] Calculate the gas collection network admittance matrix, where G 集 is the admittance matrix of the gas gathering network, Q is the initial value of the pipeline volume flow rate, and α is the flow index;

[0097] According to formula (3):

[0098] Y 集 =A×G 集 ×A T (3)

[0099] Calculate the gas collection network coefficient matrix, where Y 集is the gas gathering network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose;

[0100] According to formula (4):

[0101]

[0102] Calculate the friction resistance of the gas production network, where S 采 is the friction resistance of the gas production network, C1 is a constant, x is the gas mass content, ρ tp is the density of the mixture in the gas pipeline network, ρ g is the gas phase density, M represents the initial value of the pipeline mass flow rate, and p is the average absolute pressure of the pipeline;

[0103] According to formula (5):

[0104]

[0105] Calculate the gas production network admittance matrix, where G 采 is the admittance matrix of the gas production network, M is the initial value of the pipeline mass flow rate, and α is the flow index;

[0106] According to formula (6):

[0107] Y 采 =A×G 采 ×A T (6)

[0108] Calculate the gas production network coefficient matrix, where Y 采 is the gas production network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose.

[0109] Step 202: Obtaining the node pressure and pipeline volume flow of the gas gathering network, and the node pressure and pipeline mass flow of the gas production network according to the pipeline friction, admittance matrix, and coefficient matrix;

[0110] In this step, the step of obtaining the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction, admittance matrix and coefficient matrix includes:

[0111] According to formula (7):

[0112] Y 集 P 集 =-q (7)

[0113] Calculate the node pressure and pipeline volume flow of the gas gathering network, where P 集 is the node pressure of the gas gathering network, q is the pipeline volume flow rate;

[0114] Specifically, the volume flow of branch pipe networks is known or easily calculated, while the volume flow of ring pipe networks is often unknown and difficult to calculate. The node pressures of a gas gathering network include both known and unknown node pressures. Typically, the pressure of at least one reference node in the network must be known before the solution can be obtained. Using the linear approximation method, the unknown node pressures of the gas gathering network can be calculated using Equation (7).

[0115] According to formula (8):

[0116] Y 采 P 采 =-m (8)

[0117] Calculate the node pressure and pipeline mass flow of the gas production network, where P 采 is the pressure of the gas pipeline network node, and m is the pipeline mass flow rate. Specifically, the linear approximation method is used to solve the pressure at the unknown node and the pipeline mass flow rate of equation group (8).

[0118] Step 203: performing thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network;

[0119] In this step, the gas field gathering and transportation network is thermally calculated based on the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network further includes:

[0120] Constructing a correlation matrix between pipeline ends and nodes according to the pipeline segment node correlation matrix and the pipeline volume flow;

[0121] Specifically, if Figure 3 As shown, the process of constructing the pipeline end and node association matrix based on the pipeline segment node association matrix and the pipeline volume flow further includes:

[0122] Step 301: determining whether the pipeline volume flow rate of the middle pipeline of the gas field gathering and transportation network is less than 0;

[0123] Step 302: If yes, use the element of the pipeline in the pipeline segment node association matrix as the element of the pipeline in the pipeline end and node association matrix;

[0124] Step 303: If not, the inverse numbers of the elements of the pipeline in the pipeline segment node association matrix are used as the elements of the pipeline in the pipeline end and node association matrix;

[0125] Step 304: according to the structure of the gas field gathering and transportation pipeline network, the value of the element corresponding to the inlet pipeline in the pipeline-node association matrix is ​​set to 1.

[0126] According to the correlation matrix between the pipeline end and the node and the pipeline mass flow, formula (9) is used:

[0127]

[0128] Calculate the node temperature of the gas pipeline network, where T i,采 is the node temperature of the gas pipeline network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,采 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, T in,j,采 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,采 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,采 It is calculated based on the mass flow of the pipeline, and the outflow node flow m out,i,采 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,采 The inflow node temperature at the end of the pipeline includes a set known node temperature and an unknown node temperature, wherein the unknown node temperature is the set initial value of the node temperature;

[0129] According to the correlation matrix between the pipe end and the node and the pipe volume flow, formula (10) is used:

[0130]

[0131] Calculate the node temperature of the gas pipeline network, where T i,集 is the node temperature of the gas gathering network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,集 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, T in,j,集 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,集 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,集It is calculated based on the pipeline volume flow and the natural gas density in the pipeline, and the outflow node flow m out,i,集 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,集 The inflow node temperature at the end of the pipeline includes the set known node temperature and unknown node temperature, wherein the unknown node temperature is the gas network node temperature T i,采 .

[0132] Step 204: recalculating the natural gas physical properties according to the BWRS state equation and the node pressures of the gas gathering network, the node pressures of the gas production network, and the node temperatures;

[0133] It should be noted that the calculation method of the BWRS state equation described in this step is a common knowledge method in the art and will not be repeated in the examples of this specification.

[0134] Step 205: Determine whether the difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether the difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether the difference between the node temperature and the initial value of the node temperature is less than a third threshold;

[0135] In the embodiment of this article, if the difference between the pipeline volume flow obtained in this iteration and the pipeline volume flow obtained in the previous iteration (i.e., the initial value of the pipeline volume flow) is less than the first threshold, the difference between the pipeline mass flow obtained in this iteration and the pipeline mass flow obtained in the previous iteration (i.e., the initial value of the pipeline mass flow) is less than the second threshold, and the difference between the node temperature obtained in this iteration and the node temperature obtained in the previous iteration (i.e., the initial value of the node) is less than the third threshold, it means that the iterative calculation has converged, and the pipeline volume flow, pipeline mass flow and node pressure obtained in this iteration can be used as the single-phase flow-multiphase flow coupling calculation results of the gas field gathering and transportation pipeline network.

[0136] In this step, if not, execute step 206 , otherwise execute step 207 .

[0137] Step 206: using the pipeline volume flow rate as the pipeline volume flow initial value, the pipeline mass flow rate as the pipeline mass flow initial value, and the node temperature as the node temperature initial value, re-execute step 201 to iteratively calculate the pipeline volume flow rate of the gas gathering network, the pipeline mass flow rate of the gas production network, and the node temperature;

[0138] It should be noted that in this step, in addition to using the calculated pipeline volume flow as the initial value of the pipeline volume flow in step 201, using the calculated pipeline mass flow as the initial value of the pipeline mass flow in step 201, and using the calculated node temperature as the initial value of the node temperature in step 201, it is also necessary to replace the natural gas physical property parameters recalculated in step 204 with the natural gas physical property parameters in step 201. That is, during the iterative process, only the basic pipeline network parameters in step 201 remain unchanged, while the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, the initial value of the node temperature, and the natural gas physical property parameters all change with the iteration. Then, based on the iteratively calculated gas gathering network friction resistance, gas production network friction resistance, gas gathering network admittance matrix, gas production network admittance matrix, gas gathering network coefficient matrix, and gas production network coefficient matrix, the node pressure, pipeline volume flow rate, and pipeline mass flow rate of the gas gathering network are obtained. Based on the pipeline volume flow rate and pipeline mass flow rate, a thermal calculation is performed on the gas field gathering and transportation network to obtain the node temperature. The natural gas physical properties are recalculated based on the BWRS state equation and the node pressure and node temperature. Step 205 is then executed to determine whether the difference between the pipeline volume flow rate and the initial value of the pipeline volume flow rate is less than a first threshold, whether the difference between the pipeline mass flow rate and the initial value of the pipeline mass flow rate is less than a second threshold, and whether the difference between the node temperature and the initial value of the node temperature is less than a third threshold.

[0139] Step 207: The calculated node pressures of the gas gathering network, the node pressures of the gas production network, and the node temperatures are used as the coupled calculation results of the unidirectional flow and multiphase flow of the gas field gathering and transportation network.

[0140] In the embodiments of this article, after obtaining the coupled calculation results of the unidirectional flow and multiphase flow of the gas field gathering and transportation network, the coupled calculation results can be used to guide the production of the gas field.

[0141] Taking the implementation scenario of predicting whether there is a risk of hydrate formation in a gas field gathering and transportation network as an example, it can be determined whether there is a risk of hydrate formation in the gas field gathering and transportation network based on the obtained node pressures, node pressures, and node temperatures of the gas gathering and production network;

[0142] In this step, judging whether the gas field gathering and transportation network has a risk of hydrate generation based on the obtained node pressures, node pressures, and node temperatures of the gas gathering and production network further includes:

[0143] According to formula (11),

[0144]

[0145] Calculate the pressure at a distance x from the starting point of the pipe, where p x is the pressure at the distance x from the starting point of the pipeline, p L is the starting pressure of the pipeline, pR is the pressure at the end of the pipeline, L represents the length of the pipeline, and the pressure at the starting point of the pipeline p L , and the pressure at the end of the pipeline p R It is calculated based on the node pressure of the gas gathering network and the node pressure of the gas production network;

[0146] According to formulas (12)-(14),

[0147]

[0148]

[0149]

[0150] Calculate the critical temperature of hydrate formation at a distance x from the starting point of the pipeline, where T represents the critical temperature of hydrate formation at a distance x from the starting point of the pipeline, and f i is the fugacity of mixed gas component i, where f i The pressure p at the distance x from the starting point of the pipeline is calculated using the BWRS state equation. x The calculated value, f i 0 is the fugacity of pure basic hydrate component i, θ i is the filling rate of gas component i in the hydrate connection pores, C i is the Langmuir constant, where x i is the molar percentage of basic hydrate and mixed basic hydrate formed by gas component i, β is the hydrate structure type parameter, a w is the activity of water in the water-rich phase, λ2 represents the ratio of the number of pores in the hydrate to the number of water molecules, P represents the pressure value, α is the mobility index, and A ij is the binary interaction coefficient, A ij Obtained by regression of typical binary hydrate formation data, A, B, X, Y, Z are model parameters;

[0151] Using formula (15),

[0152] T 实 =T0+(T L -T0)e -ax (15)

[0153] Calculate the temperature at a distance x from the starting point of the pipe, where T 实 represents the temperature at the starting point x of the pipeline, T L is the starting temperature of the pipeline, where the starting temperature of the pipeline T L It is calculated based on the node temperature, T0 is the ambient temperature, Where K is the total heat transfer coefficient of the pipeline, c p represents the specific heat capacity at constant pressure, G represents the mass flow rate of the pipeline, and D represents the diameter of the pipeline;

[0154] Determine the temperature T at the distance x from the starting point of the pipeline 实 Is it less than the critical temperature T for hydrate formation at the distance x from the pipeline starting point? If so, there is a risk of hydrate formation at the distance x from the pipeline starting point.

[0155] In some other embodiments of the present invention, it is also possible to Figure 2 The method shown guides gas field production in actual formations. For example, a gas field gathering and transportation network in an actual formation is first modeled to obtain a segment node association matrix. Then, based on the initial values ​​of the pipeline volume flow rate, the initial values ​​of the pipeline mass flow rate, and the initial values ​​of the node temperature of the gas gathering network in the production plan, a coupled calculation is performed for the multiphase flow and single-phase flow gas gathering network to obtain the node pressures of the gas gathering network, the node pressures of the gas production network, and the node temperatures of the gas field gathering and transportation network. Based on these node pressures, the node pressures of the gas gathering network, and the node temperatures of the gas field gathering and transportation network, it is determined whether the gas field gathering and transportation network is at risk of hydrate formation. If so, the initial values ​​of the pipeline volume flow rate, the initial values ​​of the pipeline mass flow rate, and the initial values ​​of the node temperature are adjusted and the coupled calculation is performed again until hydrate formation is prevented in the gas field gathering and transportation network. At this point, gas field production in the actual formation is performed according to the corresponding initial values ​​of the pipeline volume flow rate, the initial values ​​of the pipeline mass flow rate, and the initial values ​​of the node temperature, thereby preventing hydrate formation in the gas field gathering and transportation network.

[0156] By the method of the embodiment of the present invention, firstly, the pipeline friction and admittance matrix of the gas field gathering and transportation pipeline network is calculated according to the basic parameters of the pipeline network in the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the physical properties of natural gas, and a coefficient matrix is ​​generated. Then, hydraulic calculations are performed on the gas gathering pipeline network and the gas production pipeline network in the gas field gathering and transportation pipeline network respectively to obtain the node pressure and pipeline volume flow of the gas gathering pipeline network, and the node pressure and pipeline mass flow of the gas production pipeline network. Then, thermal calculations are performed using the pipeline volume flow of the gas gathering pipeline network and the pipeline mass flow of the gas production pipeline network to obtain the node temperature of the gas field gathering and transportation pipeline network, and the calculated node pressure of the gas gathering pipeline network, the node pressure of the gas production pipeline network and the The natural gas physical properties are recalculated based on the node temperature. Then, it is determined whether the calculated pipeline volume flow, pipeline mass flow, and node pressure are less than their respective preset threshold values. If not, the calculated pipeline volume flow, pipeline mass flow, and node temperature are used as the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, and the initial value of the node temperature. The hydraulic calculation of the gas gathering network and the gas production network and the thermal calculation of the gas field gathering network are re-performed until the obtained pipeline volume flow, pipeline mass flow, and node pressure are greater than or equal to their respective preset threshold values. At this time, the hydraulic and thermal coupling iteration of the gas field gathering network is terminated. Finally, the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature are used as the single-phase flow-multiphase flow coupling calculation results of the gas field gathering network. A coupled solution method for the multiphase flow gas production network and the unidirectional flow gas gathering network is realized, which improves the correlation of the network calculation parameters and makes the calculation results closer to actual production. Finally, the node temperature of the gas field gathering network is calculated based on the pipeline mass flow of the multiphase flow gas production network and the pipeline volume flow of the single-phase flow gas gathering network. The thermal calculation of the gas field gathering network is carried out on the basis of the results of the hydraulic calculation of the gas field gathering network, and the coupled calculation results of the unidirectional flow and multiphase flow of the gas field gathering network are obtained, so as to guide the exploitation of the gas field according to the said coupled calculation results.

[0157] Based on the same inventive concept, the embodiment of the present invention also provides a single-phase flow-multiphase flow coupling calculation device for a gas field gathering and transportation network, such as Figure 4 As shown, it includes an iterative parameter calculation unit 401, a pipeline flow calculation unit 402, a node temperature calculation unit 403, a natural gas physical property parameter calculation unit 404 and a result verification unit 405;

[0158] The iterative parameter calculation unit 401 is used to calculate the pipeline friction and admittance matrix of the gas field gathering and transportation pipeline network according to the basic parameters of the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical properties, and generate a coefficient matrix;

[0159] The pipeline flow calculation unit 402 is used to obtain the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction resistance, admittance matrix and coefficient matrix;

[0160] The node temperature calculation unit 403 is used to perform thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network;

[0161] The natural gas physical property parameter calculation unit 404 is used to recalculate the natural gas physical property parameters according to the BWRS state equation and the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature;

[0162] The result verification unit 405 is used to determine whether the difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether the difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether the difference between the node temperature and the initial value of the node temperature is less than a second threshold; if not, the pipeline volume flow is used as the initial value of the pipeline volume flow, the pipeline mass flow is used as the initial value of the pipeline mass flow, and the node temperature is used as the initial value of the node temperature, and the pipeline friction and admittance matrix of the gas field gathering and transportation network is calculated according to the basic parameters of the pipeline network, the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, the set initial value of the node temperature, and the recalculated natural gas physical properties, and a coefficient matrix is ​​generated. If so, the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature are used as the coupled calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network, so as to guide the exploitation of the gas field according to the coupled calculation results.

[0163] Since the principle of solving the problem by the above device is similar to that of the above method, the implementation of the above device can refer to the implementation of the above method, and the repeated parts will not be repeated.

[0164] like Figure 5As shown, a computer device provided in an embodiment of this document is shown. The apparatus herein may be a computer device in this embodiment, executing the method described above. The computer device 502 may include one or more processors 504, such as one or more central processing units (CPUs), each of which may implement one or more hardware threads. The computer device 502 may also include any memory 506 for storing any type of information, such as code, settings, data, etc. For example, without limitation, the memory 506 may include any one or more combinations of the following: any type of RAM, any type of ROM, a flash memory device, a hard disk, an optical disk, etc. More generally, any memory may use any technology to store information. Furthermore, any memory may provide volatile or non-volatile retention of information. Furthermore, any memory may represent a fixed or removable component of the computer device 502. In one embodiment, when the processor 504 executes associated instructions stored in any memory or combination of memories, the computer device 502 may perform any operation of the associated instructions. The computer device 502 also includes one or more drive mechanisms 508, such as a hard disk drive mechanism, an optical disk drive mechanism, etc., for interacting with any memory.

[0165] The computer device 502 may also include an input / output module 510 (I / O) for receiving various inputs (via input devices 512) and for providing various outputs (via output devices 514). A specific output mechanism may include a presentation device 516 and an associated graphical user interface (GUI) 518. In other embodiments, the input / output module 510 (I / O), input devices 512, and output devices 514 may not be included, and the computer device 502 may simply be a computer device in a network. The computer device 502 may also include one or more network interfaces 520 for exchanging data with other devices via one or more communication links 522. One or more communication buses 524 couple the components described above together.

[0166] The communication link 522 may be implemented in any manner, for example, via a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 522 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.

[0167] Corresponding to Figure 2-Figure 3 The embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the above steps are performed.

[0168] The embodiment of the present invention also provides a computer readable instruction, wherein when the processor executes the instruction, the program causes the processor to execute the following Figure 2-Figure 3 The method shown.

[0169] It should be understood that in the various embodiments of this document, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this document.

[0170] It should also be understood that in the embodiments herein, the term "and / or" merely describes an association between associated objects, indicating that three possible relationships exist. For example, "A and / or B" could represent: A alone, A and B simultaneously, or B alone. Furthermore, the character " / " in this document generally indicates an "or" relationship between the associated objects.

[0171] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the composition and steps of each example according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this document.

[0172] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0173] In the several embodiments provided herein, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices, or units, or can be an electrical, mechanical, or other form of connection.

[0174] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments herein.

[0175] In addition, the functional units in the various embodiments herein may be integrated into a single processing unit, each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0176] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this article is essentially or the part that contributes to the existing technology, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of this article. The aforementioned storage medium includes: various media that can store program code, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0177] This article uses specific embodiments to illustrate the principles and implementation methods of this article. The description of the above embodiments is only used to help understand the methods and core ideas of this article. At the same time, for those skilled in the art, based on the ideas of this article, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation to this article.

Claims

1. A single-phase flow-multiphase flow coupling calculation method for a gas field gathering and transportation network, characterized in that: include: Calculating the pipeline friction and admittance matrix of the gas field gathering and transportation network according to the basic parameters of the gas field gathering and transportation network, the set initial value of the pipeline volume flow of the gas gathering and transportation network, the set initial value of the pipeline mass flow of the gas production and transportation network, the set initial value of the node temperature, and the physical properties of natural gas, and generating a coefficient matrix; Obtaining the node pressure and pipeline volume flow of the gas gathering network, and the node pressure and pipeline mass flow of the gas production network according to the pipeline friction resistance, admittance matrix and coefficient matrix; Performing thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network; Recalculating the natural gas physical properties according to the BWRS state equation and the node pressures of the gas gathering network, the node pressures of the gas production network, and the node temperatures; Determining whether a difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether a difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether a difference between the node temperature and the initial value of the node temperature is less than a third threshold; If not, taking the pipeline volume flow as the pipeline volume flow initial value, taking the pipeline mass flow as the pipeline mass flow initial value, taking the node temperature as the node temperature initial value, and calculating the pipeline friction resistance and admittance matrix of the gas field gathering and transportation pipeline network according to the pipeline network basic parameters, the pipeline volume flow initial value, the pipeline mass flow initial value, the set node temperature initial value, and the recalculated natural gas physical property parameters, and generating a coefficient matrix; If so, the obtained node pressure, node pressure and node temperature of the gas gathering network are used as the coupling calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network, so as to guide the production of the gas field according to the coupling calculation results.

2. The method according to claim 1, characterized in that The basic parameters of the pipe network include the gas collection pipe diameter D 集 , Gas collecting pipe length L 集 , gas gathering network friction coefficient λ 集 , Gas pipe diameter D 采 , Gas pipe length L 采 , friction coefficient of gas pipeline network λ 采 , a pipe section node association matrix, the natural gas physical property parameters include the natural gas compressibility factor Z, the natural gas relative density Δ, and the natural gas average temperature T, wherein the natural gas average temperature T is calculated based on the set node temperature initial value and the pipe section node association matrix, the pipeline friction resistance of the gas field gathering and transportation network includes the gas gathering network friction resistance and the gas production network friction resistance, the admittance matrix includes the gas gathering network admittance matrix and the gas production network admittance matrix, and the coefficient matrix includes the gas gathering network coefficient matrix and the gas production network coefficient matrix.

3. The method according to claim 2, characterized in that The steps of calculating the pipeline friction and admittance matrix of the gas field gathering and transportation pipeline network according to the basic parameters of the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical properties, and generating a coefficient matrix include: According to the formula Calculate the friction resistance of the gas gathering network, where S 集 is the friction resistance of the gas collecting pipe network, and C is a constant; According to the formula Calculate the gas collection network admittance matrix, where G 集 is the admittance matrix of the gas gathering network, Q is the initial value of the pipeline volume flow rate, and α is the flow index; According to formula Y 集 =A×G 集 ×A T Calculate the gas collection network coefficient matrix, where Y 集 is the gas gathering network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose; According to the formula Calculate the friction resistance of the gas production network, where S 采 is the friction resistance of the gas production network, C1 is a constant, x is the gas mass content, ρ tp is the density of the mixture in the gas pipeline network, ρ g is the gas phase density, M represents the initial value of the pipeline mass flow rate, and p is the average absolute pressure of the pipeline; According to the formula Calculate the gas production network admittance matrix, where G 采 is the admittance matrix of the gas production network, M is the initial value of the pipeline mass flow rate, and α is the flow index; According to formula Y 采 =A×G 采 ×A T Calculate the gas production network coefficient matrix, where Y 采 is the gas production network coefficient matrix, A is the pipe segment node association matrix, and T is the matrix transpose.

4. The method according to claim 3, characterized in that The natural gas physical property parameters also include the compression factor of the natural gas at the compressor inlet, the natural gas density, the specific heat capacity of the gas at the inflow node, and the specific heat capacity of the gas at the compressor outlet.

5. The method according to claim 4, characterized in that The steps of obtaining the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction, admittance matrix and coefficient matrix include: According to the formula: AND 集 P 集 =-q, Calculate the node pressure and pipeline volume flow of the gas gathering network, where P 集 is the node pressure of the gas gathering network, q is the pipeline volume flow rate; According to the formula: Y 采 P 采 =-m, Calculate the node pressure and pipeline mass flow of the gas production network, where P 采 is the node pressure of the gas pipeline network, and m is the pipeline mass flow rate.

6. The method according to claim 5, characterized in that The gas field gathering and transportation network is thermally calculated based on the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network, further comprising: Constructing a correlation matrix between pipeline ends and nodes according to the pipeline segment node correlation matrix and the pipeline volume flow; According to the correlation matrix between the pipeline end and the node and the pipeline mass flow, the formula is used: Calculate the node temperature of the gas pipeline network, where T i,采 is the node temperature of the gas pipeline network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,采 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, T in,j,采 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,采 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,采 It is calculated based on the mass flow of the pipeline, and the outflow node flow m out,i,采 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,采 The inflow node temperature at the end of the pipeline includes a set known node temperature and an unknown node temperature, wherein the unknown node temperature is the set initial value of the node temperature; According to the correlation matrix between the pipe end and the node and the pipe volume flow, the formula is used: Calculate the node temperature of the gas pipeline network, where T i,集 is the node temperature of the gas gathering network, b in,ij is an element of a pipeline end node in the association matrix between the pipeline end and the node, m in,j,集 is the mass flow rate from the end of the pipe into the node, c in,j is the specific heat capacity of the gas flowing into the node, T in,j,集 is the temperature of the inflow node at the end of the pipe, m g,i is the node load flow, c g,i is the specific heat capacity of the node load gas, T g,i is the node load temperature, m out,i,集 is the outflow node flow, c out,i is the specific heat capacity of the gas flowing out of the node, where the mass flow rate m at the end of the pipe flowing into the node in,j,集 It is calculated based on the pipeline volume flow and the natural gas density in the pipeline, and the outflow node flow m out,i,集 It is calculated based on the mass flow rate of the inflow node at the end of the pipeline and the number of pipelines, T in,j,集 The inflow node temperature at the end of the pipeline includes the set known node temperature and unknown node temperature, wherein the unknown node temperature is the gas network node temperature T i,采 .

7. The method according to claim 6, characterized in that The step of constructing a pipeline end-node association matrix based on the pipeline segment node association matrix and the pipeline volume flow includes: Determining that the pipeline volume flow rate of the middle pipeline of the gas field gathering and transportation pipeline network is less than 0; If so, the element of the pipeline in the pipeline segment node association matrix is ​​used as the element of the pipeline in the pipeline end and node association matrix; If not, the inverse numbers of the elements of the pipeline in the pipeline segment node association matrix are used as the elements of the pipeline in the pipeline end and node association matrix; According to the structure of the gas field gathering and transportation pipeline network, the value of the element corresponding to the inlet pipeline in the pipeline-node association matrix is ​​set to 1.

8. A single-phase flow-multiphase flow coupling calculation device for a gas field gathering and transportation network, characterized in that: include: Iterative parameter calculation unit, pipeline flow calculation unit, node temperature calculation unit, natural gas physical property parameter calculation unit, result verification unit; The iterative parameter calculation unit is used to calculate the pipeline friction resistance and admittance matrix of the gas field gathering and transportation pipeline network according to the basic parameters of the gas field gathering and transportation pipeline network, the set initial value of the pipeline volume flow of the gas gathering pipeline network, the set initial value of the pipeline mass flow of the gas production pipeline network, the set initial value of the node temperature, and the natural gas physical property parameters, and generate a coefficient matrix; The pipeline flow calculation unit is used to obtain the node pressure and pipeline volume flow of the gas gathering network, the node pressure and pipeline mass flow of the gas production network according to the pipeline friction resistance, admittance matrix and coefficient matrix; The node temperature calculation unit is used to perform thermal calculation on the gas field gathering and transportation network according to the pipeline volume flow of the gas gathering network and the pipeline mass flow of the gas production network to obtain the node temperature of the gas field gathering and transportation network; The natural gas physical property parameter calculation unit is used to recalculate the natural gas physical property parameters according to the BWRS state equation and the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature; The result verification unit is used to determine whether the difference between the pipeline volume flow and the initial value of the pipeline volume flow is less than a first threshold, whether the difference between the pipeline mass flow and the initial value of the pipeline mass flow is less than a second threshold, and whether the difference between the node temperature and the initial value of the node temperature is less than a second threshold; if not, the pipeline volume flow is used as the initial value of the pipeline volume flow, the pipeline mass flow is used as the initial value of the pipeline mass flow, and the node temperature is used as the initial value of the node temperature, and the pipeline friction and admittance matrix of the gas field gathering and transportation network is calculated according to the basic parameters of the pipeline network, the initial value of the pipeline volume flow, the initial value of the pipeline mass flow, the set initial value of the node temperature, and the recalculated natural gas physical properties, and a coefficient matrix is ​​generated. If so, the node pressure of the gas gathering network, the node pressure of the gas production network, and the node temperature are used as the coupled calculation results of the unidirectional flow-multiphase flow of the gas field gathering and transportation network, so as to guide the exploitation of the gas field according to the coupled calculation results.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images