Hydrogen-doped natural gas regulation and control system and method based on data analysis

Abstract

The invention relates to the technical field of control systems, in particular to a hydrogen-doped natural gas regulation and control system and method based on data analysis, and the system comprises a data preprocessing module which collects original hydrogen-doped natural gas process data, divides gas flow, pressure and temperature parameters according to a time sequence, and calculates an average value and maximum and minimum values in each time period; according to the method, short-term fluctuation of gas flow and pressure is extracted through a convolution kernel in the feature extraction stage, local feature vectors are generated, accurate input is provided for time sequence modeling, the time sequence modeling captures trend change based on long-term dependency, a stable time sequence model is established, and the prediction precision of the hydrogen doping proportion is improved; in the prediction optimization stage, the proportion is adjusted under the constraint of energy demand, carbon emission and other targets, multi-target optimization is carried out, deviation is corrected in real time in the dynamic adjustment stage, the adaptability and flexibility of the system are ensured, the response capability of the model is enhanced through optimization of multiple links, and the adaptability to complex environment changes is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of control systems, and in particular to a hydrogen-blended natural gas control system and method based on data analysis. Background Art

[0002] The field of control system technology focuses on how to design, implement and optimize systems for automated control. This field involves controlling the behavior of equipment or processes through mathematical models and algorithms to achieve desired performance goals. Control systems are widely used in many industries such as industrial production, aerospace, transportation, robotics, etc., aiming to improve system efficiency, stability and safety. With the improvement of computer technology and data processing capabilities, modern control systems have developed towards more intelligent, automated and real-time directions, using sensors, actuators and feedback mechanisms to precisely adjust complex systems.

[0003] Among them, the hydrogen-blended natural gas control system based on data analysis refers to the real-time monitoring and analysis of the proportion of hydrogen blended into natural gas, and the use of advanced control algorithms to adjust and optimize the hydrogen blending process. The main purpose of this system is to optimize the energy structure, improve the environmental protection and economic efficiency of natural gas use, and ensure the stability and safety of natural gas supply. Through data analysis, the system can dynamically adjust the hydrogen blending ratio according to different operating conditions to adapt to different needs and environmental conditions.

[0004] Existing technologies lack automation and real-time performance in data processing and system control. Traditional systems rely on manual intervention in abnormal data processing and lack automatic detection and supplementation mechanisms, resulting in missing data or abnormal values ​​that cannot be corrected in a timely manner, affecting prediction accuracy. The feature extraction and modeling parts fail to fully capture long-term trends and ignore the long-term dependencies between data, resulting in the system's inability to adapt to environmental and demand fluctuations in real time. Traditional methods usually rely on fixed rules to adjust the hydrogen blending ratio and are unable to make timely optimizations when the load or environment changes, resulting in the system's inability to balance energy demand and emission targets. For example, when climate conditions or demand change, the existing system is difficult to respond quickly, resulting in untimely adjustments to the hydrogen blending ratio, affecting system performance and economic benefits. Summary of the invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose a hydrogen-blended natural gas control system and method based on data analysis.

[0006] In order to achieve the above-mentioned purpose, the present invention adopts the following technical scheme: A hydrogen-blended natural gas control system based on data analysis comprises:

[0007] The data preprocessing module collects the original hydrogen-blended natural gas process data, divides the gas flow, pressure, and temperature parameters in chronological order, calculates the average value and the maximum and minimum values ​​in each time period, identifies data anomalies and interpolates and replaces the abnormal values, calculates the change range of adjacent time periods and standardizes all data, performs time series interpolation to supplement missing data, and generates standardized time series data;

[0008] The feature extraction module divides the sequence by time window based on the time series data after the standardization processing, inputs the data of multiple time windows into the convolution kernel for calculation, extracts the short-term fluctuation and change amplitude of the gas flow and pressure, completes the vector construction of the local feature of the time series, and normalizes the time window feature to generate a local feature vector;

[0009] The time series modeling module analyzes the trend in the time window based on the local feature vector, calculates the long-term dependency between time series, gradually transfers the hidden state to calculate the changes in adjacent time periods, predicts the long-term trend, captures the dependency of time series through recursion, and generates a time series dependency model;

[0010] The prediction optimization module calculates the expected value of the hydrogen blending ratio based on the time series dependency model, constrains the energy demand, carbon emissions, stability indicators and the hydrogen blending ratio, adjusts the ratio value and applies the constraint condition to the predicted value, performs iterative correction, and generates an optimized hydrogen blending ratio prediction;

[0011] The dynamic adjustment module monitors the deviation between the real-time feedback data and the predicted value based on the optimized hydrogen blending ratio prediction, analyzes the deviation trend and corrects the prediction model parameters, recalculates the hydrogen blending ratio, updates the adjusted parameters, applies them to the current prediction results, and generates a dynamically adjusted hydrogen blending ratio prediction.

[0012] The standardized time series data includes time series data, outlier processing results, and missing data filling results. The local feature vector includes gas flow characteristics, pressure change characteristics, and short-term fluctuation characteristics. The time series dependency model includes time period hidden state, long-term trend capture, and time series recursive relationship. The optimized hydrogen blending ratio prediction includes energy demand adjustment results, carbon emission optimization results, and stability balance results. The dynamically adjusted hydrogen blending ratio prediction includes feedback deviation correction results, adjusted parameters, and updated prediction results.

[0013] As a further solution of the present invention, the step of obtaining the time series data after the standardization process is specifically as follows:

[0014] Collect the original hydrogen-blended natural gas process data, and divide the gas flow, pressure and temperature parameters in the time series into multiple intervals according to time, calculate the average, maximum and minimum values ​​of the gas flow, pressure and temperature in each time interval, and obtain a statistical parameter set;

[0015] Based on the statistical parameter set, an outlier threshold is set and data points exceeding the threshold are identified, and the outlier data points are replaced by a linear interpolation method to generate a statistical parameter set after interpolation processing;

[0016] The statistical parameter set after the interpolation process is used to calculate the variation range of gas flow, pressure and temperature in adjacent time intervals, using the formula:

[0017]

[0018] Perform standardization processing and calculate and generate a standardized set of change amplitudes;

[0019] Among them, S i represents the standardized change in the ith interval, X i+1 represents the gas parameter value in the i+1th time interval, X i represents the gas parameter value in the i-th time interval, is the mean of the gas parameters in all intervals, n is the total number of intervals, ω is the weighting factor, γ is the variance adjustment coefficient, and δ is the stabilization increment;

[0020] In combination with the standardized set of variation amplitudes, a time series interpolation method is applied to complete the missing data, and the completed data is merged with the standardized variation amplitude to obtain a standardized time series data set.

[0021] As a further solution of the present invention, the step of obtaining the local feature vector is specifically:

[0022] Dividing the standardized time series data into multiple subsequences according to fixed time windows, extracting the change values ​​of gas flow and pressure for each subsequence, and obtaining a time window change amplitude set by calculating the difference of the parameters in each time window;

[0023] Using the time window change amplitude set, the data of multiple time windows are input into the convolution kernel, the local response value of the convolution kernel to each window is calculated, and the weight and relative change ratio are combined for processing, using the formula:

[0024]

[0025] Extract local change features and generate a set of local feature values;

[0026] Among them, F iis the local eigenvalue of the i-th time window, W k is the kth convolution kernel weight, V i+k and V i are the gas parameter values ​​of the i+kth and ith time windows respectively, α is the enhancement index, γ is the variance adjustment coefficient, is the mean value of the gas parameter in all time windows, m is the total number of time windows, δ is the stabilization increment, η is the maximum value adjustment coefficient, and max(V) is the maximum value of the gas parameter in the time window;

[0027] The local feature value set is integrated into a feature vector according to the time window sequence, and normalized to generate a local feature vector.

[0028] As a further solution of the present invention, the step of acquiring the timing dependency model is specifically as follows:

[0029] Based on the local feature vector, the time window is divided and the local feature values ​​of multiple windows are extracted, the trend change rate in each time window is calculated, and the trend change rate and the time window sequence are combined to obtain a trend sequence set;

[0030] Based on the trend sequence set, the trend cumulative effect between time windows is analyzed, and the hidden state value is gradually transferred using the formula:

[0031]

[0032] Calculate the long-term dependencies between adjacent time windows and generate a sequence of hidden state values;

[0033] Among them, H t is the hidden state value of the t-th time window, φ i is the weight of the i-th time window, T t-i Represents the trend change value of the ti time window, T t represents the trend change value of the t-th time window, ψ is the long-term dependence adjustment coefficient, ω k is the short-term adjustment coefficient, γ j is the variance adjustment factor, is the trend mean, δ is the stabilization parameter, and η is the offset parameter;

[0034] Combined with the hidden state value sequence, the dependency relationship between time series is recursively calculated, the long-term dependency of multiple time windows is cumulatively analyzed, and the dependency relationship of multiple windows is integrated to generate a temporal dependency model.

[0035] As a further solution of the present invention, the steps for obtaining the optimized hydrogen blending ratio prediction are specifically as follows:

[0036] Using the time series dependency model, inputting energy demand, carbon emission and stability index parameters, calculating the expected value of hydrogen blending ratio in each window in combination with the time window sequence, analyzing the constraint conditions and integrating the results, and generating a ratio constraint set;

[0037] Based on the ratio constraint set, the hydrogen blending ratio is gradually adjusted and the constraint conditions are applied, using the formula:

[0038]

[0039] Calculate and generate a corrected set of ratio values;

[0040] Among them, P t is the corrected hydrogen blending ratio in the tth time window, λ i is the weight of the i-th window, E t-i is the energy demand of the ti window, S t is the stability index of the tth window, α is the stability adjustment coefficient, C t is the carbon emission value of the tth window, γ is the carbon emission weight, β j is the variance adjustment factor, is the mean value of hydrogen blending ratio, δ is the stabilization parameter, ∈ is the energy demand adjustment coefficient, and max(E) is the maximum energy demand value;

[0041] The corrected ratio value set is used to recursively adjust the ratio value of each time window, and the correction results are integrated with the time window constraints and energy, emission, and stability requirements to generate an optimized hydrogen blending ratio prediction.

[0042] As a further solution of the present invention, the steps for obtaining the dynamically adjusted hydrogen blending ratio prediction are specifically as follows:

[0043] Based on the optimized prediction of hydrogen blending ratio, real-time feedback data is collected and the deviation between the feedback data and the predicted value is calculated, the deviation is integrated according to the time window and then the deviation change trend is calculated to generate a deviation trend set;

[0044] Based on the deviation trend set, the deviation variation law between time windows is analyzed and the model parameters are corrected using the formula:

[0045]

[0046] Adjust parameter values ​​to generate a revised set of model parameters;

[0047] Among them, θ t is the model parameter adjusted in the tth time window, θ t-1 is the model parameter of the previous time window, λ i is the weight coefficient, D t-iis the deviation value of the ti-th window, Δ t is the cumulative deviation, φ is the deviation adjustment coefficient, β j is the variance adjustment factor, is the mean deviation, δ is the stabilization parameter, ∈ is the deviation amplitude adjustment coefficient, max(D) is the maximum deviation, n represents the number of previous time windows used for deviation calculation, and m represents the number of previous time windows used for variance calculation;

[0048] The hydrogen blending ratio value is recalculated and adjusted using the revised model parameter set, and the revised ratio value is applied to the prediction result of the current time window to generate a dynamically adjusted hydrogen blending ratio prediction.

[0049] A method for regulating hydrogen-blended natural gas based on data analysis, wherein the method for regulating hydrogen-blended natural gas based on data analysis is performed based on the above-mentioned hydrogen-blended natural gas regulation system based on data analysis, and comprises the following steps:

[0050] S1: Collect the original hydrogen-blended natural gas process data, divide the time periods of gas flow, pressure, and temperature parameters in chronological order, calculate the average, maximum, and minimum values ​​of gas flow, pressure, and temperature in each time period, identify abnormal values ​​and perform interpolation replacement, calculate the change range of adjacent time periods and standardize the data to generate standardized time series data;

