Oil and gas field electric power system operation state prediction method and device and storage medium

Abstract

The invention discloses an oil and gas field power system operation state prediction method and device and a storage medium. The method comprises the steps of obtaining a topological structure, branch parameters and real-time measurement data of a target oil and gas field power system; based on the topological structure, the branch parameters and the real-time measurement data, constructing an induction motor dynamic equivalent model corresponding to the target oil and gas field power system; on the basis of the induction motor dynamic equivalent model, the topological structure, the branch parameters and the real-time measurement data, measurement functions of different operation state parameters of the target oil and gas field power system are constructed, and Jacobian matrixes corresponding to the different measurement functions are determined; and determining an initial operation state parameter, carrying out iterative solution on the initial operation state parameter based on the measurement function and the Jacobian matrix, and finally taking the operation state obtained by the final solution meeting an iterative solution convergence condition as the predicted operation state of the target oil and gas field power system.

Description

Technical Field

[0001] This invention relates to the field of smart grid technology, and in particular to a method, apparatus and storage medium for predicting the operating status of an oil and gas field power system. Background Technology

[0002] Power system state estimation is one of the core functions of energy management systems and serves as the foundation for advanced applications such as power flow analysis, security early warning, and optimized scheduling.

[0003] Currently, in the estimation of power grid operating status, the load is usually modeled as a constant power load or a simple static load model. However, the core electrical equipment in oil and gas fields (such as pumping units and water injection pumps) are typical induction motor loads, which have higher requirements for voltage control. This leads to a fundamental deviation between the model on which the state estimation is based and the actual physical characteristics of the system, resulting in a large error between the estimation results and the actual state, and low reliability. Summary of the Invention

[0004] This invention provides a method, device, and storage medium for predicting the operating status of an oil and gas field power system, which mainly improves the accuracy of predicting the operating status of the oil and gas field power system.

[0005] According to a first aspect of the present invention, a method for predicting the operating status of an oil and gas field power system is provided, comprising: In response to the predicted operating status signal of the target oil and gas field power system, the topology, branch parameters, and real-time measurement data of the target oil and gas field power system are acquired. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power. Based on the topology, the branch parameters, and the real-time measurement data, a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field is constructed. Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, measurement functions for different operating state parameters of the target oil and gas field power system are constructed, and the Jacobian matrix corresponding to different measurement functions is determined. The initial operating state parameters are determined, and the initial operating state parameters are iteratively solved based on the measurement function and the Jacobian matrix. Finally, the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the target oil and gas field power system.

[0006] Optionally, based on the topology, the branch parameters, and the real-time measurement data, a dynamic equivalent model of the induction motors corresponding to the power system of the target oil and gas field is constructed, including: Based on the aforementioned topology, branch parameters, and real-time measurement data, the equivalent admittance function of the induction motor in the target oil and gas field's power system is constructed. ,in, , The stator impedance of the induction motor, For resistors, For imaginary numbers, For the stator reactance of the induction motor, For rotor impedance, The resistance of the rotor impedance The reactance component of the rotor impedance. For the reactance of the excitation branch, s is the magnetizing reactance, and s is the slip of the induction motor; The equivalent admittance function Decomposed into equivalent conductance and equivalent susceptance ,in, ; Based on the equivalent conductivity and the equivalent susceptance Determine the node-injected active power function of the induction motor connected to node i in the power system of the target oil and gas field. and node-injected reactive power function ,in, , , Let x be the voltage amplitude at node i, and let x be the operating state variable of the target oil and gas field power system. Based on the aforementioned topology, branch parameters, and real-time measurement data, the electromagnetic torque function of the induction motor is constructed. and load torque function ,in, , , The mechanical power output by the induction motor is s, where s is the slip rate. For synchronous angular velocity, For constant loss torque, This is the load factor; When the electromagnetic torque of the induction motor and the load torque are equal, the electromagnetic torque function is... and the load torque function Performing a function transformation, the slip ratio constructor is obtained as follows: ,in, The output mechanical power of the induction motor under slip s and voltage V; The equivalent admittance function The node injects active power function The node injects reactive power function The slip rate constructor is used as the dynamic equivalent model of the induction motor.

[0007] Optionally, the step of constructing measurement functions for different operating state parameters of the target oil and gas field power system based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data includes: Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, the node voltage amplitude measurement function is determined. ,in, x is the running state vector. Let be the voltage at node i; Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, the node-injected active power measurement function is determined. and node-injected reactive power measurement function ,in, , x is the running state vector, The voltage amplitude at node i, The voltage amplitude at node j, For the real part of the branch ij admittance, For the imaginary part of the admittance of branch ij, Let be the phase angle difference between node i and node j; Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, a motor slip measurement function is constructed. ,in, x is the running state vector. The electromagnetic torque of the induction motor under slip s and voltage V, Let be the load torque of the induction motor at slip s.

[0008] Optionally, determining the Jacobian matrix corresponding to different measurement functions includes: Determine the node voltage amplitude measurement function The partial derivative matrix of the node voltage amplitude measurement function The partial derivative matrix is ​​used as the node voltage amplitude measurement function. The Jacobian matrix; Determine the active power measurement function injected at the node. The partial derivative matrix of the function is used to inject the active power measurement function into the node. The partial derivative matrix is ​​used as the active power measurement function injected at the node. The Jacobian matrix; Determine the reactive power measurement function injected at the node. The partial derivative matrix of the node is used to inject reactive power measurement function. The partial derivative matrix is ​​used as the reactive power measurement function injected at the node. The Jacobian matrix; Determine the motor slip measurement function The partial derivative matrix of the motor slip measurement function is used to... The partial derivative matrix is ​​used as the slip measurement function of the motor. The Jacobian matrix.

[0009] Optionally, the initial operating state parameters include the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip. The process of iteratively solving the initial operating state parameters based on the measurement function and the Jacobian matrix, and finally using the operating state obtained from the final solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system, includes: Based on the initial values ​​of the node voltage amplitude, the node voltage phase angle, and the induction motor slip, the initial value of the operating state vector of the target oil and gas field power system is determined, and this initial value of the operating state vector is used as the current operating state vector at the current iteration number k. ; The current running state vector Substituting the values ​​into the measurement function, we can perform measurement estimation to obtain the measurement estimate. The current running state vector Substituting the values ​​into the Jacobian matrix, we can perform Jacobian matrix estimation to obtain the Jacobian matrix estimator. And based on the actual quantity measurement of the power system of the target oil and gas field. and the measurement estimate Determine the measurement residual ,in, ; Based on the measurement residual The Jacobian matrix estimator Determine the change in the current running state corresponding to the current running state vector. ,in, ,in, The weight matrix is ​​set according to actual needs; Based on the change in the current operating status The predicted operating status of the power system of the target oil and gas field is determined.

[0010] Optionally, the change based on the current operating state Determining the predicted operating status of the target oil and gas field's power system includes: Determine the change in the current operating state. Is it less than a preset change threshold? If the current operating state changes If the change is less than the preset threshold, then based on the current running state vector... and the change in the current operating state Determine the predicted operating status of the power system of the target oil and gas field. ,in, ; If the current operating state changes If the change is greater than or equal to the preset change threshold, then the... As the new current running state vector under the new current iteration number k+1 Based on the measurement function and the Jacobian matrix, the new current running state vector is... The solution is iteratively solved, and the final solution that meets the convergence condition is taken as the predicted operating state of the target oil and gas field power system.

[0011] Optionally, after taking the final solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system, the method further includes: Each of the multiple predicted operating states is taken as a target predicted operating state. The target's predicted running state Substituting the values ​​into the measurement function for vector estimation yields the predicted operating state vector. And based on the predicted running state vector The predicted running state vector The corresponding actual running state vector Determine the target prediction running state Corresponding abnormal operation status evaluation value ,in, ,in, These are the weighting coefficients; Determine the abnormal operating status evaluation value If the value exceeds a preset evaluation threshold, then based on the actual operating state vector... and the predicted running state vector Determine the measurement residual ,in, ; Based on the measurement residual The target predicted operating state Jacobian matrix The residual covariance matrix is ​​determined, and the residual covariance matrix is ​​standardized to obtain the target prediction running state. Corresponding negative data evaluation values; The predicted operating state corresponding to the maximum adverse data evaluation value in each predicted operating state is taken as the abnormal predicted operating state, and the abnormal predicted operating states in the predicted operating states are removed to obtain the cleaned predicted operating state corresponding to the power system of the target oil and gas field.

[0012] According to a second aspect of the present invention, a device for predicting the operating status of an oil and gas field power system is provided, comprising: The acquisition unit is used to acquire the topology, branch parameters, and real-time measurement data of the target oil and gas field power system in response to the operation status prediction signal of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power. The model building unit is used to build a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field based on the topology, the branch parameters, and the real-time measurement data. The function construction unit is used to construct measurement functions for different operating state parameters of the target oil and gas field power system based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, and to determine the Jacobian matrix corresponding to different measurement functions. The solution unit is used to determine the initial operating state parameters, and iteratively solve the initial operating state parameters based on the measurement function and the Jacobian matrix. Finally, the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the target oil and gas field power system.

[0013] According to a third aspect of the present invention, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the above-described method for predicting the operating status of an oil and gas field power system.

[0014] According to a fourth aspect of the present invention, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the above-mentioned method for predicting the operating status of an oil and gas field power system.

[0015] The present invention provides a method, apparatus, and storage medium for predicting the operating state of an oil and gas field power system. Compared with the current method of estimating the operating state of the power grid by modeling the load as a constant power load or a simple static load model, the present invention constructs a dynamic equivalent model of the induction motor corresponding to the target oil and gas field power system and predicts the operating state of the power system based on the dynamic equivalent model of the induction motor. This ensures that the model on which the state estimation is based is consistent with the actual physical characteristics of the power system, improves the simulation capability of the system's dynamic behavior, and thus improves the accuracy of the power system operating state prediction. By determining the Jacobian matrix corresponding to different measurement functions, the relationship between the measurement functions and state variables can be described more accurately, which helps to converge to the true solution faster and more accurately in the iterative solution process. Using real-time measurement data for state prediction can promptly detect abnormal changes in the system and provide real-time data support for system stability analysis. Attached Figure Description

[0016] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 A flowchart of a method for predicting the operating status of an oil and gas field power system provided by an embodiment of the present invention is shown; Figure 2 A flowchart of another method for predicting the operating status of an oil and gas field power system provided by an embodiment of the present invention is shown; Figure 3 A schematic diagram of the structure of a predictive device for the operating status of an oil and gas field power system provided in an embodiment of the present invention is shown; Figure 4 A schematic diagram of another oil and gas field power system operation status prediction device provided in an embodiment of the present invention is shown; Figure 5 A schematic diagram of the physical structure of a computer device provided in an embodiment of the present invention is shown. Detailed Implementation

[0017] The present invention will be described in detail below with reference to the accompanying drawings and embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in the present application can be combined with each other.

[0018] Currently, in the process of estimating the operating status of the power grid, the load is usually modeled as a constant power load or a simple static load model. This results in a fundamental deviation between the model on which the state estimation is based and the actual physical characteristics of the system, leading to a large error between the estimation results and the actual state.

[0019] To address the aforementioned problems, embodiments of the present invention provide a method for predicting the operating status of an oil and gas field power system, such as... Figure 1 As shown, the method includes: 101. In response to the operation status prediction signal of the target oil and gas field power system, acquire the topology, branch parameters, and real-time measurement data of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix. The real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power.

[0020] In this embodiment of the invention, the topology information of the power grid is extracted from the database or design documents of the oil and gas field power system, including the connection relationships between nodes and line parameters. Branch parameters mainly include branch admittance and node admittance matrices, obtained by reading power system design data or through real-time calculation. Branch admittance reflects the line's ability to conduct current, while the node admittance matrix describes the electrical relationships between nodes. Using a SCADA system (Supervisory Control and Data Acquisition) or other real-time monitoring equipment, real-time measurement data such as node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power are acquired. This data forms the basis for state estimation. This embodiment of the invention utilizes real-time measurement data for state prediction, enabling timely detection of abnormal changes in the system and providing real-time data support for system stability analysis.

[0021] 102. Based on the topology, branch parameters, and real-time measurement data, construct a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field.

[0022] In this embodiment of the invention, a state variable vector x is defined, including node voltage magnitude, node voltage phase angle, and induction motor slip. Based on the topology, branch parameters, and real-time measurement data, an equivalent admittance function, node injected active power function, node injected reactive power function, and slip constructor for the induction motor are constructed. These constructed functions are then integrated to form a complete dynamic equivalent model of the induction motor. This model comprehensively reflects the dynamic characteristics of the induction motor in the oil and gas field power system, providing an accurate foundation for subsequent state estimation. This embodiment of the invention, by constructing a dynamic equivalent model of the induction motor corresponding to the target oil and gas field power system and predicting the power system's operating state based on this model, ensures that the model used for state estimation is consistent with the actual physical characteristics of the power system, improving the simulation capability of the system's dynamic behavior and thus enhancing the accuracy of power system operating state prediction.

