CN109145444A - Rainfall is obtained to the method and apparatus of Effect of Underground Water Level specific gravity - Google Patents
Rainfall is obtained to the method and apparatus of Effect of Underground Water Level specific gravity Download PDFInfo
- Publication number
- CN109145444A CN109145444A CN201810961180.8A CN201810961180A CN109145444A CN 109145444 A CN109145444 A CN 109145444A CN 201810961180 A CN201810961180 A CN 201810961180A CN 109145444 A CN109145444 A CN 109145444A
- Authority
- CN
- China
- Prior art keywords
- level
- rainfall
- wavelet
- ground water
- curve
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of acquisition rainfalls to the method and apparatus of Effect of Underground Water Level specific gravity, belongs to fwaater resources protection field.The described method includes: carrying out wavelet transform process to annual rainfall original series, rainfall matrix of wavelet coefficients is obtained；Wavelet transform process is carried out to level of ground water original series, obtains level of ground water matrix of wavelet coefficients.Rainfall wavelet variance curve is established according to rainfall matrix of wavelet coefficients, level of ground water wavelet variance curve is established according to level of ground water matrix of wavelet coefficients.According to rainfall wavelet variance curve and level of ground water wavelet variance curve, obtaining rainfall is more than the constant time range of predeterminable level as driving domain to Effect of Underground Water Level degree.The ratio that the sum of the level of ground water wavelet variance in driving domain accounts for level of ground water wavelet variance summation is obtained, ratio is determined as influence specific gravity of the rainfall to level of ground water.The present invention can obtain influence specific gravity of the rainfall to level of ground water, and then effectively instruct the formulation of regional science water use strategy.
Description
Technical field
The present invention relates to fwaater resources protection fields, in particular to a kind of to obtain rainfall to the method for Effect of Underground Water Level specific gravity
And device.
Background technique
Groundwater level fluctuation is cleared different factors by using rainfall as the meteorological factor of representative and mankind's activity double influence
Influence specific gravity to level of ground water, can effective guidance water.For example, if influence ratio of the somewhere rainfall to level of ground water
80% or more is focused on, illustrates that this area's level of ground water is mainly influenced by rainfall, the effect of human activity is weaker, by regulating and controlling people
Class activity intensity is unable to improve the situation of groundwater level fluctuation.If influence specific gravity of the somewhere rainfall to level of ground water 30% with
Under, illustrate that this area's level of ground water depends primarily on mankind's activity water, it can be effective by Reasonable Regulation And Control strength of human activities
Level of ground water is prevented to reduce.Therefore, it is necessary to provide a kind of acquisition rainfall to the method for Effect of Underground Water Level specific gravity.
Summary of the invention
The embodiment of the invention provides a kind of acquisition rainfalls to the method and apparatus of Effect of Underground Water Level specific gravity, can solve
State technical problem.Specific technical solution is as follows:
A kind of acquisition rainfall is provided to the method for Effect of Underground Water Level specific gravity, which comprises
Wavelet transform process is carried out to annual rainfall original series, obtains rainfall matrix of wavelet coefficients, the annual rainfall
Amount original series include the rainfall of history measurement, the rainfall matrix of wavelet coefficients be rainfall wavelet coefficient and time and
Three-dimensional relationship curved surface between period；
Wavelet transform process is carried out to level of ground water original series, obtains level of ground water matrix of wavelet coefficients, the underground
Water level original series include the level of ground water of history measurement, and the level of ground water matrix of wavelet coefficients is level of ground water wavelet coefficient
Three-dimensional relationship curved surface between time and period；
Rainfall wavelet variance curve is established according to the rainfall matrix of wavelet coefficients, the rainfall wavelet variance is bent
Line is rainfall wavelet variance and the relation curve between the period；
Level of ground water wavelet variance curve, the level of ground water small echo are established according to the level of ground water matrix of wavelet coefficients
Variance curve is level of ground water wavelet variance and the relation curve between the period；
According to the rainfall wavelet variance curve and the level of ground water wavelet variance curve, rainfall is obtained to underground
Water level influence degree is more than the constant time range of predeterminable level as driving domain；
The ratio that the sum of the level of ground water wavelet variance in the driving domain accounts for level of ground water wavelet variance summation is obtained, it will
The ratio is determined as influence specific gravity of the rainfall to level of ground water.
It is described that wavelet transform process is carried out to annual rainfall original series in a kind of possible design, obtain rainfall
Matrix of wavelet coefficients, comprising:
Each rainfall in the annual rainfall original series is standardized, rainfall is obtained and standardizes sequence
Column, the rainfall standardization sequence includes the first standard value being standardized to each rainfall；
Wavelet transformation is carried out to the rainfall standardised series, obtains the rainfall matrix of wavelet coefficients.
It is described that wavelet transform process is carried out to level of ground water original series in a kind of possible design, obtain underground water
Position matrix of wavelet coefficients, comprising:
Each level of ground water in the level of ground water original series is standardized, level of ground water standard is obtained
Change sequence, the level of ground water standardization sequence includes the second mark being standardized to each level of ground water
Quasi- value；
Wavelet transformation is carried out to the level of ground water standardised series, obtains the level of ground water matrix of wavelet coefficients.
It is described that rainfall wavelet variance curve packet is established according to rainfall matrix of wavelet coefficients in a kind of possible design
It includes:
The square value of the corresponding wavelet coefficient of the rainfall matrix of wavelet coefficients is accumulated by integration variable of the time
Point, obtain corresponding wavelet variance of each period；
According to the corresponding wavelet variance of each periodic quantity in the period, it is bent to establish the rainfall wavelet variance
Line.
It is described that level of ground water wavelet variance song is established according to level of ground water matrix of wavelet coefficients in a kind of possible design
Line includes:
The square value of the corresponding wavelet coefficient of the level of ground water matrix of wavelet coefficients is carried out by integration variable of the time
Integral, obtains corresponding wavelet variance of each period；
According to the corresponding wavelet variance of each periodic quantity in the period, the level of ground water wavelet variance is established
Curve.
It is described according to the rainfall wavelet variance curve and the level of ground water small echo side in a kind of possible design
Poor curve, obtaining rainfall is more than the constant time range of predeterminable level as driving domain to Effect of Underground Water Level degree, comprising:
N period 1 section is determined according to the rainfall wavelet variance curve, it is every in the n period 1 section
Line segment in a period 1 section is the wave crest section of the rainfall wavelet variance curve；
M second round section is determined according to the level of ground water wavelet variance curve, in the m second round section
Line segment in each second round section is the wave crest section of the level of ground water wavelet variance curve；
Period 1 section in the presence of overlapping and the is found out from n period 1 section and m second round section
The part being overlapped between the period 1 found out section and second round section is determined as resonance cycle section by two cycles section；
To the corresponding rainfall wavelet coefficient curve of periodic quantity each in the resonance cycle section and level of ground water wavelet systems
Number curve carries out correlation analysis, obtains the corresponding rainfall wavelet coefficient curve of each periodic quantity and level of ground water wavelet systems
The confidence level of the related coefficient of number curve, the rainfall wavelet coefficient curve are the rainfall matrix of wavelet coefficients last week
The corresponding curve of time value, the level of ground water wavelet coefficient curve are the corresponding song of periodic quantity on level of ground water matrix of wavelet coefficients
Line；
The confidence level of related coefficient in the resonance cycle section is chosen more than the collection cooperation that the periodic quantity of preset value forms
To drive domain.
It is described that n period 1 section, packet are determined according to the rainfall wavelet variance curve in a kind of possible design
It includes:
Each wave crest section is determined according to the rainfall wavelet variance curve, if the quantity of the wave crest section is greater than p, is chosen
For the corresponding constant time range of wave crest section of p as period 1 section, p is the first default wave crest quantity before height ranking；
If the quantity of the wave crest section is less than p, whole corresponding constant time ranges of wave crest section is chosen as period 1 section.
It is described that m second round section is determined according to the level of ground water wavelet variance curve in a kind of possible design,
Include:
Each wave crest section is determined according to the rainfall wavelet variance curve, if the quantity of the wave crest section is more than or equal to q,
For the corresponding constant time range of wave crest section of q as second round section, q is the second default wave crest quantity before selection height ranking；
If the quantity of the wave crest section is less than q, whole corresponding constant time ranges of wave crest section is chosen as second round section.
A kind of acquisition rainfall is provided to the device of Effect of Underground Water Level specific gravity, described device includes:
First conversion module obtains rainfall wavelet systems for carrying out wavelet transform process to annual rainfall original series
Matrix number, the annual rainfall original series include the rainfall of history measurement, and the rainfall matrix of wavelet coefficients is rainfall
Measure wavelet coefficient and the three-dimensional relationship curved surface between time and period；
Second conversion module obtains level of ground water small echo for carrying out wavelet transform process to level of ground water original series
Coefficient matrix, the level of ground water original series include the level of ground water of history measurement, the level of ground water matrix of wavelet coefficients
It is level of ground water wavelet coefficient and the three-dimensional relationship curved surface between time and period；
First establishes module, for establishing rainfall wavelet variance curve, institute according to the rainfall matrix of wavelet coefficients
Stating rainfall wavelet variance curve is rainfall wavelet variance and the relation curve between the period；
Second establishes module, bent for establishing level of ground water wavelet variance according to the level of ground water matrix of wavelet coefficients
Line, the level of ground water wavelet variance curve are level of ground water wavelet variance and the relation curve between the period；
First obtains module, for bent according to the rainfall wavelet variance curve and the level of ground water wavelet variance
Line, obtaining rainfall is more than the constant time range of predeterminable level as driving domain to Effect of Underground Water Level degree；
Second obtains module, for obtaining the level of ground water wavelet variance in driving domain and to account for level of ground water wavelet variance total
The ratio is determined as influence specific gravity of the rainfall to level of ground water by the ratio of sum.
In a kind of possible design, first conversion module, comprising:
First processing units, for being standardized to each rainfall in the annual rainfall original series,
It obtains rainfall and standardizes sequence, the rainfall standardization sequence includes being standardized to each rainfall
The first standard value arrived；
It is small to obtain the rainfall for carrying out wavelet transformation to the rainfall standardised series for first converter unit
Wave system matrix number.
Technical solution bring beneficial effect provided in an embodiment of the present invention includes at least:
Acquisition rainfall provided in an embodiment of the present invention can obtain rainfall to the method and apparatus of Effect of Underground Water Level specific gravity
The influence specific gravity to level of ground water is measured, and then effectively instructs to use water.For example, for rainfall to Effect of Underground Water Level specific gravity
Lesser region, mankind's activity water is larger to Effect of Underground Water Level, and underground water can effectively be prevented by changing strength of human activities
The reduction of position.For rainfall to the biggish region of Effect of Underground Water Level specific gravity, mankind's activity water to Effect of Underground Water Level compared with
It is small, improve the situation that mankind's activity can not be effectively improved level of ground water reduction.
Detailed description of the invention
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is one provided in an embodiment of the present invention and obtains rainfall to the flow chart of the method for Effect of Underground Water Level specific gravity；
Fig. 2 is that provided in an embodiment of the present invention another obtains rainfall to the process of the method for Effect of Underground Water Level specific gravity
Figure.
Specific embodiment
Unless otherwise defined, all technical terms used in the embodiment of the present invention all have usual with those skilled in the art
The identical meaning understood.To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to this hair
Bright embodiment is described in further detail.
Referring to Fig. 1, the embodiment of the invention provides a kind of method for obtaining rainfall to Effect of Underground Water Level specific gravity, this method
Include:
Step 101: wavelet transform process being carried out to annual rainfall original series, obtains rainfall matrix of wavelet coefficients, year
Rainfall original series include the rainfall year by year of history measurement, and rainfall matrix of wavelet coefficients is rainfall wavelet coefficient and year
Three-dimensional relationship curved surface between part and period.
Step 102: wavelet transform process is carried out to level of ground water original series, obtains level of ground water matrix of wavelet coefficients,
Level of ground water original series include the level of ground water of history measurement, level of ground water matrix of wavelet coefficients be rainfall wavelet coefficient with
Three-dimensional relationship curved surface between time and period.
Wherein, above-mentioned steps 101 and step 102 can carry out simultaneously, or first carry out step 101 and carry out step 102 again,
Or it first carries out step 102 and carries out step 101 again.
Step 103: rainfall wavelet variance curve being established according to rainfall matrix of wavelet coefficients, rainfall wavelet variance is bent
Line is rainfall wavelet variance and the relation curve between the period.
Step 104: level of ground water wavelet variance curve, level of ground water small echo are established according to level of ground water matrix of wavelet coefficients
Variance curve is level of ground water wavelet variance and the relation curve between the period.
Wherein, above-mentioned steps 103 and step 104 can carry out simultaneously, or first carry out step 103 and carry out step 104 again,
Or it first carries out step 104 and carries out step 103 again.
Step 105: according to rainfall wavelet variance curve and level of ground water wavelet variance curve, obtaining rainfall to underground
Water level influence degree is more than the constant time range of predeterminable level as driving domain.
Step 106: obtaining the ratio that the sum of the level of ground water wavelet variance in driving domain accounts for level of ground water wavelet variance summation
Value, is determined as influence specific gravity of the rainfall to level of ground water for ratio.
After carrying out the available rainfall matrix of wavelet coefficients of wavelet transformation to annual rainfall original series, wavelet coefficient square
The X-axis expression of years of battle array, Y-axis indicate the period, and Z axis indicates wavelet coefficient, to wavelet coefficient square using the time as integration variable
It is integrated, rainfall wavelet variance curve can be obtained, the horizontal axis of rainfall wavelet variance curve is the period, and the longitudinal axis is small echo
Variance.
Similarly, level of ground water matrix of wavelet coefficients can be obtained according to level of ground water original series, it is small according to level of ground water
Wave system matrix number obtains level of ground water wavelet variance curve.
According to rainfall wavelet variance curve and level of ground water wavelet variance curve, rainfall is selected to Effect of Underground Water Level
Biggish constant time range, as driving domain.Calculate the sum of level of ground water wavelet variance and level of ground water wavelet variance in driving domain
Summation, the sum of level of ground water wavelet variance and the ratio of level of ground water wavelet variance summation in domain will be driven to be determined as rainfall pair
The influence specific gravity of level of ground water.
As it can be seen that method provided in an embodiment of the present invention can obtain influence specific gravity of the rainfall to level of ground water, Jin Eryou
Imitate guidance water.For example, for rainfall to the lesser region of Effect of Underground Water Level specific gravity, mankind's activity water is to underground
Water level is affected, and regulation strength of human activities can effectively prevent the reduction of level of ground water.For rainfall to level of ground water
The biggish region of specific gravity is influenced, mankind's activity water is smaller to Effect of Underground Water Level, and improving mankind's activity can not be effectively improved ground
It is lauched the low situation of potential drop.
Referring to fig. 2, a kind of method the embodiment of the invention provides acquisition rainfall to Effect of Underground Water Level specific gravity, this method
Include:
Step 201: each rainfall in annual rainfall original series being standardized, rainfall standard is obtained
Change sequence, rainfall original series include the rainfall year by year of history measurement, and it includes to each rainfall that rainfall, which standardizes sequence,
Measure the first standard value being standardized.
Optionally, for the i-th annual rainfall in annual rainfall original series, the N of i=1,2,3 ...P, NPFor annual rainfall
The total number of annual rainfall in original series, that is, the length of annual rainfall original series.Wherein, annual rainfall original series
Length is greater than or equal to 30.For example, when choosing the rainfall of 40 years history measurement, NPIt is 40, to the rainfall by as follows
First formula is standardized to obtain the first standard value:
First formula:
In formula, G 'iFor the first standard value of the i-th annual rainfall, unit is millimeter；
GiFor the i-th annual rainfall in annual rainfall original series, unit is millimeter；
GmaxAnd GminMaximum value and minimum value respectively in annual rainfall original series, the unit of the two are millimeter.
After obtaining corresponding first standard value of each rainfall in annual rainfall original series, rainfall standardizes sequence
Including each first standard value.It is understood that length and the year of the rainfall standardization sequence after standardization
The length of rainfall original series is identical, is NP。
Step 202: wavelet transformation being carried out to rainfall standardised series, obtains rainfall matrix of wavelet coefficients.
Rainfall standardization sequence { P ' } is imported in MATLAB software, and as follows in the input of MATLAB software command window
Order carries out continuous wavelet transform, to obtain the matrix of wavelet coefficients of rainfall standardization sequence.
zP(xP,yP)=cwt (P ', 0.1:0.1:NP,'morl')
In formula, zP(xP,yP) it is rainfall wavelet coefficient, dimensionless, which is time xPWith period yPFunction；
NPThe length of sequence is standardized for rainfall, unit is year；
Cwt is continuous wavelet transform designator, and morl is Morlet fundamental wave symbol.
Step 203: each level of ground water in level of ground water original series being standardized, level of ground water is obtained
Sequence is standardized, level of ground water original series include the level of ground water of history measurement, and it includes to every that level of ground water, which standardizes sequence,
The second standard value that a level of ground water is standardized.
Optionally, for 1 year level of ground water in level of ground water original series, the N of i=1,2,3 ...G, NGFor underground water
The total number of level of ground water in the original series of position, that is, the length of level of ground water original series.Wherein, level of ground water original series
Length be greater than or equal to 30.For example, when choosing the level of ground water of 40 years history measurement, NGIt is 40, to the level of ground water
It is standardized to obtain the second standard value by following second formula:
Second formula:
In formula, P 'iFor the second standard value of the i-th annual rainfall, unit is rice；
PiFor 1 year level of ground water in level of ground water original series, unit is rice；
PmaxAnd PminMaximum value and minimum value respectively in level of ground water original series, the unit of the two are rice；
After obtaining corresponding second standard value of each level of ground water in level of ground water original series, level of ground water standardization
Sequence includes each second standard value.It is understood that the length of the level of ground water standardization sequence after standardization
Degree is identical as the length of level of ground water original series, is NG。
Step 204: wavelet transformation being carried out to level of ground water standardised series, obtains level of ground water matrix of wavelet coefficients.
The standardization sequence { G ' } of rainfall and level of ground water is imported in MATLAB software, and in MATLAB software command window
Mouth input is such as issued orders, and continuous wavelet transform is carried out, to obtain the matrix of wavelet coefficients of level of ground water standardization sequence.
zG(xG,yG)=cwt (G ', 0.1:0.1:NG,'morl')
In formula, zGFor the wavelet coefficient of level of ground water, dimensionless, which is time xGWith period yGFunction；
NGThe length of sequence, year are standardized for level of ground water；
Cwt is continuous wavelet transform designator；Morl is Morlet fundamental wave symbol.
Step 205, to the square value of the corresponding wavelet coefficient of rainfall matrix of wavelet coefficients using the time as integration variable into
Row integral obtains the corresponding wavelet variance of each periodic quantity in the period；
According to the corresponding wavelet variance of periodic quantity each in the period, rainfall wavelet variance curve is established.
Optionally, the square value of the corresponding wavelet coefficient of rainfall matrix of wavelet coefficients is carried out by integration variable of the time
Integral, can carry out according to third formula:
Third formula:
In formula, VP(yP) indicate rainfall wavelet variance, zP(xP,yP) indicate rainfall wavelet coefficient.
It is understood that the horizontal axis of rainfall wavelet variance curve is period yP, the longitudinal axis is rainfall wavelet variance VP
(yP)。
Step 206, to the square value of the corresponding wavelet coefficient of level of ground water matrix of wavelet coefficients using the time as integration variable
It is integrated, obtains the corresponding wavelet variance of each periodic quantity in the period；
According to the corresponding wavelet variance of periodic quantity each in the period, level of ground water wavelet variance curve is established.
Optionally, to the square value of the corresponding wavelet coefficient of level of ground water matrix of wavelet coefficients using the time as integration variable into
Row integral, can carry out according to the 4th formula:
4th formula:
In formula, VG(yG) indicate level of ground water wavelet variance, zG(xG,yG) indicate level of ground water wavelet coefficient.
It is understood that the horizontal axis of drop level of ground water wavelet variance curve is period yG, the longitudinal axis is rainfall small echo side
Poor VG(yG)。
Step 207, n period 1 section is determined according to rainfall wavelet variance curve, it is every in n period 1 section
Line segment in a period 1 section is the wave crest section of rainfall wavelet variance curve.
Optionally, determine that n period 1 section can be carried out by following manner according to rainfall wavelet variance curve:
Each wave crest section is determined according to rainfall wavelet variance curve, if the quantity of wave crest section is greater than p, chooses crest height
The corresponding constant time range of wave crest section of p is as period 1 section before ranking, at this point, n is equal to p；If the quantity of wave crest section is less than p, choosing
Take whole corresponding constant time ranges of wave crest section as period 1 section, at this point, n is equal to the quantity of wave crest section, p is the first default wave
Peak number amount.For example, p=3, when the quantity of wave crest section is greater than 3,3 wave crest section is corresponding before selection crest height ranking
Constant time range as the period 1 section.When the quantity of wave crest section is less than 3, the corresponding constant time range of whole wave crest sections is chosen as the
One constant time range.
Step 208, m second round section is determined according to level of ground water wavelet variance curve, in m second round section
Line segment in each second round section is the wave crest section of level of ground water wavelet variance curve.
Optionally, m second round section is determined according to level of ground water wavelet variance curve, comprising:
Each wave crest section is determined according to rainfall wavelet variance curve, if the quantity of wave crest section is more than or equal to q, chooses wave crest
The corresponding constant time range of wave crest section of q is as second round section before height ranking, at this point, m is equal to q；If the quantity of wave crest section is less than
Q chooses whole corresponding constant time ranges of wave crest section as period 1 section, at this point, m is equal to the quantity of wave crest section, q is second pre-
If wave crest quantity.For example, q=3 chooses before crest height ranking 3 wave crest section when the quantity of wave crest section is greater than 3
Corresponding constant time range is as second round section.When the quantity of wave crest section is less than 3, chooses the corresponding constant time range of whole wave crest sections and make
For second round section.
Step 209, period 1 section in the presence of overlapping and the is found out from n period 1 section and m second round section
The part being overlapped between the period 1 found out section and second round section is determined as resonance cycle section by two cycles section.
Period 1 section and second round section are the section of horizontal axis on coordinate system, find out the lap in two sections,
As resonance cycle section.
Step 210, small to the corresponding rainfall wavelet coefficient curve of periodic quantity each in resonance cycle section and level of ground water
Wave system number curve carries out correlation analysis, obtains the corresponding rainfall wavelet coefficient curve of each periodic quantity and level of ground water wavelet systems
The confidence level of the related coefficient of number curve, the confidence level for choosing related coefficient are gathered i.e. more than the periodic quantity composition of preset value
To drive domain.Rainfall wavelet coefficient curve is the corresponding curve of periodic quantity on rainfall matrix of wavelet coefficients, and level of ground water is small
Wave system number curve is the corresponding curve of level of ground water matrix of wavelet coefficients periodic quantity.
Wherein, correlation analysis is carried out, the confidence level for obtaining related coefficient can be by the correlation analysis work of SPSS software
Tool carries out.
Since resonance cycle section is the lap of period 1 section and second round section, in the period 1 section, Ge Gezhou
Time value has corresponding rainfall wavelet coefficient curve, has corresponding level of ground water wavelet coefficient curve in second round section, therefore,
In resonance cycle section, same period value has rainfall and level of ground water wavelet coefficient curve simultaneously.To in resonance cycle section, often
The corresponding rainfall wavelet coefficient curve of a periodic quantity and level of ground water wavelet coefficient curve carry out correlation analysis, available every
The confidence level of the related coefficient of the corresponding rainfall wavelet coefficient curve of a periodic quantity and level of ground water wavelet coefficient curve, choosing
The set for taking the confidence level of related coefficient to form more than the periodic quantity of preset value is to drive domain.
Wherein, preset value can be 0.01,0.05 etc..
Step 211, the ratio that the sum of the level of ground water wavelet variance in driving domain accounts for level of ground water wavelet variance summation is obtained
Value, is determined as influence specific gravity of the rainfall to level of ground water for ratio.
Wherein, obtaining the ratio for driving the sum of the level of ground water wavelet variance in domain to account for level of ground water wavelet variance summation can
To be carried out by following 5th formula:
5th formula:
In formula,Influence specific gravity for rainfall to level of ground water, %；To drive domain.
As it can be seen that method provided in an embodiment of the present invention can obtain influence specific gravity of the rainfall to level of ground water, Jin Eryou
Imitate guidance water.For example, for rainfall to the lesser region of Effect of Underground Water Level specific gravity, mankind's activity water is to underground
Water level is affected, and groundwater level fluctuation can effectively be prevented by regulating and controlling mankind's activity with water intensity.For rainfall to underground water
The biggish region of position influence specific gravity, mankind's activity water is smaller to Effect of Underground Water Level, and improving mankind's activity can not be effectively improved
The situation that level of ground water reduces.
The embodiment of the invention provides a kind of acquisition rainfalls to the device of Effect of Underground Water Level specific gravity, which includes:
First conversion module obtains rainfall wavelet systems for carrying out wavelet transform process to annual rainfall original series
Matrix number, annual rainfall original series include the rainfall year by year of history measurement, and rainfall matrix of wavelet coefficients is that rainfall is small
Wave system number and the three-dimensional relationship curved surface between time and period.
Second conversion module obtains level of ground water small echo for carrying out wavelet transform process to level of ground water original series
Coefficient matrix, level of ground water original series include the level of ground water of history measurement, and level of ground water matrix of wavelet coefficients is underground water
Position wavelet coefficient and the three-dimensional relationship curved surface between time and period；
First establishes module, for establishing rainfall wavelet variance curve, rainfall according to rainfall matrix of wavelet coefficients
Wavelet variance curve is rainfall wavelet variance and the relation curve between the period；
Second establishes module, for establishing level of ground water wavelet variance curve, ground according to level of ground water matrix of wavelet coefficients
Lower water level wavelet variance curve is level of ground water wavelet variance and the relation curve between the period；
First obtains module, for obtaining drop according to rainfall wavelet variance curve and level of ground water wavelet variance curve
Rainfall is more than the constant time range of predeterminable level as driving domain to Effect of Underground Water Level degree；
Second obtains module, and the sum of the level of ground water wavelet variance for obtaining in driving domain accounts for level of ground water wavelet variance
Ratio is determined as influence specific gravity of the rainfall to level of ground water by the ratio of summation.
Optionally, the first conversion module, comprising:
First processing units are obtained for being standardized to each rainfall in annual rainfall original series
Rainfall standardizes sequence, and it includes the first standard being standardized to each rainfall that rainfall, which standardizes sequence,
Value；
First converter unit obtains rainfall wavelet coefficient square for carrying out wavelet transformation to rainfall standardised series
Battle array.
Optionally, the second conversion module, comprising:
The second processing unit is obtained for being standardized to each level of ground water in level of ground water original series
It obtains level of ground water and standardizes sequence, it includes being standardized to each level of ground water that level of ground water, which standardizes sequence,
Second standard value；
Second converter unit obtains level of ground water wavelet systems for carrying out wavelet transformation to level of ground water standardised series
Matrix number.
Optionally, it first establishes module and includes:
First integral unit, for the square value to the corresponding wavelet coefficient of rainfall matrix of wavelet coefficients with the time be product
Variation per minute is integrated, and the corresponding wavelet variance of each periodic quantity in the period is obtained；
First establishing unit, for it is small to establish rainfall according to the corresponding wavelet variance of periodic quantity each in the period
Wave variance curve.
Optionally, it second establishes module and includes:
Second integral unit, for establishing level of ground water wavelet variance curve packet according to level of ground water matrix of wavelet coefficients
It includes:
The square value of the corresponding wavelet coefficient of level of ground water matrix of wavelet coefficients is integrated by integration variable of the time,
Obtain the corresponding wavelet variance of each periodic quantity in the period；
Second establishes unit, for establishing level of ground water according to the corresponding wavelet variance of periodic quantity each in the period
Wavelet variance curve.
Optionally, the first acquisition module includes:
First determination unit, it is first week a in n for determining n period 1 section according to rainfall wavelet variance curve
The line segment in section of each period 1 in phase section is the wave crest section of rainfall wavelet variance curve；
Second determination unit, for determining m second round section according to level of ground water wavelet variance curve, at m second
The line segment in section of each second round in constant time range is the wave crest section of level of ground water wavelet variance curve；
Third determination unit, for finding out first week in the presence of overlapping from n period 1 section and m second round section
The part being overlapped between the period 1 found out section and second round section is determined as resonance cycle by phase section and second round section
Section；
Analytical unit, for the corresponding rainfall wavelet coefficient curve of periodic quantity each in resonance cycle section and underground water
Position wavelet coefficient curve carries out correlation analysis, obtains the corresponding rainfall wavelet coefficient curve of each periodic quantity and level of ground water is small
The confidence level of the related coefficient of wave system number curve, the confidence level for choosing related coefficient is more than the periodic quantity composition collection of preset value
Closing is driving domain.Rainfall wavelet coefficient curve is the corresponding curve of periodic quantity on rainfall matrix of wavelet coefficients, underground water
Position wavelet coefficient curve is the corresponding curve of periodic quantity on level of ground water matrix of wavelet coefficients.
Selection unit, the set formed for choosing the periodic quantity that related coefficient reliability in resonance cycle section is more than preset value
As driving domain.
The above is only illustrative embodiments of the invention, the protection scope being not intended to limit the invention is all in the present invention
Spirit and principle within, any modification, equivalent replacement, improvement and so on, should be included in protection scope of the present invention it
It is interior.
Claims (10)
1. a kind of obtain rainfall to the method for Effect of Underground Water Level specific gravity, which is characterized in that the described method includes:
Wavelet transform process is carried out to annual rainfall original series, obtains rainfall matrix of wavelet coefficients, the annual rainfall is former
Beginning sequence includes the rainfall year by year of history measurement, the rainfall matrix of wavelet coefficients be rainfall wavelet coefficient and time and
Three-dimensional relationship curved surface between period；
Wavelet transform process is carried out to level of ground water original series, obtains level of ground water matrix of wavelet coefficients, the level of ground water
Original series include the level of ground water of history measurement, and the level of ground water matrix of wavelet coefficients is level of ground water wavelet coefficient and year
Three-dimensional relationship curved surface between part and period；
Rainfall wavelet variance curve is established according to the rainfall matrix of wavelet coefficients, the rainfall wavelet variance curve is
Rainfall wavelet variance and the relation curve between the period；
Level of ground water wavelet variance curve, the level of ground water wavelet variance are established according to the level of ground water matrix of wavelet coefficients
Curve is level of ground water wavelet variance and the relation curve between the period；
According to the rainfall wavelet variance curve and the level of ground water wavelet variance curve, rainfall is obtained to level of ground water
Influence degree is more than the constant time range of predeterminable level as driving domain；
The ratio that the sum of the level of ground water wavelet variance in the driving domain accounts for level of ground water wavelet variance summation is obtained, it will be described
Ratio is determined as influence specific gravity of the rainfall to level of ground water.
2. the method according to claim 1, wherein described carry out at wavelet transformation annual rainfall original series
Reason, obtains rainfall matrix of wavelet coefficients, comprising:
Each rainfall in the annual rainfall original series is standardized, rainfall is obtained and standardizes sequence,
The rainfall standardization sequence includes the first standard value being standardized to each rainfall；
Wavelet transformation is carried out to the rainfall standardised series, obtains the rainfall matrix of wavelet coefficients.
3. the method according to claim 1, wherein described carry out at wavelet transformation level of ground water original series
Reason, obtains level of ground water matrix of wavelet coefficients, comprising:
Each level of ground water in the level of ground water original series is standardized, level of ground water is obtained and standardizes sequence
Column, the level of ground water standardization sequence includes the second standard being standardized to each level of ground water
Value；
Wavelet transformation is carried out to the level of ground water standardised series, obtains the level of ground water matrix of wavelet coefficients.
4. the method according to claim 1, wherein described establish rainfall according to rainfall matrix of wavelet coefficients
Wavelet variance curve includes:
The square value of the corresponding wavelet coefficient of the rainfall matrix of wavelet coefficients is integrated by integration variable of the time, is obtained
To corresponding wavelet variance of each period；
According to each period corresponding wavelet variance, the rainfall wavelet variance curve is established.
5. the method according to claim 1, wherein described establish underground according to level of ground water matrix of wavelet coefficients
Water level wavelet variance curve includes:
The square value of the corresponding wavelet coefficient of the level of ground water matrix of wavelet coefficients is integrated by integration variable of the time,
Obtain corresponding wavelet variance of each period；
According to each period corresponding wavelet variance, the level of ground water wavelet variance curve is established.
6. the method according to claim 1, wherein according to the rainfall wavelet variance curve and the underground
Water level wavelet variance curve obtains rainfall and is used as driving domain to the constant time range that Effect of Underground Water Level degree is more than predeterminable level,
Include:
N period 1 section is determined according to the rainfall wavelet variance curve, in each of described n period 1 section the
Line segment in one constant time range is the wave crest section of the rainfall wavelet variance curve；
M second round section is determined according to the level of ground water wavelet variance curve, in each of described m second round section
Line segment in second round section is the wave crest section of the level of ground water wavelet variance curve；
The period 1 section and second week in the presence of overlapping are found out from the n period 1 section and m second round section
The part being overlapped between the period 1 found out section and second round section is determined as resonance cycle section by phase section；
It is bent to the corresponding rainfall wavelet coefficient curve of periodic quantity each in the resonance cycle section and level of ground water wavelet coefficient
Line carries out correlation analysis, obtains the corresponding rainfall wavelet coefficient curve of each periodic quantity and level of ground water wavelet coefficient is bent
The confidence level of the related coefficient of line, the rainfall wavelet coefficient curve are periodic quantity on the rainfall matrix of wavelet coefficients
Corresponding curve, the level of ground water wavelet coefficient curve are the corresponding curve of periodic quantity on level of ground water matrix of wavelet coefficients；
The confidence level for choosing related coefficient in the resonance cycle section is used as more than the set that the periodic quantity of preset value forms drives
Dynamic domain.
7. according to the method described in claim 6, it is characterized in that, described determine n according to the rainfall wavelet variance curve
A period 1 section, comprising:
Each wave crest section is determined according to the rainfall wavelet variance curve, if the quantity of the wave crest section is greater than p, chooses wave crest
For the corresponding constant time range of wave crest section of p as period 1 section, p is the first default wave crest quantity before height ranking；
If the quantity of the wave crest section is less than p, whole corresponding constant time ranges of wave crest section is chosen as period 1 section.
8. according to the method described in claim 6, it is characterized in that, described determine according to the level of ground water wavelet variance curve
M second round section, comprising:
Each wave crest section is determined according to the rainfall wavelet variance curve, if the quantity of the wave crest section is more than or equal to q, is chosen
For the corresponding constant time range of wave crest section of q as second round section, q is the second default wave crest quantity before crest height ranking；
If the quantity of the wave crest section is less than q, whole corresponding constant time ranges of wave crest section is chosen as second round section.
9. a kind of obtain rainfall to the device of Effect of Underground Water Level specific gravity, which is characterized in that described device includes:
First conversion module obtains rainfall wavelet coefficient square for carrying out wavelet transform process to annual rainfall original series
Battle array, the annual rainfall original series include the rainfall year by year of history measurement, and the rainfall matrix of wavelet coefficients is rainfall
Measure wavelet coefficient and the three-dimensional relationship curved surface between time and period；
Second conversion module obtains level of ground water wavelet coefficient for carrying out wavelet transform process to level of ground water original series
Matrix, the level of ground water original series include the level of ground water of history measurement, and the level of ground water matrix of wavelet coefficients is ground
Lower water level wavelet coefficient and the three-dimensional relationship curved surface between time and period；
First establishes module, for establishing rainfall wavelet variance curve, the drop according to the rainfall matrix of wavelet coefficients
Rainfall wavelet variance curve is annual rainfall wavelet variance and the relation curve between the period；
Second establishes module, for establishing level of ground water wavelet variance curve, institute according to the level of ground water matrix of wavelet coefficients
Stating level of ground water wavelet variance curve is level of ground water wavelet variance and the relation curve between the period；
First obtains module, for obtaining according to the rainfall wavelet variance curve and the level of ground water wavelet variance curve
Taking rainfall is more than the constant time range of predeterminable level as driving domain to Effect of Underground Water Level degree；
Second obtains module, and the sum of the level of ground water wavelet variance for obtaining in driving domain accounts for level of ground water wavelet variance summation
Ratio, the ratio is determined as influence specific gravity of the rainfall to level of ground water.
10. device according to claim 9, which is characterized in that first conversion module, comprising:
First processing units are obtained for being standardized to each rainfall in the annual rainfall original series
Rainfall standardizes sequence, and the rainfall standardization sequence includes being standardized to each rainfall
First standard value；
First converter unit obtains the rainfall wavelet systems for carrying out wavelet transformation to the rainfall standardised series
Matrix number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810961180.8A CN109145444B (en) | 2018-08-22 | 2018-08-22 | Method and device for acquiring influence proportion of rainfall on underground water level |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810961180.8A CN109145444B (en) | 2018-08-22 | 2018-08-22 | Method and device for acquiring influence proportion of rainfall on underground water level |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109145444A true CN109145444A (en) | 2019-01-04 |
| CN109145444B CN109145444B (en) | 2021-04-23 |
Family
ID=64790826
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810961180.8A Active CN109145444B (en) | 2018-08-22 | 2018-08-22 | Method and device for acquiring influence proportion of rainfall on underground water level |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109145444B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10177076A (en) * | 1996-12-19 | 1998-06-30 | Toshiba Corp | Amount-of-rainfall predicting device |
| CN101824808A (en) * | 2010-05-26 | 2010-09-08 | 北京师范大学 | Ecological water supplement technology for lakes based on multi-scale wavelet transform |
| CN102033851A (en) * | 2010-09-30 | 2011-04-27 | 南京大学 | Hydrologic time series wavelet correlation analysis method |
-
2018
- 2018-08-22 CN CN201810961180.8A patent/CN109145444B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10177076A (en) * | 1996-12-19 | 1998-06-30 | Toshiba Corp | Amount-of-rainfall predicting device |
| CN101824808A (en) * | 2010-05-26 | 2010-09-08 | 北京师范大学 | Ecological water supplement technology for lakes based on multi-scale wavelet transform |
| CN102033851A (en) * | 2010-09-30 | 2011-04-27 | 南京大学 | Hydrologic time series wavelet correlation analysis method |
Non-Patent Citations (1)
| Title |
|---|
| 迟光耀: "基于小波分析与Mann_Kendall法的岩溶大泉动态研究", 《中国岩溶》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109145444B (en) | 2021-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20060116792A1 (en) | Irrigation controller | |
| CN108416468B (en) | Flood advanced early warning and forecasting method | |
| CN107169276B (en) | Method for determination of amount and device are adopted in a kind of underground water control | |
| CN109765916A (en) | A kind of unmanned surface vehicle path following control device design method | |
| Kamal et al. | GIS-based irrigation water management for precision farming of rice | |
| CN102526879A (en) | Device and system for adaptively adjusting pacing frequency | |
| CN103808380A (en) | Flow rapid tracking method for ultrasonic flow metering device | |
| CN104897245A (en) | Zigbee-based multi-point acquisition self-correction water level detector and measuring method | |
| CN108008151A (en) | A kind of moving state identification method and system based on 3-axis acceleration sensor | |
| CN109145444A (en) | Rainfall is obtained to the method and apparatus of Effect of Underground Water Level specific gravity | |
| Siva et al. | Smart watering of plants | |
| CN206637942U (en) | Refrigerator image capturing system and refrigerator | |
| CN104298806A (en) | Hydropower station dynamic property computer-assisted testing method | |
| CN116592958A (en) | Method and device for synchronously monitoring water quality and water quantity | |
| CN109117457A (en) | Bodily form curvature big data calculation method | |
| CN206497347U (en) | A kind of fertigation control system controlled based on PLC | |
| CN115363606A (en) | Pelvic floor muscle and abdominal muscle coordination analysis method, device and equipment | |
| Norsworthy et al. | Productivity and costs in the private economy, 1973 | |
| CN104977556B (en) | The average waveform capture rate method of testing measured based on dead time | |
| Jobbágy et al. | Evaluation of efficiency of precision irrigation for potatoes | |
| CN105699981A (en) | Hardware and software cooperation method for measuring cloud layer height and thickness | |
| CN107808304A (en) | A kind of utility model patent based on big data is worth pre- evaluation method | |
| CN207301758U (en) | A kind of earth source heat pump unit control system controlled based on fuzzy control and PLC | |
| EP1203936B1 (en) | Device for evaluating, with acoustic waves, the temperature of a room | |
| CN110675045A (en) | Environment assessment method and device based on winter easy-to-maintain |
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 |