Natural gas pipeline network gas source tracking calculation method and system

Abstract

The invention discloses a gas source tracking calculation method and system for a natural gas pipeline network, which greatly reduces the number of nodes in a network by reforming and optimizing a constructed original pipeline network topology network structure to obtain a new pipeline network topology network structure, so that the matrix scale is remarkably reduced, the calculation complexity is effectively reduced, and meanwhile, the calculation efficiency is improved. And judging whether each row of data in the inverse matrix meets the condition that the total flow of the nodes is equal to the sum of the flow which flows to other nodes from the nodes and is downloaded by other nodes and the sum of the downloading quantity of the nodes, so that input errors and calculation errors are prevented, the accuracy of subsequent calculation results is improved, and the gas source proportion in the natural gas pipeline network is accurately tracked.

Description

Technical Field

[0001] The present invention belongs to the technical field of natural gas pipeline network gas source analysis, and relates to a natural gas pipeline network gas source tracking calculation method and system. Background Art

[0002] Tracking the gas source proportion in the natural gas pipeline network involves obtaining the gas source proportion distribution of user nodes and calculating the gas source proportion distribution at any node and pipeline position inside the pipeline network. When there is a gas quality problem at a certain user node, with the help of gas source proportion tracking, the gas source that may have problems can be quickly located. Similarly, if the quality of a certain gas source changes, using the gas source proportion tracking technology, it is also possible to intuitively analyze which nodes and pipelines may be affected by this gas source. In addition, since the wellhead prices of different gas sources are different, therefore, the benchmark wellhead price of natural gas corresponding to each user must be based on the results of gas source tracking for accounting.

[0003] In view of the fact that natural gas pipeline networks usually have complex characteristics of multiple gas sources and multiple users intertwined, and there are diverse gas supply paths and numerous transportation restrictions, therefore, in the daily production scheduling process, it is quite difficult to accurately define the supply and demand relationship between each gas source and its corresponding specific users. At the same time, in 2019, the national pipeline network was established, and the parallel operation of the main pipelines has further increased the complexity of the pipeline network structure.

[0004] One of the technical difficulties in achieving rapid calculation of gas source proportion tracking in natural gas pipeline networks lies in simplifying the complex network structure of the pipeline network. The actual natural gas pipeline network is huge in scale, and may contain tens of thousands of nodes, hundreds of gas sources, and the total mileage can reach tens of thousands of kilometers. In this case, relying on manual means to carry out gas source proportion tracking is almost an impossible task. If common hydraulic simulation software or simulation models are used for simulation, it has extremely high requirements for computer performance, and faces challenges such as low operation efficiency and non-convergent results. At the same time, in some simulation models, the way of directly assuming the gas source as the pipeline boundary is unreasonable. In addition, when using traditional complex network theory for calculation, the operation complexity increases rapidly with the increase of the order of the matrix. Therefore, on the premise of ensuring the accuracy of the calculation results, further simplifying the calculation method has become a major problem in the natural gas pipeline network gas source proportion tracking technology. Summary of the Invention

[0005] The purpose of the present invention is to solve the problems in the prior art that when tracking and calculating the gas source of the network management, the calculation is complex, the calculation accuracy is low, and it is difficult to accurately track the gas source in the natural gas pipeline network, and to provide a natural gas pipeline network gas source tracking calculation method and system.

[0006] To achieve the above object, the present invention adopts the following technical solutions: A method for calculating the gas source tracking of a natural gas pipeline network, comprising the following steps: Obtain network management parameters and construct an original pipeline network topology network structure based on the network management parameters; Reorganize and optimize the constructed original pipeline network topology network structure to obtain a new pipeline network topology network structure; Extract the distribution matrix of the new pipeline network topology network structure, calculate the inverse matrix of the distribution matrix, and determine whether each row of data in the inverse matrix satisfies: the total node flow is equal to the sum of the flows flowing from the node to other nodes and downloaded by other nodes plus the sum of the node's own downloads; If so, calculate the gas source proportion of each user based on the inverse matrix to obtain the gas source tracking result; If not, re-obtain the network management parameters for calculation.

[0007] Preferably, the network management parameters include the upload amount of the gas source station, the download amount of the user station, and the flow and flow direction in each pipeline in the pipeline network.

[0008] Preferably, the reorganizing and optimizing the constructed original pipeline network topology network structure to obtain a new pipeline network topology network structure includes: Obtain the original pipeline network topology network structure according to the flow direction of natural gas in the pipeline ; Among them, contains from to of original nodes; is the set of original network edges; Based on the original pipeline network topology network structure, traverse all nodes. When a certain node has only a single upstream pipe segment, combine the node, the upstream node, and the pipe segment between them as a whole into a super node, and obtain a new pipeline network topology network structure based on the integrated nodes.

[0009] Preferably, when calculating the inverse matrix of the distribution matrix, it further includes: Determine whether the distribution matrix of the new pipeline network topology network structure complies with the law of conservation of flow; If so, calculate the inverse matrix of the distribution matrix; If not, re-obtain the network management parameters for calculation.

[0010] Preferably, the extracting the distribution matrix of the new pipeline network topology network structure includes: The distribution matrix of the new pipeline network topology network structure is , and the element of the matrix is calculated by the following formula:

[0011] Wherein: is the Dirac function, that is, when then ; when then ; is the element value of the adjacency matrix of the directed network ; is the flow corresponding to the directed edge ; is the total flow of the node ;

[0012] 6. A method for calculating the gas source tracking of a natural gas pipeline network according to claim 5, characterized in that, calculating the inverse matrix of the distribution matrix, and determining whether each row of data in the inverse matrix satisfies: the total flow of the node is equal to the sum of the flow flowing from the node to other nodes and downloaded by other nodes and the sum of the self-download amount of the node, including: Calculating the inverse matrix of the calculation matrix , and the elements in the inverse matrix are denoted as ; Checking whether each row of data in the inverse matrix satisfies the following formula:

[0013] Wherein: is the total flow of the node , and the total flow is equal to the sum of the flows of all adjacent downstream pipe segments of the node and the download amount of the node ; is the download amount when the node is a user station, represents the flow flowing from the node to the node and downloaded by the node ;

[0014] Preferably, calculating the gas source proportion of each user based on the inverse matrix to obtain the gas source tracking result, including: Given the gas source station The proportion of the gas supply volume to the corresponding node for the user is: For:

[0015] Wherein, is the upload amount when the node is a gas source station, is the total flow of the node . The total flow is equal to the node The sum of the flows of all adjacent downstream pipe segments and the node download volume.

[0016] A natural gas pipeline network gas source tracking calculation system, comprising: A parameter acquisition module, configured to acquire network management parameters and construct an original pipeline network topology network structure based on the network management parameters; A new pipeline network topology network structure acquisition module, configured to reorganize and optimize the constructed original pipeline network topology network structure to obtain a new pipeline network topology network structure; A gas source ratio acquisition module, configured to extract the distribution matrix of the new pipeline network topology network structure, calculate the inverse matrix of the distribution matrix, and determine whether each row of data in the inverse matrix satisfies: the total node flow is equal to the sum of the flows flowing from the node to other nodes and downloaded by other nodes and the sum of the node's own download volume. If so, calculate the gas source proportion of each user based on the inverse matrix to obtain the gas source tracking result. If not, re-acquire the network management parameters for calculation.

[0017] A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the method are implemented.

[0018] A computer-readable storage medium storing a computer program, characterized in that when the computer program is executed by a processor, the steps of the method are implemented.

[0019] Compared with the prior art, the present invention has the following beneficial effects: The present invention discloses a natural gas pipeline network gas source tracking calculation method. By reorganizing and optimizing the constructed original pipeline network topology network structure, a new pipeline network topology network structure is obtained, significantly reducing the number of nodes in the network, thereby significantly reducing the matrix scale, effectively reducing the calculation complexity. At the same time, it is judged whether each row of data in the inverse matrix satisfies that the total node flow is equal to the sum of the flows flowing from the node to other nodes and downloaded by other nodes and the sum of the node's own download volume, preventing input errors and calculation errors, improving the accuracy of subsequent calculation results, and accurately tracking the gas source proportion in the natural gas pipeline network. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following will briefly introduce the drawings required in the embodiments. It should be understood that the following drawings only show some embodiments of the present invention, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts.

[0021] Figure 1 This is the application flow chart of the rapid tracking calculation method for the gas source proportion in the natural gas pipeline network described in the present invention; Figure 2 This is the schematic diagram of the network structure of the original pipeline network in the embodiment of the present invention; Figure 3 This is the schematic diagram of the network structure after renormalization improvement in the embodiment of the present invention. Detailed implementation manners

[0022] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Components of the embodiments of the present invention described and illustrated herein usually can be arranged and designed in various different configurations.

[0023] Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed present invention, but merely represents selected embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0024] It should be noted that: similar reference numerals and letters denote similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings.

[0025] In the description of the embodiments of the present invention, it should be noted that if terms such as "upper", "lower", "horizontal", "inner", etc. are used to indicate the orientation or positional relationship, it is based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship when the product of the present invention is usually placed. It is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation to the present invention. In addition, terms such as "first", "second", etc. are only used for distinguishing descriptions, and cannot be understood as indicating or implying relative importance.

[0026] In addition, if the term "horizontal" appears, it does not mean that the component is required to be absolutely horizontal, but it can be slightly inclined. For example, "horizontal" only means that its direction is more horizontal relative to "vertical", and does not mean that the structure must be completely horizontal, but it can be slightly inclined.

[0027] In the description of the embodiments of the present invention, it should also be noted that unless otherwise clearly specified and limited, if the terms "set", "installed", "connected", "connected" appear, they should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected, or indirectly connected through an intermediate medium, and it can be the communication inside two components. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

[0028] The following further describes the present invention in detail with reference to the drawings: See Figure 1 , the present invention discloses a method and system for tracking and calculating the gas source of a natural gas pipeline network. By using the means of network renormalization to reduce the number of nodes in the network, while ensuring the accuracy of the calculation results, the calculation complexity is greatly reduced, and the problem of difficult gas source tracking in large and complex natural gas pipeline networks is successfully solved.

[0029] Step 1: Obtain pipeline network parameters such as the upload volume of the gas source station, the download volume of the user station, and the flow rate and flow direction in each pipeline in the pipeline network.

[0030] Specifically, it includes the following steps: The pipeline network parameters such as the upload volume of the gas source station, the download volume of the user station, and the flow rate and flow direction in each pipeline in the pipeline network can be directly obtained from the pipeline network operator. The pipeline network operator plans the pipeline transportation plan according to the daily specified requirements of the natural gas shipper, so as to determine the specific transportation volume of each pipeline section.

[0031] Step 2: Obtain the directed network topological network structure of the original pipeline network.

[0032] Specifically, it includes the following steps: Regard the gas source station and the user station as nodes, and regard each pipeline as an edge. According to the flow direction of the natural gas in the pipeline, obtain a directed complex network .

[0033] Among them, contains from to of original nodes (stations), is the set of original edges of the network.

[0034] Step 3: Obtain the new directed network after renormalization.

[0035] Specifically, it includes the following steps: For the original network structure, network renormalization is performed based on the principle of full mixing in the station yard and the principle of fair distribution. Traverse all nodes. When a certain node has only a single upstream pipe segment, the node, together with the upstream node and the pipe segment between them, is merged as a whole into a super node. The download volume of the super node is the sum of the download volumes of all the merged original nodes. The gas source station yard cannot be processed for renormalization.

[0036] Step 4: Calculate the distribution matrix of the new network .

[0037] Specifically, it includes the following steps: Calculate each element of the matrix through the following formula: :

[0038] In the formula: is the Dirac function, that is, when , ; when , ; is the element value of the adjacency matrix of the directed network . That is, when there is a directed edge from node to node , ; when there is no directed edge from node to node , . In the formula: is the flow corresponding to the directed edge . In the formula: is the total flow of node . The total flow is equal to the sum of the flows of all adjacent downstream pipe segments of node and the download volume of node . Check the sum of each column of data in the matrix to determine whether it satisfies the law of conservation of flow.

[0039] Step 5: Calculate the inverse matrix of the matrix .

[0040] Specifically, it includes the following steps: You can choose to use the inv() function in the MATLAB software to calculate the inverse matrix of the matrix . Check whether each row of data in the matrix satisfies the following formula:

[0041] In the formula: is the total flow of node The total flow is equal to the sum of the flows of all adjacent downstream pipe segments of node and the download volume of node ; is the download volume when node is a user station.

[0042] Check whether the sum of the data in each column of matrix meets the requirement of observing the law of conservation of flow and whether each row of data in the check matrix in step 5 satisfies formula (2). This process significantly prevents input errors and calculation errors from occurring.

[0043] Step 6: Calculate the gas source proportion distribution of user nodes.

[0044] Specifically, it includes the following steps: Calculate the gas supply proportion from a given gas source station to a given node for the corresponding user through the following formula:

[0045] In the formula: is the upload volume when node is a gas source station; is the total flow of node ; The total flow is equal to the sum of the flows of all adjacent downstream pipe segments of node and the download volume of node ; is an element of matrix , which is calculated in step 5.

[0046] The present invention also discloses a specific embodiment. Refer to Figures 2 to 3 : As Figure 2 shown, in this embodiment, a set of natural gas pipeline network models is pre-constructed, and the flow configuration of each pipe segment and the daily specified upload and download volumes of each station are listed in detail.

[0047] Among them, the gas source station is represented by a red hollow ring, and its upload volume is marked in the upper left corner of the ring and displayed in red numbers. The ordinary station is presented as a dark blue solid circle, and its download volume is marked in the lower right corner of the circle and displayed in blue numbers. The side arrow indicates the direction of the gas flow, and the corresponding pipeline flow is marked in black numbers on the arrow. The flow unit is .

[0048] As Figure 3As shown, after traversing all the original nodes (excluding the gas source node), the following renormalization strategy is determined: Integrate nodes 11, 12, 13, 14, 15, and 16 into a super node, named S11; Merge nodes 3, 4, and 5 into a super node, labeled S3; Merge nodes 7 and 8 into super node S7; Merge nodes 17 and 19 into super node S17.

[0049] Furthermore, through Figure 3 and the formula , calculate each element of the matrix , as follows:

[0050] After checking the sum of each column of the matrix , it conforms to the law of conservation of flow, that is, the sum of the column data of the user nodes must be equal to zero.

[0051] Furthermore, calculate the inverse matrix of the matrix , , as follows:

[0052] After checking, each row of data in the matrix satisfies formula (2).

[0053] Furthermore, through the formula , calculate the proportion of the gas supply volume from a given gas source station to a given node for the corresponding user , see Table 1: Table 1 Proportion of gas supply volume for users

[0054] The method disclosed in this embodiment realizes significant reduction of the matrix scale and effectively reduces the computational complexity through the optimization based on network renormalization in step 3. Compared with the prior art, the present invention has the following beneficial effects: The method for quickly tracking the proportion of gas sources in a natural gas pipeline network improved based on network renormalization in the present invention can effectively reduce the number of nodes in the network, and while ensuring the accuracy of the calculation results, greatly reduce the computational complexity.

[0055] This embodiment also discloses a natural gas pipeline network gas source tracking calculation system, including: A parameter acquisition module, configured to acquire network management parameters and construct an original pipe network topology network structure based on the network management parameters; A new pipe network topology network structure acquisition module, configured to reorganize and optimize the constructed original pipe network topology network structure to obtain a new pipe network topology network structure; A gas source ratio acquisition module, configured to extract the distribution matrix of the new pipe network topology network structure, calculate the inverse matrix of the distribution matrix, and determine whether each row of data in the inverse matrix satisfies: the total node flow is equal to the sum of the flows flowing from the node to other nodes and downloaded by other nodes plus the self-download volume of the node. If so, calculate the gas source ratio of each user based on the inverse matrix to obtain the gas source tracking result. If not, re-acquire the network management parameters for calculation.

[0056] The schematic diagram of the terminal device provided by an embodiment of the present invention. The terminal device in this embodiment includes: a processor, a memory, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps in the above-mentioned various method embodiments are implemented. Alternatively, when the processor executes the computer program, the functions of each module / unit in the above-mentioned various device embodiments are implemented.

[0057] The computer program can be divided into one or more modules / units, and the one or more modules / units are stored in the memory and executed by the processor to complete the present invention.

[0058] The terminal device can be a computing device such as a desktop computer, a notebook, a palm computer, and a cloud server. The terminal device may include, but is not limited to, a processor and a memory.

[0059] The processor may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.

[0060] The memory can be used to store the computer program and / or module. The processor realizes various functions of the terminal device by running or executing the computer program and / or module stored in the memory, and calling the data stored in the memory.

[0061] If the modules / units integrated in the terminal device are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, to implement all or part of the processes in the above-described embodiment methods of the present invention, it can also be completed by instructing relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above-described various method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file, or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0062] The above are only the preferred embodiments of the present invention and are not used to limit the present invention. For those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A natural gas pipeline network gas source tracking calculation method, characterized in that: The following steps are involved: Obtain network management parameters, and build the original network topology structure based on the network management parameters; Reorganize and optimize the constructed original pipe network topology to obtain a new pipe network topology; Extract the allocation matrix of the new pipe network topology, calculate the inverse matrix of the allocation matrix, and determine whether each row of data in the inverse matrix satisfies: the total flow of the node is equal to the sum of the flow from the node to other nodes and downloaded by other nodes plus the sum of the node's own downloads; If yes, the gas source proportion of each user is calculated based on the inverse matrix to obtain the gas source tracking result; If not, re-obtain the network management parameters for calculation.

2. A natural gas pipeline network gas source tracing calculation method according to claim 1, characterized in that: The network management parameters include the upload volume of the gas source station, the download volume of the user station, and the flow rate and flow direction in each pipeline in the pipeline network.

3. A natural gas pipeline network gas source tracing calculation method according to claim 1, characterized in that: The reorganization and optimization of the constructed original pipe network topology structure to obtain a new pipe network topology structure includes: Obtain the original pipeline network topology structure based on the flow direction of natural gas in the pipeline ; in, Contains from arrive of Original nodes; is the set of original edges of the network; Based on the original pipe network topology, all nodes are traversed. When a node has only a single upstream pipe segment, the node, the upstream node and the pipe segment between them are merged into a super node as a whole. Based on the integrated node, the new pipe network topology is obtained.

4. A natural gas pipeline network gas source tracing calculation method according to claim 1, characterized in that: The calculation of the inverse matrix of the allocation matrix also includes: Determine whether the distribution matrix of the new pipe network topology structure complies with the law of flow conservation; If so, calculate the inverse matrix of the allocation matrix; If not, re-obtain the network management parameters for calculation.

5. A natural gas pipeline network gas source tracing calculation method according to claim 1, characterized in that: The method of extracting the allocation matrix of the new pipe network topology structure includes: The allocation matrix of the new pipe network topology is: ,matrix Elements Calculated by the following formula: Where: Dirac Function, that is, when hour, ;when hour, ; For a directed network The element values ​​of the adjacency matrix; For directed edges The corresponding flow rate; For Node The total flow rate.

6. A natural gas pipeline network gas source tracing calculation method according to claim 5, characterized in that: Calculate the inverse matrix of the allocation matrix and determine whether each row of data in the inverse matrix satisfies the following conditions: the total traffic of a node is equal to the sum of the traffic from the node to other nodes and downloaded by other nodes plus the sum of the traffic downloaded by the node itself, including: Calculate the calculation matrix The inverse matrix , the inverse matrix The elements in are denoted as ; Check the inverse matrix Does each row of data satisfy the following formula: Where: For Node The total flow is equal to the total flow of the node The flow rate of all adjacent downstream pipe sections and the node The sum of the downloads of For Node The number of downloads when the user is on the site. Represents a slave node Flow to Node and is node Download traffic.

7. A natural gas pipeline network gas source tracing calculation method according to claim 1, characterized in that: The gas source proportion of each user is calculated based on the inverse matrix to obtain the gas source tracking result, including: Given gas source station To a given node Gas supply ratio for corresponding users for: In the formula, For Node is the loading amount at the gas source station, For Node The total flow of the node The flow rate of all adjacent downstream pipe sections and the node The sum of downloads.

8. A natural gas pipeline network gas source tracking calculation system, characterized in that: include: A parameter acquisition module is used to acquire network management parameters and build the original network topology structure based on the network management parameters; A new pipe network topology network structure acquisition module is used to reorganize and optimize the constructed original pipe network topology network structure to obtain a new pipe network topology network structure; The gas source ratio acquisition module is used to extract the distribution matrix of the new pipe network topology network structure, calculate the inverse matrix of the distribution matrix, and determine whether each row of data in the inverse matrix satisfies: the total flow of the node is equal to the sum of the flow from the node to other nodes and downloaded by other nodes plus the sum of the node's own downloads. If so, the gas source ratio of each user is calculated based on the inverse matrix to obtain the gas source tracking results. If not, the network management parameters are re-obtained for calculation.

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

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

Images