CN106802962A - A kind of method for calculating and correcting ocean essential average - Google Patents

Abstract

The invention discloses a kind of method for calculating and correcting ocean essential average, the method is based on oceanography Period Process and the regular understanding of ocean essential urban agglomeration, with wavelet analysis technology, the period of change of ocean essential multiple different time scales is decomposed, the mechanical periodicity component (periodic component) in multiple time scales is extracted；There is the property of statistically uniformity with the cycle based on ocean essential with position phase spatial and temporal distributions situation, space-time correction is carried out to partially different marine site ocean essential survey data, it is used to calculate the spatial domain average of target marine site each periodic component of ocean essential, using the time change of each periodic component spatial domain average of Model fitting, and then by the time domain average of time integral calculating ocean essential periodic component, then by the multifrequency period average of target marine site ocean essential in each periodic component time domain mean value computation minimum positive period.

Description

A kind of method for calculating and correcting ocean essential average

Technical field

The present invention relates to a kind of method for calculating and correcting target marine site ocean essential average, more particularly to survey data There is a kind of method of the partially different Marine Chemistry key element mean value computation of space-time and correction, specifically using the method for wavelet analysis, Period of change on analysis ocean essential Multiple Time Scales, decomposes periodic component of the ocean essential in multiple time scales, profit With the change of ocean essential periodic component in Model fitting each time scale, by partially different Sea Area Survey data estimation target Marine site ocean essential multifrequency period average.

Background technology

China has carried out a series of marine environmental qualities investigation or has monitored, its medium-and-large-sized sea since the late 1950s Foreign environmental survey mainly including late 1950s first time the whole nation comprehensive marine investigation, the mid-1970s first time it is complete The whole nation that state marine pollution baseline visit, the national littoral zone at the beginning of the eighties and tidal flat resources comprehensive survey, the beginning of the nineties carry out Second National marine pollution baseline visit and China of development in 2006 that sea island resources comprehensive survey, the end of the nineties carry out Coastal ocean comprehensive survey and evaluation etc..Marine site is polluted for emphasis, National Bureau of Oceanography and other relevant departments also successively organize Multiple ad hoc survey is carried out and great research is special, such as the pollution surveys of East Sea coastal waters, 20th century 80 at the end of 20th century 70 The estuary area and its Adjacent Sea Area pollutant of mid-nineties 90 and the biogeochemistry of nutritive salt, in the nineties in 20th century The yellow Bohai environment quality of phase and Ecological Changes research, migration, the conversion process of Bohai Sea typical environment load at the beginning of 21 century And environmental carrying capacity research etc..Additionally, National Bureau of Oceanography also successively carries out in typical marine site since middle 1960s The long-term monitoring plans such as ocean standard section is investigated, the monitoring of marine pollution tendency；National environmental protection is total since 20 century 70s Main coastal cities environmental administration has also carried out the monitoring of stretch of coastal water marine environmental quality in succession belonging to office.By these investigate or Monitoring, obtains valuable data in terms of the conventional Marine Environmental Elements of the subjects such as the hydrology, chemistry, biology, geology.However, from tune Look into from the point of view of coverage, the only Bohai Sea has and few in number is related to the ocean essential distributed data data in whole marine site, the Huanghai Sea, east Sea, the ocean essential Distribution Data overwhelming majority at the South Sea are only confined in some importance or typical marine site；Additionally, these investigation or The task of monitoring and executable unit are different, and most of arrangement of control time and erect-position lacks systematicness and comparativity, thus difficult Valuable regularity cognition is obtained on a larger time scale with using these data.Under the situation of whole world change, utilize Historical summary, obtains ocean essential numerical value change trend from historical summary, and the differentiation for studying marine environment has important Meaning.Effective ways of the network analysis using these valuable survey data are highly desirable at present.