[0051] S2: Based on the standardized processing time series data, the data is divided according to the set time window, the data of each time window is sent to the convolution kernel for calculation, the short-term fluctuation and change amplitude of the gas flow and pressure are extracted, and the local feature vector is constructed;

[0052] S3: Based on the local feature vector, analyze the trend in the time window, calculate the long-term dependency between time series, gradually transfer the hidden state, calculate the changes in adjacent time periods, capture the time series dependency, and generate a time series dependency model;

[0053] S4: Based on the time series dependency model, the expected value of the hydrogen blending ratio is calculated, and the hydrogen blending ratio is adjusted in combination with the constraints including energy demand, carbon emissions, and stability, and the constraints are applied to the predicted value for iterative correction to generate an optimized hydrogen blending ratio prediction;

[0054] S5: Based on the optimized hydrogen blending ratio prediction, real-time monitoring of the deviation between the feedback data and the predicted value, analyzing the deviation trend, correcting the prediction model parameters, recalculating the hydrogen blending ratio, updating the adjusted parameters, and generating a dynamically adjusted hydrogen blending ratio prediction;

[0055] S6: Based on the dynamically adjusted hydrogen blending ratio prediction, combined with real-time feedback data, difference analysis and trend adjustment are performed to monitor the real-time changes of the hydrogen blending ratio, and the hydrogen blending ratio is iteratively optimized according to demand and stability requirements to generate an optimized value of the hydrogen blending ratio.

[0056] Compared with the prior art, the advantages and positive effects of the present invention are:

[0057] In the present invention, the time series data processing and control in the process of hydrogen-blended natural gas are optimized through the steps of data preprocessing, local feature extraction, time series modeling and dynamic adjustment. In the feature extraction stage, the short-term fluctuations of gas flow and pressure are extracted by convolution kernels to generate local feature vectors, providing accurate input for time series modeling. Time series modeling captures trend changes based on long-term dependencies, establishes a stable time series model, and improves the prediction accuracy of hydrogen blending ratio. In the prediction optimization stage, the ratio is adjusted under the constraints of energy demand, carbon emissions and other targets, and multi-objective optimization is performed. In the dynamic adjustment stage, deviations are corrected in real time to ensure the adaptability and flexibility of the system. The optimization of multiple links enhances the responsiveness of the model and improves its adaptability to complex environmental changes. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 is a system flow chart of the present invention;

[0059] Figure 2 This is a flow chart of the steps for obtaining time series data after standardization processing of the present invention;

[0060] Figure 3 Flow chart of the steps of obtaining the local feature vector of the present invention;

[0061] Figure 4 A flow chart of the steps for obtaining the timing dependency model of the present invention;

[0062] Figure 5 A flow chart of the steps for obtaining the prediction of the optimized hydrogen blending ratio of the present invention;

[0063] Figure 6 The flowchart is a step chart of obtaining the prediction of the hydrogen blending ratio after dynamic adjustment according to the present invention. DETAILED DESCRIPTION

[0064] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0065] In the description of the present invention, it should be understood that the terms "length", "width", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate positions or positional relationships based on the positions or positional relationships shown in the drawings, and are 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 limiting the present invention. In addition, in the description of the present invention, "multiple" means two or more, unless otherwise clearly and specifically defined.

[0066] Embodiment 1

[0067] See also Figure 1 The present invention provides a technical solution: a hydrogen-blended natural gas control system based on data analysis includes:

[0068] The data preprocessing module collects the original hydrogen-blended natural gas process data, divides the gas flow, pressure, and temperature parameters in chronological order, calculates the average value and the maximum and minimum values ​​in each time period, identifies data anomalies and interpolates and replaces the abnormal values, calculates the change range of adjacent time periods and standardizes all data, performs time series interpolation to supplement missing data, and generates standardized time series data;

[0069] The feature extraction module divides the sequence by time window based on the standardized time series data, inputs the data of multiple time windows into the convolution kernel for calculation, extracts the short-term fluctuation and change amplitude of gas flow and pressure, completes the vector construction of local time series features, and normalizes the time window features to generate local feature vectors;

[0070] The time series modeling module analyzes the trend within the time window based on the local feature vector, calculates the long-term dependency between time series, gradually transfers the hidden state to calculate the changes in adjacent time periods, predicts the long-term trend, captures the dependency of time series through recursion, and generates a time series dependency model;

[0071] The prediction optimization module calculates the expected value of hydrogen blending ratio based on the time-series dependency model, constrains energy demand, carbon emissions, stability indicators and hydrogen blending ratio, adjusts the ratio value and applies the constraint conditions to the predicted value, performs iterative correction, and generates an optimized hydrogen blending ratio prediction;

[0072] The dynamic adjustment module monitors the deviation between the real-time feedback data and the predicted value based on the optimized hydrogen blending ratio prediction, analyzes the deviation trend and corrects the prediction model parameters, recalculates the hydrogen blending ratio, updates the adjusted parameters, applies them to the current prediction results, and generates a dynamically adjusted hydrogen blending ratio prediction.

[0073] The time series data after standardization include time series data, outlier processing results, and missing data filling results. The local feature vectors include gas flow characteristics, pressure change characteristics, and short-term fluctuation characteristics. The time series dependency model includes time period hidden state, long-term trend capture, and time series recursive relationship. The optimized hydrogen blending ratio prediction includes energy demand adjustment results, carbon emission optimization results, and stability balance results. The dynamically adjusted hydrogen blending ratio prediction includes feedback deviation correction results, adjusted parameters, and updated prediction results.

[0074] See also Figure 2 , the specific steps for obtaining the standardized time series data are as follows:

[0075] Collect the original hydrogen-blended natural gas process data, and divide the gas flow, pressure and temperature parameters in the time series into multiple intervals according to time, calculate the average, maximum and minimum values ​​of the gas flow, pressure and temperature in each time interval, and obtain a statistical parameter set;

[0076] The raw data collected includes multiple parameters such as gas flow, pressure and temperature. The changes of these data at different time points can reflect the operating status of the system and possible problems. The time segmentation of these data is to be able to observe the changes of parameters in each time interval more carefully. The statistical data of each time interval, such as the average, maximum and minimum values, are calculated to summarize the data characteristics in each interval and provide basic information for subsequent data processing. Through this step, we can obtain a data set containing all key statistical information, which will be used to identify and process abnormal data, as well as for further data analysis and model building.

[0077] Based on the statistical parameter set, an outlier threshold is set and data points exceeding the threshold are identified, and the outlier data points are replaced by a linear interpolation method to generate a statistical parameter set after interpolation processing;

[0078] The outliers in the data are further analyzed using a set of statistical parameters. This analysis is based on the setting of outlier thresholds, which are based on previous data and statistical knowledge to determine what level of data deviation should be considered anomalies. By identifying and removing these outliers, subsequent analysis can be prevented from being misled by these data. Linear interpolation of these abnormal data points is to fill the gaps left by the removed data points and ensure the integrity of the data, which is crucial for accurate time series analysis. The interpolated statistical parameter set will be used to calculate the amplitude of changes in adjacent time intervals in order to explore the dynamic changes of time series.

[0079] The statistical parameter set after interpolation is used to calculate the change range of gas flow, pressure and temperature in adjacent time intervals using the formula:

[0080]

[0081] Perform standardization processing and calculate and generate a standardized set of change amplitudes;

[0082] Among them, S i represents the standardized change in the ith interval, X i+1 represents the gas parameter value in the i+1th time interval, X i represents the gas parameter value in the i-th time interval, is the mean of the gas parameters in all intervals, n is the total number of intervals, ω is the weighting factor, γ is the variance adjustment coefficient, and δ is the stabilization increment;

[0083] formula:

[0084]

[0085] The benefit of the formula is that by weighting the absolute differences and the adjusted variances, the formula effectively normalizes the volatility in the time series data so that the results are not overly affected by a single extreme data point.

[0086] Detailed explanation of the formula and the process of formula calculation and derivation:

[0087] Set X i+1 =105,X i =100, n=10, ω=1.5, γ=0.8, δ=1.

[0088] First calculate the absolute value difference:

[0089] |X i+1 -X i |=|105-100|=5

[0090] Then calculate the adjusted variance:

[0091]

[0092] Assume X j Evenly distributed between 95 and 105, we get:

[0093]

[0094] Add δ=1 to get 95.4, and finally calculate the standardized change:

[0095]

[0096] This result shows that after considering the weights and tuning parameters, i To X i+1The standardized change in is approximately 0.768, which reflects the standardized fluctuation from one time point to another.

[0097] Combined with the standardized change amplitude set, the time series interpolation method is applied to complete the missing data. At the same time, the completed data is merged with the standardized change amplitude to obtain a standardized time series data set.

[0098] Through the standardized change amplitude, the application of time series interpolation method is to deal with those discontinuous points caused by missing data. Such completion is to predict unknown data points through the time series pattern of existing data. The completed data is merged with the standardized change amplitude in order to obtain a complete and unified standard data set, which provides support for further data analysis and decision-making. The obtained standardized time series data set provides a prepared data foundation for subsequent data analysis and application.

[0099] See also Figure 3 , the steps to obtain the local feature vector are as follows:

[0100] The standardized time series data is divided into multiple subsequences according to fixed time windows. The change values ​​of gas flow and pressure are extracted for each subsequence, and the difference of parameters in each time window is calculated to obtain the time window change amplitude set.

[0101] Time windows are divided based on the standardized time series data. The division of time windows needs to be set in combination with the frequency of data collection and the dynamic characteristics of gas flow. Each time window should contain enough data points to ensure the effectiveness of statistical analysis. For example, a time window is selected to cover data collection points within 10 minutes to ensure that the changes in flow and pressure within this time period can be effectively captured. At the same time, the number of data points and window length in each window need to be adjusted according to the specific situation to adapt to different monitoring environments and equipment performance. Next, the maximum, minimum and average values ​​of gas flow and pressure in each time window are calculated. These statistical data will provide a basis for the next step of calculating the change amplitude and generate a time window change amplitude set.

[0102] Using the time window change amplitude set, the data of multiple time windows are input into the convolution kernel, the local response value of the convolution kernel to each window is calculated, and the weight and relative change ratio are combined for processing, using the formula:

[0103]

[0104] Extract local change features and generate a set of local feature values;

[0105] Among them, F iis the local eigenvalue of the i-th time window, W k is the kth convolution kernel weight, V i+k and V i are the gas parameter values ​​of the i+kth and ith time windows respectively, α is the enhancement index, γ is the variance adjustment coefficient, is the mean value of the gas parameter in all time windows, m is the total number of time windows, δ is the stabilization increment, η is the maximum value adjustment coefficient, and max(V) is the maximum value of the gas parameter in the time window;

[0106] formula:

[0107]

[0108] The benefit of the formula is that by adjusting the weight W k , the enhancement index α, the variance adjustment coefficient γ, the stabilization increment δ, and the maximum adjustment coefficient η can flexibly adjust the sensitivity to local fluctuations and the normalization process, thereby improving the accuracy and robustness of feature extraction.

[0109] Detailed explanation of the formula and the process of formula calculation and derivation:

[0110] Assume that within a specific time window, the observed values ​​of V are as follows: V i =100, V i+1 =105, V i+2 =110, weight W is 1, enhancement index α=2, variance adjustment coefficient γ=1, average value of all time windows The stabilization increment δ = 1, the maximum adjustment coefficient η = 0.1, the maximum value max(V) = 110, and the calculation process is as follows:

[0111]

[0112]

[0113] This result shows that the local eigenvalue F obtained by convolution kernel processing i It can effectively reflect the fluctuation of gas parameters within the time window, and this fluctuation feature will be used to evaluate the stability of gas flow and detect anomalies in subsequent analysis.

[0114] The local eigenvalue set is integrated into a feature vector according to the time window order, and normalized to generate a local feature vector.

[0115] The process of integrating the local eigenvalue set into the feature vector in the order of time windows needs to pay attention to the alignment and scaling of the data to ensure that the features in each time window are properly normalized to avoid the deviation of the comparison between different time windows due to different scales. The normalization of the feature vector can adopt the minimum-maximum normalization, that is, by subtracting the minimum value from the eigenvalue and dividing it by the difference between the maximum and minimum values. In this way, the processed feature vector will be between 0 and 1 and can be directly used for subsequent machine learning model training, thereby improving the model's adaptability to data and prediction accuracy, and generating local feature vectors.

[0116] See also Figure 4 , the specific steps for obtaining the timing dependency model are:

[0117] Based on the local feature vector, the time window is divided and the local feature values ​​of multiple windows are extracted. The trend change rate in each time window is calculated. The trend change rate and the time window sequence are combined to obtain the trend sequence set.

[0118] Based on the call of local feature vectors, this step starts with detailed time window division. First, the system collects gas flow and pressure data through a series of sensors. These data are obtained by real-time monitoring by the sensors in each predetermined time window. After that, they are aggregated through the data acquisition system to form the original data set. Next, the original data is preprocessed, including removing noise and outliers to ensure the accuracy and reliability of the data. Then, the data processing algorithm is used to calculate the trend change rate in each time window. This calculation process includes mathematical modeling and trend analysis of each window data. Trend analysis can identify the change pattern of each time window, laying the foundation for the creation of the next trend sequence set. The data processing process strictly follows the standards of data integrity and time series analysis to ensure the accuracy of the results and the effectiveness of implementation. The generated trend sequence set will serve as the basis for model establishment for subsequent dependency analysis and prediction model construction.

[0119] Based on the trend sequence set, the trend cumulative effect between time windows is analyzed, and the hidden state value is gradually transferred using the formula:

[0120]

[0121] Calculate the long-term dependencies between adjacent time windows and generate a sequence of hidden state values;

[0122] Among them, H t is the hidden state value of the t-th time window, φ i is the weight of the i-th time window, T t-i Represents the trend change value of the ti time window, T trepresents the trend change value of the t-th time window, ψ is the long-term dependence adjustment coefficient, ω k is the short-term adjustment coefficient, γ j is the variance adjustment factor, is the trend mean, δ is the stabilization parameter, and η is the offset parameter;

[0123] formula:

[0124]

[0125] The benefit of the formula is that by dynamically adjusting the weight parameter φ i , ψ, and ω k To adapt to the long-term and short-term changes of trends within the time window, it can simulate and predict long-term dependencies more accurately, and enhance the model's adaptability and prediction accuracy to dynamic changes in time series data.

[0126] Detailed explanation of the formula and the process of formula calculation and derivation:

[0127] Assume that in a specific time window t, a trend value T is observed t = 0.5 and the values ​​T of the previous three time windows t-1 =0.45, T t-2 =0.4, T t-3 =0.35. Select weight parameters φ1=0.8, φ2=0.15, φ3=0.05, long-term dependence adjustment coefficient ψ=1.2, short-term adjustment coefficient ω1=0.9. Variance adjustment factor γ1=γ2=γ3=1, trend mean The stabilization parameter δ=0.02 and the offset parameter η=0.1.

[0128]

[0129] This result shows that the hidden state value H calculated by the current model t When it is close to 1, it means that in a given time window, the long-term dependence is very strong and the trend change is very stable. The prediction model's prediction of the future state will be very close to the actual observation value, which improves the reliability and accuracy of the prediction.

[0130] Combined with the hidden state value sequence, the dependency relationship between time series is recursively calculated, the long-term dependency of multiple time windows is cumulatively analyzed, and the dependency relationship of multiple windows is integrated to generate a temporal dependency model.

[0131] When recursively calculating the dependencies between time series, the process first involves integrating the hidden state value sequence. In this process, the system analyzes the hidden state values ​​of each time window one by one through advanced algorithms to identify and integrate long-term trend dependencies. The algorithm takes into account the time sequence and data integrity of each state value to ensure the accuracy of the analysis. Through continuous data input, the model can more accurately capture the dependencies of time series and then effectively predict long-term trends. As enough time window data is analyzed, the model's predictive ability gradually increases. Through the cumulative analysis of these dependencies, a comprehensive time series dependency model can eventually be generated. This model can not only predict future states, but also provide a scientific basis for decision-making. The resulting time series dependency model provides strong support for real-time dynamic predictions, ensuring the efficiency and reliability of the model in practical applications.

[0132] See also Figure 5 The specific steps for obtaining the optimized hydrogen blending ratio prediction are as follows:

[0133] Using the time-series dependency model, input energy demand, carbon emissions and stability index parameters, calculate the expected value of hydrogen blending ratio in each window based on the time window sequence, analyze the constraints and integrate the results to generate a ratio constraint set;

[0134] When calling the time-dependent model, you first need to ensure that the input parameters include energy demand, carbon emissions and stability indicators. These data are obtained through actual energy consumption records, carbon emission reports and stability assessment reports to ensure the authenticity and accuracy of the data. After preprocessing, the input data is input into the model in the form of a time window series. The model processes these data to evaluate the expected value of the hydrogen blending ratio in each time window. These ratio values ​​are calculated based on current and expected energy demand, environmental protection standards and system stability requirements, ensuring that the predicted hydrogen blending ratio meets both energy demand and carbon emission standards without disrupting the stable operation of the system.

[0135] Based on the ratio constraint set, the hydrogen blending ratio is gradually adjusted and the constraints are applied, using the formula:

[0136]

[0137] Calculate and generate a corrected set of ratio values;

[0138] Among them, P t is the corrected hydrogen blending ratio in the tth time window, λ i is the weight of the i-th window, E t-i is the energy demand of the ti window, S t is the stability index of the tth window, α is the stability adjustment coefficient, C tis the carbon emission value of the tth window, γ is the carbon emission weight, β j is the variance adjustment factor, is the mean value of hydrogen blending ratio, δ is the stabilization parameter, ∈ is the energy demand adjustment coefficient, and max(E) is the maximum energy demand value;

[0139] formula:

[0140]

[0141] The benefit of the formula is that by integrating energy demand, stability indicators and carbon emission data, it provides a method for adjusting the hydrogen blending ratio that takes multiple factors into consideration, thereby increasing the accuracy and practicality of the model predictions.