[0023] 103. Based on the dynamic equivalent model of the induction motor, topology, branch parameters, and real-time measurement data, construct measurement functions for different operating state parameters of the target oil and gas field power system, and determine the Jacobian matrix corresponding to different measurement functions.

[0024] In this embodiment of the invention, key parameters (such as equivalent admittance, slip constructor, etc.) in the dynamic equivalent model of the induction motor are integrated with topology (node ​​connection relationships, line parameters, etc.), branch parameters (branch admittance, node admittance matrix, etc.), and real-time measurement data (node ​​voltage amplitude, line active / reactive power, etc.) to form a complete dataset. The integrated data is cleaned to remove outliers and missing values, and then normalized or standardized to improve data quality and the accuracy of subsequent calculations. Based on the integrated and preprocessed data, measurement functions for different operating state parameters of the target oil and gas field power system are constructed. These measurement functions describe the relationship between measured values ​​and state variables. Specifically, the node voltage amplitude measurement function directly uses the node voltage amplitude from the real-time measurement data as the measured value, constructing a direct mapping relationship between the node voltage amplitude and the node voltage amplitude state variable. This function reflects the real-time changes in node voltage. Node-injected power measurement functions: Active power: Based on the node-injected active power function in the dynamic equivalent model of the induction motor, and combined with topology and branch parameters, a measurement function for node-injected active power is constructed. This function considers the influence of the motor's dynamic characteristics on the node-injected active power. Reactive power: Similarly, a measurement function for node-injected reactive power is constructed to reflect the real-time changes in node-injected reactive power. Line power measurement function: Based on the topology and branch parameters, and combined with real-time measurement data of line active / reactive power, a measurement function for line power is constructed. This function describes the relationship between the active / reactive power transmitted on the line and the state variables. Motor slip measurement function: Using the slip constructor in the dynamic equivalent model of the induction motor, a measurement function for motor slip is constructed. This function reflects the real-time changes in motor slip and is an important indicator for evaluating the motor's operating status. The Jacobian matrix is ​​the partial derivative matrix of the measurement function with respect to the state variables, used to describe the linear relationship between the measurement function and the state variables.

[0025] 104. Determine the initial operating state parameters, and iteratively solve the initial operating state parameters based on the measurement function and Jacobian matrix. Finally, take the operating state obtained by the last solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system.

[0026] In this embodiment of the invention, initial operating state parameters such as the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip of the target oil and gas field power system are randomly set according to actual needs. For example, operating state data of the target oil and gas field power system over a past period, including node voltage amplitude, node voltage phase angle, and induction motor slip, are used as the initial parameter setting basis. Further, the initial operating state parameters are substituted into the measurement function to calculate the estimated measurement values ​​for each measurement point. The calculated estimated measurement values ​​are compared with the actual real-time measurement data to calculate the difference between the two, i.e., the residual. Using the Jacobian matrix and the residual, the correction amount for the current state parameters is determined. The correction amount is applied to the current state parameters to update the initial operating state parameters or the state parameters obtained in the previous iteration, resulting in new state parameters. The calculated residual is then checked to see if it meets the pre-set convergence conditions. The convergence conditions can be set according to actual needs, such as the absolute value of the residual being less than a small threshold, or the rate of change of the residual being less than a given value. If the convergence condition is met, it indicates that the current state parameters are sufficiently close to the actual operating state of the system, and the iterative solution process ends. If the convergence condition is not met, the process returns to continue with steps such as measurement estimate calculation, residual calculation, correction calculation, and parameter update until the convergence condition is met. Once the iterative solution process meets the convergence condition, the state parameters obtained from the final solution are used as the predicted operating state of the target oil and gas field power system. These state parameters include node voltage amplitude, node voltage phase angle, and slip of induction motors, which can comprehensively reflect the operating state of the system at the current moment or at a future moment. By determining the Jacobian matrix corresponding to different measurement functions, this embodiment of the invention can more accurately describe the relationship between measurement functions and state variables, which helps to converge to the true solution faster and more accurately during the iterative solution process.

[0027] The present invention provides a method for predicting the operating state of an oil and gas field power system. Compared with the current method of estimating the operating state of the power grid by modeling the load as a constant power load or a simple static load model, this invention constructs a dynamic equivalent model of the induction motor corresponding to the target oil and gas field power system and predicts the operating state of the power system based on the dynamic equivalent model of the induction motor. This ensures that the model on which the state estimation is based is consistent with the actual physical characteristics of the power system, improves the simulation capability of the system's dynamic behavior, and thus improves the accuracy of the power system operating state prediction. By determining the Jacobian matrix corresponding to different measurement functions, the relationship between the measurement functions and state variables can be described more accurately, which helps to converge to the true solution faster and more accurately in the iterative solution process. Using real-time measurement data for state prediction can promptly detect abnormal changes in the system and provide real-time data support for system stability analysis.

[0028] Furthermore, to better illustrate the above process of predicting the operating status of oil and gas field power systems, as a refinement and extension of the above embodiments, this invention provides another method for predicting the operating status of oil and gas field power systems, such as... Figure 2 As shown, the method includes: 201. In response to the operation status prediction signal of the target oil and gas field power system, acquire the topology, branch parameters, and real-time measurement data of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix. The real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power.

[0029] Specifically, data such as the topology and branch parameters of the power system in the oil and gas field are obtained from the database, and measurement data such as node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power are measured in real time through measuring devices.

[0030] 202. Based on the topology, branch parameters, and real-time measurement data, construct a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field.

[0031] For embodiments of the present invention, in order to accurately reflect the operating status of the oil and gas field power system, it is first necessary to construct a dynamic equivalent model of the induction motor corresponding to the oil and gas field power system. Based on this, step 202 specifically includes: constructing the equivalent admittance function of the induction motor in the target oil and gas field power system based on the topology, the branch parameters, and the real-time measurement data. ,in, , The stator impedance of the induction motor, For resistors, For imaginary numbers, For the stator reactance of the induction motor, For rotor impedance, The resistance of the rotor impedance The reactance component of the rotor impedance. For the reactance of the excitation branch, Let be the magnetizing reactance, and s be the slip of the induction motor; the equivalent admittance function Decomposed into equivalent conductance and equivalent susceptance ,in, Based on the equivalent conductivity and the equivalent susceptance Determine the node-injected active power function of the induction motor connected to node i in the power system of the target oil and gas field. and node-injected reactive power function ,in, , , Let x be the voltage amplitude at node i, and x be the operating state variable of the target oil and gas field power system; based on the topology, the branch parameters, and the real-time measurement data, construct the electromagnetic torque function of the induction motor. and load torque function ,in, , , The mechanical power output by the induction motor is s, where s is the slip rate. For synchronous angular velocity, For constant loss torque, The load factor is used; when the electromagnetic torque of the induction motor and the load torque are equal, the electromagnetic torque function is... and the load torque function Performing a function transformation, the slip ratio constructor is obtained as follows: ,in, The output mechanical power of the induction motor under slip s and voltage V; the equivalent admittance function The node injects active power function The node injects reactive power function The slip rate constructor is used as the dynamic equivalent model of the induction motor.

[0032] Specifically, firstly, based on the topology, branch parameters, and real-time measurement data of the target oil and gas field's power system, the relevant parameters of the induction motor are determined. These parameters include the stator resistance. Stator Reactor Rotor resistance Rotor reactance Reactance of the excitation branch Parameters such as slip (s) can be obtained through analysis of power system design documents, measurement of equipment parameters, and acquisition of real-time monitoring data. Based on these parameters, the equivalent admittance function of the induction motor is constructed according to the following formula. :

[0033] It should be noted that in determining the equivalent admittance function, it is essential to ensure that the units of all parameters are consistent and that calculations are performed according to the rules of complex number operations. For example, for addition, subtraction, multiplication, and division of complex numbers, the real and imaginary parts must be operated on separately to ensure the accuracy of the calculation results. Furthermore, mathematical computation tools (such as MATLAB software) are used to calculate the equivalent admittance function. The input is decomposed into the equivalent conductance of the real part using the software's internal complex number processing function. Equivalent susceptance of the imaginary part Furthermore, based on the equivalent conductance obtained from the decomposition... and equivalent susceptance Combined with the voltage amplitude V of node i in the power system of the target oil and gas field i This operating state variable determines the nodal injected active power function of the induction motor connected to node i. and node-injected reactive power function Furthermore, based on the topology, branch parameters, and real-time measurement data, the mechanical power output of the induction motor is determined. Constant loss torque Load factor and synchronous angular velocity Based on the parameters mentioned above, the load torque function is constructed according to the following formula. :

[0034] It should be noted that when constructing the load torque function, attention should be paid to the range of values ​​for the slip s. Furthermore, it is essential to ensure that the units of all parameters are consistent with the formula requirements to guarantee the accuracy of the function calculation. Further, when the electromagnetic torque of the induction motor and the load torque are equal, i.e. The electromagnetic torque function and load torque function A function transformation is performed, and the slip constructor is obtained based on the transformation relationship. Finally, the equivalent admittance function, nodal injection power function, and slip constructor constructed above are integrated to form a complete dynamic equivalent model of the induction motor.

[0035] 203. Based on the dynamic equivalent model of the induction motor, topology, branch parameters, and real-time measurement data, construct measurement functions for different operating state parameters of the target oil and gas field power system, and determine the Jacobian matrix corresponding to different measurement functions.

[0036] Among them, different operating state parameters include node voltage amplitude, node injected power, motor slip, etc.

[0037] In this embodiment of the invention, in order to determine the operating status of the oil and gas field power system, it is also necessary to construct a measurement function for the operating status parameters. Based on this, step 203 specifically includes: determining the node voltage amplitude measurement function based on the induction motor dynamic equivalent model, the topology, the branch parameters, and the real-time measurement data. ,in, x is the running state vector. The voltage at node i; based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, determine the active power measurement function injected into the node. and node-injected reactive power measurement function ,in, , x is the running state vector, The voltage amplitude at node i, The voltage amplitude at node j, For the real part of the branch ij admittance, For the imaginary part of the admittance of branch ij, The phase angle difference between node i and node j; based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, a motor slip measurement function is constructed. ,in, x is the running state vector. The electromagnetic torque of the induction motor under slip s and voltage V, Let be the load torque of the induction motor at slip s.

[0038] Specifically, the node voltage amplitude measurement function is constructed directly based on the physical meaning of the node voltage. Let the voltage at node i be... Then the node voltage amplitude measurement function The function value corresponding to node i is This means that the measurement function directly uses the voltage value of node i in the operating state vector as the measurement result to reflect the node voltage amplitude, an operating state parameter. For example, in a power system of an oil and gas field with multiple nodes, each node has a corresponding voltage amplitude measurement function. Through these functions, the voltage amplitude information of each node can be obtained in real time. Simultaneously, a node-injected active power measurement function is constructed. and node-injected reactive power measurement function Node-injected active power measurement function and node-injected reactive power measurement function The injected active power at node i is accurately calculated by comprehensively considering factors such as voltage magnitude, branch admittance, and phase angle difference between node i and its connected node j. Simultaneously, a measurement function for the injected active power at motor nodes is constructed. and reactive power measurement function As shown below:

[0039]

[0040] Where i and j represent the node numbers in the power system, V i Let V be the voltage magnitude at node i. j For the voltage magnitude at node j, G ij B is the real part of the branch admittance (branch conductance). ij Let be the imaginary part of the branch ij admittance (branch susceptance). Let be the phase angle difference between node i and node j. Simultaneously, construct the active power measurement function for injecting power into the motor nodes. and reactive power measurement function :

[0041]

[0042] Among them, V i Let be the voltage magnitude at node i. This is the equivalent conductance of the induction motor. Let be the equivalent susceptance of the induction motor, and s be the slip of the motor. Simultaneously, determine the active power measurement function for branch ij. and reactive power measurement function :

[0043]

[0044] Among them, V i and V j G represents the voltage magnitudes at nodes i and j at the two ends of branch ij, respectively. ij b is the real part of the branch ij admittance. ij It is the imaginary part of the branch ij admittance. This represents the phase angle difference between nodes i and j at both ends of branch ij. Simultaneously, the motor slip measurement function is determined. .

[0045] Furthermore, after determining the measurement function, it is also necessary to determine the Jacobian matrix. Based on this, step 203 specifically includes: determining the node voltage amplitude measurement function. The partial derivative matrix of the node voltage amplitude measurement function The partial derivative matrix is ​​used as the node voltage amplitude measurement function. The Jacobian matrix; determine the node-injected active power measurement function. The partial derivative matrix of the function is used to inject the active power measurement function into the node. The partial derivative matrix is ​​used as the active power measurement function injected at the node. The Jacobian matrix; determine the nodal injected reactive power measurement function. The partial derivative matrix of the node is used to inject reactive power measurement function. The partial derivative matrix is ​​used as the reactive power measurement function injected at the node. The Jacobian matrix; determine the slip measurement function of the motor. The partial derivative matrix of the motor slip measurement function is used to... The partial derivative matrix is ​​used as the slip measurement function of the motor. The Jacobian matrix.

