Disclosure of Invention
In view of the above, the invention aims to provide a ground source heat pump system simulation method and device based on a shallow soil source g-DTM model, so as to solve the technical problem that the ground source heat pump system is difficult to optimize due to the fact that the ground source heat pump system is difficult to simulate.
In a first aspect, an embodiment of the present invention provides a ground source heat pump system simulation method based on a shallow soil source g-DTM model, including:
An acquisition step of acquiring input parameters, wherein the input parameters comprise: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, ground source heat pump unit configuration parameters and ground source heat pump system quasi-operation parameters; iterative step, carrying out iterative calculation based on the input parameters and a target model, and determining current temperature data of the shallow soil source meeting a first preset iteration termination condition, wherein the target model comprises: a shallow soil source g-DTM model and a shallow soil source response model, the current temperature data comprising: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model for determining current temperature data under the conditions of specific mass flow and heat extraction quantity; executing, namely obtaining target parameters based on the current temperature data of the shallow soil source, the simulated operation parameters of the ground source heat pump system and a ground source heat pump unit model, wherein the target parameters comprise: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit; determining, if the target parameter does not meet a second preset iteration termination condition, repeating the iteration step and the execution step by using the current temperature data and the target parameter until the target parameter meets the second preset iteration termination condition, and determining the target parameter meeting the second preset iteration termination condition as a simulation result; performing iterative computation based on the input parameters and a target model, determining current temperature data of the shallow soil source meeting a first preset iteration termination condition, including: calculating rated total power of the ground source heat pump unit and initial heat extraction quantity of the shallow soil source based on the configuration parameters of the ground source heat pump unit; and carrying out iterative computation by using the initial heat extraction quantity, the input parameters, the target model and a numerical gradient descent algorithm, and determining the corresponding current temperature data of the shallow soil source meeting the first preset iteration termination condition.
Further, performing iterative computation by using the initial heat extraction amount, the input parameter, the target model and a numerical gradient descent algorithm, and determining the corresponding current temperature data of the shallow soil source meeting the first preset iteration termination condition, where the current temperature data includes: determining initial temperature data of the shallow soil source based on the input parameters of the shallow soil source and the shallow soil source g-DTM model; performing iterative computation based on the initial temperature data, the initial heat extraction amount, the input parameters of the shallow soil source, the shallow soil source response model and a numerical gradient descent algorithm, and calculating an absolute value of a numerical gradient corresponding to the current iteration; if the absolute value of the numerical gradient corresponding to the current iteration is larger than or equal to the first convergence value, calculating the step length corresponding to the numerical gradient corresponding to the current iteration, and calculating the inlet water temperature corresponding to the current iteration by utilizing the step length; calculating an absolute value of a difference value between an objective function value corresponding to the current iteration and an objective function value corresponding to a previous iteration of the current iteration, or calculating an absolute value of a difference value between inlet water temperature corresponding to the current iteration and inlet water temperature corresponding to a previous iteration of the current iteration, wherein the objective function value is an absolute value of a difference value between an heat extraction amount corresponding to the current iteration and the initial heat extraction amount; if the absolute value of the difference value between the objective function value corresponding to the current iteration and the objective function value corresponding to the previous iteration of the current iteration is smaller than or equal to a second convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source, or if the absolute value of the difference value between the inlet water temperature corresponding to the current iteration and the inlet water temperature corresponding to the previous iteration of the current iteration is smaller than or equal to a third convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source; if the absolute value of the numerical gradient corresponding to the current iteration is smaller than the first convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source; and calculating the current outlet water temperature of the shallow soil source based on the current inlet water temperature of the shallow soil source and the initial heat extraction amount.
Further, the second preset iteration termination condition is that an absolute value of a difference between the heat extraction amount corresponding to the current iteration and the heat extraction amount corresponding to the previous iteration of the current iteration is smaller than or equal to a fourth convergence value, and an absolute value of a difference between the total power corresponding to the current iteration and the total power corresponding to the previous iteration of the current iteration is smaller than or equal to a fifth convergence value.
In a second aspect, an embodiment of the present invention further provides a ground source heat pump system simulation device based on a shallow soil source g-DTM model, including: an obtaining unit, configured to obtain an input parameter, where the input parameter includes: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, ground source heat pump unit configuration parameters and ground source heat pump system quasi-operation parameters; the iteration unit is used for carrying out iterative computation based on the input parameters and a target model, and determining the current temperature data of the shallow soil source meeting a first preset iteration termination condition, wherein the target model comprises: a shallow soil source g-DTM model and a shallow soil source response model, the current temperature data comprising: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model for determining current temperature data under the conditions of specific mass flow and heat extraction quantity; the execution unit is used for obtaining target parameters based on the current temperature data of the shallow soil source, the quasi-operation parameters of the ground source heat pump system and a ground source heat pump unit model, wherein the target parameters comprise: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit; the determining unit is used for controlling the iteration unit and the execution unit to work repeatedly by utilizing the current temperature data and the target parameter if the target parameter does not meet a second preset iteration termination condition until the target parameter meets the second preset iteration termination condition, and determining the target parameter meeting the second preset iteration termination condition as a simulation result; wherein, the iteration unit is used for: calculating rated total power of the ground source heat pump unit and initial heat extraction quantity of the shallow soil source based on the configuration parameters of the ground source heat pump unit; and carrying out iterative computation by using the initial heat extraction quantity, the input parameters, the target model and a numerical gradient descent algorithm, and determining the corresponding current temperature data of the shallow soil source meeting the first preset iteration termination condition.
Further, the iteration unit is used for determining initial temperature data of the shallow soil source based on the input parameters of the shallow soil source and the shallow soil source g-DTM model; performing iterative computation based on the initial temperature data, the initial heat extraction amount, the input parameters of the shallow soil source, the shallow soil source response model and a numerical gradient descent algorithm, and calculating an absolute value of a numerical gradient corresponding to the current iteration; if the absolute value of the numerical gradient corresponding to the current iteration is larger than or equal to the first convergence value, calculating the step length corresponding to the numerical gradient corresponding to the current iteration, and calculating the inlet water temperature corresponding to the current iteration by utilizing the step length; calculating an absolute value of a difference value between an objective function value corresponding to the current iteration and an objective function value corresponding to a previous iteration of the current iteration, or calculating an absolute value of a difference value between inlet water temperature corresponding to the current iteration and inlet water temperature corresponding to a previous iteration of the current iteration, wherein the objective function value is an absolute value of a difference value between an heat extraction amount corresponding to the current iteration and the initial heat extraction amount; if the absolute value of the difference value between the objective function value corresponding to the current iteration and the objective function value corresponding to the previous iteration of the current iteration is smaller than or equal to a second convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source, or if the absolute value of the difference value between the inlet water temperature corresponding to the current iteration and the inlet water temperature corresponding to the previous iteration of the current iteration is smaller than or equal to a third convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source; if the absolute value of the numerical gradient corresponding to the current iteration is smaller than the first convergence value, determining the inlet water temperature corresponding to the current iteration as the current inlet water temperature of the shallow soil source; and calculating the current outlet water temperature of the shallow soil source based on the current inlet water temperature of the shallow soil source and the initial heat extraction amount.
Further, the second preset iteration termination condition is that an absolute value of a difference between the heat extraction amount corresponding to the current iteration and the heat extraction amount corresponding to the previous iteration of the current iteration is smaller than or equal to a fourth convergence value, and an absolute value of a difference between the total power corresponding to the current iteration and the total power corresponding to the previous iteration of the current iteration is smaller than or equal to a fifth convergence value.
In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is configured to store a program for supporting the processor to execute the method described in the first aspect, and the processor is configured to execute the program stored in the memory.
In a fourth aspect, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon.
In the embodiment of the present invention, through the obtaining step, input parameters are obtained, where the input parameters include: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, configuration parameters of a ground source heat pump unit and quasi-operation parameters of a ground source heat pump system; iterative step, carrying out iterative calculation based on input parameters and a target model, and determining current temperature data of a shallow soil source meeting a first preset iteration termination condition, wherein the target model comprises: shallow soil source g-DTM model and shallow soil source response model, current temperature data includes: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model used for determining current temperature data under the conditions of specific mass flow and heat extraction; the method comprises the steps of executing, based on current temperature data of a shallow soil source, obtaining target parameters by a ground source heat pump system quasi-operation parameter and a ground source heat pump unit model, wherein the target parameters comprise: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit; determining, if the target parameter does not meet the second preset iteration termination condition, repeating the iteration step and the execution step by using the current temperature data and the target parameter until the target parameter meets the second preset iteration termination condition, determining the target parameter meeting the second preset iteration termination condition as a simulation result, and realizing more accurate simulation of the ground source heat pump system, thereby solving the technical problem that the ground source heat pump system is difficult to simulate and further difficult to perform accurate optimization control on the ground source heat pump system, and laying a foundation for design optimization, performance prediction or evaluation, operation decision control, fault diagnosis and the like of the ground source heat pump system; meanwhile, the efficiency and the reliability of the ground source heat pump system can be improved by matching with other optimization algorithms, and the technical effects of energy consumption and operation cost are reduced.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the above objects, features and advantages of the present invention more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Embodiment one:
According to an embodiment of the present invention, there is provided an embodiment of a ground source heat pump system simulation method based on a shallow soil source g-DTM model, it being noted that the steps shown in the flowchart of the drawings may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is shown in the flowchart, in some cases the steps shown or described may be performed in an order different from that herein.
Fig. 1 is a flowchart of a simulation method of a ground source heat pump system based on a shallow soil source g-DTM model according to an embodiment of the present invention, as shown in fig. 1, the method includes the following steps:
Step S102, an acquisition step, namely acquiring input parameters, wherein the input parameters comprise: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, ground source heat pump unit configuration parameters and ground source heat pump system quasi-operation parameters;
step S104, an iteration step, namely carrying out iterative computation based on the input parameters and a target model, and determining current temperature data of the shallow soil source meeting a first preset iteration termination condition, wherein the target model comprises: a shallow soil source g-DTM model and a shallow soil source response model, the current temperature data comprising: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model for determining current temperature data under the conditions of specific mass flow and heat extraction quantity;
Step S106, executing the step, and obtaining target parameters based on the current temperature data of the shallow soil source, the to-be-operated parameters of the ground source heat pump system and a ground source heat pump unit model, wherein the target parameters comprise: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit;
And step S108, determining, namely repeatedly executing the iteration step and the execution step by using the current temperature data and the target parameter if the target parameter does not meet a second preset iteration termination condition until the target parameter meets the second preset iteration termination condition, and determining the target parameter meeting the second preset iteration termination condition as a simulation result.
In the embodiment of the present invention, through the obtaining step, input parameters are obtained, where the input parameters include: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, configuration parameters of a ground source heat pump unit and quasi-operation parameters of a ground source heat pump system; iterative step, carrying out iterative calculation based on input parameters and a target model, and determining current temperature data of a shallow soil source meeting a first preset iteration termination condition, wherein the target model comprises: shallow soil source g-DTM model and shallow soil source response model, current temperature data includes: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model used for determining current temperature data under the conditions of specific mass flow and heat extraction; the method comprises the steps of executing, based on current temperature data of a shallow soil source, obtaining target parameters by a ground source heat pump system quasi-operation parameter and a ground source heat pump unit model, wherein the target parameters comprise: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit; determining, if the target parameter does not meet the second preset iteration termination condition, repeating the iteration step and the execution step by using the current temperature data and the target parameter until the target parameter meets the second preset iteration termination condition, determining the target parameter meeting the second preset iteration termination condition as a simulation result, and realizing more accurate simulation of the ground source heat pump system, thereby solving the technical problem that the ground source heat pump system is difficult to simulate and further difficult to perform accurate optimization control on the ground source heat pump system, and laying a foundation for design optimization, performance prediction or evaluation, operation decision control, fault diagnosis and the like of the ground source heat pump system; meanwhile, the efficiency and the reliability of the ground source heat pump system can be improved by matching with other optimization algorithms, and the technical effects of energy consumption and operation cost are reduced.
Configuration parameters of shallow soil source: specific heat capacity c p, drilling depth H, internal thermal resistance R b of drilling and comprehensive heat conductivity coefficient of rock-soil medium of fluid working mediumZero point T 0 of excess temperature (rock-soil medium surface temperature), within the current time interval, under the action of unit rectangular pulse heat flow, the dimensionless temperature response G of the drilling holes, the minimum value T 1,min,summer of the inlet water temperature of the shallow soil source in summer, the maximum value T 1,max,summer of the inlet water temperature of the shallow soil source in summer, the minimum value T 1,min,winter of the inlet water temperature of the shallow soil source in winter and the maximum value T 1,max,winter of the inlet water temperature of the shallow soil source in winter;
configuration parameters of the ground source heat pump unit: summer performance parameters a i, i epsilon [0,9], rated cold load Q GSHP,rated,cooling, rated refrigeration power P rated,cooling, winter performance parameters b i, i epsilon [0,9], rated cold load Q GSHP,rated,heating, rated heating power P rated,cooling and total number n GSHP;
shallow soil source side historical operating data: mass flow rate of borehole j ground heat exchanger in kth time interval In the kth time interval, inlet temperature/>, of borehole j borehole heat exchangerIn the kth time interval, outlet temperature/>, of borehole j borehole heat exchangerA time discrete step delta tau;
the ground source heat pump system is currently intended to operate and set parameters: the total mass flow M of the shallow soil source, the running load Q of the ground source heat pump system, the water temperature T in,load at the inlet of the load side and the water temperature T out,load at the outlet of the load side.
According to the embodiment of the invention, the ground source heat pump system is provided with the buried pipe heat exchanger in the shallow soil, so that heat in the shallow soil is extracted in winter, and indoor heat is released into the shallow soil in summer, thereby realizing cooling and heating of the ground source heat pump system.
The shallow soil source g-DTM model is a g-DTM (g Function Based DISCRETE TRANSFER Matrix) model Based on a Green Function method and a linear superposition principle, updates the temperature state of the shallow soil source through historical operation data (namely historical heat flow sequences of all drilling holes of the shallow soil source) at the ground source side, rapidly calculates the water outlet temperature response of the shallow soil source under the specific water inlet temperature of specific mass flow, and has the expression For drilling depth,Is the comprehensive heat conductivity coefficient of the rock-soil medium, c p is the specific heat capacity of the fluid working medium, R b is the internal thermal resistance of the drill hole, T 0 is the zero point of excessive temperature, the surface temperature of the rock-soil medium, ζ is an intermediate variable, γ is an intermediate variable, I is the flow, T 2 is the current moment, the water outlet temperature of the shallow soil source, T 1 is the inlet temperature column vector of each drill hole buried pipe heat exchanger in the current time interval,T i,1 is the inlet temperature of the borehole i buried pipe heat exchanger in the current time interval, DEG C, T 1 is the water inlet temperature of the shallow soil source at the current time, N is the number of boreholes, T 2 is the inlet temperature column vector of each borehole buried pipe heat exchanger in the current time interval,T i,2 is the outlet temperature of the borehole i borehole heat exchanger during the current time interval,MK is the mass flow diagonal matrix of each borehole buried pipe heat exchanger in the current time interval, M i is the mass flow of each borehole i buried pipe heat exchanger in the current time interval, the mass flow of each borehole is the mass flow of the total mass flow M of the shallow soil source averaged to each borehole, G is the dimensionless temperature response of each borehole under the action of unit rectangular pulse heat flow in the current time interval, G i,j is the dimensionless temperature response of the j-th borehole buried pipe heat exchanger at the outer wall surface of the i-th borehole under the action of unit rectangular pulse heat flow in the current time interval, S is the shallow soil source temperature state in the current time interval, namely the excessive temperature of the middle position of each borehole wall surface of the shallow soil source caused by the stepped heat flow of each borehole buried pipe heat exchanger based on the historical moment,S i is the excess temperature at the borehole wall of the i-th borehole caused by the step heat flow of each borehole heat exchanger at the historical time,K is the current moment, S i is the excess temperature at the wall of the drill hole of the ith drill hole caused by the step heat flow of each drill hole buried pipe heat exchanger group at the historical moment,For the k-th time interval, mass flow rate of borehole j ground borehole heat exchanger,For the kth time interval, the inlet temperature of the borehole j borehole heat exchanger,For the k-th time interval, the outlet temperature of the borehole j ground borehole heat exchanger,In order to realize the dimensionless temperature response of the j-th borehole buried heat exchanger at the outer wall surface of the i-th borehole under the action of unit rectangular pulse heat flow in the time interval ([ (k-1) delta tau, k delta tau ]).
The shallow soil source response model is described below.
Under the condition of specific mass flow M and heat extraction quantity Q soil,target, the shallow soil source response model determines the inlet water temperature T 1 and the outlet water temperature T 2 of the shallow soil source through a numerical gradient descent method.
The calculation formula of the heat extraction amount of the shallow soil source is as follows:
Qsoil=cpM(T2-T1)
wherein:
q soil -shallow soil source heat extraction quantity, kW;
c p -specific heat capacity of fluid working medium, kJ/(kg.K);
M, total mass flow of shallow soil source inlet, kg/s;
T 2 -shallow soil source water outlet temperature, DEG C;
T 1 -shallow soil source water inlet temperature, DEG C;
The calculation formula of the heat extraction quantity of the shallow soil source is combined with a shallow soil source model, and under the condition that the specific mass flow M and the shallow soil source temperature state S are known, namely the inlet water temperature T 1 of the shallow soil source is given, the outlet water temperature T 2 of the shallow soil source can be obtained, and then the heat extraction quantity Q soil of the shallow soil source can be determined, namely the functional relation Q soil=h(T1 exists.
1) Based on the heat extraction quantity Q soil,target, judging the shallow soil source operation mode:
2) Constructing an objective function f (T 1)=|h(T1)-Qsoil,target |
3) T=0, initialize
Wherein:
t 1,min, the minimum inlet water temperature of the shallow soil source, and the temperature of the shallow soil source;
T 1,max, the maximum value of the inlet water temperature of the shallow soil source, and the temperature is lower than the temperature;
T 1,min,summer, the minimum inlet water temperature of the shallow soil source in summer, and the temperature of the shallow soil source in summer is lower than the temperature;
t 1,max,summer, the maximum value of the water temperature imported from the shallow soil source in summer, and the temperature is lower than the temperature;
T 1,min,winter, the minimum water temperature of the inlet of the shallow soil source in winter, and the temperature is lower than the temperature;
T 1,max,winter, the maximum value of the inlet water temperature of the shallow soil source in winter, and the temperature is lower than the temperature;
4) For small increment DeltaT 1, calculating the numerical gradient When the distance of the numerical gradientStopping iteration, letEpsilon 1 is the first convergence value,The small increment delta T 1 is the inlet water temperature of the shallow soil source at the T moment and can be set by a worker according to the actual situation;
otherwise, the step size lambda t,λt needs to satisfy the following formula:
5) Order the And calculateWhen Stopping iteration, letTurning (7), ε 2 is the second convergence value and ε 3 is the third convergence value;
6) Let t=t+1 and return to (4);
7) Outputting shallow soil source inlet temperature Shallow soil Source Outlet temperature
The ground source heat pump unit model will be described below.
In summer, the ground source heat pump unit power P cooling is described as a function of the load side outlet water temperature T load,out, the ground source side inlet water temperature T source,in and the unit load rate PLR:
wherein:
P cooling -ground source heat pump unit power, kW;
a i, i epsilon [0,9] -summer performance parameters of the ground source heat pump unit;
T load,out -the outlet water temperature of the load side of the ground source heat pump unit, namely the outlet water temperature of chilled water, at DEG C;
T source,in -the ground source side water inlet temperature of the ground source heat pump unit, namely the return water temperature of cooling water, at DEG C;
PLR-unit load factor, ratio of ground source heat pump unit operating load Q GSHP,cooling to ground source heat pump unit rated cooling load Q GSHP,rated,cooling;
P rated,cooling -rated refrigeration power of ground source heat pump unit, kW;
In winter, the ground source heat pump unit power P heating is described as a function of the load side water inlet temperature T load,in, the ground source side water inlet temperature T source,out, and the unit load rate PLR:
wherein:
p heating -ground source heat pump unit power, kW;
b i, i epsilon [0,9] -a winter performance parameter of the ground source heat pump unit;
T load,in -the water inlet temperature of the load side of the ground source heat pump unit, namely the water return temperature of hot water, at DEG C;
T source,out -the ground source side water outlet temperature of the ground source heat pump unit, namely the shallow soil source water inlet temperature, DEG C;
PLR-unit load factor, ratio of ground source heat pump unit operating load Q GSHP,heating to ground source heat pump unit rated heat load Q GSHP,rated,heating;
the calculation formula of the running number of the ground source heat pump unit in the cooling season is as follows:
wherein:
q cooling -ground source heat pump system cold load, kW;
q GSHP,rated,cooling -rated cold load of the ground source heat pump unit, kW;
n GSHP -total number of ground source heat pump units;
n GSHP,operation -the running number of the ground source heat pump units;
i-i ground source heat pump units.
In the heating season, the calculation formula of the running number of the ground source heat pump units is shown as follows.
Wherein:
Q heating -the heat load of the ground source heat pump system, kW;
Q GSHP,rated,heating -rated heat load of the ground source heat pump unit, kW;
n GSHP -total number of ground source heat pump units;
n GSHP,operation -the running number of the ground source heat pump units;
i-i ground source heat pump units.
Embodiments of the present invention will be described below with reference to the above models.
And calculating the heat extraction quantity of the shallow soil source under the current operation working condition based on the rated working condition performance of the ground source heat pump unit.
All the ground source heat pump units of the ground source heat pump system have the same model number and performance parameters, and the running ground source heat pump units are uniformly distributed with the current running load. When the cycle number t=0, the load of the ground source heat pump system in the cooling season is Q cooling =q, the load of the ground source heat pump system in the heating season is Q heating =q, if the inlet and outlet water temperatures of the shallow soil source are unknown, the power of the ground source heat pump unit is determined by utilizing the inlet and outlet water temperatures of the design of the shallow soil source in the heating season or the cooling season, and the total power P GSHP of the ground source heat pump unit and the heat extraction quantity Q soil of the shallow soil source are shown in the following calculation formulas:
and (5) cooling season:
Heating season:
determining specific heat extraction amount based on shallow soil source response model Inlet water temperature/>, of lower shallow soil sourceAnd outlet water temperature
Inlet water temperature of shallow soil sourceAnd outlet water temperatureCalculating total power/>, of the ground source heat pump unit at the time t+1 by combining current operation condition parametersAnd determining the heat extraction quantity/>, of the shallow soil source at the t+1 time
Inlet water temperature obtained through shallow soil source response modelAnd outlet water temperatureTotal power/>, of ground source heat pump unitHeat extraction amount of shallow soil sourceThe calculation formula is as follows:
and (5) cooling season:
Heating season:
Let t=t+1;
Iteration termination criteria: if it is AndMeet the convergence criterion, letOtherwise, the convergence criterion is not satisfied, ε 4 is a fourth convergence value, ε 5 is a fifth convergence value.
Outputting the performance parameters of the ground source heat pump system, and outputting the inlet temperature of the shallow soil sourceOutlet temperatureHeat taking amountTotal power
According to the ground source heat pump system high-performance simulation method based on the shallow soil source g-DTM model, the shallow soil source response model and the ground source heat pump unit model are coupled and considered, and the performance of the ground source heat pump system is predicted by utilizing the technology of iterative successive approximation of the fixed point. Not only lays a foundation for design optimization, performance prediction or evaluation, operation decision control, fault diagnosis and the like of the ground source heat pump system; meanwhile, the efficiency and the reliability of the ground source heat pump system can be improved by matching with other optimization algorithms, and the energy consumption and the running cost are reduced.
Embodiment two:
The embodiment of the invention also provides a ground source heat pump system simulation device based on the shallow soil source g-DTM model, which is used for carrying out the ground source heat pump system simulation method based on the shallow soil source g-DTM model provided by the embodiment of the invention, and the following is a specific introduction of the device provided by the embodiment of the invention.
As shown in fig. 2, fig. 2 is a schematic diagram of the above ground source heat pump system simulation device based on the shallow soil source g-DTM model, where the device includes:
An obtaining unit 10, configured to obtain input parameters, where the input parameters include: configuration parameters of a shallow soil source, historical operation data of the shallow soil source, ground source heat pump unit configuration parameters and ground source heat pump system quasi-operation parameters;
The iteration unit 20 is configured to perform iterative computation based on the input parameters and a target model, and determine current temperature data of the shallow soil source that meets a first preset iteration termination condition, where the target model includes: a shallow soil source g-DTM model and a shallow soil source response model, the current temperature data comprising: the method comprises the steps that the current inlet water temperature of a shallow soil source and the current outlet water temperature of the shallow soil source are obtained, a shallow soil source g-DTM model is a model constructed based on a Green function method and a linear superposition principle, and a shallow soil source response model is a model for determining current temperature data under the conditions of specific mass flow and heat extraction quantity;
the execution unit 30 is configured to obtain target parameters based on the current temperature data of the shallow soil source, the to-be-operated parameter of the ground source heat pump system, and a ground source heat pump unit model, where the target parameters include: the current heat extraction amount of the shallow soil source and the current total power of the ground source heat pump unit;
And the determining unit 40 is configured to control the iteration unit and the execution unit to work repeatedly by using the current temperature data and the target parameter if the target parameter does not meet a second preset iteration termination condition, until the target parameter meets the second preset iteration termination condition, and determine the target parameter meeting the second preset iteration termination condition as a simulation result.
Embodiment III:
An embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is configured to store a program that supports the processor to execute the method described in the first embodiment, and the processor is configured to execute the program stored in the memory.
Embodiment four:
The embodiment of the invention also provides a computer readable storage medium, and a computer program is stored on the computer readable storage medium, and when the computer program is executed by a processor, the steps of the method in the first embodiment are executed.
In addition, in the description of embodiments of the present invention, unless explicitly stated and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
In the description of the present invention, it should be noted that the directions or positional relationships indicated by the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. are based on the directions or positional relationships shown in the drawings, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the several embodiments provided by the present application, it should be understood that the disclosed systems, devices, and methods may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
Finally, it should be noted that: the above examples are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, but it should be understood by those skilled in the art that the present invention is not limited thereto, and that the present invention is described in detail with reference to the foregoing examples: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.