[0142] Detailed explanation of the formula and the process of formula calculation and derivation:

[0143] Set energy requirement E t is 500, the stability index S t is 0.8, carbon emission C t is 300, weight λ i is 0.5, the adjustment coefficient α is 1.2, the carbon emission weight γ is 1.1, and the variance adjustment factor β j The average hydrogen ratio is 0.03. is 0.2, the stabilization parameter δ is 0.05, the energy demand adjustment coefficient ∈ is 0.02, and the maximum energy demand max(E) is 600. The calculation results are:

[0144]

[0145] The results show that under given parameters, the adjusted hydrogen blending ratio is 46.8%, which indicates that by applying this formula, the hydrogen blending ratio can be effectively adjusted according to current environmental and operational data, ensuring energy supply while reducing environmental impact.

[0146] Using the corrected ratio value set, the ratio value of each time window is recursively adjusted. The correction results are integrated with the time window constraints and energy, emission, and stability requirements to generate an optimized hydrogen blending ratio prediction.

[0147] Based on the corrected ratio value set, the ratio value of each time window is recursively adjusted. This process involves comparing the current ratio value of each time window with the ideal ratio value calculated by the model, and then gradually adjusting the actual hydrogen blending ratio according to the actual situation to approach the ideal ratio value. This method can ensure that the hydrogen blending ratio of each time window can be optimized and adjusted while meeting energy demand, reducing carbon emissions and ensuring system stability. The optimized hydrogen blending ratio prediction will be directly used to adjust the hydrogen blending ratio in actual operation.

[0148] See also Figure 6 The specific steps for obtaining the dynamically adjusted hydrogen blending ratio prediction are as follows:

[0149] Based on the optimized prediction of hydrogen blending ratio, real-time feedback data is collected and the deviation between the feedback data and the predicted value is calculated. The deviation is integrated according to the time window and the deviation change trend is calculated to generate a deviation trend set.

[0150] Monitor real-time feedback data, analyze the deviation between the predicted value and the actual data, and modify the model parameters according to the deviation trend. The real-time feedback data collected includes data streams obtained from sensors and monitoring systems. These data reflect the actual performance of the hydrogen blending ratio. Through the real-time data processing system, the deviation value is calculated at fixed time intervals. This process involves comparing the collected data with the model predicted value. The deviation analysis includes calculating the mean and standard deviation of the deviation, evaluating the persistence and volatility of the deviation to determine whether the model parameters need to be adjusted, and dynamically adjusting the model according to the deviation trend, specifically including the use of statistical and machine learning methods to identify the pattern and cause of the deviation. The process of modifying parameters is based on the nature of the deviation and the initial settings of the model, including adjusting the learning rate, retraining part or all of the model, or modifying the structure of the model to better adapt to data changes. The updated parameters will be used to optimize the prediction algorithm of the hydrogen blending ratio to ensure the accuracy and reliability of the prediction results.

[0151] Based on the deviation trend set, analyze the deviation change rules between time windows and correct the model parameters using the formula:

[0152]

[0153] Adjust parameter values ​​to generate a revised set of model parameters;

[0154] Among them, θ t is the model parameter adjusted in the tth time window, θ t-1 is the model parameter of the previous time window, λ i is the weight coefficient, D t-i is the deviation value of the ti-th window, Δ t is the cumulative deviation, φ is the deviation adjustment coefficient, β j is the variance adjustment factor, is the mean deviation, δ is the stabilization parameter, ∈ is the deviation amplitude adjustment coefficient, max(D) is the maximum deviation, n represents the number of previous time windows used for deviation calculation, and m represents the number of previous time windows used for variance calculation;

[0155] formula:

[0156]

[0157] The benefit of the formula is that it can dynamically adjust model parameters and respond to deviation changes in real time, thereby improving the accuracy and adaptability of predictions. Detailed explanation of the formula and the process of formula calculation and derivation: Use simulated data to illustrate. Assume that in five consecutive time windows, the deviation data D is 2, -1, 3, -2, 1, and the weight coefficient λ is i Take 0.1, 0.2, 0.3, 0.4, the adjustment coefficient φ is 0.5, and the variance adjustment factor β j Take 0.1, 0.2, 0.3, 0.4, the stabilization parameter δ is 1, the deviation amplitude adjustment coefficient ∈ is 0.1, and the calculation steps are:

[0158] Calculating weight adjustments =-0.6

[0159] 2. Calculate the deviation adjustment term φ·Δ5 and get -1.5

[0160] 3. Calculate the variance adjustment part in the denominator =2.9

[0161] 4. Maximum deviation max(D) is 3

[0162] 5. Substitute the above values ​​into the main formula to calculate

[0163]

[0164] This means that the value of θ5 is about 0.808 smaller than that of θ4. This result shows that the adjusted parameter θ5 can more accurately reflect the actual data changes and improve the response speed and accuracy of the prediction model.

[0165] The hydrogen blending ratio value is recalculated and adjusted using the revised model parameter set, and the revised ratio value is applied to the prediction result of the current time window to generate a dynamically adjusted hydrogen blending ratio prediction.

[0166] Using the revised model parameters, the hydrogen blending ratio is recalculated and adjusted to achieve optimal energy utilization and carbon emission efficiency by adjusting the ratio. This process includes using the new model parameters to re-evaluate the hydrogen blending ratio under current operating conditions, and updating the predicted value for each time window to ensure that each update is based on the latest model parameters and actual operating data, so as to achieve dynamic adjustment of the hydrogen blending ratio prediction. During the adjustment process, special attention is paid to the feedback of real-time data, as well as changes in environmental factors and market demand, to ensure that the adjustment of the hydrogen blending ratio can maximize the adaptation to changing operating needs and external conditions while ensuring safety and efficiency. The prediction results will be used for further energy management and scheduling strategy optimization.

[0167] A method for regulating hydrogen-blended natural gas based on data analysis, the method for regulating hydrogen-blended natural gas based on data analysis is performed based on the above-mentioned hydrogen-blended natural gas regulation system based on data analysis, and comprises the following steps:

