CN109255158A - Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium - Google Patents

Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium Download PDF

Info

Publication number
CN109255158A
CN109255158A CN201810929370.1A CN201810929370A CN109255158A CN 109255158 A CN109255158 A CN 109255158A CN 201810929370 A CN201810929370 A CN 201810929370A CN 109255158 A CN109255158 A CN 109255158A
Authority
CN
China
Prior art keywords
water discharge
warm water
mathematical model
waters
boundary condition
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810929370.1A
Other languages
Chinese (zh)
Other versions
CN109255158B (en
Inventor
陈小莉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201810929370.1A priority Critical patent/CN109255158B/en
Publication of CN109255158A publication Critical patent/CN109255158A/en
Application granted granted Critical
Publication of CN109255158B publication Critical patent/CN109255158B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

This application involves a kind of warm water discharge prototype measurement waters background water temperature calculation method, device, computer equipment and storage mediums.Method includes: to utilize the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, hydrologic parameter, warm water discharge water intaking and draining source sink term and waters model meshes, determine the first input boundary condition of the mathematical model after engineering operation during warm water discharge prototype measurement, according to the first input boundary condition, actual measurement hydrodynamic force and water temperature result after engineering operation during warm water discharge prototype measurement, calibration is carried out to mathematical model, obtain mathematical model Reasonable Parameters value, according to second of the mathematical model during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of mathematical model the input boundary condition, calculate the background water temperature in waters to be detected, the available more accurate reasonable background water temperature after calibration.

Description

Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium
Technical field
This application involves waters warm water discharge detection technique fields, more particularly to a kind of warm water discharge prototype measurement waters background Water temperature calculates method, apparatus, computer equipment and storage medium.
Background technique
With the development of the times, influence of the warm water discharge to water environment is gradually concerned by people, and commonly uses remote sensing technique Temperature field observation is carried out to warm water discharge discharge receiving waters, and then determines temperature rise area coverage caused by warm water discharge.And according to Remote sensing technique determines that an important link of temperature rise area coverage caused by warm water discharge is to determine without the water under the influence of warm water discharge Domain background water temperature.
It is previous to determine that one is artificial pre-segmentations there are mainly two types of the methods without the waters background water temperature under the influence of warm water discharge Warm water discharge temperature rise area method, if this method be do temperature according to temperature rise area periphery and obtain temperature rise area benchmark distribution of water temperature because Temperature rise area, which chooses, has randomness, and temperature rise area peripheral point temperature may be influenced by warm water discharge, peripheral point temperature with do not arranged by temperature The water temperature that water influences not necessarily is consistent, and what is determined with this is inaccurate without the waters background water temperature under the influence of warm water discharge.Second It kind is using some region of average temperature outside warm water discharge influence area as the waters background water temperature under the influence of no warm water discharge, This method has ignored the spatial variations of waters natural temperature, is also without the waters background water temperature under the influence of warm water discharge with what this was determined Inaccuracy.
Summary of the invention
Based on this, it is necessary to which in view of the above technical problems, providing one kind can accurately determine without the temperature under the influence of warm water discharge Drain prototype measurement waters background water temperature calculation method, device, computer equipment and storage medium.
A kind of warm water discharge prototype measurement waters background water temperature calculation method, which comprises
Utilize the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, hydrologic parameter, temperature Draining water intaking and draining source sink term and preset waters model meshes determine the number after engineering operation during warm water discharge prototype measurement Learn the first input boundary condition of model；
According to the actual measurement hydrodynamic force and water during warm water discharge prototype measurement after the first input boundary condition, engineering operation Temperature obtains the Reasonable Parameters value of the mathematical model as a result, to mathematical model progress calibration；Wherein, the mathematical modulo Type includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating the hydrodynamic force and water temperature point in waters Cloth rule；
According to the mathematical modulo during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of the mathematical model Second input boundary condition of type, calculates the background water temperature in waters to be detected；Wherein, the second input boundary condition is basis Meteorologic parameter, hydrologic parameter and the waters model meshes institute after the engineering operation that satellite obtains during warm water discharge prototype measurement It determines.
A kind of warm water discharge prototype measurement waters background water temperature computing device, described device include:
Determining module, for utilizing the meteorological ginseng during warm water discharge prototype measurement after the engineering operation obtained by satellite Number, hydrologic parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes, warm water discharge is former after determining engineering operation First input boundary condition of the mathematical model during type observation；
Mathematical model calibration module, for being seen according to warm water discharge prototype after the first input boundary condition, engineering operation Actual measurement hydrodynamic force and water temperature during survey obtain the reasonable ginseng of the mathematical model as a result, to mathematical model progress calibration Number value；Wherein, the mathematical model includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, is used for mould The hydrodynamic force and distribution of water temperature rule in quasi- waters；
Computing module, for warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation according to the mathematical model Second input boundary condition of the mathematical model of period, calculates the background water temperature in waters to be detected；Wherein, second input side Boundary's condition is meteorologic parameter, hydrologic parameter and the water after the engineering operation obtained according to satellite during warm water discharge prototype measurement Domain model grid is determined.
A kind of computer equipment can be run on a memory and on a processor including memory, processor and storage Computer program, the processor perform the steps of when executing the computer program
Utilize the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, hydrologic parameter, temperature Draining water intaking and draining source sink term and preset waters model meshes determine the number after engineering operation during warm water discharge prototype measurement Learn the first input boundary condition of model；
According to the actual measurement hydrodynamic force and water during warm water discharge prototype measurement after the first input boundary condition, engineering operation Temperature obtains the Reasonable Parameters value of the mathematical model as a result, to mathematical model progress calibration；Wherein, the mathematical modulo Type includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating the hydrodynamic force and water temperature point in waters Cloth rule；
According to the mathematical modulo during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of the mathematical model Second input boundary condition of type, calculates the background water temperature in waters to be detected；Wherein, the second input boundary condition is basis Meteorologic parameter, hydrologic parameter and the waters model meshes institute after the engineering operation that satellite obtains during warm water discharge prototype measurement It determines.
A kind of computer readable storage medium, is stored thereon with computer program, and the computer program is held by processor It is performed the steps of when row
Utilize the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, hydrologic parameter, temperature Draining water intaking and draining source sink term and preset waters model meshes determine the number after engineering operation during warm water discharge prototype measurement Learn the first input boundary condition of model；
According to the actual measurement hydrodynamic force and water during warm water discharge prototype measurement after the first input boundary condition, engineering operation Temperature obtains the Reasonable Parameters value of the mathematical model as a result, to mathematical model progress calibration；Wherein, the mathematical modulo Type includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating the hydrodynamic force and water temperature point in waters Cloth rule；
According to the mathematical modulo during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of the mathematical model Second input boundary condition of type, calculates the background water temperature in waters to be detected；Wherein, the second input boundary condition is basis Meteorologic parameter, hydrologic parameter and the waters model meshes institute after the engineering operation that satellite obtains during warm water discharge prototype measurement It determines.
Above-mentioned warm water discharge prototype measurement waters background water temperature calculation method, by for simulating the distribution of waters natural temperature The mathematical model of rule carries out the hydrodynamic force and water temperature contrast verification in engineering waters during warm water discharge prototype measurement after engineering operation, Calibration is carried out to mathematical model parameter, the Reasonable Parameters value of mathematical model is obtained, according to the Reasonable Parameters value of mathematical model And the second input boundary condition of the mathematical model after engineering operation during warm water discharge prototype measurement, calculate temperature after engineering operation The background water temperature during prototype measurement is drained, because it is contemplated that the waters natural temperature regularity of distribution under physical action mechanism, thus More accurate true background distribution of water temperature can be obtained.In addition, overcoming temperature rise area because the above method is without estimating temperature rise area The shortcomings that choosing has randomness, and temperature rise area peripheral point temperature may be influenced by warm water discharge, it is available more accurate reasonable Background water temperature.
Detailed description of the invention
Fig. 1 is a kind of applied environment figure of warm water discharge prototype measurement waters background water temperature calculation method in one embodiment；
Fig. 2 is a kind of flow diagram of warm water discharge prototype measurement waters background water temperature calculation method in one embodiment；
Fig. 3 is the flow diagram of the refinement step of step 201 in one embodiment；
Fig. 4 is the flow diagram of the refinement step of step 301 in one embodiment；
Fig. 5 is the flow diagram of the refinement step of step 202 in one embodiment；
Fig. 6 is the flow diagram of the refinement step of step 203 in one embodiment；
Fig. 7 is a kind of structural block diagram of warm water discharge prototype measurement waters background water temperature computing device in one embodiment；
Fig. 8 is waters model meshes schematic diagram in one embodiment；
Fig. 9 is hydrologic survey flow measurement point and tidal level measuring point distribution map；
Figure 10 is tidal level verification result schematic diagram；
Figure 11 is flow field verification result schematic diagram；
Figure 12 is the comparison schematic diagram of flood slack the model calculation and remote sensing temperature field；
Figure 13 is the comparison schematic diagram of low slack tide the model calculation and remote sensing temperature field；
Figure 14 is the statistical distribution of the surface temperature difference of surface temperature difference and actual measurement that the model of V1-V9 measuring point calculates Schematic diagram；
Figure 15 is each measuring point surface temperature contrast schematic diagram of each feature tide state of large, medium and small tide；
Figure 16 is that the background distribution of water temperature that the mathematical model after utilizing works are run during warm water discharge prototype measurement obtains is shown It is intended to；
Figure 17 is the absolute temperature field distribution schematic diagram in the case of having warm water discharge using satellite remote sensing；
Figure 18 is warm water discharge temperature rise field distribution schematic diagram；
Figure 19 is the internal structure chart of computer equipment in one embodiment.
Specific embodiment
It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not For limiting the application.
Background water temperature calculation method in waters provided by the present application, can be applied in application environment as shown in Figure 1.Its In, data acquisition device 10 is connected with computer equipment 20.Optionally, which is satellite, can using satellite To obtain a wide range of and accurate meteorological model remotely-sensed data.Optionally, which includes at least one processor 21 with memory 22.Optionally, the processor 21 can for CPU (Central Processing Unit, central processing unit), IPU (Intelligence Processing Unit, intelligent processor) etc..Optionally, which can be monokaryon Processor or multi-core processor.
Wherein, including the mathematics during warm water discharge prototype measurement after engineering operation in the memory 22 of the computer equipment 20 Model, the mathematical model after the engineering operation during warm water discharge prototype measurement can choose different equation form and software and carry out Simulation.Optionally, the governing equation of the mathematical model in the application include continuity equation, the equation of momentum, matter transportation equation and Heat exchange equation.
In one embodiment, as shown in Fig. 2, providing a kind of warm water discharge prototype measurement waters background water temperature calculation method Flow diagram, in this way be applied to Fig. 1 in computer equipment 20 for be illustrated, comprising the following steps:
Step 201, the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, the hydrology are utilized Parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes, determine warm water discharge prototype measurement after engineering operation First input boundary condition of the mathematical model of period；
Wherein, the meteorologic parameter after engineering operation during warm water discharge prototype measurement, hydrologic parameter are got using satellite Satellite data, satellite data using National Weather environmental forecasting center (NCEP) Meteorology Forecast System analysis of data again (CSFR), the minimum time resolution ratio of data is 1 hour, and horizontal resolution is up to 0.2 degree.
Wherein, the meteorologic parameter after engineering operation during warm water discharge prototype measurement include: gas epidemic disaster, wind speed, wind direction, The spatial and temporal distributions such as solar radiation, atmospheric radiation, cloud amount, air pressure.Hydrologic parameter after engineering operation during warm water discharge prototype measurement It include: spatial and temporal distributions such as water temperature, run-off etc..Wherein, waters model meshes use a kind of multi-level far and near area nested grid, Waters is divided into multiple regions by the grid, as shown in figure 8, waters to be detected is divided into 4 regions, respectively region 1., region 2., region 3. and region 4..
Wherein, the first input boundary condition includes: flow field boundary condition and temperature boundaries condition.Flow field boundary condition packet Include: water level opens boundary, flow field initial fields, bottom friction, water surface friction, water intaking and draining source sink term etc..Optionally, water level opens side Boundary uses tidal boundary；Flow field initial fields are calculated using stationary stream condition；Bottom friction is the bottom type according to waters to be detected It determines；Water surface friction is calculated using water surface wind-stress, and water surface wind-stress indicates are as follows:
In formula (2-1), ρaIndicate atmospheric density, CdIndicate wind drag stress coefficient, U10Indicate the wind of the above 10m of Free Surface Speed.
Wherein, temperature boundaries condition includes: bank circle, bottom boundary, opens boundary, Free Surface, temperature initial fields, temperature Draining water intaking and draining source sink term etc..Optionally, bank circle and bottom boundary are indicated using adiabatic boundary；Boundary is opened using environment Water temperature indicates；Free Surface is calculated using marine ecosystem model model and is indicated；Temperature initial fields are using ambient water thermal field as initial ?；Warm water discharge is fetched water and is drained source sink term and determines heat flux according to given flow and temperature value, wherein water intaking temperature is to fetch water Place calculates the calculating temperature value of grid, and temperature at discharging condensate is to fetch water temperature plus given taking and discharging temperature difference.Wherein, water intaking and Draining source sink term is given flow.
It should be noted that above-mentioned waters includes sea area etc., refer to river, river, lake, sea from the water surface to water-bed a certain range.
Step 202, according to the actual measurement hydrodynamic force during warm water discharge prototype measurement after the first input boundary condition, engineering operation With water temperature as a result, carrying out calibration to mathematical model, the Reasonable Parameters value of mathematical model is obtained；Wherein, mathematical model includes connecting Continuous equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating the hydrodynamic force and distribution of water temperature rule in waters；
Wherein, the actual measurement hydrodynamic force after engineering operation during warm water discharge prototype measurement and water temperature result include: engineering waters Observed current data flow velocity, flow direction, water level and actual measurement water temperature time series value on fixed point.
Wherein, the simulation process of the mathematical model are as follows:
Waters Hydrodynamic field is simulated first with continuity equation and the equation of momentum, then based on Hydrodynamic field as a result, utilizing object Matter Boltzman transport equation Boltzmann and heat exchange equation can simulate water temperature field.The Hydrodynamic field and water temperature field that the modeling obtains may not Accurately, it needs to carry out parameter calibration to mathematical model using the hydrodynamic force and water temperature value of actual measurement, obtaining can accurate simulation engineering water The Rational Parameters value of the mathematical model of domain hydrodynamic force and distribution of water temperature.
Wherein, mathematical model is assumed based on Boussinesq and vertical hydrostatic pressure, using three-dimensional in simulation process Curvilinear coordinate system, the vertical coordinate representation in three-dimensional curve coordinate system are as follows:
In formula (2-2), d is reference planes water depth value below, and H is the full depth of water, and z is the vertical coordinate of physical space, and ζ is certainly By surface elevation, in bottom σ=- 1, surface layer σ=0.
Wherein, continuity equation are as follows:
In formula (2-3), t is the time, and ξ and η are the coordinate of the horizontal direction under three-dimensional curve coordinate system, Gξξ、GηηFor three-dimensional Curvilinear coordinate system is converted to the conversion coefficient of rectangular coordinate system, and U is the vertical mean flow rate in the direction ξ, and V is the vertical flat of the direction η Equal flow velocity, Q indicate source sink term.
The equation of momentum are as follows:
In formula (2-4) and formula (2-5), u, v, ω are the velocity components under curvilinear orthogonal coordinate system, and f is Coriolis ginseng Number, FξFor the turbulent fluctuation momentum flux in the direction ξ, FηFor the turbulent fluctuation momentum flux in the direction η, PξFor the hydraulic pressure force gradient in the direction ξ, PηFor η The hydraulic pressure force gradient in direction, νVFor the vertical coefficient of viscosity, MξFor the source sink term of momentum on the direction ξ, MηFor the source of momentum on the direction η Remittance item.
Under shallow water assumption, the vertical equation of momentum is reduced to hydrostatic pressing equation:
Wherein, P indicates pressure, and g indicates that constant, ρ indicate water body density.
Hydrostatic pressing indicates after integral are as follows:
Wherein, PatmIndicate atmospheric pressure.
It should be noted that in simulations, using baroclinic model, water body density p is the function of temperature and salinity, use UNESCO state equation indicates that each solution procedure seeks density contrast caused by the distribution of temperature and salinity by state equation, Reflected in pressure term P.
Vertical flow velocity ω is solved by continuity equation under σ coordinate system:
In formula (2-8), vertical flow velocity ω is defined in equal σ plane, indicates the catenary motion relative to σ plane.
Wherein, F is being solvedξAnd FηWhen, it is related to horizontal coefficient of viscosity νH, expression formula are as follows:
ν H=ν 2D+ν3D+νmol (2-9)
Wherein, in formula (2-9), νmolFor molecular viscosity coefficient, ν2DThe horizontal turbulent fluctuation phase of sub- grid not solved for horizontal direction The turbulence viscosity coefficient of pass can pass through GUI UI Preferences ν in Delft3D2DFor whole audience steady state value or with spatial variations Value, can also be solved by HLES method.ν3DFor the related turbulence viscosity coefficient of three-dimensional turbulent fluctuation, can be asked by Turbulent Model Solution, wherein Delft3D provides Prandtl and mixes long model, and model is k-L model and k- ε model.
Wherein, vertical coefficient of viscosity νVExpression formula are as follows:
νV=νmol+max(νback, ν3D) (2-10)
Wherein, in formula (2-10), νbackFor vertical viscosity background value, mainly reflect other blending effects not solved (the interior wave action in such as stratified flow).
Above-mentioned equation is the component part of hydrodynamic force mathematical model, utilizes above-mentioned equation and matter transportation equation and heat exchange Equation can simulate waters distribution of water temperature.
Wherein, matter transportation equation are as follows:
In formula (2-11), C is water temperature, DHFor horizontal proliferation coefficient, DVFor vertical diffusion coefficient, λdFor substance decay coefficient (/day), λdThe determination method analogy water flow coefficient of viscosity, S is source sink term
When solving matter transportation equation, need to consider the heat exchange between water body Free Surface and atmosphere.Water body absorbs Origin of heat in solar radiation and atmospheric radiation, water body transmits heat mainly by water body radiation, evaporation and to wandering outward Heat.Heat exchange equation is needed using heat flux calculation formula, and the heat flux calculation formula constant value that different models use slightly has It is different, it is desirable that the meteorologic parameter of input difference can pass through tentative calculation reasonable selection.
The total heat flux Q of Free SurfacetotIt is calculated using following formula:
Qtot=Qsn+Qan-Qbr-Qev-Qco (2-12)
In formula (2-12), QsnFor solar radiation flux, QanFor atmospheric radiation flux, QbrFor the outside radiation flux of water body, Qev For evaporative heat loss flux, QcoFor heat loss through convection flux.Each subitem heat flux calculation formula is as follows:
(1) net solar radiation heat flux Qsn:
Qsn=(1- α) Qsc(1.0-0.4Fc-0.38Fc 2) (2-13)
In formula (2-13), α is reflection coefficient, takes 0.06, FcFor cloud amount coefficient, QscIt is logical for solar radiation under the conditions of clear sky Amount, QscCalculation formula are as follows:
In formula (2-14), S is solar constant, and value 1368 [J/m2s], γ is solar elevation, with geographic latitude and when Between it is related.
(2) net long-wave radiation Qeb:Qeb=Qbr-Qan (2-15)
In marine ecosystem model model, radiation from atmosphere to water body and water body by think it is all long-wave radiation to external radiation A part, the two merge into net averaged long wave radiation flux Qeb。
QebIt is calculated using following formula:
In formula (2-16), ε is emission ratio, and σ is Stefan-Boltsman constant, 5.67 × 10-8 of value [J/ M2sK4], T is bucket temperature [K], eaFor vapour pressure, FcFor cloud amount coefficient.
(3) evaporative heat loss Qev:
Qev=Qev,forced+Qev,free (2-17)
Wherein, Qev,forcedTo force evaporative heat loss flux, driven by wind, Qev,freeFor free evaporative heat loss flux, by big Buoyancy caused by density contrast drives in gas boundary layer.
Wherein, evaporative heat loss flux is forced to calculate using following formula:
Qev,forced=LVρaceU10{qs(Ts)-qa(Ta)} (2-18)
L in formulaVFor evaporation latent heat [J/kg], LV=2.5*106-2.3*103Ts, ρaFor atmospheric density, qsFor saturated air Humidity ratio, qaFor 10m high air humidity ratio, ceFor Dalton number, value 0.0015, U10For function of wind speed.
Wherein, free evaporative heat loss flux is calculated using following formula:
In formula (2-19), ksFor the coefficient of heat conduction,For sea to the high-altitude 10m average air density, remaining parameter is the same.
(4) heat loss through convection
Heat loss through convection is also divided into forced convection heat dissipation Q in marine ecosystem model modelco,forcedIt radiates with free convection Qco,free, heat loss through convection total amount is the sum of two parts:
Qco=Qco,forced+Qco,free (2-20)
Forced convection heat dissipation flux Qco,forcedIt is calculated using following formula:
Qco,forced=ρacpcHU10(Ts-Ta) (2-21)
ρ in formulaaFor atmospheric density, cpFor air specific heat, value 1004.0 [J/kgK], cHFor Stanton number, value 0.00145, TsFor water surface temperature, TaFor air themperature.
Free convection heat dissipation flux Qco,freeIt is calculated using following formula:
In formula (2-22), ksFor the coefficient of heat conduction,For sea to the high-altitude 10m average air density, remaining parameter is the same.
Waters hydrodynamic force and distribution of water temperature are simulated using above-mentioned formula (2-1) to formula (2-22).
Step 203, according to the number during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of mathematical model The the second input boundary condition for learning model, calculates the background water temperature in waters to be detected；Wherein, the second input boundary condition is basis Meteorologic parameter, hydrologic parameter and waters model meshes institute after the engineering operation that satellite obtains during warm water discharge prototype measurement are really It is fixed.
In embodiments of the present invention, when needing to calculate the background water temperature of waters period to be detected, it is necessary first to right Mathematical model carries out calibration, obtains the Reasonable Parameters value of mathematical model.After obtaining mathematical model Reasonable Parameters value, again The second input boundary condition is determined using meteorologic parameter, hydrologic parameter and the waters model meshes that satellite obtains, and utilizes second The mathematical model for inputting boundary condition and aforementioned calibration, calculates the background water temperature in waters to be detected.
Wherein, meteorologic parameter, hydrologic parameter indicate the remotely-sensed data of satellite acquisition in above-mentioned certain period.Wherein, on Stating the period can be a few hours, one day, a couple of days etc..Preferably, in the present embodiment using by hour data.
Above-mentioned waters background water temperature calculation method method, by the mathematical model for simulating waters distribution of water temperature rule Calibration is carried out, obtains mathematical model Reasonable Parameters value, on this basis calculated background water temperature, it is contemplated that physical action machine The lower waters natural temperature regularity of distribution is made, it is hereby achieved that more accurate true background water temperature.In addition, because the above method without It need to estimate temperature rise area, overcoming temperature rise area and choosing has randomness, and temperature rise area peripheral point temperature may be influenced by warm water discharge Disadvantage, available more accurate true background water temperature.
As an alternative embodiment, as shown in figure 3, being the flow diagram of the refinement step of step 201, specifically Include:
Step 301, the meteorologic parameter during warm water discharge prototype measurement after engineering operation, hydrologic parameter are interpolated into waters mould Type grid obtains the third input boundary condition of mathematical model；
In embodiments of the present invention, grid dividing is carried out to waters to be detected, obtains waters model meshes.Later, by gas As parameter, hydrologic parameter are interpolated into waters model meshes, the third input boundary condition of mathematical model is obtained.
Wherein, because the real-time change of waters background distribution of water temperature is a wide range of meteorological and the coupling of hydrology field as a result, It is limited by off-lying sea water temperature data and accurately obtains difficult problem, in order to eliminate as much as boundary condition assignment to engineering waters temperature The influence of analog result, simulating boundary should be as far as possible far from Project Areas, and significance level height is simulated in engineering near region, and water temperature simulates resolution ratio It is required that high, near region water front precision is affected to analog result.Therefore, grid dividing is carried out to waters to be detected, such as Fig. 8 institute Show, engineering far field model takes relatively crude grid, and engineering near region uses thinner grid.Engineering far field is used to engineering near region Two layers of transition region, grid gradually refine.Mathematical model uses four layers of nested grid in total in embodiment, wherein engineering far region modeling Off-lying sea range tapers into power plant periphery Engineering Zone 4. area as shown in the area Tu8Zhong①, through 2. area and 3. area's transition.① The grid plan scale average out to 5000m × 5000m in area, 2. the plane grid scale in area is 1100m × 1100m, and 3. area is flat Surface grids scale is 220m × 220m, and 4. the plane grid scale in area is 44m × 44m.It is in fig. 8, vertical to use σ grid, It is analyzed through tentative calculation, in order to balance simulation precision and computational efficiency, it is vertical upper to be finally divided into 6 layers.
Step 302, warm water discharge is fetched water and is drained source sink term assignment to waters model meshes, obtain the temperature row of mathematical model The input of water source sink term；
Step 303, the third of mathematical model is inputted into boundary condition and warm water discharge source sink term input the first input of composition side Boundary's condition.
In embodiments of the present invention, it is inputted using third input boundary condition, the warm water discharge source sink term of mathematical model, is existing Flow field boundary condition algorithm and existing temperature boundaries condition algorithm, determine the first input boundary condition.Because of flow field side Boundary's condition algorithm, temperature boundaries condition algorithm are all the prior art, are no longer repeated in detail herein.
Above-mentioned warm water discharge prototype measurement waters background water temperature calculation method, by the meteorologic parameter acquired by satellite, water Literary parameter interpolation has carried out calibration into waters model meshes, to mathematical model parameter, utilizes the mathematical model calculated Bottom water temperature, it is contemplated that the natural temperature regularity of distribution in waters under physical action mechanism, it is hereby achieved that more accurate true sheet Bottom water temperature.In addition, carrying out grid dividing to entire waters to be detected, estimation temperature rise area is not needed, temperature rise area is overcome and chooses tool There are the shortcomings that randomness, temperature rise area peripheral point temperature may be influenced by warm water discharge, available more accurate reasonable background Water temperature.
As an alternative embodiment, as shown in figure 4, being the flow diagram of the refinement step of step 301, specifically Include:
Step 401, the meteorologic parameter during warm water discharge prototype measurement after engineering operation, hydrologic parameter are subjected to file format Conversion, text file needed for obtaining mathematical model input；
Step 402, text file is interpolated into waters model meshes, obtains the third input boundary condition of mathematical model.
It in embodiments of the present invention, is that can not directly be interpolated into waters using the meteorologic parameter of satellite acquisition, hydrologic parameter In model meshes, need to format meteorologic parameter, hydrologic parameter, it is preferred that using wgrib2 software and Meteorologic parameter, the format conversion of hydrologic parameter are text formatting required for mathematical model, and are interpolated into water by Matlab software Domain model grid obtains the third input boundary condition of mathematical model.
It should be noted that before carrying out document format conversion, needing will be by meteorologic parameter that satellite obtains, the hydrology Parameter and Project Areas bank Weather Station Data compare verification, confirm meteorologic parameter, hydrologic parameter reasonability.If passing through satellite It is not much different between the parameter and Project Areas bank Weather Station Data of acquisition, then can directly uses.Larger difference if it exists, can Use after suitably correcting to meteorologic parameter, hydrologic parameter using the methods of interpolation fitting.
Above-mentioned waters background water temperature calculation method, it is available by being formatted to meteorologic parameter, hydrologic parameter The text file that mode input needs is interpolated into waters model meshes, available number by the text file that mode input needs The third for learning model inputs boundary condition.
As an alternative embodiment, as shown in figure 5, being the flow diagram of the refinement step of step 202, specifically Include:
Step 501, the first input boundary condition is input to mathematical model and solved, obtain calculating hydrodynamic force distribution and Calculate distribution of water temperature；
In embodiments of the present invention, the first input boundary condition is input in mathematical model, moving model is calculated Hydrodynamic force distribution and calculating distribution of water temperature.Wherein, calculating hydrodynamic force distribution and calculating distribution of water temperature is the hydrodynamic(al) after moving model The calculated result of power and water temperature.
Step 502, whether the distribution of contrast conting hydrodynamic force and actual measurement hydrodynamic force distribution are consistent, and contrast conting distribution of water temperature It is whether consistent with actual measurement distribution of water temperature；
In embodiments of the present invention, the distribution of actual measurement hydrodynamic force and actual measurement distribution of water temperature are the true of the waters to be detected of measurement Data distribution.
Step 503, consistent with actual measurement hydrodynamic force distribution if calculating hydrodynamic force distribution, and calculate distribution of water temperature and actual measurement water temperature Distribution is consistent, then the parameter value of the mathematical model during warm water discharge prototype measurement after engineering operation is determined as mathematical model Reasonable Parameters value.
In embodiments of the present invention, using existing comparison method, contrast conting data distribution and measured data is gone to be distributed Whether it is consistent, optionally, existing comparison method is to compare root mean square, standard deviation etc..Because of above-mentioned data distribution comparison method It is the prior art, will not be repeated here herein.
Wherein, the calculated result that data distribution is exactly the hydrodynamic force and water temperature after moving model is calculated.If calculating hydrodynamic force Distribution is consistent with actual measurement hydrodynamic force distribution, and it is consistent with actual measurement distribution of water temperature to calculate distribution of water temperature, then it represents that mathematical model can Accurate rationally reflection hydrodynamic force distribution and distribution of water temperature.If there is any one distribution inconsistent in the two, adjusting parameter is needed Value, and solution and verifying are re-started, it is consistent with actual measurement hydrodynamic force distribution until calculating hydrodynamic force distribution, and calculate water temperature point Cloth is consistent with actual measurement distribution of water temperature.
Wherein, hydrodynamic force distribution verifying is to choose Project Areas nearby tidal level station and stream gaging station progress tidal level process and flow velocity mistake Journey verifying, precision will meet relevant code requirement, if not satisfied, need to adjust the parameters such as roughness, rerun model until Meet relevant code requirement.Specifically, hydrologic survey flow measurement point and tidal level measuring point distribution map are as shown in figure 9, tidal level verifying knot Fruit is as shown in Figure 10, and flow field verifying is as shown in figure 11, wherein the main verifying flow velocity of flow field verifying and flow direction.It can be with from Figure 10 Find out, is coincide using the calculated tidal level process phase of model and actual measurement good, except part peak, valley deviation are somewhat larger, meter The tidal level value of calculation is most of to coincide well with actual measurement.It is verified from the tidal current speed of the trend survey station in Figure 11 with vertical layered is flowed to As a result as can be seen that the flow direction and actual measurement good agreement simulated, flow velocity size and actual measurement major part are preferably coincide, illustrate model energy The hydrodynamic force distribution in enough rationally reflection engineering waters.
Distribution of water temperature verifying is divided into two steps, and the first step verifies plane distribution, and plane distribution uses satellite remote sensing surface temperature Field is verified, and mainly verifying the temperature field being calculated from offshore to the distribution in offshore direction and along the distribution in bank direction is The no remote sensing distribution with actual measurement is consistent.Second step verifies the absolute precision of water temperature, and the absolute precision of water temperature is using before engineering construction The measuring point water temperature of hydrologic survey is verified, and verifying calculates the absolute precision of water temperature, and standard deviation is within 0.5 DEG C.If meter The distribution of water temperature of calculation and actual measurement remote sensing distribution of water temperature are inconsistent or standard deviation is excessive, then need to adjust horizontal proliferation coefficient Etc. parameters, rerun model until meeting relevant code requirement.
Above-mentioned waters background water temperature calculation method, be distributed by contrast conting hydrodynamic force and survey hydrodynamic force distribution whether one It causes, and whether contrast conting distribution of water temperature and actual measurement distribution of water temperature are consistent, obtain the water that can rationally accurately reflect engineering waters The mathematical model parameter reasonable value of power distribution and distribution of water temperature, can calculate reasonable background water temperature on this basis.
As an alternative embodiment, as shown in fig. 6, being the flow diagram of the refinement step of step 203, specifically Include:
Step 601, the warm water discharge source sink term input for deleting the mathematical model in the first input boundary condition, obtains engineering fortune Second input boundary condition of the mathematical model after row during warm water discharge prototype measurement；
In embodiments of the present invention, mathematical model when having the waters water temperature under warm water discharge emission behaviour after model engineering operation It inputs comprising a warm water discharge source sink term input in boundary condition, warm water discharge source sink term input is a constraint condition.It deletes After warm water discharge source sink term input, obtain simulating the background water temperature computation model input perimeter strip that the former sight phase discharge without warm water discharge Part, i.e., the second input boundary condition.
Step 602, by the second input boundary condition and the Reasonable Parameters value of mathematical model, it is input to mathematical model progress It solves, obtains the background water temperature after engineering operation during warm water discharge prototype measurement.
It will be solved without thermal discharge model using obtaining the second input boundary condition in step 203 and be input to the former sight phase, The background water temperature in available waters to be detected.
Above-mentioned warm water discharge prototype measurement waters background water temperature calculation method, by the number for simulating waters distribution of water temperature It learns model and carries out calibration, the mathematical model parameter value of energy feasible simulation waters distribution of water temperature rule is obtained, using this calibration mistake The calculated background water temperature of mathematical model, it is contemplated that the natural temperature regularity of distribution in waters under physical action mechanism.In addition, because upper Method is stated without estimating temperature rise area, overcoming temperature rise area and choosing has randomness, and temperature rise area peripheral point temperature may be arranged by temperature The shortcomings that water influences, therefore using the available more accurate true background water temperature of this calculation method.
Wherein, can also utilize meteorologic parameter without warm water discharge before the engineering operation obtained by satellite, hydrologic parameter and Waters model meshes determine the 4th boundary condition；4th boundary condition is input to mathematical model to solve, obtains engineering fortune Hydrodynamic force and water temperature calculated result before row；Wherein, the hydrodynamic force before engineering operation and water temperature calculated result include engineering operation The preceding distribution of calculating hydrodynamic force and calculating distribution of water temperature.Hydrodynamic force and water temperature calculated result and engineering fortune before comparing engineering operation Whether the hydrodynamic force before row is consistent with water temperature measured result；If hydrodynamic force and water temperature calculated result and engineering fortune before engineering operation Hydrodynamic force before row is consistent with water temperature measured result, then will be determined as before purpose project operation without warm water discharge without temperature before engineering operation Draining.Wherein, if the calculating hydrodynamic force distribution before engineering operation is consistent with the actual measurement hydrodynamic force distribution before engineering operation, and engineering Calculating distribution of water temperature before operation is consistent with the actual measurement distribution of water temperature before engineering operation, then it represents that the mathematical model can be rationally quasi- The really hydrodynamic force distribution and distribution of water temperature in the reflection preceding waters of engineering.It is tied using being calculated before the engineering operation without warm water discharge mathematical model The comparison of fruit and measured result can verify reasonability of the mathematical model simulation without the natural temperature distribution under warm water discharge, into one Step illustrates the reliability of mathematical model simulation background water temperature.
During water temperature observed result carries out model verifying before utilizing works, engineering waters before power plant runs is chosen Remote sensing image data and the model calculation have carried out contrast verification.Remote sensing image derives from Landsat7 satellite data, satellite It is in the course of transit larger by the probability of sky cloud effect, therefore power plant's lesser satellite image of periphery cloud amount can only be chosen as far as possible and carried out pair Than verifying.The comparison of the model calculation and remote sensing temperature field that July in 2005 (flood slack) on the 28th is given in Figure 12 is illustrated Scheme, give the model calculation on May 28th, 2006 (low slack tide) and the comparison schematic diagram in remote sensing temperature field in Figure 13, The data on July 28th, 2005 are influenced by cloud, and having the part of cloud in remote sensing temperature image is blank value, can be with from remote sensing image Find out that remote sensing temperature is more uneven, the jump of adjacent pixel temperature value is larger, and has certain band to influence.However model calculates knot Fruit is then more continuous Temperature Distribution.But from the whole audience absolute temperature Data-Statistics of Figure 12 and Figure 13, the model calculation with Remote sensing temperature difference is generally within 2 DEG C, it is contemplated that remote sensing absolute temperature accuracy is limited, but temperature Relative distribution is able to reflect water The trend of domain bulk temperature variation, therefore can be verified from temperature Relative distribution angle and utilize the calculated temperature field of model (the model calculation) opposite variation is reasonable.Shown using the calculated temperature field of model and the comparison of remote sensing Temperature Distribution:
(1) in terms of from offshore to offshore direction temperature by the high trend being lower, offshore high temperature bandwidth and the remote sensing of calculating become Gesture is consistent.
(2) from along in terms of the Temperature Distribution of bank direction, remote sensing temperature field is shallowly curved with two sides to the north of killing of dog mouthpiece on the south the angle Fu Zhou Always there is higher temperature band in area, the factory site warm water discharge influence area angle Fu Zhou to one band of killing of dog mouthpiece is by south and north flanks height temperate zone shadow It rings, on July 28th, 2005 flood slack temperature field, the angle Fu Zhou north and south side height temperate zone forms a piece of, and killing of dog mouthpiece north side then exists One relative low temperature area, be since northwards flood current causes, on May 28th, 2006 low slack tide temperature field, the angle Fu Zhou southern side high temperature Band does not bypass the angle Fu Zhou, and it is to utilize mould since flood current to the south causes that development trend to the south, which is presented, in killing of dog mouthpiece north side height temperate zone Similar distribution trend is also presented in the temperature field that type calculates.
It can be seen that the background temperature field offshore variation simulated using mathematical model from above-mentioned comparing result and change along bank Trend and remote sensing temperature results are almost the same, show that the temperature field space Relative distribution of simulation is reasonable.
Figure 14 gives all measurement moment of all measuring points of the V1-V9 in Fig. 9, the superficial water calculated using mathematical model The statistical distribution of temperature difference and the surface temperature difference of actual measurement, it can be seen that temperature difference largely falls in ± 0.5 DEG C of section, statistics It is 0.44 DEG C that analysis, which obtains calculating with the surface temperature root-mean-square error of actual measurement,.
Figure 15 gives each measuring point surface temperature comparison of each feature tide state of large, medium and small tide, it can be seen that model calculated Water temperature variation tendency is overall consistent with the water temperature variation tendency of actual measurement, and above-mentioned verification result shows that model energy accurate simulation goes out temperature Field distribution situation.
Above-mentioned warm water discharge prototype measurement waters background water temperature calculation method, using before the engineering operation without warm water discharge mathematical modulo Type, can verify can be gone out without the hydrodynamic force distribution and distribution of water temperature under warm water discharge using the mathematical model with accurate simulation.Because of work Calculated before Cheng Yunhang use mathematical model and engineering operation after the mathematical model of use is calculated during warm water discharge prototype measurement Governing equation is identical with model parameter value, by improving engineering fortune to the verifying without warm water discharge mathematical model before engineering operation The reliability that mathematical model used by calculating during warm water discharge prototype measurement after row calculates purpose project waters.Using can be with Rationally the mathematical model of the reflection waters natural temperature regularity of distribution calculates background water temperature, because it is contemplated that physical action mechanism is lauched The domain natural temperature regularity of distribution, it is hereby achieved that more accurate true background water temperature.
It in embodiments of the present invention, is as shown in figure 16 to be obtained using the mathematical model that water temperature data calibration is crossed after actual measurement Background water temperature distribution schematic diagram is the absolute temperature field distribution in the case of having warm water discharge for utilizing satellite remote sensing as shown in figure 17 Schematic diagram deducts the background distribution of water temperature of Figure 16 using the absolute temperature field distribution of Figure 17, it can obtain as shown in figure 18 Warm water discharge temperature rise field distribution.
It should be understood that although each step in the flow chart of Fig. 2-6 is successively shown according to the instruction of arrow, These steps are not that the inevitable sequence according to arrow instruction successively executes.Unless expressly stating otherwise herein, these steps Execution there is no stringent sequences to limit, these steps can execute in other order.Moreover, at least one in Fig. 2-6 Part steps may include that perhaps these sub-steps of multiple stages or stage are not necessarily in synchronization to multiple sub-steps Completion is executed, but can be executed at different times, the execution sequence in these sub-steps or stage is also not necessarily successively It carries out, but can be at least part of the sub-step or stage of other steps or other steps in turn or alternately It executes.
In one embodiment, as shown in fig. 7, providing a kind of warm water discharge prototype measurement waters background water temperature calculates dress It sets, comprising: determining module 701, mathematical model calibration module 702 and computing module 703, in which:
Determining module 701, for utilizing the meteorology during warm water discharge prototype measurement after the engineering operation obtained by satellite Parameter, hydrologic parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes, determine warm water discharge after engineering operation First input boundary condition of the mathematical model during prototype measurement；
Mathematical model calibration module 702, for according to warm water discharge prototype measurement after the first input boundary condition, engineering operation The actual measurement hydrodynamic force and water temperature of period obtains the Reasonable Parameters value of mathematical model as a result, to mathematical model progress calibration；Its In, mathematical model includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating the hydrodynamic(al) in waters Power and distribution of water temperature rule；
Computing module 703, for warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation according to mathematical model Second input boundary condition of the mathematical model of period, calculates the background water temperature in waters to be detected；Wherein, the second input perimeter strip Part is meteorologic parameter, hydrologic parameter and the waters model net after the engineering operation obtained according to satellite during warm water discharge prototype measurement Lattice are determined.
As an alternative embodiment, determining module 701 is specifically used for: by warm water discharge prototype measurement after engineering operation The meteorologic parameter of period, hydrologic parameter are interpolated into waters model meshes, obtain the third input boundary condition of mathematical model；It will be warm Draining water intaking and draining source sink term assignment obtain the warm water discharge source sink term input of mathematical model to waters model meshes；By mathematics The third input boundary condition and warm water discharge source sink term input composition the first input boundary condition of model.
As an alternative embodiment, further including grid dividing module, grid dividing module is used for water to be detected Domain carries out grid dividing, obtains waters model meshes.
As an alternative embodiment, determining module 701 also particularly useful for: warm water discharge prototype after engineering operation is seen Meteorologic parameter, hydrologic parameter during survey carry out document format conversion, text file needed for obtaining mathematical model input；It will be literary This document is interpolated into waters model meshes, obtains the third input boundary condition of mathematical model.
As an alternative embodiment, mathematical model calibration module 702 is specifically used for: inputting boundary condition for first It is input to mathematical model to be solved, obtains calculating hydrodynamic force distribution and calculates distribution of water temperature；Contrast conting hydrodynamic force distribution and Whether consistent survey hydrodynamic force distribution, and whether contrast conting distribution of water temperature and actual measurement distribution of water temperature are consistent；If calculating hydrodynamic force Distribution is consistent with actual measurement hydrodynamic force distribution, and calculating distribution of water temperature is consistent with actual measurement distribution of water temperature, then arranges temperature after engineering operation The parameter value of mathematical model during water prototype measurement is determined as the Reasonable Parameters value of mathematical model.
As an alternative embodiment, computing module 703 is specifically used for: deleting the number in the first input boundary condition The warm water discharge source sink term input for learning model obtains the second input of the mathematical model after engineering operation during warm water discharge prototype measurement Boundary condition；By the second input boundary condition and the Reasonable Parameters value of mathematical model, it is input to mathematical model and is solved, obtained Background water temperature after to engineering operation during warm water discharge prototype measurement.
Specific restriction about warm water discharge prototype measurement waters background water temperature computing device may refer to above for temperature The restriction of prototype measurement waters background water temperature calculation method is drained, details are not described herein.Above-mentioned warm water discharge prototype measurement waters sheet Modules in bottom water temperature computing device can be realized fully or partially through software, hardware and combinations thereof.Above-mentioned each module It can be embedded in the form of hardware or independently of in the processor in computer equipment, computer can also be stored in a software form In memory in equipment, the corresponding operation of the above modules is executed in order to which processor calls.
In one embodiment, a kind of computer equipment is provided, which can be server, internal junction Composition can be as shown in figure 19.The computer equipment include by system bus connect processor, memory, network interface and Database.Wherein, the processor of the computer equipment is for providing calculating and control ability.The memory packet of the computer equipment Include non-volatile memory medium, built-in storage.The non-volatile memory medium is stored with operating system, computer program and data Library.The built-in storage provides environment for the operation of operating system and computer program in non-volatile memory medium.The calculating The database of machine equipment is for storing meteorological data and hydrographic data.The network interface of the computer equipment is used for and external end End passes through network connection communication.To realize a kind of waters background water temperature calculation method when the computer program is executed by processor.
It will be understood by those skilled in the art that structure shown in Figure 19, only part relevant to application scheme The block diagram of structure, does not constitute the restriction for the computer equipment being applied thereon to application scheme, and specific computer is set Standby may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
In one embodiment, a kind of computer equipment is provided, including memory, processor and storage are on a memory And the computer program that can be run on a processor, processor are performed the steps of when executing computer program using by defending Meteorologic parameter, hydrologic parameter, warm water discharge water intaking and draining source after the engineering operation that star obtains during warm water discharge prototype measurement are converged And preset waters model meshes, determine the first input side of the mathematical model after engineering operation during warm water discharge prototype measurement Boundary's condition；According to the actual measurement hydrodynamic force and water temperature knot during warm water discharge prototype measurement after the first input boundary condition, engineering operation Fruit carries out calibration to mathematical model, obtains the Reasonable Parameters value of mathematical model；Wherein, mathematical model includes continuity equation, moves Equation, transport of substances equation and heat exchange equation are measured, for simulating the hydrodynamic force and distribution of water temperature rule in waters；According to mathematical modulo Second input boundary condition of the mathematical model after the Reasonable Parameters value and engineering operation of type during warm water discharge prototype measurement, meter Calculate the background water temperature in waters to be detected；Wherein, the second input boundary condition is warm water discharge after the engineering operation obtained according to satellite Meteorologic parameter, hydrologic parameter and waters model meshes during prototype measurement are determined.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated Machine program performs the steps of when being executed by processor utilizes the warm water discharge prototype measurement phase after the engineering operation obtained by satellite Between meteorologic parameter, hydrologic parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes, determine engineering operation First input boundary condition of the mathematical model during warm water discharge prototype measurement afterwards；According to the first input boundary condition, engineering fortune Actual measurement hydrodynamic force and water temperature after row during warm water discharge prototype measurement obtain mathematical model as a result, to mathematical model progress calibration Reasonable Parameters value；Wherein, mathematical model includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, is used Hydrodynamic force and distribution of water temperature rule in simulation waters；According to warm water discharge after the Reasonable Parameters value and engineering operation of mathematical model Second input boundary condition of the mathematical model during prototype measurement, calculates the background water temperature in waters to be detected；Wherein, second is defeated Entering boundary condition is meteorologic parameter, hydrologic parameter and the water after the engineering operation obtained according to satellite during warm water discharge prototype measurement Domain model grid is determined.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein, To any reference of memory, storage, database or other media used in each embodiment provided herein, Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM (PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms, Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance Shield all should be considered as described in this specification.
The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.

Claims (9)

1. a kind of warm water discharge prototype measurement waters background water temperature calculation method, which is characterized in that the described method includes:
Utilize the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, hydrologic parameter, warm water discharge Water intaking and draining source sink term and preset waters model meshes determine the mathematical modulo after engineering operation during warm water discharge prototype measurement First input boundary condition of type；
According to the actual measurement hydrodynamic force and water temperature knot during warm water discharge prototype measurement after the first input boundary condition, engineering operation Fruit carries out calibration to the mathematical model, obtains the Reasonable Parameters value of the mathematical model；Wherein, the mathematical model packet Continuity equation, the equation of momentum, transport of substances equation and heat exchange equation are included, for simulating hydrodynamic force and the distribution of water temperature rule in waters Rule；
According to the mathematical model during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation of the mathematical model Second input boundary condition, calculates the background water temperature in waters to be detected；Wherein, the second input boundary condition is according to satellite Meteorologic parameter, hydrologic parameter and waters model meshes institute after the engineering operation of acquisition during warm water discharge prototype measurement are really It is fixed.
2. method according to claim 1, which is characterized in that described to utilize warm water discharge after the engineering operation obtained by satellite Meteorologic parameter, hydrologic parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes during prototype measurement, really Determine the first input boundary condition of the mathematical model after engineering operation during warm water discharge prototype measurement, comprising:
Meteorologic parameter during warm water discharge prototype measurement after the engineering operation, hydrologic parameter are interpolated into the waters model net Lattice obtain the third input boundary condition of the mathematical model；
Source sink term assignment is fetched water and drained to the warm water discharge to the waters model meshes, obtains the temperature row of the mathematical model The input of water source sink term；
The third of the mathematical model is inputted into boundary condition and warm water discharge source sink term input composition first input side Boundary's condition.
3. method according to claim 1 or claim 2, which is characterized in that described to utilize temperature after the engineering operation obtained by satellite Drain meteorologic parameter, hydrologic parameter, warm water discharge water intaking and the draining source sink term and preset waters model net during prototype measurement Lattice determine the first input boundary condition of the mathematical model after engineering operation during warm water discharge prototype measurement, before further include:
Grid dividing is carried out to the waters to be detected, obtains the waters model meshes.
4. method according to claim 2, which is characterized in that it is described will be during warm water discharge prototype measurement after the engineering operation Meteorologic parameter, hydrologic parameter be interpolated into the waters model meshes, obtain the mathematical model third input boundary condition, Include:
Meteorologic parameter during warm water discharge prototype measurement after the engineering operation, hydrologic parameter are subjected to document format conversion, obtained To text file needed for mathematical model input；
The text file is interpolated into the waters model meshes, obtains the third input boundary condition of the mathematical model.
5. method according to claim 1, which is characterized in that described according to the first input boundary condition, engineering operation The actual measurement hydrodynamic force during warm water discharge prototype measurement and water temperature obtain the number as a result, to mathematical model progress calibration afterwards Learn the Reasonable Parameters value of model, comprising:
The first input boundary condition is input to the mathematical model to solve, obtain calculating hydrodynamic force distribution and is calculated Distribution of water temperature；
It compares the calculating hydrodynamic force distribution and whether actual measurement hydrodynamic force distribution is consistent, and compare the calculating distribution of water temperature and reality Whether consistent survey distribution of water temperature；
If the calculating hydrodynamic force distribution is consistent with the actual measurement hydrodynamic force distribution, and the calculating distribution of water temperature and the actual measurement Distribution of water temperature is consistent, then the parameter value of the mathematical model during warm water discharge prototype measurement after the engineering operation is determined as institute State the Reasonable Parameters value of mathematical model.
6. method according to claim 2, which is characterized in that the Reasonable Parameters value and work according to the mathematical model Second input boundary condition of the mathematical model after Cheng Yunhang during warm water discharge prototype measurement, calculates the background water in waters to be detected Temperature, comprising:
The warm water discharge source sink term input for deleting the mathematical model in the first input boundary condition obtains the engineering fortune Second input boundary condition of the mathematical model after row during warm water discharge prototype measurement；
By the second input boundary condition and the Reasonable Parameters value of the mathematical model, it is input to the mathematical model and carries out It solves, obtains the background water temperature after the engineering operation during warm water discharge prototype measurement.
7. a kind of warm water discharge prototype measurement waters background water temperature computing device, which is characterized in that described device includes:
Determining module, for utilizing the meteorologic parameter during warm water discharge prototype measurement after the engineering operation obtained by satellite, water Literary parameter, warm water discharge water intaking and draining source sink term and preset waters model meshes, warm water discharge prototype is seen after determining engineering operation First input boundary condition of the mathematical model during survey；
Mathematical model calibration module, for according to the warm water discharge prototype measurement phase after the first input boundary condition, engineering operation Between actual measurement hydrodynamic force and water temperature as a result, carrying out calibration to the mathematical model, the Reasonable Parameters for obtaining the mathematical model take Value；Wherein, the mathematical model includes continuity equation, the equation of momentum, transport of substances equation and heat exchange equation, for simulating water Hydrodynamic force and the distribution of water temperature rule in domain；
Computing module, for during warm water discharge prototype measurement after the Reasonable Parameters value and engineering operation according to the mathematical model Mathematical model second input boundary condition, calculate the background water temperature in waters to be detected；Wherein, the second input perimeter strip Part is meteorologic parameter, hydrologic parameter and the waters mould after the engineering operation obtained according to satellite during warm water discharge prototype measurement Type grid is determined.
8. a kind of computer equipment, including memory and processor, be stored on the memory to run on a processor Computer program, which is characterized in that the processor realizes any one of claims 1 to 6 when executing the computer program The step of the method.
9. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program quilt The step of processor realizes method described in any one of claims 1 to 6 when executing.
CN201810929370.1A 2018-08-15 2018-08-15 Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium Active CN109255158B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810929370.1A CN109255158B (en) 2018-08-15 2018-08-15 Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810929370.1A CN109255158B (en) 2018-08-15 2018-08-15 Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium

Publications (2)

Publication Number Publication Date
CN109255158A true CN109255158A (en) 2019-01-22
CN109255158B CN109255158B (en) 2019-09-17

Family

ID=65049856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810929370.1A Active CN109255158B (en) 2018-08-15 2018-08-15 Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium

Country Status (1)

Country Link
CN (1) CN109255158B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110298080A (en) * 2019-05-30 2019-10-01 中国船舶重工集团公司第七一九研究所 Floating nuclear power plant warm water discharge thermal diffusion method for numerical simulation based on CFD
CN111914496A (en) * 2020-08-15 2020-11-10 四川大学 Frozen reservoir water temperature-ice condition inversion and prediction method based on thermal coupling model
CN113740371A (en) * 2021-09-15 2021-12-03 中国科学院空天信息创新研究院 Method, device and computer readable storage medium for determining water body thermal pollution
CN113790822A (en) * 2021-09-15 2021-12-14 中国科学院空天信息创新研究院 Method and device for detecting abnormity of ground measured temperature data and readable storage medium
CN115147239A (en) * 2022-06-10 2022-10-04 自然资源部第一海洋研究所 Calculation method for environmental background water temperature of newly-built coastal power plant temperature rise area and temperature rise calculation method
CN115144095A (en) * 2022-09-06 2022-10-04 自然资源部第一海洋研究所 Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing
CN115796076A (en) * 2022-11-30 2023-03-14 生态环境部卫星环境应用中心 Nuclear power warm water drainage reference temperature determination method and device fused with hydrodynamic model

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104751252A (en) * 2015-04-18 2015-07-01 国家***第一海洋研究所 Method for coastal power plant thermal discharge prediction by utilizing actually-measured surface heat transfer coefficient

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104751252A (en) * 2015-04-18 2015-07-01 国家***第一海洋研究所 Method for coastal power plant thermal discharge prediction by utilizing actually-measured surface heat transfer coefficient

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PING ZENG等: "Transport of waste heat from a nuclear power plant into coastal water", 《COASTAL ENGINEERING》 *
刘思等: "某滨海核电厂原型观测环境水体温升分布研究", 《给水排水》 *
刘思等: "温排水数值预报模型校验分析研究", 《给水排水》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110298080A (en) * 2019-05-30 2019-10-01 中国船舶重工集团公司第七一九研究所 Floating nuclear power plant warm water discharge thermal diffusion method for numerical simulation based on CFD
CN111914496A (en) * 2020-08-15 2020-11-10 四川大学 Frozen reservoir water temperature-ice condition inversion and prediction method based on thermal coupling model
CN113740371A (en) * 2021-09-15 2021-12-03 中国科学院空天信息创新研究院 Method, device and computer readable storage medium for determining water body thermal pollution
CN113790822A (en) * 2021-09-15 2021-12-14 中国科学院空天信息创新研究院 Method and device for detecting abnormity of ground measured temperature data and readable storage medium
CN113740371B (en) * 2021-09-15 2022-11-18 中国科学院空天信息创新研究院 Method, device and computer readable storage medium for determining water body thermal pollution
CN115147239A (en) * 2022-06-10 2022-10-04 自然资源部第一海洋研究所 Calculation method for environmental background water temperature of newly-built coastal power plant temperature rise area and temperature rise calculation method
WO2023237134A1 (en) * 2022-06-10 2023-12-14 自然资源部第一海洋研究所 Environmental background water temperature estimating method and temperature rise calculation method for temperature rise area of newly built coastal power plant
CN115144095A (en) * 2022-09-06 2022-10-04 自然资源部第一海洋研究所 Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing
CN115144095B (en) * 2022-09-06 2022-11-11 自然资源部第一海洋研究所 Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing
US11830635B1 (en) 2022-09-06 2023-11-28 First Institute of Oceanography, Ministry of Natural Resources Method and system for determining background water temperature of thermal discharge from operating nuclear power plants based on remote sensing
CN115796076A (en) * 2022-11-30 2023-03-14 生态环境部卫星环境应用中心 Nuclear power warm water drainage reference temperature determination method and device fused with hydrodynamic model
CN115796076B (en) * 2022-11-30 2023-07-21 生态环境部卫星环境应用中心 Nuclear power temperature drainage reference temperature determining method and device fused with hydrodynamic model

Also Published As

Publication number Publication date
CN109255158B (en) 2019-09-17

Similar Documents

Publication Publication Date Title
CN109255158B (en) Warm water discharge prototype measurement waters background water temperature calculation method, device and storage medium
Schenkel et al. An examination of tropical cyclone position, intensity, and intensity life cycle within atmospheric reanalysis datasets
Kuiry et al. Coupled 1D–Quasi-2D flood inundation model with unstructured grids
Piroozmand et al. Coupled CFD framework with mesoscale urban climate model: Application to microscale urban flows with weak synoptic forcing
Garcia-Dorado et al. Fast weather simulation for inverse procedural design of 3d urban models
Wei et al. Coupling of a regional atmospheric model (RegCM3) and a regional oceanic model (FVCOM) over the maritime continent
Bryan et al. Isohaline salinity budget of the North Atlantic salinity maximum
Mirzabozorg et al. Mathematical modeling and numerical analysis of thermal distribution in arch dams considering solar radiation effect
Twigt et al. Coupled 1D–3D hydrodynamic modelling, with application to the Pearl River Delta
Lauritzen et al. CAM-SE–CSLAM: Consistent coupling of a conservative semi-Lagrangian finite-volume method with spectral element dynamics
Diansky et al. Numerical simulation of the Caspian sea circulation using the marine and atmospheric research system
Al Mohit et al. A numerical estimate of water level elevation due to a cyclone associated with a different landfall angle
Telišman Prtenjak et al. Impact of mesoscale meteorological processes on anomalous radar propagation conditions over the northern Adriatic area
Zhang et al. A two-way nested coupled tide-surge model for the Taiwan Strait
Jiang et al. Estimating smoothly varying open boundary conditions for a 3D internal tidal model with an improved independent point scheme
Zhang et al. A method on estimating time-varying vertical eddy viscosity for an Ekman layer model with data assimilation
Vashist et al. Minimisation of overestimation of river flows in 1D-hydrodynamic modeling
Zhang et al. Variational estimation of wave-affected parameters in a two-equation turbulence model
Fujisaki et al. Improvement of short-term sea ice forecast in the southern Okhotsk Sea
Olsen et al. Mesoscale and microscale downscaling for the Wind Atlas of Mexico (WAM) project
Wong et al. A compressible nonhydrostatic cell-integrated semi-Lagrangian semi-implicit solver (CSLAM-NH) with consistent and conservative transport
Rasheed et al. Wind farm modeling in a realistic environment using a multiscale approach
Sainte-Rose et al. Monitoring and Performance Evaluation of Plastic Cleanup Systems: Part II—Results and Analysis
Kiss et al. Wind profile and shear stress at reed-open water interface–recent research achievements in Lake Fertő
Mori et al. Future changes in extreme storm surge based on a maximum potential storm surge model for East Asia

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Chen Xiaoli

Inventor after: Ji Ping

Inventor after: Zhao Yijun

Inventor after: Yuan Jue

Inventor after: Qin Xiao

Inventor after: Zhang Qiang

Inventor after: Liu Zanqiang

Inventor after: Liu Yan

Inventor before: Chen Xiaoli

CB03 Change of inventor or designer information