Marine environmental quality management needs the average value of hard objectives marine site ocean essential.Target marine site ocean will in theory Plain average value should be equal to the ocean essential space integral amount of storing divided by seawater cumulative volume, i.e. spatial domain average.Using numerical solution by When the ocean essential space integral amount of storing calculates target marine site ocean essential spatial domain average, investigation station's bit quantity is not required nothing more than enough It is many, and require that its space layout is relatively uniform.But in actual marine environmental quality is investigated or is monitored, investigation erect-position is not only Limited amount, space layout also tends to uneven, it is necessary to application space Discretization Procedure estimation target marine site ocean essential Spatial domain average.However, different historical summaries not only investigate erect-position sets inconsistent, the field of investigation is also not quite similar, and is and target Marine site is distinguished, referred to as partially different marine site, and the ocean essential spatial domain average in target marine site and partially different marine site lacks comparativity, thus needs Carry out free-air correction.

Marine environmental quality management not only needs the spatial domain average of target marine site ocean essential under clear and definite special time, also needs The time trend of GPRS spatial domain average, is normally applied the mathematic(al) mean of multiple survey data in certain hour yardstick at present Value characterizes the average in the range of ocean essential certain hour.However, the time scale of selection is different, the average value of gained is often deposited In larger difference, the comparativity of the time scale average value of different time domain is reduced.

Due to oceanography processes such as tidal motion, ocean circulation, atmospheric sedimentation, runoff input and biogeochemical processes Periodicity, ocean essential show to a certain extent the multifrequency periods such as diurnal periodicity, the cycle moon, annual period, climatic cycle change Change, ocean essential numerical value also tends to periodically reproduction.If according to the time scale reasonable selection of ocean essential mechanical periodicity The time scale of its time domain average value is calculated, the time domain average of gained should have bigger stability and comparativity.The time domain Average should be spatially the average value of certain marine site ocean essential space integral amount of storing, and should be in time certain hour model Enclose the average value of the interior ocean essential amount of storing time integral.By ocean in target marine site (covering multiple cycles) minimum positive period The average value of the key element amount of storing mesh integrator referred to as multifrequency period average.It is equal that target marine site ocean essential multifrequency is calculated in theory Under conditions of value is, it is necessary to time interval is consistent, with the certain frequency METHOD FOR CONTINUOUS DETERMINATION ocean essential instantaneous value, then integration is asked Average.Target marine site ocean essential multifrequency average is calculated using numerical solution, it is sufficiently high not require nothing more than investigation frequency, and It is required that control time layout is relatively uniform.However, in actual marine environmental quality is investigated or is monitored, not only investigating frequency not to the utmost It is identical, control time interval it is also inconsistent, or even discontinuously, it is difficult to apply numerical solution.

The inventive method can control basic law according to ocean essential spatial and temporal distributions, make full use of history survey data, To being corrected analysis in the presence of a series of partially different historical summaries on investigation marine site, on control time, the ocean essential of gained is more Frequency Periodic Mean can more accurately reflect the real conditions of target marine site ocean essential, and this becomes for disclosing target marine site ocean essential Law, the especially variation tendency of China's ocean essential in recent decades, deepen the understanding of ocean essential time space distribution, With important scientific meaning and application value.

The content of the invention

The present invention provides a kind of method for calculating and correcting ocean essential average, based on ocean essential spatial and temporal distributions situation Regularity, with wavelet analysis technology, the means such as binding model fitting, in time domain and the different Sea Area Survey money of spatial domain colonel's positively biased Material, the multifrequency period changed according to ocean essential selects the time scale of time domain mean value computation, calculates target marine site ocean essential Multifrequency period average, improves the stability and comparativity of ocean essential time domain average.

The technical solution adopted by the present invention is：It is a kind of calculate and correction ocean essential average method, the method include with Lower step：

6. target marine site ocean essential space lattice discretization

For the target marine site that ocean essential survey data has high-spatial and temporal resolution, based on Kriging optimal interpolation methods, Using Matlab software programming three-dimensional interpolation programs, by target marine site ocean essential erect-position survey data space lattice discretization, The target marine site ocean essential gridded data data of certain hour sequence is obtained, target marine site sea under each control time is calculated Foreign essential factors space average value, obtains the time series data of target marine site ocean essential spatial domain average.

7. ocean essential multifrequency period analysis

With wavelet analysis technology analysis target marine site ocean essential gridded data time series data, decomposing ocean will Periodic component (P1, P2, P3, P4 ...) on the different time scales of element change, obtains the period of change of multiple time scales (T1, T2, T3, T4 ...), such as diurnal periodicity, the cycle moon, annual period, climatic cycle.

8. the correction of partially different marine site ocean essential periodic component grid values and spatial domain average

For the partially different marine site for only being partially overlapped with target marine site, it is partially different to there is space-time in its survey data and target marine site, Need correction.Ocean essential standard point of reference is selected in target marine site first, by target marine site ocean essential periodic component grid Change numerical value all divided by the numerical value (formula 1) of standard point of reference, obtain the normalization grid values of ocean essential periodic component：

Wherein,It is target marine site ocean essential periodic component normalization grid value set,It is target sea Domain ocean essential periodic component grid value set,It is ocean essential periodic component standard basis point value.

Equally, by identical control time, the partially different marine site ocean essential periodic component gridding of same period, select it empty Between in the range of central gridding as partially different datum mark, by the data after gridding all divided by partially different benchmark point value (formula 2), obtain The normalization grid values of partially different marine site ocean essential periodic component：

Wherein,It is partially different marine site ocean essential periodic component normalization grid value set,It is partially different sea Domain ocean essential periodic component grid value set,It is the partially different benchmark point value of ocean essential periodic component.

In the range of the overlapping region in target marine site and partially different marine site, by the ocean essential periodic component normalization of target marine site The average of grid valuesThe average of grid values is normalized divided by partially different marine site ocean essential periodic component(formula 3), Obtain normalization coefficient (γ)：

Then by the normalization grid value set of entirely partially different marine site ocean essential periodic componentIt is multiplied by correction Coefficient (γ) (formula 4), the partially different marine site ocean essential periodic component normalization grid value set put on the basis of partially different datum mark Translate into the ocean essential periodic component normalization grid set of correction values put on the basis of standard point of reference

By target marine site ocean essential periodic component normalization grid value setPartially different marine site ocean essential Periodic component normalizes grid set of correction valuesIt is referred to as " ocean essential periodic component field intensity " (ψ (E^BZ)).Such as The partially different marine site of fruit and target marine site do not exist overlapping region, can first be corrected to target marine site have overlapping region other are partially different The normalization grid values in marine site, are converted into " ocean essential periodic component field intensity " multiplied by with correction coefficient.So pass through space school " the ocean essential periodic component field intensity " of factual survey marine site maximum coverage range can just set up.

By the average of partially different marine site ocean essential periodic component grid valuesDivided by " the ocean in the range of the same space Key element periodic component field intensity " average(formula 6), just obtains periodic component correction coefficient average (λ)：

By " the ocean essential periodic component field intensity " grid value set (ψ (E in target marine site^BZ)) and averageIt is multiplied by school Positive coefficient average (λ) (formula 7, formula 8), just obtains the grid values of the target marine site periodic component by partially different marine site periodic component estimation Set (ψ (C)) and spatial domain average

ψ (C)=λ ψ (E^BZ) (7)

For there was only average value, without the ocean essential historical summary data for specifically investigating erect-position, as long as there is definite sea Area's field of investigation, under conditions of with " ocean essential field intensity ", it is also possible to obtain the target marine site ocean essential corrected by it Spatial domain average.

9. the models fitting of ocean essential periodic component time change

According to the rule that target marine site ocean essential periodic component is changed over time, the appropriate Model fitting cycle is selected Component spatial domain averageChange with time (formula 9), obtains model parameter, time integral is carried out to it, and seek its average (formula 10), obtains target marine site ocean essential periodic component time domain average