[0168] S1: Collect the original hydrogen-blended natural gas process data, divide the time periods of gas flow, pressure, and temperature parameters in chronological order, calculate the average, maximum, and minimum values ​​of gas flow, pressure, and temperature in each time period, identify abnormal values ​​and perform interpolation replacement, calculate the change range of adjacent time periods and standardize the data to generate standardized time series data;

[0169] S2: Based on the standardized processing of time series data, the data is divided according to the set time window, and the data of each time window is sent to the convolution kernel for calculation, the short-term fluctuation and change amplitude of gas flow and pressure are extracted, and the local feature vector is constructed;

[0170] S3: Based on the local feature vector, analyze the trend in the time window, calculate the long-term dependency between time series, gradually transfer the hidden state, calculate the changes in adjacent time periods, capture the time series dependency, and generate a time series dependency model;

[0171] S4: Based on the time-series dependency model, the expected value of the hydrogen blending ratio is calculated, and the hydrogen blending ratio is adjusted in combination with the constraints including energy demand, carbon emissions, and stability. The constraints are applied to the predicted value for iterative correction to generate an optimized hydrogen blending ratio prediction;

[0172] S5: Based on the optimized hydrogen blending ratio prediction, real-time monitoring of the deviation between the feedback data and the predicted value, analysis of the deviation trend, correction of the prediction model parameters, recalculation of the hydrogen blending ratio, update of the adjusted parameters, and generation of a dynamically adjusted hydrogen blending ratio prediction;

[0173] S6: Based on the dynamically adjusted hydrogen blending ratio prediction, combined with real-time feedback data, difference analysis and trend adjustment are performed to monitor the real-time changes in the hydrogen blending ratio. The hydrogen blending ratio is iteratively optimized according to demand and stability requirements to generate an optimized value for the hydrogen blending ratio.

[0174] The above are only preferred embodiments of the present invention and are not intended to limit the present invention in other forms. Any technician familiar with the profession may use the technical contents disclosed above to change or modify them into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. A hydrogen-blended natural gas control system based on data analysis, characterized in that: The system comprises: The data preprocessing module collects the original hydrogen-blended natural gas process data, divides the gas flow, pressure, and temperature parameters in chronological order, calculates the average value and the maximum and minimum values ​​in each time period, identifies data anomalies and interpolates and replaces the abnormal values, calculates the change range of adjacent time periods and standardizes all data, performs time series interpolation to supplement missing data, and generates standardized time series data; The feature extraction module divides the sequence by time window based on the time series data after the standardization processing, inputs the data of multiple time windows into the convolution kernel for calculation, extracts the short-term fluctuation and change amplitude of the gas flow and pressure, completes the vector construction of the local feature of the time series, and normalizes the time window feature to generate a local feature vector; The time series modeling module analyzes the trend in the time window based on the local feature vector, calculates the long-term dependency between time series, gradually transfers the hidden state to calculate the changes in adjacent time periods, predicts the long-term trend, captures the dependency of time series through recursion, and generates a time series dependency model; The prediction optimization module calculates the expected value of the hydrogen blending ratio based on the time series dependency model, constrains the energy demand, carbon emissions, stability indicators and the hydrogen blending ratio, adjusts the ratio value and applies the constraint condition to the predicted value, performs iterative correction, and generates an optimized hydrogen blending ratio prediction; The dynamic adjustment module monitors the deviation between the real-time feedback data and the predicted value based on the optimized hydrogen blending ratio prediction, analyzes the deviation trend and corrects the prediction model parameters, recalculates the hydrogen blending ratio, updates the adjusted parameters, applies them to the current prediction results, and generates a dynamically adjusted hydrogen blending ratio prediction.

2. The hydrogen-blended natural gas control system based on data analysis according to claim 1 is characterized in that: The standardized time series data includes time series data, outlier processing results, and missing data filling results. The local feature vector includes gas flow characteristics, pressure change characteristics, and short-term fluctuation characteristics. The time series dependency model includes time period hidden state, long-term trend capture, and time series recursive relationship. The optimized hydrogen blending ratio prediction includes energy demand adjustment results, carbon emission optimization results, and stability balance results. The dynamically adjusted hydrogen blending ratio prediction includes feedback deviation correction results, adjusted parameters, and updated prediction results.

3. The hydrogen-blended natural gas control system based on data analysis according to claim 2 is characterized in that: The specific steps for obtaining the standardized time series data are as follows: Collect the original hydrogen-blended natural gas process data, and divide the gas flow, pressure and temperature parameters in the time series into multiple intervals according to time, calculate the average, maximum and minimum values ​​of the gas flow, pressure and temperature in each time interval, and obtain a statistical parameter set; Based on the statistical parameter set, an outlier threshold is set and data points exceeding the threshold are identified, and the outlier data points are replaced by a linear interpolation method to generate a statistical parameter set after interpolation processing; The statistical parameter set after the interpolation process is used to calculate the variation range of gas flow, pressure and temperature in adjacent time intervals, using the formula: Perform standardization processing and calculate and generate a standardized set of change amplitudes; Among them, S i represents the standardized change in the ith interval, X i+1 represents the gas parameter value in the i+1th time interval, X i represents the gas parameter value in the i-th time interval, is the mean of the gas parameters in all intervals, n is the total number of intervals, ω is the weighting factor, γ is the variance adjustment coefficient, and δ is the stabilization increment; In combination with the standardized set of variation amplitudes, a time series interpolation method is applied to complete the missing data, and the completed data is merged with the standardized variation amplitude to obtain a standardized time series data set.

4. The hydrogen-blended natural gas control system based on data analysis according to claim 3 is characterized in that: The steps of obtaining the local feature vector are specifically as follows: Dividing the standardized time series data into multiple subsequences according to fixed time windows, extracting the change values ​​of gas flow and pressure for each subsequence, and obtaining a time window change amplitude set by calculating the difference of the parameters in each time window; Using the time window change amplitude set, the data of multiple time windows are input into the convolution kernel, the local response value of the convolution kernel to each window is calculated, and the weight and relative change ratio are combined for processing, using the formula: Extract local change features and generate a set of local feature values; Among them, F i is the local eigenvalue of the i-th time window, W k is the kth convolution kernel weight, V i+k and V i are the gas parameter values ​​of the i+kth and ith time windows respectively, α is the enhancement index, γ is the variance adjustment coefficient, is the mean value of the gas parameter in all time windows, m is the total number of time windows, δ is the stabilization increment, η is the maximum value adjustment coefficient, and max(V) is the maximum value of the gas parameter in the time window; The local feature value set is integrated into a feature vector according to the time window sequence, and normalized to generate a local feature vector.

