CN106597546B - A method of detection river bed surface layer silt physical property - Google Patents
A method of detection river bed surface layer silt physical property Download PDFInfo
- Publication number
- CN106597546B CN106597546B CN201611121801.9A CN201611121801A CN106597546B CN 106597546 B CN106597546 B CN 106597546B CN 201611121801 A CN201611121801 A CN 201611121801A CN 106597546 B CN106597546 B CN 106597546B
- Authority
- CN
- China
- Prior art keywords
- silt
- signal
- frequency
- reflection
- formula
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000002344 surface layer Substances 0.000 title claims abstract description 23
- 230000000704 physical effect Effects 0.000 title claims abstract description 11
- 238000001514 detection method Methods 0.000 title abstract description 11
- 230000035699 permeability Effects 0.000 claims abstract description 18
- 239000013049 sediment Substances 0.000 claims abstract description 11
- 239000002245 particle Substances 0.000 claims description 17
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000001052 transient effect Effects 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000001228 spectrum Methods 0.000 claims description 9
- 239000011148 porous material Substances 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 239000004576 sand Substances 0.000 claims description 5
- 230000006854 communication Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims description 2
- 239000010410 layer Substances 0.000 claims 1
- 239000000463 material Substances 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6246—Permeability
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to acoustics and geophysics field, and in particular to a method of detection river bed surface layer silt physical property.This method analyzes shallow seismic profile data, calculates attenuation coefficient α, search porosity according to porosity-reflection coefficient curve, calculate permeability κ, calculate silt velocity of wave c by extracting shallow seismic profile sonic data；Permeability-attenuation coefficient curve is calculated, according to step attenuation coefficient α, inquires permeability；It is calculated by successive ignition, as porosity βnWith previous iteration value βn‑1When difference is smallerStop iteration；Silt density p is calculated, is calculated reflection R (f)；Export permeability κ, silt velocity of wave c, reflection R (f), attenuation coefficient α, silt density p.This method can obtain a variety of silt physical parameters, abundant shoal materials information, underwater sediment physical parameter can be obtained merely with shallow bottom section plotter, it saves human and material resources, improve working efficiency, this method can obtain more accurate silt parameter information, provide data supporting for researchs such as river dynamics.
Description
Technical field
The invention belongs to acoustics and geophysics field, and in particular to a method of detection river bed surface layer silt physical property.
Background technique
Currently, the acquisition of river bed siltation surface layer silt physical parameter (density, permeability, porosity etc.) is mainly taken by machinery
Sample obtains, and the perturbation of mechanical sampling causes sample tests to generate deviation.In addition, mechanical sampling spends a large amount of manpower objects
Power, efficiency are lower.Acoustic method is also commonly used to detection river sediment characteristic, presently mainly judges that silt becomes silted up according to sonar image
Product structure, still lacks effective ways to the extraction of silt physical parameter.
The present invention is a kind of new method using sub-bottom profiler detection river bed surface layer silt parameter.Based on pore media
Acoustic theory receives signal by analysis sub-bottom profiler, and acoustic reflection coefficient and the decaying for extracting sub-bottom profiler are special
Sign obtains the parameters such as river sediment density, velocity of wave, porosity and permeability.The present invention does not need mechanical sampling test result work
For foundation, receiving signal merely with sub-bottom profiler can be obtained related silt parameter.Compared with mechanical sampling, have quick
Efficient feature.
Summary of the invention
In conclusion for overcome the deficiencies in the prior art, the present invention provides a kind of detection river bed surface layer silt physical property
Method, related formula and parameter is as follows:
K is wave number,For angular frequency, f is frequency, and c is pore media velocity of wave, ρfAnd ρsRespectively pore-fluid density and
Sand grain density, and be known constant.ρ is silt volume density.η indicates pore water viscosity, and κ is permeability, KγIt is silt
The bulk modulus of grain, μ is the modulus of shearing of silt frame, KbAnd KfRespectively indicate the volume of silt skeleton modulus and pore-fluid
Modulus, F are the deviation that pore-fluid is generated with frequency increase.The velocity of wave c of waterω=1500kg/m3, Kr=3.2 × 1010Pa, Kf
=2.395 × 109Pa, η=0.001Kg/m*s.
ρ=β ρf+(1-β)ρs (2)
τx=K0τz (13)
τy=K0τz (14)
τz=(1- β) (ρs-ρf)gz (15)
The value range of δ is 0.1~0.3, K0=0.5, g are acceleration of gravity, and z is water-sediment interface depth below
Degree, δb=0.1, δs=0.1.
F is determined by (16)-(19) formula:
Attenuation coefficient α is determined by (20)-(24) formula:
A=C2-HM (20)
Attenuation coefficient α=Im { k1Im expression imaginary part (24)
Reflection R (β, κ, f) is determined by (25) formula:
Wherein,B1=ρf
ω2, B2=Hk1 2-Ck1 2G1, B3=Hk2 2-Ck2 2G2C2=Mk1 2G1-Ck1 2, C3=Mk2 2G2-Ck2 2；
2.1 from shallow seismic profile data extract reflection R (f) calculation method it is as follows: transmitting and receive signal, from
Sub-bottom profiler, which receives, extracts water-surface layer sediment interface reflection coefficient in signal.Emit signal S (f)=E (f) Hs
(f), received spectrum can indicate are as follows:
r1For energy converter to water-bed distance, E (f) is the frequency spectrum of source signal, HsIt (f) is the transfer function of Sonar system, E*
(f) complex conjugate function for being E (f), R (f) are reflection coefficient.Reflection coefficient is calculated using the received secondary wave signal of sonar, two
Subwave may be expressed as: by compressed spectrum
r2For the sum of the distance and the depth of water that energy converter distance is water-bed.Water-bed reflection signal is reflected by air-water interface,
Then it is received by underwater reflection by energy converter, reception signal is secondary wave.Reflection coefficient is -1.Formula (26) is obtained divided by formula (27)
To reflection coefficient:
2.2 relaxation time trWith the method for building up of reflection signal frequency shift Δ f relation curve
The calculation method of signal transient frequency:
S (t) is signal,For Hilbert transformation, E (t) is function envelope, and φ (t) is phase, φ (t)=- jloge
[z (t)/E (t)], the instantaneous frequency of s (t) are
The signal of sub-bottom profiler transmitting is linear FM signal S (f), is calculated by the calculation method of signal transient frequency
The instantaneous frequency f of S (f)1。
The transfer function of silt:
Wherein f is frequency, and x is propagation distance of the sound wave in silt,V is propagation of the sound wave in silt
Speed.Only consider that frequency shift effect of the sound wave in communication process, reflection signal are expressed as B'(f)=S (f) Ha(f).X=1m is taken,
By the calculation method of signal transient frequency, B'(f is calculated in numerical value) instantaneous frequency f2, Δ f=f1-f2.Obtain trWith frequency displacement
The relation curve of Δ f；
2.3 methods for calculating attenuation coefficient α:
According to 2.2, signal transient frequency f is reflected from river water-surface layer silt is extracted in A'(f)1With surface layer silt-bottom mud
The instantaneous frequency f of husky interface reflection signal2, Δ f=f1-f2, by trRelaxation time t is obtained with the relation curve of frequency displacement Δ fr,
By formula (30), attenuation coefficient of the sound wave in silt:
α=k'f2 (31)
2.4 methods that porosity is searched according to porosity-reflection coefficient curve:
Porosity β and median particle diameter φ (unit phi, φ=- log2dmm, dmmFor the partial size of sand grain) experience it is public
Formula
β=0.208+0.0943 φ -0.00334 φ2 (32)
The empirical equation of permeability κ and porosity β, i.e. Kozeny-Carman formula
D is silt median particle diameter, and K is empirical estimation, when particle is circular tube shaped, K=2；When particle is spherical, K
=5,
By formula (25) and formula (33), the relation curve template of median particle diameter d, porosity β, frequency f and reflection coefficient is obtained.
A method of detection river bed surface layer silt physical property, its step are as follows:
(1) shallow seismic profile sonic data is extracted, calculates acoustic reflection coefficient R (f) by 2.1；
(2) by 2.2, shallow seismic profile data are analyzed, frequency displacement Δ f is calculated, finds relaxation time tr；
(3) FM signal centre frequency is taken by 2.3, f, calculates attenuation coefficient α；
(4) by 2.4, porosity is searched according to porosity-reflection coefficient curve, convolution (32), (33) calculate permeability
κ, κ is only used as initial value here；
(5) silt velocity of wave c is calculated by formula (1) and formula (33)；
(6) permeability-attenuation coefficient curve is calculated, according to the attenuation coefficient α of step (2), inquires permeability；
(7) step (4)-(7) are repeated, are calculated by successive ignition, as porosity βnWith previous iteration value βn-1Difference is smaller
When (| βn-1-βn|/βn< 3%), stop iteration；
(8) silt density p is calculated by formula (2), calculates reflection R (f) by formula (25)；
(9) permeability κ, silt velocity of wave c, reflection R (f), attenuation coefficient α, silt density p are exported.
The utility model has the advantages that
(1) this method can obtain a variety of silt physical parameters, enrich shoal materials information.
(2) this method can obtain underwater sediment physical parameter merely with shallow bottom section plotter, save human and material resources, and improve
Working efficiency.
(3) this method can obtain more accurate silt parameter information, provide data supporting for researchs such as river dynamics.
Detailed description of the invention
Fig. 1 is sampling process schematic diagram of the present invention；
Fig. 2 is calculation flow chart of the present invention；
Fig. 3 is relation curve (f=3700Hz) figure of reflection R (f) and porosity β；
Fig. 4 is unit apart from frequency displacement Δ f and relaxation time trRelation curve；
Fig. 5 is that sub-bottom profiler receives 47 μ s of signal sampling interval.
Specific embodiment
Below with reference to embodiment, technical solution of the present invention is further described in detail.
Embodiment 1
A method of detection river bed surface layer silt physical property, by taking Xiaolangdi reservoir area as an example:
36 section of Xiaolangdi reservoir area, 1 measuring point, sub-bottom profiler transmitted signal bandwidth 500-7000Hz, time width
20ms.The frequency displacement Δ f at the upper and lower interface of surface layer silt is 88Hz, frequency displacement the Δ f=18Hz, attenuation coefficient α of unit distance=
0.04, reflection R=0.17.By calculating, permeability κ=4.46 × 10-13, porosity β=67.66%, silt velocity of wave c
=1370m/s, silt density p=1549.4kg/m3。
Embodiment 2
A method of detection river bed surface layer silt physical property, by taking 4 section of Reservoir Area of Sanmenxia as an example:
4 section of Reservoir Area of Sanmenxia, 1 measuring point, sub-bottom profiler transmitted signal bandwidth 500-7000Hz, time width 20ms.It is single
The frequency displacement Δ f of position distance is 71Hz, attenuation coefficient α=0.094, reflection R=0.16.By calculating, permeability κ=3.89
×10-11, porosity β=0., silt velocity of wave c=1503m/s, silt density p=1412kg/m3。
Embodiment 3
A method of detection river bed surface layer silt physical property, by taking 8 section of Reservoir Area of Sanmenxia as an example:
8 section of Reservoir Area of Sanmenxia, 1 measuring point, sub-bottom profiler transmitted signal bandwidth 500-7000Hz, time width 20ms.It is single
The frequency displacement Δ f=104Hz of position distance, attenuation coefficient α=0.136, reflection R=0.34.By calculating, permeability κ=
6.49×10-12, porosity β=0.46, silt velocity of wave c=1623m/s, silt density p=1921kg/m3。
Although above the present invention is described in detail with a general description of the specific embodiments,
On the basis of the present invention, it can be made some modifications or improvements, this will be apparent to those skilled in the art.Cause
This, these modifications or improvements, fall within the scope of the claimed invention without departing from theon the basis of the spirit of the present invention.
Claims (2)
1. a kind of method for detecting river bed surface layer silt physical property, it is characterised in that: related formula and parameter is as follows:
K is wave number,For angular frequency, f is frequency, and c is pore media velocity of wave, ρfAnd ρsRespectively pore-fluid density and silt
Grain density, and be known constant, ρ is silt volume density, and η indicates pore water viscosity, and κ is permeability, KrIt is sand grain
Bulk modulus, μ are the modulus of shearing of silt frame, KbAnd KfThe bulk modulus of silt skeleton modulus and pore-fluid is respectively indicated,
F is the deviation that pore-fluid is generated with frequency increase, the velocity of wave c of waterω=1500kg/m3, Kr=3.2 × 1010Pa, Kf=
2.395×109Pa, η=0.001Kg/m*s；
ρ=β ρf+(1-β)ρs (2)
τx=K0τz (13)
τy=K0τz (14)
τz=(1- β) (ρs-ρf)gz (15)
The value range of δ is 0.1~0.3, K0=0.5, g are acceleration of gravity, and z is water-sediment interface depth below, δb
=0.1, δs=0.1；
F(ε)It is determined by (16)-(19) formula:
Attenuation coefficient is determined by (20)-(24) formula:
A=C2-HM (20)
Attenuation coefficient=Im { k1Im expression imaginary part (24)
Reflection R (β, κ, f) is determined by (25) formula:
Wherein,B1=ρf
ω2, B2=Hk1 2-Ck1 2G1, B3=Hk2 2-Ck2 2G2, C2=Mk1 2G1-Ck1 2, C3=Mk2 2G2-Ck2 2；
2.1 from shallow seismic profile data extract reflection R (f) calculation method it is as follows: transmitting and receive signal, from shallowly
Layer section plotter, which receives, extracts river water-surface layer sediment interface reflection coefficient in signal, emit signal S (f)=E (f) Hs(f),
Received spectrum can indicate are as follows:
r1For energy converter to water-bed distance, E (f) is the frequency spectrum of source signal, HsIt (f) is the transfer function of Sonar system, E*(f)
For the complex conjugate function of E (f), R (f) is acoustic reflection coefficient；Sound wave reflection system is calculated using the received secondary wave signal of sonar
Number R (f), secondary wave may be expressed as: by compressed spectrum
r2For the sum of the distance and the depth of water that energy converter distance is water-bed, water-bed reflection signal is reflected by air-water interface, then
It is received by underwater reflection by energy converter, reception signal is secondary wave, and reflection coefficient is negative, obtains sound wave reflection by formula (28)
Coefficients R (f):
2.2 relaxation time trWith the method for building up of reflection signal frequency shift Δ f relation curve
The calculation method of signal transient frequency:
S (t) is signal,For Hilbert transformation, E (t) is function envelope, and φ (t) is phase, φ (t)=- jloge[z
(t)/E (t)], the instantaneous frequency of s (t) is
The signal of sub-bottom profiler transmitting is linear FM signal S (f), calculates S (f) by the calculation method of signal transient frequency
Instantaneous frequency f1,
The transfer function of silt:
Wherein f is frequency, and x is propagation distance of the sound wave in silt,V is spread speed of the sound wave in silt,
Only consider that frequency shift effect of the sound wave in communication process, reflection signal are expressed as B'(f)=S (f) Ha(f), x=1m is taken, by believing
B'(f is calculated in the calculation method of number instantaneous frequency, numerical value) instantaneous frequency f2, Δ f=f1-f2；Obtain trWith frequency displacement Δ f
Relation curve；
2.3 methods for calculating attenuation coefficient α of the sound wave in silt:
According to 2.2, signal transient frequency f is reflected from river water-surface layer silt is extracted in A'(f)1It is handed over surface layer silt-bottom sediment
The instantaneous frequency f of interface reflection signal2, Δ f=f1-f2, by trRelaxation time t is obtained with the relation curve of frequency displacement Δ fr,
Attenuation coefficient α=k'f of the sound wave in silt2 (31)
2.4 methods that porosity is searched according to porosity-reflection coefficient curve:
The unit of median particle diameter φ is phi, φ=- log2dmm, dmmFor the partial size of sand grain, porosity β and median particle diameter φ
Empirical equation
β=0.208+0.0943 φ -0.00334 φ2 (32)
The empirical equation of permeability κ and porosity β, i.e. Kozeny-Carman formula
D is silt median particle diameter, and K is empirical estimation, when particle is circular tube shaped, K=2；When particle is spherical, K=5,
By formula (25) and formula (33), the relation curve template of median particle diameter d, porosity β, frequency f and reflection coefficient is obtained.
2. a kind of method for detecting river bed surface layer silt physical property, its step are as follows:
(1) shallow seismic profile sonic data is extracted, by the calculation method for extracting reflection R (f) from shallow seismic profile data
As follows: transmitting and reception signal receive from sub-bottom profiler and extract river water-surface layer sediment interface reflection system in signal
Number emits signal S (f)=E (f) Hs(f), received spectrum can indicate are as follows:
r1For energy converter to water-bed distance, E (f) is the frequency spectrum of source signal, HsIt (f) is the transfer function of Sonar system, E*(f)
For the complex conjugate function of E (f), R (f) is reflection coefficient；Reflection R (f) is calculated using the received secondary wave signal of sonar, two
Subwave may be expressed as: by compressed spectrum
r2For the sum of the distance and the depth of water that energy converter distance is water-bed, water-bed reflection signal is reflected by air-water interface, then
It is received by underwater reflection by energy converter, reception signal is secondary wave, and reflection coefficient is negative, obtains sound wave reflection by formula (28)
Coefficients R (f):
It calculates acoustic reflection coefficient R (f)；
(2) shallow seismic profile data are analyzed, frequency displacement Δ f is calculated, finds relaxation time tr
The calculation method of signal transient frequency:
S (t) is signal,For Hilbert transformation, E (t) is function envelope, and φ (t) is phase, φ (t)=- jloge[z(t)/
E (t)], the instantaneous frequency of s (t) is
The signal of sub-bottom profiler transmitting is linear FM signal S (f), calculates S (f) by the calculation method of signal transient frequency
Instantaneous frequency f1,
The transfer function of silt:
Wherein f is frequency, and x is propagation distance of the sound wave in silt,V is spread speed of the sound wave in silt,
Only consider that frequency shift effect of the sound wave in communication process, reflection signal are expressed as B'(f)=S (f) Ha(f), x=1m is taken, by believing
B'(f is calculated in the calculation method of number instantaneous frequency, numerical value) instantaneous frequency f2, Δ f=f1-f2；Obtain trWith frequency displacement Δ f
Relation curve, calculate frequency displacement Δ f, find relaxation time tr；
(3) f takes FM signal centre frequency, calculates attenuation coefficient α
Signal transient frequency f is reflected from river water-surface layer silt is extracted in A'(f)1With surface layer silt-bottom sediment interface reflection
The instantaneous frequency f of signal2, Δ f=f1-f2, by trRelaxation time t is obtained with the relation curve of frequency displacement Δ fr,
α=k'f2(31)；
(4) method that porosity is searched according to porosity-reflection coefficient curve:
The unit of median particle diameter φ is phi, φ=- log2dmm, dmmFor the partial size of sand grain, porosity β and median particle diameter φ
Empirical equation
β=0.208+0.0943 φ -0.00334 φ2 (32)
The empirical equation of permeability κ and porosity β, i.e. Kozeny-Carman formula
D is silt median particle diameter, and K is empirical estimation, when particle is circular tube shaped, K=2；When particle is spherical, K=5,
By formula (25) Obtain the relation curve template of median particle diameter d, porosity β, frequency f and reflection coefficient；
(5) by formula (1)Calculate pore media velocity of wave c；
(6) permeability-attenuation coefficient curve is calculated, according to the attenuation coefficient α of step (3), inquires permeability；
(7) step (4)-(6) are repeated, are calculated by successive ignition, as porosity βnWith previous iteration value βn-1When difference is smaller, |
βn-1-βn|βn< 3% stops iteration；
(8) by formula (2) ρ=β ρf+(1-β)ρsSilt volume density ρ is calculated, by formula (25)It calculates reflection R (f)；
(9) permeability κ, pore media velocity of wave c, attenuation coefficient α, silt volume density ρ are exported,
By formula (25)What is calculated is anti-
Penetrate coefficients R (f).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611121801.9A CN106597546B (en) | 2016-12-08 | 2016-12-08 | A method of detection river bed surface layer silt physical property |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611121801.9A CN106597546B (en) | 2016-12-08 | 2016-12-08 | A method of detection river bed surface layer silt physical property |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106597546A CN106597546A (en) | 2017-04-26 |
| CN106597546B true CN106597546B (en) | 2019-04-02 |
Family
ID=58597446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611121801.9A Expired - Fee Related CN106597546B (en) | 2016-12-08 | 2016-12-08 | A method of detection river bed surface layer silt physical property |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106597546B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109440721B (en) * | 2018-12-21 | 2019-08-06 | 中国水利水电科学研究院 | A kind of identification of river channel sedimentation object and preprocess method |
| CN109738344B (en) * | 2019-01-15 | 2021-06-15 | 厦门大学 | Real-time measurement method for suspended sediment concentration at sea entrance based on sound attenuation |
| CN113504159B (en) * | 2021-06-16 | 2023-01-03 | 中国农业大学 | Method and device for detecting and analyzing opaque particulate matter and electronic equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5815465A (en) * | 1997-04-11 | 1998-09-29 | Turgut; Altan | Method and apparatus of classifying marine sediment |
| CN1900743A (en) * | 2005-07-01 | 2007-01-24 | 中国科学院声学研究所 | Method and system for measuring sea bed deposite sediment property by shallow stratigraphic section instrument |
| CN101545974A (en) * | 2009-04-23 | 2009-09-30 | 杭州瑞声海洋仪器有限公司 | Stratigraphic profile acoustic detection method for sub-bottom profiler |
| US8009512B2 (en) * | 2006-07-17 | 2011-08-30 | Biosonics, Inc. | Networked sonar observation of selected seabed environments |
| CN103884625A (en) * | 2014-04-15 | 2014-06-25 | 黄河水利委员会黄河水利科学研究院 | Data-fusion-based reservoir bottom mud density detection method and system |
-
2016
- 2016-12-08 CN CN201611121801.9A patent/CN106597546B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5815465A (en) * | 1997-04-11 | 1998-09-29 | Turgut; Altan | Method and apparatus of classifying marine sediment |
| CN1900743A (en) * | 2005-07-01 | 2007-01-24 | 中国科学院声学研究所 | Method and system for measuring sea bed deposite sediment property by shallow stratigraphic section instrument |
| US8009512B2 (en) * | 2006-07-17 | 2011-08-30 | Biosonics, Inc. | Networked sonar observation of selected seabed environments |
| CN101545974A (en) * | 2009-04-23 | 2009-09-30 | 杭州瑞声海洋仪器有限公司 | Stratigraphic profile acoustic detection method for sub-bottom profiler |
| CN103884625A (en) * | 2014-04-15 | 2014-06-25 | 黄河水利委员会黄河水利科学研究院 | Data-fusion-based reservoir bottom mud density detection method and system |
Non-Patent Citations (6)
| Title |
|---|
| A Method for Estimating the Physical and Acoustic Properties of the Sea Bed Using Chirp Sonar Data;Steven G. Schock;《IEEE JOURNAL OF OCEANIC ENGINEERING》;20041031;第29卷(第4期);第1200-1217页 |
| An effective density fluid model for acoustic propagation in sediments derived from Biot theory;Kevin L. Williams;《The Journal of the Acoustical Society of America》;20011130;第110卷(第5期);第2276-2281页 |
| Estimating Geoacoustic Properties of Surficial Sediments in the North Mien-Hua Canyon Region With a Chirp Sonar Profiler;Linus Y. S. Chiu et al.;《IEEE JOURNAL OF OCEANIC ENGINEERING》;20150131;第40卷(第1期);第222-236页 |
| Sonar attenuation modeling for classification of marine sediments;Lester R. LeBlanc et al.;《The Journal of the Acoustical Society of America》;19920131;第91卷(第1期);第116-126页 |
| 小浪底库区深层淤积泥沙物理特性分析;郑军等;《人民黄河》;20141031;第36卷(第10期);第23-25页 |
| 水声勘探中的浅地层剖面技术应用;谢向文等;《港工技术》;20091031;第46卷(第5期);第53-55页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106597546A (en) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Correlations between the in situ acoustic properties and geotechnical parameters of sediments in the Yellow Sea, China | |
| CN102393196B (en) | Ocean internal wave acoustic detection method | |
| CN104749568B (en) | A kind of sorting technique of the Shallow Water Target depth based on hydrophone array | |
| Xu | Measuring currents in submarine canyons: technological and scientific progress in the past 30 years | |
| CN106597546B (en) | A method of detection river bed surface layer silt physical property | |
| Leifer et al. | The acoustic signature of marine seep bubbles | |
| CN104820218B (en) | Shallow sea bottom single-parameter inversion method based on frequency domain autocorrelation | |
| CN104007176A (en) | Full-wave field detection system and method of complex geotechnical engineering medium | |
| CN106646661A (en) | Comprehensive survey system for hydrogeology of mineral deposit | |
| Dolgikh et al. | A hydroacoustic system that radiates at frequencies of 19− 26 Hz | |
| Bjørnestad et al. | Lagrangian measurements of orbital velocities in the surf zone | |
| CN106959152A (en) | A kind of method of SVEL in measurement seawater | |
| CN2852153Y (en) | Device for measuring acoustical attenuation and frequency relation of marine bottom sediment | |
| CN103884625B (en) | Data-fusion-based reservoir bottom mud density detection method and system | |
| CN107389786B (en) | A kind of gas incursion check signal processing method based on time domain | |
| CN107942316A (en) | Concentrate suspension movement velocity method of estimation in a kind of water based on multibeam sonar beamformer output signal | |
| CN206057253U (en) | A kind of acoustic characteristic measurement apparatus | |
| CN110780340B (en) | Inversion method for transverse wave acoustic velocity or thickness of shallow sea elastic soft sedimentary layer | |
| CN112540403A (en) | Stress calculation and dynamic response stage division method based on simulated earthquake weak layer | |
| KR101349562B1 (en) | Estimation method of gas content in the bubbly sediment layer with the difference frequency wave | |
| Li et al. | Acoustic parameters inversion and sediment properties in the Yellow River reservoir | |
| CN204286983U (en) | A kind of Reservoir Sediment density detection system based on data fusion | |
| Rui et al. | Application of Sub-Bottom Profiler to Study Riverbed Structure and Sediment Density | |
| Soloway et al. | Modeling explosion generated Scholte waves in sandy sediments with power law dependent shear wave speed | |
| CN100456046C (en) | Acoustic method and system for measuring multi-metal nodule ore in sea bottom |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190402 Termination date: 20191208 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |