CN105718666B - A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea - Google Patents
A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea Download PDFInfo
- Publication number
- CN105718666B CN105718666B CN201610041809.8A CN201610041809A CN105718666B CN 105718666 B CN105718666 B CN 105718666B CN 201610041809 A CN201610041809 A CN 201610041809A CN 105718666 B CN105718666 B CN 105718666B
- Authority
- CN
- China
- Prior art keywords
- wave
- radar
- sea
- parameter
- stream
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention proposes a kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea, comprising steps of needing sea condition to be simulated and used radar sensor type setting Ocean surface simulation parameter and radar parameter to establish wave-stream couple electromagnetic scattering model according to user；Obtain wave-stream coupling sea evolutionary process and characterization transient state sea；Further using Huygen's principle, gram western Hough is approximate and dirac Green's function, according to wave-stream couple electromagnetic scattering model inverting Radar backscattering coefficients γ and derive after normalization to radar cross section σ0.The present invention can obtain wave-ocean current coupling sea and electromagnetic wave phase interaction mechanism, promote the research of wave-stream couple electromagnetic scattering signatures, can be used as the technological means of ocean current radar inverting research and development.
Description
Technical field
The invention belongs to marine technology fields, and in particular to a kind of normalized radar backscatter cross section number on wave-stream coupling sea
It is worth analogy method.
Background technique
Ocean current has very important influence to oceanic environment, it is the conveying tie of energy and nutritive salt, while can also
With by ocean rubbish, various polluted-waters are distributed to each corner of World Oceans, to influence global water body, Air-sea fluxes.
Therefore to ocean surface ocean current carry out it is round-the-clock, on a large scale monitor it is most important, although people can use acoustic Doppler sea
Section plotter (ADCP) accurately measurement sea and underwater ocean current profile is flowed, but similar direct measuring method is quasi- real far from meeting
When, Large Area Synchronous measurement demand, satellite remote sensing means then become the only resource for meeting the demand.Regrettably, at present
Until, although there is fragmentary limitation algorithm to inverting sea flow velocity, such as have tried to utilize synthetic aperture radar (SAR)
Doppler frequency center deviates and carries out ocean current inverting along mark interference (ATI) technology, but without the biography being directly monitored to ocean current
Sensor is broadcast TV programs by satellite, and main cause is that wave-ocean current coupling sea is not understood completely with electromagnetic wave interaction mechanism, to lead
Cause the shortcoming of ocean current remote sensing inversion method.
Summary of the invention
Goal of the invention: in order to solve the above technical problems, the present invention proposes a kind of radar raster-displaying on wave-stream coupling sea
Section method for numerical simulation.
A kind of technical solution: technical solution proposed by the present invention are as follows: normalized radar backscatter cross section number on wave-stream coupling sea
It is worth analogy method, comprising the following steps:
Step 1: Ocean surface simulation parameter being set and radar parameter, Ocean surface simulation parameter include: wave significant wave height hs, wave
Dominant wavelength λ0, amplitude and dimensions in frequency factor b, Fractal Sea roughness value s and the main wave wave number K of wave0；Radar parameter includes:
Radar illumination range scale L, radar incidence angle θ, radar wavelength λ, radar incidence wave wave number k；
Wherein, b > 1, K0=2 π/λ0, k=2 π/λ, 1 < s < 2；
Step 2: extra large apparent height being calculated according to Ocean surface simulation parameter and radar parameter, if extra large apparent height are as follows:
Wherein, N is wave sinusoidal component number, and value meets bN-bN-1≥10；σ is the standard of sea apparent height ξ (x, t)
Variance, σ and significant wave height hsRelationship are as follows: hs=4 σ；C is normaliztion constant,
ωnFor the harmonic wave angular frequency of n-th wave component,bnFor the amplitude and frequency of n-th wave
Scale factor, g are acceleration of gravity；It is the wave phase of n-th wave；C0For ocean current flow velocity；
Step 3: using Huygen's principle, gram western Hough approximation and dirac Green's function, according to extra large apparent height ξ
The expression formula of (x, t) inverting Radar backscattering coefficients γ, γ are as follows:
Wherein, JpnFor Bessel function, p=(p0, p1... pn…pN-1)T, T is transposition oeprator；For time phase
Bit function:
Wherein
Step 4: being derived according to Radar backscattering coefficients γ after normalizing to radar cross section σ0, σ0Expression formula
Are as follows:
σ0=20log10γ (3)
Step 5: radar parameter used by needing sea condition to be simulated and analog sea state according to user determines
The input parameter value of step 1- step 3 is calculated the corresponding time from formula (1) and formula (2) and is irradiated after the normalization in space to thunder
Up to scattering section expression formula.
The utility model has the advantages that the present invention needs sea condition to be simulated and used radar sensor type to be arranged according to user
Ocean surface simulation parameter and radar parameter establish wave-stream couple electromagnetic scattering model, obtain wave-stream coupling sea evolutionary process and table
Levy transient state sea；Further using Huygen's principle, gram western Hough approximation and dirac Green's function, coupled according to wave-stream
It EM scatter model inverting Radar backscattering coefficients γ and derives after normalization to radar cross section σ0.Energy of the present invention
Wave-ocean current coupling sea and electromagnetic wave phase interaction mechanism is enough obtained, wave-stream couple electromagnetic scattering signatures is promoted to grind
Study carefully, can be used as the technological means of ocean current radar inverting research and development.
Detailed description of the invention
When Fig. 1 is different action of ocean current pre- same sea, radar scattering interface is normalized with radar incidence angle variation relation
Figure.
Specific embodiment
The present invention will be further explained with reference to the accompanying drawing.
The technical problems to be solved by the invention are as follows: obtain wave-ocean current coupling sea and electromagnetic wave phase interaction pass
System, to solve this technical problem, the present invention propose following technical scheme:
The following steps are included:
Step 1: Ocean surface simulation parameter being set and radar parameter, Ocean surface simulation parameter include: wave significant wave height hs, wave
Dominant wavelength λ0, amplitude and dimensions in frequency factor b, Fractal Sea roughness value s and the main wave wave number K of wave0=2 π/λ0；Radar ginseng
Number includes: radar illumination range scale L, radar incidence angle θ, radar wavelength λ, radar incidence wave wave number k=2 π/λ；
Step 2: extra large apparent height being calculated according to Ocean surface simulation parameter and radar parameter, if extra large apparent height is ξ (x, t)；
Step 3: using Huygen's principle, gram western Hough approximation and dirac Green's function, according to extra large apparent height ξ
The expression formula of (x, t) formula inverting Radar backscattering coefficients γ, γ are as follows:
In formula (1), b > 1, b are rational, then sea apparent height ξ (x, t) is periodic function, if b is irrational number,
Then sea apparent height ξ (x, t) is acyclic function, has randomness；1 < s < 2, to metric space sea figure
Degree of roughness；N is wave sinusoidal component number, and value meets bN-bN-1≥10；σ is the standard side of sea apparent height ξ (x, t)
Difference, σ and significant wave height hsRelationship are as follows: hs=4 σ；C is normaliztion constant,
ωnIt is the harmonic wave angular frequency of n-th wave component,G is acceleration of gravity；It is n-th sea
The wave phase of wave；C0It with propagation of ocean waves direction is then negative when being positive, is opposite for ocean current flow velocity；JpnFor Bessel function,
P=(p0, p1... pn…pN-1)T, T is transposition oeprator；
Wherein
Step 4: being derived according to Radar backscattering coefficients γ after normalizing to radar cross section σ0, σ0Expression formula
Are as follows:
σ0=20log10γ (3)
Step 5: radar parameter used by needing sea condition to be simulated and analog sea state according to user determines
The input parameter value of step 1- step 3 is calculated the corresponding time from formula (2) and formula (3) and is irradiated after the normalization in space to thunder
Up to scattering section expression formula.
Fig. 1 is different action of ocean current when same sea, the mould that normalization radar cross section changes with radar incidence angle
Quasi- result.Can be clearly seen that from figure normalization radar cross section and radar incidence angle relationship coincidence formula (2) and
Formula (3).
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (1)
1. a kind of wave-stream coupling sea normalized radar backscatter cross section method for numerical simulation, it is characterised in that the following steps are included:
Step 1: Ocean surface simulation parameter being set and radar parameter, Ocean surface simulation parameter include: wave significant wave height hs, the main wave of wave
Long λ0, amplitude and dimensions in frequency factor b, Fractal Sea roughness value s and the main wave wave number K of wave0；Radar parameter includes: radar
Irradiate range scale L, radar incidence angle θ, radar wavelength λ, radar incidence wave wave number k；
Wherein, b > 1, K0=2 π/λ0, k=2 π/λ, 1 < s < 2；
Step 2: extra large apparent height is calculated according to Ocean surface simulation parameter and radar parameter, if extra large apparent height is ξ (x, t), ξ (x,
T) expression formula are as follows:
Wherein, N is wave sinusoidal component number, and value meets bN-bN-1≥10；σ is the standard variance of sea apparent height ξ (x, t),
σ and significant wave height hsRelationship are as follows: hs=4 σ；C is normaliztion constant,
ωnFor the harmonic wave angular frequency of n-th wave component,bnFor n-th wave amplitude and dimensions in frequency because
Son, g are acceleration of gravity；It is the wave phase of n-th wave；C0For ocean current flow velocity；
Step 3: using Huygen's principle, gram western Hough approximation and dirac Green's function, according to extra large apparent height ξ (x, t)
The expression formula of inverting Radar backscattering coefficients γ, γ are as follows:
Wherein,For Bessel function, p=(p0, p1... pn…pN-1)T, T is transposition oeprator；For time phase
Function:
Wherein
Step 4: being derived according to Radar backscattering coefficients γ after normalizing to radar cross section σ0, σ0Expression formula are as follows:
σ0=20log10γ (3)
Step 5: radar parameter used by needing sea condition to be simulated and analog sea state according to user determines step
The input parameter value of 1- step 3 dissipates after calculating the normalization in corresponding time and irradiation space from formula (2) and formula (3) to radar
Penetrate section expression.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610041809.8A CN105718666B (en) | 2016-01-21 | 2016-01-21 | A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610041809.8A CN105718666B (en) | 2016-01-21 | 2016-01-21 | A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105718666A CN105718666A (en) | 2016-06-29 |
| CN105718666B true CN105718666B (en) | 2019-02-12 |
Family
ID=56154911
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610041809.8A Active CN105718666B (en) | 2016-01-21 | 2016-01-21 | A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105718666B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107144826B (en) * | 2017-04-28 | 2019-07-23 | 华中科技大学 | A method of establishing improved spectrometer modulated signal model |
| CN112240907B (en) * | 2019-07-18 | 2021-09-10 | 中国科学院声学研究所 | Acoustic propagation calculation method considering influence of bubble mixing layer under shallow-sea fluctuating sea surface |
| CN112213708B (en) * | 2020-10-14 | 2022-04-22 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Shore-based sea clutter observation method |
| CN113126053B (en) * | 2021-04-08 | 2024-02-13 | 中国人民解放军92942部队 | Method and device for evaluating radar cross section of ship |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101493522A (en) * | 2009-02-27 | 2009-07-29 | 武汉理工大学 | Anomalous wave detecting method and device on sea surface |
| CN101697009A (en) * | 2009-10-27 | 2010-04-21 | 武汉理工大学 | Sea wave surface reduction method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060149513A1 (en) * | 2005-01-06 | 2006-07-06 | Sangpil Yoon | Method and system for simulating a surface acoustic wave on a modeled structure |
| US9594857B2 (en) * | 2013-07-15 | 2017-03-14 | California Institute Of Technology | Methods for designing quadruple-ridged flared horn antennas |
-
2016
- 2016-01-21 CN CN201610041809.8A patent/CN105718666B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101493522A (en) * | 2009-02-27 | 2009-07-29 | 武汉理工大学 | Anomalous wave detecting method and device on sea surface |
| CN101697009A (en) * | 2009-10-27 | 2010-04-21 | 武汉理工大学 | Sea wave surface reduction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105718666A (en) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105718666B (en) | A kind of normalized radar backscatter cross section method for numerical simulation on wave-stream coupling sea | |
| Fer | Near-inertial mixing in the central Arctic Ocean | |
| CN105181145B (en) | A kind of underwater polarization field analogy method of rough water transmitted light | |
| Wang et al. | Modeling of sediment and heavy metal transport in Taihu Lake, China | |
| Cosoli et al. | Surface circulation in the Gulf of Trieste (northern Adriatic Sea) from radar, model, and ADCP comparisons | |
| CN106768179B (en) | The measurement method of tidal level based on the station continuous operation GNSS signal-to-noise ratio data | |
| Wang et al. | Modeling tides in Monterey Bay, California | |
| Noble et al. | On the longshelf structure and dynamics of subtidal currents on the eastern United States continental shelf | |
| CN109696234A (en) | The determination method and system of horizontal distance between a kind of launch point and receiving point | |
| CN104091048A (en) | Underwater target depth measuring method adopting passive vertical double-vector hydrophones | |
| Poulain | Tidal currents in the Adriatic as measured by surface drifters | |
| CN104062645A (en) | Method for measuring phase difference between parametric array difference-frequency wave and same-frequency small-amplitude wave signal | |
| Chen et al. | Estimate of tidal constituents in nearshore waters using X-band marine radar image sequences | |
| Barclay et al. | On the spatial properties of ambient noise in the Tonga Trench, including effects of bathymetric shadowing | |
| Soloviev et al. | A near-surface microstructure sensor system used during TOGA COARE. Part II: Turbulence measurements | |
| Hsu et al. | Numerical study of monsoon effect on Green Island wake | |
| Yu et al. | Measurements of midfrequency acoustic backscattering from a sandy bottom in the South Yellow Sea of China | |
| Brown et al. | Intercomparison of the Charnock and COARE bulk wind stress formulations for coastal ocean modelling | |
| Lara et al. | Design and application of a passive acoustic monitoring system in the Spanish implementation of the marine strategy framework directive | |
| Wijesekera et al. | Surface wave effects on high-frequency currents over a shelf edge bank | |
| Mobarek | Directional spectra of laboratory wind waves | |
| Hisaki | Short‐wave directional properties in the vicinity of atmospheric and oceanic fronts | |
| RU2625094C1 (en) | Method for determining directions and distance to signal source | |
| Es-sadaoui et al. | Design and performance measurement of a digital embedded system implementing a I/Q demodulator for AoA range-based localization in underwater environment | |
| Deal et al. | Modeling acoustic vector fields for inverse problems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: 211800 No. 219 Ningliu Road, Jiangbei New District, Nanjing City, Jiangsu Province Patentee after: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY Address before: The Olympic Avenue in Jianye District of Nanjing city of Jiangsu Province, No. 69 210019 Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY |
|
| CP02 | Change in the address of a patent holder |