5. The hydrogen-blended natural gas control system based on data analysis according to claim 4 is characterized in that: The steps for obtaining the timing dependency model are specifically as follows: Based on the local feature vector, the time window is divided and the local feature values ​​of multiple windows are extracted, the trend change rate in each time window is calculated, and the trend change rate and the time window sequence are combined to obtain a trend sequence set; Based on the trend sequence set, the trend cumulative effect between time windows is analyzed, and the hidden state value is gradually transferred using the formula: Calculate the long-term dependencies between adjacent time windows and generate a sequence of hidden state values; Among them, H t is the hidden state value of the t-th time window, φ i is the weight of the i-th time window, T t-i Represents the trend change value of the ti time window, T t represents the trend change value of the t-th time window, ψ is the long-term dependence adjustment coefficient, ω k is the short-term adjustment coefficient, γ j is the variance adjustment factor, is the trend mean, δ is the stabilization parameter, and η is the offset parameter; Combined with the hidden state value sequence, the dependency relationship between time series is recursively calculated, the long-term dependency of multiple time windows is cumulatively analyzed, and the dependency relationship of multiple windows is integrated to generate a temporal dependency model.

6. The hydrogen-blended natural gas control system based on data analysis according to claim 5 is characterized in that: The steps for obtaining the optimized hydrogen blending ratio prediction are specifically as follows: Using the time series dependency model, inputting energy demand, carbon emission and stability index parameters, calculating the expected value of hydrogen blending ratio in each window in combination with the time window sequence, analyzing the constraint conditions and integrating the results, and generating a ratio constraint set; Based on the ratio constraint set, the hydrogen blending ratio is gradually adjusted and the constraint conditions are applied, using the formula: Calculate and generate a corrected set of ratio values; Among them, P t is the corrected hydrogen blending ratio in the tth time window, λ i is the weight of the i-th window, E t-i is the energy demand of the ti window, S t is the stability index of the tth window, α is the stability adjustment coefficient, C t is the carbon emission value of the tth window, γ is the carbon emission weight, β j is the variance adjustment factor, is the mean value of hydrogen blending ratio, δ is the stabilization parameter, ∈ is the energy demand adjustment coefficient, and max(E) is the maximum energy demand value; The corrected ratio value set is used to recursively adjust the ratio value of each time window, and the correction results are integrated with the time window constraints and energy, emission, and stability requirements to generate an optimized hydrogen blending ratio prediction.

7. The hydrogen-blended natural gas control system based on data analysis according to claim 6 is characterized in that: The steps for obtaining the dynamically adjusted hydrogen blending ratio prediction are specifically as follows: Based on the optimized prediction of hydrogen blending ratio, real-time feedback data is collected and the deviation between the feedback data and the predicted value is calculated, the deviation is integrated according to the time window and then the deviation change trend is calculated to generate a deviation trend set; Based on the deviation trend set, the deviation variation law between time windows is analyzed and the model parameters are corrected using the formula: Adjust parameter values ​​to generate a revised set of model parameters; Among them, θ t is the model parameter adjusted in the tth time window, θ t-1 is the model parameter of the previous time window, λ i is the weight coefficient, D t-i is the deviation value of the ti-th window, Δ t is the cumulative deviation, φ is the deviation adjustment coefficient, β j is the variance adjustment factor, is the mean deviation, δ is the stabilization parameter, ∈ is the deviation amplitude adjustment coefficient, max(D) is the maximum deviation, n represents the number of previous time windows used for deviation calculation, and m represents the number of previous time windows used for variance calculation; The hydrogen blending ratio value is recalculated and adjusted using the revised model parameter set, and the revised ratio value is applied to the prediction result of the current time window to generate a dynamically adjusted hydrogen blending ratio prediction.

8. A method for regulating hydrogen-blended natural gas based on data analysis, characterized in that: The hydrogen-blended natural gas control system based on data analysis according to any one of claims 1 to 7 comprises the following steps: Collect the original hydrogen-blended natural gas process data, divide the time periods of gas flow, pressure, and temperature parameters in chronological order, calculate the average, maximum, and minimum values ​​of gas flow, pressure, and temperature in each time period, identify abnormal values ​​and perform interpolation replacement, calculate the change range of adjacent time periods and standardize the data to generate standardized time series data; Based on the standardized processing time series data, the data is divided according to the set time window, the data of each time window is sent to the convolution kernel for calculation, the short-term fluctuation and change amplitude of the gas flow and pressure are extracted, and the local feature vector is constructed; Based on the local feature vector, the trend in the time window is analyzed, the long-term dependency between time series is calculated, the hidden state is gradually transferred, the changes in adjacent time periods are calculated, the time series dependency is captured, and a time series dependency model is generated; Based on the time series dependency model, the expected value of the hydrogen blending ratio is calculated, and the hydrogen blending ratio is adjusted in combination with the constraints including energy demand, carbon emissions, and stability, and the constraints are applied to the predicted value for iterative correction to generate an optimized hydrogen blending ratio prediction; Based on the optimized hydrogen blending ratio prediction, real-time monitoring of the deviation between feedback data and the predicted value, analyzing the deviation trend, correcting the prediction model parameters, recalculating the hydrogen blending ratio, updating the adjusted parameters, and generating a dynamically adjusted hydrogen blending ratio prediction; Based on the dynamically adjusted hydrogen blending ratio prediction, combined with real-time feedback data, difference analysis and trend adjustment are performed to monitor real-time changes in the hydrogen blending ratio, and the hydrogen blending ratio is iteratively optimized according to demand and stability requirements to generate an optimized hydrogen blending ratio value.

Images