Wherein, (t, α, β, are m) periodic component spatial domain average-time change model to f, and t is the time, and α, β and m are models Parameter, (α, β m) are f (t, α, β, m) the definite integral function in the range of T time to F.

For the ocean essential periodic component on different time domain same time yardstick, its mechanical periodicity has common trait, Model f (t, α, β, parameter m) can be divided into variable element (α) and invariant parameter (β, m).α is by time domain periodic component size system About, may be varied from different time domain；The shape and relative position of β and m characterization model curves, it is basic in different time domain Do not change.The time domain relatively low for the investigation frequency, can only be fitted mould using adjacent Model in Time Domain parameter beta and the average of m Shape parameter α.

The history survey data in the partially different marine site for not exclusively being overlapped with target marine site, corrects its different time chi first The spatial domain average of periodic component on degreeThen summarizingOn the basis of-t Changing Patterns, appropriate model, fitting are selected Model parameter, is integrated and asks the average value in cycle to the time, obtains each periodic component time domain average in target marine site.

10. the calculating of ocean essential multifrequency period average

The minimum for calculating target marine site ocean essential change different time scales some cycles (T1, T2, T3, T4 ...) is public Multiple, obtains the minimum positive period (T) of ocean essential time domain change, then calculate minimum by the time domain mean value weighting of each periodic component Ocean essential multifrequency period average in positive period

Wherein,The respectively periodic component time domain average of cycle T 1, T2, T3, T4, The amount of cycles of i, j, m, l for cycle T 1, T2, T3, T4 in the range of minimum positive period (T).

The inventive method is based on oceanography Period Process and the regular understanding of ocean essential urban agglomeration, uses Wavelet analysis technology, decomposes the period of change of ocean essential multiple different time scales, extracts the cycle in multiple time scales Change component (periodic component)；There is the property of statistically uniformity with the cycle based on ocean essential with position phase spatial and temporal distributions situation Matter, space-time correction is carried out to partially different marine site ocean essential survey data, is used to calculate target marine site each periodic component of ocean essential Spatial domain average, using the time change of each periodic component spatial domain average of Model fitting, so by time integral calculate The time domain average of ocean essential periodic component, then by target marine site ocean in each periodic component time domain mean value computation minimum positive period The multifrequency period average of key element.

Beneficial effect：The present invention is recognized using the periodicity of oceanography process and the regularity of ocean essential urban agglomeration Know, with wavelet analysis and mathematics model-fitting technique, calculate ocean essential multifrequency period average.The inventive method considers ocean The periodicity of factor change, specify that the foundation of mean value computation selection time yardstick on ocean essential time dimension, so as to improve The accuracy and comparativity of ocean essential mean value computation.The inventive method can be advised according to ocean essential time-space resalved technique Rule, to being corrected analysis, the ocean essential of gained in the presence of a series of partially different historical summaries on investigation marine site, on control time Multifrequency period average can more be played and make full use of history survey data research sea close to the real conditions of target marine site ocean essential The effect of foreign factor change rule.Computational methods of the present invention and mathematical tool are all the technical field common method and work Tool, is easy to be promoted the use of in vast environmental monitoring and marine management department.

Brief description of the drawings

Fig. 1 is ocean essential multifrequency period mean value computation alignment technique route schematic diagram；

Fig. 2 is target marine site and partially different marine site relation schematic diagram --- by taking the Huanghai Sea as an example.

Specific embodiment

The present invention will be further described with reference to the accompanying drawings and detailed description.

A kind of method for calculating and correcting ocean essential average, multiple weeks are decomposed into by the time series variation of ocean essential The superposition of phase component, with reference to ocean essential periodic component with position phase spatial and temporal distributions situation statistically consistent property Matter, can be by partially different Sea Area Survey material computation correction target marine site ocean essential periodic component spatial domain average, using models fitting The change of periodic component spatial domain average on ocean essential different time scales, can correct each cycle point of target marine site ocean essential The time domain average of amount, and then calculate ocean essential multifrequency period average (such as Fig. 1).Target marine site ocean essential multifrequency period average The reliability of bearing calibration is calculated, can be checked by the reliability of its result of calculation, and the reliability of result of calculation can be with There is the survey data of complete time sequence by target marine site to verify.

The inventive method specifically includes following steps：

1st, target marine site ocean essential space lattice discretization：

Determine 120 ° -126 ° of Huanghai Sea E, N 33-35 ° scopes marine site is target marine site, with 1981-2010 day interval salt Used as research data, survey data has high-spatial and temporal resolution to degree time series data.Based on Kriging optimal interpolation methods, should Three-dimensional slotting spatial value program is worked out with griddata3 functions in Matlab software toolkits, by target marine site salinity website data By mesh scale (0.1 longitude long, 0.1 latitude wide, deep 15m) interpolation discretization, target under 1981-2008 year sequences is obtained Marine site salinity spatial discretization data, calculate daily target marine site salinity spatial domain average, obtain target marine site salinity spatial domain average Time series data.

2nd, ocean essential multifrequency period analysis：

With wavelet analysis technology, using target marine site salinity under Morlet wavelet decomposition 1981-2008 year sequences Gridded data data, obtains cycle moon of salinity altercation, annual period, climatic cycle (1mon, 12mon, 108mon) and each The amplitude (0.2psu, 1.15psu, 1.1psu) of periodic component.

Interpretation of result：By earth rotation, moon effect, act on to day quasi-periodic motion and controlled, oceanography process bag Include tidal motion, ocean circulation, atmospheric sedimentation, runoff input and biogeochemical process etc. and show the diurnal to a certain extent Phase, the cycle moon, annual period, the change of climatic cycle quasi-periodic, the periodicity of these oceanography key element spatial and temporal distributions restraining factors, So that there is ocean essential multifrequency period to change, analyzed by ocean essential multifrequency period, it can be found that ocean essential changes Multifrequency period and each periodic component amplitude.

3rd, partially different marine site ocean essential periodic component grid values and spatial domain correction for mean：

Determine 121.5 ° -126 ° of Huanghai Sea E, 34.7 ° of -39.2 ° of scope marine sites of N are partially different marine site A (such as Fig. 2), with 1981- Weekly interval salinity time series data in 1991 is used as partially different marine site A survey data.Around investigation marine site central gridding Survey data density highest, interpolation result is more reliable, thus target marine site and partially different marine site A central gridding respectively as Standard point of reference and partially different datum mark, by survey data (Grid data) normalization after discretization, obtain target marine site salinity Periodic component normalizes grid value setPartially different marine site salinity periodic component normalization grid value set

The salinity periodic component to the target marine site in the range of the A of coincidence marine site and partially different marine site A normalizes net respectively Lattice value is averagedBe divided by obtain normalization coefficient (γ), then will whole partially different marine site A salinity cycles point Amount normalization grid value setCorrection coefficient (γ) is multiplied by, the ocean for obtaining being put on the basis of standard point of reference will Plain periodic component normalizes grid set of correction valuesBy target marine site salinity periodic component normalization grid values collection ClosePartially different marine site salinity periodic component normalization grid set of correction valuesComposition " the salinity cycle point Amount field intensity " (ψ (E^BZ)).By partially different marine site A salinity periodic component grid values averageDivided by " the salt in the range of the same space Degree periodic component field intensity " averageObtain periodic component correction coefficient average (λ).Again by " the salinity cycle in target marine site Component field intensity " grid value set (ψ (E^BZ)) and averageCorrection coefficient average (λ) is multiplied by, is just obtained by partially different marine site A salt The grid value set (ψ (C)) and spatial domain average of the target marine site periodic component of degree periodic component estimation

Interpretation of result：With the cycle with position phase " salinity periodic component field intensity " be salinity altercation to orographic condition, tidal motion, Ocean circulation, atmospheric sedimentation, runoff input and the response of biogeochemical process quasi-periodic active force, embody salinity empty Between be distributed the stability and regularity of situation, also cause that periodic component correction coefficient has uniformity spatially, therefore just have May be by " the salinity periodic component field intensity " grid value set (ψ (E in target marine site^BZ)) and averageIt is multiplied by correction coefficient equal Value (λ), obtains the grid value set (ψ (C)) and sky of the target marine site periodic component by partially different marine site A salinity periodic component estimation Domain average

4th, the models fitting of ocean essential periodic component time change：

The cycle moon, annual period, climatic cycle periodic component spatial domain average according to target marine site salinityThe rule for changing over time, selection SIN function is fitted change of each periodic component spatial domain average with the time Change, obtain model parameter：

Mean time change function in periodic component spatial domain is integrated, and seeks its average, calculate target marine site salinity month Cycle, annual period, the periodic component time domain average of climatic cycle

Interpretation of result：(the cycle moon, annual period, climatic cycle) periodic component grid value set of target marine site salinity, body Existing is particular moment in lower cycle active force (tidal motion, ocean circulation, atmospheric sedimentation, runoff input and Biogeochemistry Process etc.) the non-equal attribute of spatial distribution situation and its space, and be based on result that partially different marine site A survey data correction gets with The result being directly calculated based on target Sea Area Survey data is compared, and error range shows the cycle between 0.5%~15% Component correction coefficient has preferable Space Consistency, and period effects power urban agglomeration is basicly stable.Therefore, using spatial domain Average can preferably characterize salinity periodic component mean state spatially.

Salinity periodic component spatial domain average changes with time, then what is embodied is the time domain of period effects power in certain space Distribution situation and its non-equal attribute, but salinity periodic component is presented the good periodicity such as moon cycle, annual period, climatic cycle Change, shows that the time domain distribution situation of period effects power is basicly stable.Therefore, salinity week can preferably be characterized using time domain average Mean state in phase component time domain.

5th, the calculating and checking of ocean essential multifrequency period average：

Based on partially different marine site A salinity 1981-1991 space lattice discretization data informations, mesh is obtained by free-air correction Mark the grid set of correction values (ψ (C)) and spatial domain average of marine site each periodic component of 1981-1991 salinityUsing sine The time series and target marine site salinity of the periodic component spatial domain average 1981-1991 that the correction of Function Fitting target marine site is obtained The time series of periodic component spatial domain average 1992-2008, obtain a moon cycle, annual period, climatic cycle (1 month, 12 months, 108 months) periodic component time domain average, calculate minimum positive period for 108 months, and then calculate multifrequency period average and be 32.6±0.2psu。

The Salinity Data data of direct application target marine site 1981-2008 high-spatial and temporal resolutions, with wavelet analysis and Model-fitting technique obtains the cycle moon, annual period, the cycle of climatic cycle (1 month, 12 months, 108 months) of salinity altercation Component time domain average, calculates multifrequency period averageIt is 32.4 ± 0.2psu, acquired results are provided with based on partially different marine site salinity The error of material correction result is less than 1%.

Interpretation of result：To correct the multifrequency period average of calculating and be based on based on partially different marine site A salinity time series data The multifrequency period average that target marine site salinity time series data is directly calculated is contrasted, and error is less than 1%, and the two does not show Sex differernce is write, so as to the reliability for demonstrating result of calculation is come.

Using the method for the present invention, based on the multifrequency period average that partially different marine site A salinity time series data correction is calculated It is basically identical with the multifrequency period average result directly calculated based on target marine site salinity time series data.This result enters one Step demonstrate by the correction of different marine site ocean essential time series data partially calculate target marine site multifrequency period average have it is operable Property, result of calculation have reliability higher, while also show the inventive method in ocean essential mean value computation have compared with Application value high.

It should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention, Some improvements and modifications can also be made, these improvements and modifications also should be regarded as protection scope of the present invention.In the present embodiment not Clear and definite each part can use prior art to be realized.

Claims (1)

1. it is a kind of calculate and correction ocean essential average method, it is characterised in that：The method is comprised the following steps：

1. target marine site ocean essential space lattice discretization

For the target marine site that ocean essential survey data has high-spatial and temporal resolution, based on Kriging optimal interpolation methods, application Matlab software programming three-dimensional interpolation programs, target marine site ocean essential erect-position survey data space lattice discretization is obtained The target marine site ocean essential gridded data data of certain hour sequence, target marine site ocean will under calculating each control time Plain spatial averaging, obtains the time series data of target marine site ocean essential spatial domain average；

2. ocean essential multifrequency period analysis

With wavelet analysis technology analysis target marine site ocean essential gridded data time series data, decompose ocean essential and become Periodic component P1, P2, P3, P4 ... on the different time scales of change, obtain period of change T1, the T2 of multiple time scales, T3, T4 ..., such as diurnal periodicity, the cycle moon, annual period, climatic cycle；

3. the correction of partially different marine site ocean essential periodic component grid values and spatial domain average

For the partially different marine site for only being partially overlapped with target marine site, its survey data and target marine site exist space-time it is partially different, it is necessary to Correction.Ocean essential standard point of reference is selected in target marine site first, by target marine site ocean essential periodic component gridding number Value is shown in formula 1 all divided by the numerical value of standard point of reference, obtains the normalization grid values of ocean essential periodic component：

$\mathrm{\ψ}\left(E_M^{BZ}\right)=\frac{\mathrm{\ψ}\left(C_i^M\right)}{C_Z^M}---\left(1\right)$

Wherein,It is target marine site ocean essential periodic component normalization grid value set,It is target marine site sea Foreign key element periodic component grid value set,It is ocean essential periodic component standard basis point value；

Equally, by identical control time, the partially different marine site ocean essential periodic component gridding of same period, its space model is selected Central gridding in enclosing, by the data after gridding all divided by partially different benchmark point value, is shown in formula 2 as partially different datum mark, obtains partially The normalization grid values of different marine site ocean essential periodic component：

$\mathrm{\ψ}\left(E_P^{PZ}\right)=\frac{\mathrm{\ψ}\left(C_i^P\right)}{C_Z^P}---\left(2\right)$

Wherein,It is partially different marine site ocean essential periodic component normalization grid value set,It is partially different marine site sea Foreign key element periodic component grid value set,It is the partially different benchmark point value of ocean essential periodic component；

In the range of the overlapping region in target marine site and partially different marine site, by target marine site ocean essential periodic component normalization grid The average of valueThe average of grid values is normalized divided by partially different marine site ocean essential periodic componentSee formula 3, obtain normalizing Change correction coefficient γ：

$\mathrm{\γ}=\frac{{\stackrel{\‾}{E}}_{MC}^{BZ}}{{\stackrel{\‾}{E}}_{PC}^{PZ}}---\left(3\right)$

Then by the normalization grid value set of entirely partially different marine site ocean essential periodic componentCorrection coefficient γ is multiplied by, See formula 4, the partially different marine site ocean essential periodic component normalization grid value set put on the basis of partially different datum mark is translated into The ocean essential periodic component normalization grid set of correction values put on the basis of standard point of reference

$\mathrm{\ψ}\left(E_P^{BZ}\right)=\mathrm{\γ}\·\mathrm{\ψ}\left(E_P^{PZ}\right)---\left(4\right)$

By target marine site ocean essential periodic component normalization grid value setThe partially different marine site ocean essential cycle point Amount normalization grid set of correction valuesIt is referred to as " ocean essential periodic component field intensity " ψ (E^BZ)；If partially different marine site and Target marine site does not exist overlapping region, can first be corrected to the normalization in other the partially different marine sites for having overlapping region with target marine site Grid values, are converted into " ocean essential periodic component field intensity " multiplied by with correction coefficient；Can so be set up by free-air correction " the ocean essential periodic component field intensity " of factual survey marine site maximum coverage range；