[0046] Specifically, by taking the partial derivatives of each of the above measurement functions with respect to x, we obtain the partial derivative matrices corresponding to each measurement function. Finally, we use the partial derivative matrices as the Jacobian matrices of the corresponding measurement functions. .

[0047] 204. Determine the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip of the target oil and gas field power system. Based on the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip, determine the initial value of the operating state vector of the target oil and gas field power system, and use the initial value of the operating state vector as the current operating state vector under the current iteration number k.

[0048] Specifically, operational data of the target oil and gas field's power system over a past period is collected and analyzed. This historical data, considering factors such as seasonal variations and daily load patterns, is used to determine initial values ​​for node voltage amplitudes. The topology of the target oil and gas field's power system is analyzed to understand the connections and electrical distances between nodes. For radial or simple ring-shaped power systems, the voltage phase angles of each node relative to a reference node can be preliminarily estimated based on the location of the power source and line parameters. These data are then used as a reference, adjusted according to the specific characteristics of the target system, to determine the initial values ​​for node voltage phase angles. Historical operational experience of induction motors in the target oil and gas field's power system is summarized. The slip variation of the motors under different loads and operating times is understood, and combined with the current operating status of the motors, such as load size and operating time, an initial value for slip is reasonably determined. The determined initial values ​​for node voltage amplitudes, node voltage phase angles, and induction motor slip are combined in a specific order to construct the initial values ​​for the operating state vector of the target oil and gas field's power system. The determined initial value of the operating state vector is used as the current operating state vector for the current iteration number. In subsequent iterations, the operating state vector will be continuously updated and optimized based on this vector and relevant measurement functions and algorithms to gradually approximate the actual operating state of the power system.

[0049] 205. Substitute the current operating state vector into the measurement function to perform measurement estimation and obtain the measurement estimate. Substitute the current operating state vector into the Jacobian matrix to perform Jacobian matrix estimation and obtain the Jacobian matrix estimate. Based on the actual measurement and measurement estimate of the target oil and gas field power system, determine the measurement residual.

[0050] Specifically, the current operating state vector is substituted into the node voltage amplitude measurement function. Node-injected active power measurement function Node-injected reactive power measurement function Motor slip measurement function In this process, the corresponding measurement estimates are obtained; simultaneously, the current operating state vector is substituted into the node voltage amplitude measurement function. The corresponding Jacobian matrix and nodal injected active power measurement function The corresponding Jacobian matrix and nodal injected reactive power measurement function The corresponding Jacobian matrix and motor slip measurement function The corresponding Jacobian matrix yields the Jacobian matrix estimate. Simultaneously, real-time measurements of the power system (actual node voltage amplitude, actual node voltage phase angle, actual induction motor slip) are acquired using various measuring devices installed in the power system. The difference between the estimated measurements and the corresponding actual measurements is used as the measurement residual.

[0051] 206. Based on the measurement residual and the Jacobian matrix estimate, determine the change in the current operating state corresponding to the current operating state vector.

[0052] Specifically, based on measurement residuals Jacobian matrix estimator Solving the following formula yields the change in the current running state corresponding to the current running state vector. :

[0053] 207. Based on the changes in the current operating status, determine the predicted operating status of the power system of the target oil and gas field.

[0054] In this embodiment of the invention, after determining the change in the current operating state, it is necessary to determine the predicted operating state of the power system (predicted node voltage amplitude, predicted node voltage phase angle, and predicted induction motor slip rate) based on the change in the current operating state. Therefore, step 207 specifically includes: determining the change in the current operating state. Whether it is less than a preset change threshold; if the change in the current operating state is less than the preset change threshold. If the change is less than the preset threshold, then based on the current running state vector... and the change in the current operating state Determine the predicted operating status of the power system of the target oil and gas field. ,in, If the current operating state changes If the change is greater than or equal to the preset change threshold, then the... As the new current running state vector under the new current iteration number k+1 Based on the measurement function and the Jacobian matrix, the new current running state vector is... The solution is iteratively solved, and the final solution that meets the convergence condition is taken as the predicted operating state of the target oil and gas field power system.

[0055] Specifically, the preset change threshold is determined based on actual needs. ,like Then As a predicted operating state of the oil and gas field power system, otherwise, assume Continue iterating until the change in running state is obtained after n iterations. Less than the preset change threshold So far, the operating states obtained through n iterations are used as the predicted operating states of the target oil and gas field power system. It should be noted that the iteration process solves for the operating state vector, and the operating states of the power system, such as the node voltage magnitude, node voltage phase angle, and induction motor slip, are solved from the operating state vector.

[0056] Furthermore, to ensure the effectiveness of the predicted operating state, it is necessary to identify and remove abnormal data in the predicted operating states. Based on this, the method includes: taking any one of the multiple predicted operating states as a target predicted operating state. The target's predicted running state Substituting the values ​​into the measurement function for vector estimation yields the predicted operating state vector. And based on the predicted running state vector The predicted running state vector The corresponding actual running state vector Determine the target prediction running state Corresponding abnormal operation status evaluation value ,in, ,in, The weighting coefficient is used to determine the evaluation value of the abnormal operating state. If the value exceeds a preset evaluation threshold, then based on the actual operating state vector... and the predicted running state vector Determine the measurement residual ,in, Based on the measurement residual The target predicted operating state Jacobian matrix The residual covariance matrix is ​​determined, and the residual covariance matrix is ​​standardized to obtain the target prediction running state. The corresponding adverse data evaluation value; the predicted operating state corresponding to the largest adverse data evaluation value in each predicted operating state is taken as the abnormal predicted operating state, and the abnormal predicted operating state in the predicted operating state is removed to obtain the cleaned predicted operating state corresponding to the power system of the target oil and gas field.

[0057] Specifically, monitoring of adverse data involves calculating the objective function, i.e., the evaluation value of abnormal operating conditions. ,in, Set preset evaluation thresholds according to actual needs. At the significance level If If so, it is determined that there is bad data (abnormal data). At this time, the measurement residual is calculated. And determine the covariance matrix R of the measurement residual r, and determine the target prediction running state of the measurement function. The partial derivatives, i.e., the Jacobian matrix H, are then used to determine the residual covariance matrix of branch ij according to the following formula based on the covariance matrix R and the Jacobian matrix H. :

[0058] Furthermore, based on the residual covariance matrix The residual covariance matrix is ​​standardized according to the following formula to obtain the standardized residual covariance matrix. (Poor Data Evaluation Value)

[0059] Finally, find Largest measurement data This is used as abnormal data (abnormal predicted operating states). Abnormal data in the predicted operating states is then removed to obtain a high-quality operating state in the power system. In this embodiment of the invention, the high-quality operating state can also be used for subsequent state estimation of the power system.

[0060] The present invention provides another method for predicting the operating state of an oil and gas field power system. Compared with the current method of estimating the operating state of the power grid by modeling the load as a constant power load or a simple static load model, this invention constructs a dynamic equivalent model of the induction motor corresponding to the target oil and gas field power system and predicts the operating state of the power system based on the dynamic equivalent model of the induction motor. This ensures that the model on which the state estimation is based is consistent with the actual physical characteristics of the power system, improves the simulation capability of the system's dynamic behavior, and thus improves the accuracy of the power system operating state prediction. By determining the Jacobian matrix corresponding to different measurement functions, the relationship between the measurement functions and state variables can be described more accurately, which helps to converge to the true solution faster and more accurately in the iterative solution process. Using real-time measurement data for state prediction can promptly detect abnormal changes in the system and provide real-time data support for system stability analysis.

[0061] Furthermore, as Figure 1 In specific implementation, embodiments of the present invention provide a device for predicting the operating status of an oil and gas field power system, such as... Figure 3 As shown, the device includes: an acquisition unit 31, a model building unit 32, a function building unit 33, and a solution unit 34.

[0062] The acquisition unit 31 can be used to acquire the topology, branch parameters, and real-time measurement data of the target oil and gas field power system in response to the operation status prediction signal of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power.

[0063] The model building unit 32 can be used to build a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field based on the topology, the branch parameters, and the real-time measurement data.

[0064] The function construction unit 33 can be used to construct measurement functions for different operating state parameters of the target oil and gas field power system based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, and determine the Jacobian matrix corresponding to different measurement functions.

[0065] The solution unit 34 can be used to determine the initial operating state parameters, and based on the measurement function and the Jacobian matrix, iteratively solve the initial operating state parameters, and finally use the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system.

[0066] In specific application scenarios, in order to construct a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field, such as Figure 4 As shown, the model building unit 32 includes a building module 321, a decomposition module 322, a first determining module 323, and a transformation module 324.

[0067] The construction module 321 can be used to construct the equivalent admittance function of the induction motor in the target oil and gas field power system based on the topology, the branch parameters, and the real-time measurement data. ,in, , The stator impedance of the induction motor, For resistors, For imaginary numbers, For the stator reactance of the induction motor, For rotor impedance, The resistance of the rotor impedance The reactance component of the rotor impedance. For the reactance of the excitation branch, s is the excitation reactance, and s is the slip of the induction motor.

[0068] The decomposition module 322 can be used to decompose the equivalent admittance function. Decomposed into equivalent conductance and equivalent susceptance ,in, .

[0069] The first determining module 323 can be used to determine the equivalent conductivity. and the equivalent susceptance Determine the node-injected active power function of the induction motor connected to node i in the power system of the target oil and gas field. and node-injected reactive power function ,in, , , Let x be the voltage amplitude of node i, and let x be the operating state variable of the target oil and gas field power system.

[0070] The construction module 321 can also be used to construct the electromagnetic torque function of the induction motor based on the topology, the branch parameters, and the real-time measurement data. and load torque function ,in, , , The mechanical power output by the induction motor is s, where s is the slip rate. For synchronous angular velocity, For constant loss torque, This is the load factor.

[0071] The conversion module 324 can be used to convert the electromagnetic torque function when the electromagnetic torque of the induction motor and the load torque are equal. and the load torque function Performing a function transformation, the slip ratio constructor is obtained as follows: ,in, The output mechanical power of the induction motor is given by the slip s and voltage V.

[0072] The first determining module 323 can be used to determine the equivalent admittance function. The node injects active power function The node injects reactive power function The slip rate constructor is used as the dynamic equivalent model of the induction motor.

[0073] In specific application scenarios, in order to construct measurement functions for different operating state parameters of the target oil and gas field power system, the function construction unit 33 includes a second determination module 331 and a function construction module 332.

[0074] The second determining module 331 can be used to determine the node voltage amplitude measurement function based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data. ,in, x is the running state vector. Let be the voltage at node i.

[0075] The second determining module 331 can also be used to determine the node-injected active power measurement function based on the induction motor dynamic equivalent model, the topology, the branch parameters, and the real-time measurement data. and node-injected reactive power measurement function ,in, , x is the running state vector, The voltage amplitude at node i, The voltage amplitude at node j, For the real part of the branch ij admittance, For the imaginary part of the admittance of branch ij, Let be the phase angle difference between node i and node j.

[0076] The function construction module 332 can be used to construct a motor slip measurement function based on the induction motor dynamic equivalent model, the topology, the branch parameters, and the real-time measurement data. ,in, x is the running state vector. The electromagnetic torque of the induction motor under slip s and voltage V, Let be the load torque of the induction motor at slip s.

[0077] In specific application scenarios, in order to determine the Jacobian matrix corresponding to different measurement functions, the second determining module 331 can also be used to determine the node voltage amplitude measurement function. The partial derivative matrix of the node voltage amplitude measurement function The partial derivative matrix is ​​used as the node voltage amplitude measurement function. The Jacobian matrix; determine the node-injected active power measurement function. The partial derivative matrix of the function is used to inject the active power measurement function into the node. The partial derivative matrix is ​​used as the active power measurement function injected at the node. The Jacobian matrix; determine the nodal injected reactive power measurement function. The partial derivative matrix of the node is used to inject reactive power measurement function. The partial derivative matrix is ​​used as the reactive power measurement function injected at the node. The Jacobian matrix; determine the slip measurement function of the motor. The partial derivative matrix of the motor slip measurement function is used to... The partial derivative matrix is ​​used as the slip measurement function of the motor. The Jacobian matrix.

[0078] In specific application scenarios, the initial operating state parameters include the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip. In order to determine the predicted operating state of the target oil and gas field power system, the solution unit 34 includes a third determination module 341 and an estimation module 342.

[0079] The third determining module 341 can be used to determine the initial value of the operating state vector of the target oil and gas field power system based on the initial value of the node voltage amplitude, the initial value of the node voltage phase angle, and the initial value of the induction motor slip, and use the initial value of the operating state vector as the current operating state vector under the current iteration number k. .

[0080] The estimation module 342 can be used to estimate the current running state vector. Substituting the values ​​into the measurement function, we can perform measurement estimation to obtain the measurement estimate. The current running state vector Substituting the values ​​into the Jacobian matrix, we can perform Jacobian matrix estimation to obtain the Jacobian matrix estimator. And based on the actual quantity measurement of the power system of the target oil and gas field. and the measurement estimate Determine the measurement residual ,in, .

[0081] The third determining module 341 can also be used to determine the measurement residual. The Jacobian matrix estimator Determine the change in the current running state corresponding to the current running state vector. ,in, ,in, The weight matrix is ​​set according to actual needs.

[0082] The third determining module 341 can also be used to determine the amount of change in the current operating state. The predicted operating status of the power system of the target oil and gas field is determined.

[0083] In specific application scenarios, in order to determine the predicted operating status of the power system of a target oil and gas field, the third determining module 341 can be used to determine the change in the current operating status. Whether it is less than a preset change threshold; if the change in the current operating state is less than the preset change threshold. If the change is less than the preset threshold, then based on the current running state vector... and the change in the current operating state Determine the predicted operating status of the power system of the target oil and gas field. ,in, If the current operating state changes If the change is greater than or equal to the preset change threshold, then the... As the new current running state vector under the new current iteration number k+1 Based on the measurement function and the Jacobian matrix, the new current running state vector is... The solution is iteratively solved, and the final solution that meets the convergence condition is taken as the predicted operating state of the target oil and gas field power system.

[0084] In specific application scenarios, in order to identify and remove abnormal states in the predicted operating state, the device also includes an anomaly removal unit 35.

[0085] The de-identification unit 35 can be used to treat any one of the multiple predicted operating states as a target predicted operating state. The target's predicted running state Substituting the values ​​into the measurement function for vector estimation yields the predicted operating state vector. And based on the predicted running state vector The predicted running state vector The corresponding actual running state vector Determine the target prediction running state Corresponding abnormal operation status evaluation value ,in, ,in, The weighting coefficient is used to determine the evaluation value of the abnormal operating state. If the value exceeds a preset evaluation threshold, then based on the actual operating state vector... and the predicted running state vector Determine the measurement residual ,in, Based on the measurement residual The target predicted operating state Jacobian matrix The residual covariance matrix is ​​determined, and the residual covariance matrix is ​​standardized to obtain the target prediction running state. The corresponding adverse data evaluation value; the predicted operating state corresponding to the largest adverse data evaluation value in each predicted operating state is taken as the abnormal predicted operating state, and the abnormal predicted operating state in the predicted operating state is removed to obtain the cleaned predicted operating state corresponding to the power system of the target oil and gas field.

[0086] It should be noted that other corresponding descriptions of the functional modules involved in the oil and gas field power system operation status prediction device provided in this embodiment of the invention can be found in [reference needed]. Figure 1 The corresponding description of the method shown will not be repeated here.

[0087] Based on the above, Figure 1The method shown in the invention, correspondingly, also provides a computer-readable storage medium storing a computer program thereon. When executed by a processor, this program performs the following steps: In response to a prediction signal of the operating status of a target oil and gas field power system, acquiring the topology, branch parameters, and real-time measurement data of the target oil and gas field power system, wherein the branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power; based on the topology, the branch parameters, and the real-time measurement data... A dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field is constructed. Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, measurement functions for different operating state parameters of the power system of the target oil and gas field are constructed, and the Jacobian matrix corresponding to different measurement functions is determined. Initial operating state parameters are determined, and the initial operating state parameters are iteratively solved based on the measurement functions and the Jacobian matrix. Finally, the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the power system of the target oil and gas field.

[0088] Based on the above, Figure 1 The method shown and as Figure 3 The embodiment of the device shown in the invention also provides a physical structure diagram of a computer device, such as... Figure 5 As shown, the computer device includes: a processor 41, a memory 42, and a computer program stored in the memory 42 and executable on the processor. Both the memory 42 and the processor 41 are mounted on a bus 43. When the processor 41 executes the program, it performs the following steps: In response to a prediction signal of the operating status of the target oil and gas field power system, it acquires the topology, branch parameters, and real-time measurement data of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power. Based on the topology... Based on the branch parameters and the real-time measurement data, a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field is constructed. Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, measurement functions for different operating state parameters of the power system of the target oil and gas field are constructed, and the Jacobian matrix corresponding to different measurement functions is determined. Initial operating state parameters are determined, and based on the measurement functions and the Jacobian matrix, the initial operating state parameters are iteratively solved. Finally, the operating state obtained from the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the power system of the target oil and gas field.

[0089] Through the technical solution of this invention, by constructing a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field, and predicting the operating state of the power system based on the dynamic equivalent model of the induction motor, this invention can ensure that the model on which the state estimation is based is consistent with the actual physical characteristics of the power system, improve the simulation capability of the system's dynamic behavior, and thus improve the accuracy of the prediction of the power system's operating state. By determining the Jacobian matrix corresponding to different measurement functions, the relationship between the measurement functions and state variables can be described more accurately, which helps to converge to the true solution faster and more accurately in the iterative solution process. Using real-time measurement data for state prediction can promptly detect abnormal changes in the system and provide real-time data support for system stability analysis.

[0090] It is obvious to those skilled in the art that the modules or steps of the present invention described above can be implemented using general-purpose computing devices. They can be centralized on a single computing device or distributed across a network of multiple computing devices. Optionally, they can be implemented using computer-executable program code, thereby storing them in a storage device for execution by a computing device. In some cases, the steps shown or described can be performed in a different order than those presented herein, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any particular combination of hardware and software.

[0091] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for predicting the operating status of an oil and gas field power system, characterized in that, include: In response to the predicted operating status signal of the target oil and gas field power system, the topology, branch parameters, and real-time measurement data of the target oil and gas field power system are acquired. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power. Based on the topology, the branch parameters, and the real-time measurement data, a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field is constructed. Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, measurement functions for different operating state parameters of the target oil and gas field power system are constructed, and the Jacobian matrix corresponding to different measurement functions is determined. The initial operating state parameters are determined, and the initial operating state parameters are iteratively solved based on the measurement function and the Jacobian matrix. Finally, the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the target oil and gas field power system.

2. The method according to claim 1, characterized in that, Based on the aforementioned topology, branch parameters, and real-time measurement data, a dynamic equivalent model of the induction motors corresponding to the power system of the target oil and gas field is constructed, including: Based on the aforementioned topology, branch parameters, and real-time measurement data, the equivalent admittance function of the induction motor in the target oil and gas field's power system is constructed. ,in, , The stator impedance of the induction motor, For resistors, For imaginary numbers, For the stator reactance of the induction motor, For rotor impedance, The resistance of the rotor impedance The reactance component of the rotor impedance. For the reactance of the excitation branch, s is the magnetizing reactance, and s is the slip of the induction motor; The equivalent admittance function Decomposed into equivalent conductance and equivalent susceptance ,in, ; Based on the equivalent conductivity and the equivalent susceptance Determine the node-injected active power function of the induction motor connected to node i in the power system of the target oil and gas field. and node-injected reactive power function ,in, , , Let x be the voltage amplitude at node i, and let x be the operating state variable of the target oil and gas field power system. Based on the aforementioned topology, branch parameters, and real-time measurement data, the electromagnetic torque function of the induction motor is constructed. and load torque function ,in, , , The mechanical power output by the induction motor is s, where s is the slip rate. For synchronous angular velocity, For constant loss torque, This is the load factor; When the electromagnetic torque of the induction motor and the load torque are equal, the electromagnetic torque function is... and the load torque function Performing a function transformation, the slip ratio constructor is obtained as follows: ,in, The output mechanical power of the induction motor under slip s and voltage V; The equivalent admittance function The node injects active power function The node injects reactive power function The slip rate constructor is used as the dynamic equivalent model of the induction motor.

3. The method according to claim 1, characterized in that, The measurement functions for constructing different operating state parameters of the target oil and gas field power system based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data include: Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, the node voltage amplitude measurement function is determined. ,in, x is the running state vector. Let be the voltage at node i; Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, the node-injected active power measurement function is determined. and node-injected reactive power measurement function ,in, , x is the running state vector, The voltage amplitude at node i, The voltage amplitude at node j, For the real part of the branch ij admittance, For the imaginary part of the admittance of branch ij, Let be the phase angle difference between node i and node j; Based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, a motor slip measurement function is constructed. ,in, x is the running state vector. The electromagnetic torque of the induction motor under slip s and voltage V, Let be the load torque of the induction motor at slip s.

4. The method according to claim 3, characterized in that, Determining the Jacobian matrix corresponding to different measurement functions includes: Determine the node voltage amplitude measurement function The partial derivative matrix of the node voltage amplitude measurement function The partial derivative matrix is ​​used as the node voltage amplitude measurement function. The Jacobian matrix; Determine the active power measurement function injected at the node. The partial derivative matrix of the function is used to inject the active power measurement function into the node. The partial derivative matrix is ​​used as the active power measurement function injected at the node. The Jacobian matrix; Determine the reactive power measurement function injected at the node. The partial derivative matrix of the node is used to inject reactive power measurement function. The partial derivative matrix is ​​used as the reactive power measurement function injected at the node. The Jacobian matrix; Determine the motor slip measurement function The partial derivative matrix of the motor slip measurement function is used to... The partial derivative matrix is ​​used as the slip measurement function of the motor. The Jacobian matrix.

5. The method according to claim 1, characterized in that, The initial operating state parameters include the initial values ​​of node voltage amplitude, node voltage phase angle, and induction motor slip. The process of iteratively solving the initial operating state parameters based on the measurement function and the Jacobian matrix, and finally using the operating state obtained from the final solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system, includes: Based on the initial values ​​of the node voltage amplitude, the node voltage phase angle, and the induction motor slip, the initial value of the operating state vector of the target oil and gas field power system is determined, and this initial value of the operating state vector is used as the current operating state vector at the current iteration number k. ; The current running state vector Substituting the values ​​into the measurement function, we can perform measurement estimation to obtain the measurement estimate. The current running state vector Substituting the values ​​into the Jacobian matrix, we can perform Jacobian matrix estimation to obtain the Jacobian matrix estimator. And based on the actual quantity measurement of the power system of the target oil and gas field. and the measurement estimate Determine the measurement residual ,in, ; Based on the measurement residual The Jacobian matrix estimator Determine the change in the current running state corresponding to the current running state vector. ,in, ,in, The weight matrix is ​​set according to actual needs; Based on the change in the current operating status The predicted operating status of the power system of the target oil and gas field is determined.

6. The method according to claim 5, characterized in that, The amount of change based on the current operating state Determining the predicted operating status of the target oil and gas field's power system includes: Determine the change in the current operating state. Is it less than a preset change threshold? If the current operating state changes If the change is less than the preset threshold, then based on the current running state vector... and the change in the current operating state Determine the predicted operating status of the power system of the target oil and gas field. ,in, ; If the current operating state changes If the change is greater than or equal to the preset change threshold, then the... As the new current running state vector under the new current iteration number k+1 Based on the measurement function and the Jacobian matrix, the new current running state vector is... The solution is iteratively solved, and the final solution that meets the convergence condition is taken as the predicted operating state of the target oil and gas field power system.

7. The method according to claim 1, characterized in that, Finally, after taking the final solution that satisfies the convergence condition of the iterative solution as the predicted operating state of the target oil and gas field power system, the method further includes: Each of the multiple predicted operating states is taken as a target predicted operating state. The target's predicted running state Substituting the values ​​into the measurement function for vector estimation yields the predicted operating state vector. And based on the predicted running state vector The predicted running state vector The corresponding actual running state vector Determine the target prediction running state Corresponding abnormal operation status evaluation value ,in, ,in, These are the weighting coefficients; Determine the abnormal operating status evaluation value If the value exceeds a preset evaluation threshold, then based on the actual operating state vector... and the predicted running state vector Determine the measurement residual ,in, ; Based on the measurement residual The target predicted operating state Jacobian matrix The residual covariance matrix is ​​determined, and the residual covariance matrix is ​​standardized to obtain the target prediction running state. Corresponding negative data evaluation values; The predicted operating state corresponding to the maximum adverse data evaluation value in each predicted operating state is taken as the abnormal predicted operating state, and the abnormal predicted operating states in the predicted operating states are removed to obtain the cleaned predicted operating state corresponding to the power system of the target oil and gas field.

8. A device for predicting the operating status of an oil and gas field power system, characterized in that, include: The acquisition unit is used to acquire the topology, branch parameters, and real-time measurement data of the target oil and gas field power system in response to the operation status prediction signal of the target oil and gas field power system. The branch parameters include branch admittance and node admittance matrix, and the real-time measurement data includes node voltage amplitude, line active power, line reactive power, node injected active power, and node injected reactive power. The model building unit is used to build a dynamic equivalent model of the induction motor corresponding to the power system of the target oil and gas field based on the topology, the branch parameters, and the real-time measurement data. The function construction unit is used to construct measurement functions for different operating state parameters of the target oil and gas field power system based on the dynamic equivalent model of the induction motor, the topology, the branch parameters, and the real-time measurement data, and to determine the Jacobian matrix corresponding to different measurement functions. The solution unit is used to determine the initial operating state parameters, and iteratively solve the initial operating state parameters based on the measurement function and the Jacobian matrix. Finally, the operating state obtained by the final solution that satisfies the convergence condition of the iterative solution is taken as the predicted operating state of the target oil and gas field power system.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images