By the average of partially different marine site ocean essential periodic component grid valuesDivided by " the ocean essential cycle in the range of the same space Component field intensity " averageSee formula 6, just obtain periodic component correction coefficient average λ：

$\mathrm{\λ}=\frac{{\stackrel{\‾}{C}}_P}{{\stackrel{\‾}{E}}_P^{BZ}}---\left(6\right)$

By " ocean essential periodic component field intensity " the grid value set ψ (E in target marine site^BZ) and averageIt is multiplied by correction coefficient equal Value λ, is shown in formula 7, formula 8, just obtains grid value set ψ (C) of the target marine site periodic component by partially different marine site periodic component estimation And spatial domain average

ψ (C)=λ ψ (E^BZ) (7)

${\stackrel{\‾}{C}}_K=\mathrm{\λ}\·{\stackrel{\‾}{E}}^{BZ}---\left(8\right)$

For there was only average value, without the ocean essential historical summary data of erect-position are specifically investigated, as long as there is definite sea area to adjust Scope is looked into, under conditions of with " ocean essential field intensity ", it is also possible to obtain the target marine site ocean essential spatial domain corrected by it Average；

4. the models fitting of ocean essential periodic component time change

According to the rule that target marine site ocean essential periodic component is changed over time, appropriate Model fitting periodic component is selected Spatial domain averageChange with time, see formula 9, obtain model parameter, time integral is carried out to it, and seek its average, see formula 10, obtain target marine site ocean essential periodic component time domain average

${\stackrel{\‾}{C}}_K=f(t,\mathrm{\α},\mathrm{\β},m)---\left(9\right)$

${\stackrel{\‾}{C}}_S=\frac{1}{T}\underset{T}{\∫}f(t,\mathrm{\α},\mathrm{\β},m)\·dt=\frac{1}{T}F(\mathrm{\α},\mathrm{\β},m)---\left(10\right)$

Wherein, f (t, α, β, m) be periodic component spatial domain average-time change model, t is the time, and α, β and m are model parameters, (α, β m) are f (t, α, β, m) the definite integral function in the range of T time to F；

For the ocean essential periodic component on different time domain same time yardstick, its mechanical periodicity has common trait, model f (t, α, β, parameter m) are divided into variable element α and invariant parameter β, m；α is restricted by time domain periodic component size, when different Domain may be varied from；The shape and relative position of β and m characterization model curves, do not change substantially in different time domain； The time domain relatively low for the investigation frequency, using adjacent Model in Time Domain parameter beta and the average of m, only fitted model parameters α；

The history survey data in the partially different marine site for not exclusively being overlapped with target marine site, corrects on its different time scales first The spatial domain average of periodic componentThen summarizingOn the basis of-t Changing Patterns, appropriate model, model of fit are selected Parameter, is integrated and asks the average value in cycle to the time, obtains each periodic component time domain average in target marine site；

5. the calculating of ocean essential multifrequency period average

The calculating target marine site ocean essential change some cycle Ts 1, T2 of different time scales, the least common multiple of T3, T4 ..., The minimum positive period T of ocean essential time domain change is obtained, then minimum positive period is calculated by the time domain mean value weighting of each periodic component Interior ocean essential multifrequency period average

${\stackrel{\‾}{C}}_T=\frac{i\·{\stackrel{\‾}{C}}_{S1}+j\·{\stackrel{\‾}{C}}_{S2}+k\·{\stackrel{\‾}{C}}_{S3}+l\·{\stackrel{\‾}{C}}_{S4}+...}{i+j+k+l+...}---\left(11\right)$

Wherein,The respectively periodic component time domain average of cycle T 1, T2, T3, T4, i, j, m, Amount of cycles of the l for cycle T 1, T2, T3, T4 in minimum positive period T range.