CN106546523A - Detection means and method based on the porosity of porous medium of stream solid Interface Wave - Google Patents
Detection means and method based on the porosity of porous medium of stream solid Interface Wave Download PDFInfo
- Publication number
- CN106546523A CN106546523A CN201610924504.1A CN201610924504A CN106546523A CN 106546523 A CN106546523 A CN 106546523A CN 201610924504 A CN201610924504 A CN 201610924504A CN 106546523 A CN106546523 A CN 106546523A
- Authority
- CN
- China
- Prior art keywords
- wave
- porosity
- solid
- fluid
- boundary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000007787 solid Substances 0.000 title claims abstract description 70
- 238000001514 detection method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000005284 excitation Effects 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 5
- 230000000644 propagated effect Effects 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims description 52
- 239000011148 porous material Substances 0.000 claims description 34
- 238000006073 displacement reaction Methods 0.000 claims description 29
- 239000007791 liquid phase Substances 0.000 claims description 25
- 239000007790 solid phase Substances 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 11
- 239000007788 liquid Substances 0.000 claims description 10
- 230000003028 elevating effect Effects 0.000 claims description 8
- 230000033001 locomotion Effects 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 241001274660 Modulus Species 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000013480 data collection Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 description 13
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 6
- 238000001739 density measurement Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 229910052753 mercury Inorganic materials 0.000 description 3
- 239000002904 solvent Substances 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000000326 densiometry Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 230000005428 wave function Effects 0.000 description 2
- 238000009736 wetting Methods 0.000 description 2
- 206010003694 Atrophy Diseases 0.000 description 1
- 229920002749 Bacterial cellulose Polymers 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000037444 atrophy Effects 0.000 description 1
- 239000005016 bacterial cellulose Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000009123 feedback regulation Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 230000005606 hygroscopic expansion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
- 239000011850 water-based material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Disclosed by the invention is a kind of porosity of porous medium detection means based on stream solid Interface Wave and detection method.Boundary wave is produced by interface wave excitation transducer, boundary wave is propagated along interface, angle is adjusted by slide unit is easy to hydrophone to receive boundary wave on the inclined-plane of differing tilt angles, the boundary wave for receiving is amplified to the process of PC upper computer modules by signal amplifier signal, the relation between True boundary waves velocity of wave and porosity of porous medium is extracted using algorithm disclosed by the invention, and shows that inverting obtains porous media porosity.Detection method proposed by the present invention is greatly improved the precision and reliability of detection due to the True boundary waves characteristic sensitive to porous media porosity.
Description
Technical field
The present invention relates to a kind of porosity of porous medium detection means and detection method based on stream solid Interface Wave, belongs to super
Sound detection and analysis technical field.
Background technology
With the development of ultrasonic technique, ultrasound detection is more and more frequent in the application of engineering field.Porosity is that porous is situated between
The ratio of the cumulative volume of the slight void that is interconnected in matter and the apparent volume of the porous media.Porosity is that material has in itself
An important physical, for characterizing the porous and consistency of material, it directly affects the physics and mechanicalness of material
Energy, such as strength and toughness, ventilative and water absorption etc..So that detection of the ultrasound detection to porosity of porous medium has actual meaning
Justice.
Detection porosity of porous medium method mainly has at present:
(1) density measurement drainage:What density measurement drainage was determined is support overall porosity, density measurement drainage behaviour
Make simple, be especially suitable for by non-wetting solvent of water easily penetrating in material hole, and do not cause material expand and atrophy
Frame overall porosity is determined.It is but densitometry has its limitation, it is especially useful in during hydrophilic porous material porosity measurement, such as close
The hygroscopic expansion of water-based material Bacterial cellulose so that the hydrophilic porous material solid volume and material of densitometry measurement
Cumulative volume is big with the intrinsic difference of material, causes to be calculated hydrophilic material apparent density, density of material and follow-up calculating
Material porosity inaccuracy.
(2) mercury injection apparatus method:Hydrargyrum under pressure, is clamp-oned porous material with hydrargyrum as non-wetting solvent by mercury injection apparatus method
In the hole of material, volume is immersed so as to obtain total pore size volume, by calculating brace aperture rate.The non-profit of mercury injection apparatus method
Wet solvent is fixed, and is also not suitable for flexibly determining the porosity of particularity material, and hydrargyrum is also toxic heavy metal, to material and
Experiment people's secure context is all a kind of test.
(3) scanning electron microscope sem imaging method:Support is imaged with scanning electron microscope sem, chooses five images, pass through
Matlab instruments carry out image procossing and are calculated the total void area of image and the ratio of total image area, obtain final product the hole of support
Gap rate.
As stated above, the porosity of porous media can be measured, but it is little to there is measurement range, limit many,
And do not reach permissible accuracy, it is impossible to meet the requirement of engineering detecting.
The content of the invention
The defect existing to solve the measurement of above porosity of porous medium, it is to provide one kind more that the purpose of the present invention is
For the detection means and detection method of more extensive, the safer porosity of porous medium of accurate, measurement range, utilize
True- boundary waves greatly improve the range of the precision and scope of detection to porosity of porous medium sensitivity characteristic.
In order to achieve the above object, the present invention provides a kind of porosity of porous medium detection side based on stream solid Interface Wave
Method, its method is:
(1) according to displacement and the relation and stress of potential function and the relation of potential function and the boundary condition of stream liquid/solid interface
Derive the characteristic equation of stream solid Interface Wave;
(2) it is analyzed by flowing the characteristic equation of solid Interface Wave, draws the True boundary waves velocity of sound and porous media hole
The relation of degree, and then the oscillogram of the theory T rue boundary wave velocity of sound and porosity is drawn, the oscillogram of porosity and decay;
(3) and then by actual measurement, it is finally inversed by solid porosity;The data obtained by actual device measurement, are added
Porosity obtains porous media porosity to the analysis of oscillogram inverting for decaying to the True boundary wave velocities of sound, porosity.
In the step (), the derivation for flowing the characteristic equation of solid Interface Wave is as follows:
(1) double infinitely great fluid-saturated porous medium, introduces displacement potential function phis、φf、ψs、ψfRespectively solid phase is indulged
Wave function, liquid phase compressional wave potential function potential function, solid phase shear wave potential function, liquid phase shear wave potential function, and compressional wave potential function, shear wave gesture
Function meets formula (1), (2) respectively;
Wherein, ω is angular frequency, P, Q, and R is the pore media elastic constant of Biot definition, and P=A+2N, A, N are respectively hole
Gap medium Lame constants;In the case that coefficients R is to maintain constant total volume, when the fluid of certain volume enters skeleton, convection cell is produced
Raw pressure size;Coefficient Q is the coefficient of coup that solid and liquid volume change;It is Johnson's improvement respectively
Rear solid effective density, fluid effective density and the biphase inertia coupling density of solid-liquid;
(2) to formula (1), (2) two-dimensional Fourier transform, the solid phase compressional wave of pore media on its frequency wavenumber domain can be tried to achieveWith liquid phase compressional waveSolid phase shear wave potential functionLiquid phase shear wave potential functionAs shown in formula:
B in formulas_flFor the fast compressional wave coefficient of solid phase, Cs_slFor solid phase Slow P-wave coefficient, Bf_flFor the fast compressional wave coefficient of liquid phase,
Cf_slFor liquid phase Slow P-wave coefficient, Ds_tFor solid phase shear wave coefficient, Df_tFor liquid phase shear wave coefficient;
(3) and bring formula (3), (4) into formula (1), (2) respectively to obtain:
Bf_fl=n1Bs_fl, Cf_sl=n2Cs_sl, Df_t=n3Ds_t
n1、n2And n3It is the liquid phase participation coefficient of fast compressional wave, Slow P-wave and shear wave;
cfl、csl、ctRespectively pore media is indulged soon
Ripple, Slow P-wave, transverse wave speed, and cfl=sfl -1、csl=ssl -1、ct=st -1, wherein sfl、ssl、stIt is fast compressional wave respectively, indulges slowly
The multiple slowness of ripple, shear wave;
In formula,The inertia force of the seepage flow motion equivalent density related to resistance is shown in,ρaIt is the additional mass of fluid in unit volume pore media, ρa=(a∞-1)ρfφ, η are
The coefficient of viscosity, K0For static permeability;a∞It is mass coupling coefficient, is the perfect fluid accelerated relative to solid skeleton
Apparent mass coefficient, its approximate formula are a∞=(1+ φ-1)/2;H, M, C can use pore media modulus of shearing N, solid matrix volume
Modulus Ks, fluid modulus Kf, skeleton bulk moduluses KbRepresent with porosity φ:C=α M,
Ks、Kf、Kb, N be determined by experiment, Kb, N calculated from being in harmony formula by the equivalent elastic modulus of complex media.
In the step (two), show that the True boundary waves velocity of sound is as follows with the relation of porous media porosity:
When plane wave is when fluid-porous media solid interface is that y=0 is propagated, should meet following
Boundary condition:
1) interface normal stress is continuous;Represent stress and pore-fluid stress sum of solid skeletal in porous media etc.
In interface upper flow acoustic pressure and power source sum.
In formula,Expression pore media solid skeleton normal stress,For pore-fluid normal stress,For stream
Side normal stress, its expression formula see below formula:
2) interface tangential stress is continuous;As the shear stress that fluid is surveyed is zero, Gu the condition can be reduced to solid bone
The shear stress of frame is zero;
Represent that pore media solid state shear stress, its expression formula are
3) interface medium volume conservation;Represent solid skeletal displacement and pore-fluid displacement in porous media and
Equal to interface upper flow displacement
β be porous media porosity, μys2For pore media solid phase Normal Displacement, μyf2For pore media liquid phase Normal Displacement,
μyL1For fluid side Normal Displacement, its expression formula is respectively
4) interface fluid volume conservation;What is characterized is acoustic pressure conversion and the porous that interface upper and lower fluid mutually moves product
The displacement relative equilibrium of Fluid-solid action in medium;
pL1-pf2=T β (μf2-μs2)
Represent that the relative flow rates of interface fluid are exchanged by two-way pressure to cause;Wherein β is that saturation fluid is more
Hole medium porosity, T represent flow impedance;
According to two kinds of limit modes:T=0, i.e. pL1=pf2, saturation fluid porous media perforate situation is correspond to, in space
Fluid can freedom swap with interface upper flow;T=∞, correspond to saturation fluid porous media closed pore situation, due to pL1、
pf2For infinite, then μf2=μs2;
Dielectric displacement, stress are tried to achieve using potential function relation, with reference to boundary condition, arrangement is obtained following equation:
Det (M)=0
Wherein M is:
A kind of porosity of porous medium detection means based on stream solid Interface Wave, including:
Boundary wave experimental provision, for providing boundary wave experimental situation, completes multimetering or the array received of boundary wave,
Data acquisition;
Boundary wave excitation module, for exciting the boundary wave in boundary wave experimental provision;
Signal receiving module, for the directly reception interface ripple signal in water, while measure acoustic pressure and quality in sound field shaking
Each quadrature component of speed, can completely reflect sound field information;
PC upper computer modules, for receiving signal that signal receiving module transmission comes and being analyzed and process;
Display module is calculated, the signal of PC upper computer modules feedback is calculated, and result is with data or graphic form
Shown;
The boundary wave excitation module and signal receiving module are arranged on boundary wave experimental provision, and the boundary wave is excited
Module and signal receiving module are connected with PC upper computer modules, and the PC upper computer modules connection calculates display module.
The boundary wave experimental provision includes experimental trough, the slide bar being arranged on experimental trough, slide bar angular adjustment frame
With the receiver module sliding bar being arranged on slide bar and can slide on slide bar;Described slide bar one end is provided with for measuring slide bar angle
The encoder of degree, the encoder connect PC upper computer modules, and described slide bar one end is movably arranged on experiment water by encoder
On groove, its other end is slidably connected with the slide bar angular adjustment frame, in the experimental trough is placed with solid dielectric, described
Described boundary wave excitation module is set on solid dielectric;Described receiver module sliding bar one end is slidably connected on slide bar, and which is another
One end is installed by described signal receiving module.
The slide bar angular adjustment frame includes base, the support bar arranged on base, the liter that can be slided up and down in support bar
Drop bar and sliding sleeve on upper elevating lever is set, the sliding sleeve is slidably connected slide bar, the elevating lever passes through slide block and support bar
It is slidably connected.
The boundary wave excitation module uses transducer, and the transducer is arranged on solid dielectric by wedge
On.The signal receiving module uses vector hydrophone.
Described receiver module sliding bar one end is installed on slide bar by sliding shoe, and the sliding shoe is provided with control slide block fortune
Dynamic motor, the motor connect host computer, and the slide bar uses screw thread slide bar, arrange on the sliding shoe
Have and screw thread meshed gears on screw thread slide bar, the gear connects motor.
The display module that calculates includes the computing module calculated to signal and shows result of calculation
Display module, the computing module is comprising the data initialization module for initiation parameter and for calculating porous media hole
The algoritic module of rate.
The present invention extracts the relation between True boundary waves velocity of wave and porosity of porous medium, and it is more to show that inverting is obtained
Hole medium porosity, more accurately, more extensive, the safer detection means of measurement range and detection method, using True-
Boundary wave to porosity of porous medium sensitivity characteristic, greatly improve detection method that the range of the precision and scope of detection proposes by
In the True- boundary waves characteristic sensitive to porous media porosity, the precision and reliability of detection are greatly improved.
Description of the drawings
Fig. 1 is porosity attenuation curve figure;
Fig. 2 is porosity speed curve diagram;
Fig. 3 is the schematic diagram of the detection means of the present invention;
Fig. 4 is the structural representation of the detection means of the present invention.
Specific embodiment
Below in conjunction with the accompanying drawings the device in the present invention is described further.Following examples are only used for clearly saying
Bright technical scheme, and can not be limited the scope of the invention with this.
Referring to Fig. 1 to Fig. 4, it is a kind of porosity of porous medium detection based on stream liquid/solid interface boundary wave that the present embodiment is provided
Device and detection method.
The detection method of the present invention is followed the steps below:
(1) double infinitely great fluid-saturated porous medium, introduces displacement potential function phis、φf、ψs、ψfRespectively solid phase is indulged
Wave function, liquid phase compressional wave potential function potential function, solid phase shear wave potential function, liquid phase shear wave potential function, and compressional wave potential function, shear wave gesture
Function meets formula (1), (2) respectively;
Wherein, ω is angular frequency, P, Q, and R is the pore media elastic constant of Biot definition, and P=A+2N, A, N are respectively hole
Gap medium Lame constants;In the case that coefficients R is to maintain constant total volume, when the fluid of certain volume enters skeleton, convection cell is produced
Raw pressure size;Coefficient Q is the coefficient of coup that solid and liquid volume change;It is Johnson's improvement respectively
Rear solid effective density, fluid effective density and the biphase inertia coupling density of solid-liquid.
(2) to formula (1), (2) two-dimensional Fourier transform, pore media compressional wave on its frequency wavenumber domain can be tried to achieveShear wave potential functionAs shown in formula
B in formulas_flFor the fast compressional wave coefficient of solid phase, Cs_slFor solid phase Slow P-wave coefficient, Bf_flFor the fast compressional wave coefficient of liquid phase,
Cf_slFor liquid phase Slow P-wave coefficient, Ds_tFor solid phase shear wave coefficient, Df_tFor liquid phase shear wave coefficient.
(3) and bring formula (3), (4) into formula (1), (2) respectively to obtain
Bf_fl=n1Bs_fl, Cf_sl=n2Cs_sl, Df_t=n3Ds_t
n1、n2And n3It is the liquid phase participation coefficient of fast compressional wave, Slow P-wave and shear wave
cfl、csl、ctRespectively pore media is indulged soon
Ripple,
Slow P-wave, transverse wave speed, and cfl=sfl -1、csl=ssl -1、ct=st -1, wherein sfl、ssl、stBe respectively fast compressional wave,
The multiple slowness of Slow P-wave, shear wave
In formula,The inertia force of the seepage flow motion equivalent density related to resistance is shown in,ρaIt is the additional mass of fluid in unit volume pore media, ρa=(a∞-1)ρfφ, η are
The coefficient of viscosity, K0For static permeability;a∞It is mass coupling coefficient, is the perfect fluid accelerated relative to solid skeleton
Apparent mass coefficient, its approximate formula are a∞=(1+ φ-1)/2.H, M, C can use pore media modulus of shearing N, solid matrix volume
Modulus Ks, fluid modulus Kf, skeleton bulk moduluses KbRepresent with porosity φ:C=α M,Ks、Kf、Kb, N can be determined by experiment, Kb, N also can be by the equivalent bullet of complex media
Property modulus calculate from being in harmony formula.
(4) when plane wave is propagated in fluid-porous media solid interface (y=0), should meet with downstream condition
5) interface normal stress is continuous.Represent stress and pore-fluid stress sum of solid skeletal in porous media etc.
In interface upper flow acoustic pressure and power source sum.
In formula,Expression pore media solid skeleton normal stress,For pore-fluid normal stress,For stream
Side normal stress, its expression formula see below formula:
6) interface tangential stress is continuous.As the shear stress that fluid is surveyed is zero, Gu the condition can be reduced to solid bone
The shear stress of frame is zero.
Represent that pore media solid state shear stress, its expression formula are
7) interface medium volume conservation.Represent solid skeletal displacement and pore-fluid displacement in porous media and
Equal to interface upper flow displacement
β be porous media porosity, μys2For pore media solid phase Normal Displacement, μyf2For pore media liquid phase Normal Displacement,
μyL1For fluid side Normal Displacement, its expression formula is respectively
8) interface fluid volume conservation.What is characterized is acoustic pressure conversion and the porous that interface upper and lower fluid mutually moves product
The displacement relative equilibrium of Fluid-solid action in medium.
pL1-pf2=T β (μf2-μs2)
Represent that the relative flow rates of interface fluid are exchanged by two-way pressure to cause.Wherein β is that saturation fluid is more
Hole medium porosity, T represent flow impedance.Generally, it is contemplated that two kinds of limit modes:T=0, i.e. pL1=pf2, correspondence
Saturation fluid porous media perforate situation, in space fluid can freedom swap with interface upper flow;T=∞, correspondence
Saturation fluid porous media closed pore situation, due to pL1、pf2For infinite, then μf2=μs2。
Dielectric displacement, stress are tried to achieve using potential function relation, with reference to boundary condition, arrangement is obtained following equation:
Det (M)=0
Wherein M is:
Above-mentioned equation be one with regard to True- boundary waves wave number and the relational expression of porosity of porous medium, using above-mentioned
Liquid-solid interface True- interfaces wave equation can be drawing the porosity rate curve of True boundary waves, porosity attenuation curve
Figure, such as Fig. 1 and Fig. 2:
The data for obtaining are measured by device is actual again, you can com-parison and analysis are finally inversed by porous media porosity.
It is of the invention main using the relation for flowing solid True boundary waves velocity of wave and porosity of porous medium, draw theory T rue circle
The porosity rate curve of face ripple, porosity attenuation curve, by actual measurement, are finally inversed by solid porosity.
As shown in figure 3, detection means of the present invention based on the porosity of porous medium of stream liquid/solid interface True- boundary waves is main
By boundary wave excitation module 1, boundary wave experimental provision 2, signal receiving module 3, PC upper computer modules 4, display module 5 five is calculated
Part constitutes.
In the present embodiment, boundary wave excitation module 1 and signal receiving module 3 are arranged on boundary wave experimental provision 3, described
Boundary wave excitation module 1 and signal receiving module 3 are connected with PC upper computer modules 4, the connection of PC upper computer modules 4 meter
Calculate display module 5.Wherein, boundary wave experimental provision 2 includes experimental trough 10, the slide bar 11 being arranged on experimental trough 10, slides
Bar angular adjustment frame and the receiver module sliding bar 131 that be arranged on slide bar 11 and can slide on slide bar 11;The slide bar 11
One end is provided with encoder 12, and the encoder 12 connects PC upper computer modules 4, and 11 one end of the slide bar is movable by encoder 12
On experimental trough 10, it is easy to measure slide bar angle, more than 2000 lines, error is less than 0.1 degree to the encoder;Slide bar 11
The other end is slidably connected with the slide bar angular adjustment frame, in the experimental trough 10 is placed with solid dielectric 20, described solid
Described boundary wave excitation module 1 is set on body medium 20;Described 131 one end of receiver module sliding bar is connected by sliding shoe 132
Connect on slide bar 11, its other end installs described signal receiving module 2.The step of control slide block motion is provided with sliding shoe 132
Stepper motor, motor connection host computer, by its operation of PC control, from control operation of the sliding shoe on slide bar.
Above-mentioned slide bar angular adjustment frame include base 141, arrange base 141 on support bar 142, can be in support bar 142
The elevating lever 144 that slides up and down and sliding sleeve 145 on upper elevating lever 144 is set, the sliding sleeve 145 is slidably connected slide bar 11,
The elevating lever 144 is slidably connected with support bar 142 by slide block 143.The signal receiving module 2 of the present embodiment uses arrow
Amount hydrophone, boundary wave excitation module 1 use transducer, and the transducer is arranged on solid dielectric 20 by wedge 21
On.The present embodiment is capable of achieving the directly reception interface ripple signal in water by adopting vector hydrophone for signal receive transducer,
Installation receiving device need not be carried out at Liquid-solid interface, operational approach is relatively simple flexibly, and can using vector hydrophone
Each quadrature component of acoustic pressure and quality vibration velocity in sound field is measured simultaneously, can completely be reflected sound field information, functionally be connect than single
Receive transducer more powerful;Again signal is received by being analyzed in PC upper computer modules, process.
Support bar 142 and slide bar 11 in the present embodiment uses screw thread slide bar, in sliding shoe 132 and slide block 143
On be provided with and screw thread meshed gears on screw thread slide bar, the gear connects motor, step motor control gear
Operation, motor connection PC upper computer modules, controls the operation of motor by PC upper computer modules.
The present invention is manually put into solid dielectric in experimental trough with arbitrarily angled, PC upper computer module input angles
Degree, you can it is arbitrarily angled that control sliding shoe and slide block movement tilt to the slide bar on tank, and encoder measures the angle of slide bar,
Feedback regulation makes angle more accurate.By the movement of step motor control silk slide bar platform, the movement of vector hydrophone is driven, can be real
Existing multimetering or array received, data acquisition amount are big, and high degree of automation, signals collecting efficiency are greatly improved.
Additionally, the display module that calculates of the present embodiment includes the computing module that calculated to signal and by result of calculation
The display module for being shown, the computing module is comprising the data initialization module for initiation parameter and for calculating many
The algoritic module of hole medium porosity.The present invention is calculated to signal by the algorithm proposed by computing module, and result is passed
Enter microprocessor, display module is received the data or figure of microprocessor transmission, shown by display screen.
The present invention produces boundary wave by interface wave excitation transducer, and boundary wave is propagated along interface, adjusts angle just by slide unit
Boundary wave on the inclined-plane of differing tilt angles is received in hydrophone, the boundary wave for receiving is amplified to PC by signal amplifier signal
Upper computer module process, is extracted between True boundary waves velocity of wave and porosity of porous medium using algorithm disclosed by the invention
Relation, and show that inverting obtains porous media porosity.Detection method proposed by the present invention is situated between to porous due to True- boundary waves
The sensitive characteristic of matter porosity, is greatly improved the precision and reliability of detection.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, on the premise of without departing from the technology of the present invention principle, some improvement and deformation can also be made, these improve and deform
Also should be regarded as protection scope of the present invention.
Claims (10)
1. a kind of porosity of porous medium detection method based on stream solid Interface Wave, is characterized in that, its method is:
(1) boundary condition according to displacement with the relation and stress of potential function with the relation of potential function and stream liquid/solid interface is derived
Go out to flow the characteristic equation of solid Interface Wave;
(2) it is analyzed by flowing the characteristic equation of solid Interface Wave, draws the True boundary waves velocity of sound with porous media porosity
Relation, and then the oscillogram of the theory T rue boundary wave velocity of sound and porosity is drawn, the oscillogram of porosity and decay;
(3) and then by actual measurement, it is finally inversed by solid porosity;The data obtained by actual device measurement, add hole
Degree obtains porous media porosity to the analysis of oscillogram inverting for decaying to the True boundary wave velocities of sound, porosity.
2. it is according to claim 1 based on the porosity of porous medium detection method for flowing solid Interface Wave, it is characterised in that institute
State in step (), the derivation for flowing the characteristic equation of solid Interface Wave is as follows:
(1) double infinitely great fluid-saturated porous medium, introduces displacement potential function phis、φf、ψs、ψfRespectively solid phase compressional wave letter
Number, liquid phase compressional wave potential function potential function, solid phase shear wave potential function, liquid phase shear wave potential function, and compressional wave potential function, shear wave potential function
Formula (1), (2) are met respectively;
Wherein, ω is angular frequency, and P, Q, R are the pore media elastic constant of Biot definition, and P=A+2N, A, N are respectively hole and are situated between
Matter Lame constants;In the case that coefficients R is to maintain constant total volume, when the fluid of certain volume enters skeleton, convection cell is produced
Pressure size;Coefficient Q is the coefficient of coup that solid and liquid volume change;It is after Johnson improves respectively
The biphase inertia coupling density of solid effective density, fluid effective density and solid-liquid;
(2) to formula (1), (2) two-dimensional Fourier transform, the solid phase compressional wave of pore media on its frequency wavenumber domain can be tried to achieveWith
Liquid phase compressional waveSolid phase shear wave potential functionLiquid phase shear wave potential functionAs shown in formula:
B in formulas_flFor the fast compressional wave coefficient of solid phase, Cs_slFor solid phase Slow P-wave coefficient, Bf_flFor the fast compressional wave coefficient of liquid phase, Cf_slFor liquid
Phase Slow P-wave coefficient, Ds_tFor solid phase shear wave coefficient, Df_tFor liquid phase shear wave coefficient;
(3) and bring formula (3), (4) into formula (1), (2) respectively to obtain:
Bf_fl=n1Bs_fl, Cf_sl=n2Cs_sl, Df_t=n3Ds_t
n1、n2And n3It is the liquid phase participation coefficient of fast compressional wave, Slow P-wave and shear wave;
cfl、csl、ctRespectively the fast compressional wave of pore media,
Slow P-wave, transverse wave speed, and cfl=sfl -1、csl=ssl -1、ct=st -1, wherein sfl、ssl、stBe respectively fast compressional wave, Slow P-wave,
The multiple slowness of shear wave;
In formula, The inertia force of the seepage flow motion equivalent density related to resistance is shown in,ρaIt is the additional mass of fluid in unit volume pore media, ρa=(a∞-1)ρfφ, η are
The coefficient of viscosity, K0For static permeability;a∞It is mass coupling coefficient, is the perfect fluid accelerated relative to solid skeleton
Apparent mass coefficient, its approximate formula are a∞=(1+ φ-1)/2;H, M, C can use pore media modulus of shearing N, solid matrix volume
Modulus Ks, fluid modulus Kf, skeleton bulk moduluses KbRepresent with porosity φ:C=α M,
Ks、Kf、Kb, N be determined by experiment, Kb, N calculated from being in harmony formula by the equivalent elastic modulus of complex media.
3. it is according to claim 2 based on the porosity of porous medium detection method for flowing solid Interface Wave, it is characterised in that institute
State in step (two), show that the True boundary waves velocity of sound is as follows with the relation of porous media porosity:
When plane wave is when fluid-porous media solid interface is that y=0 is propagated, should meet with downstream condition:
1) interface normal stress is continuous;The stress and pore-fluid stress sum for representing solid skeletal in porous media is equal to boundary
Face upper flow acoustic pressure and power source sum.
In formula,Expression pore media solid skeleton normal stress,For pore-fluid normal stress,For fluid side
Normal stress, its expression formula see below formula:
2) interface tangential stress is continuous;As the shear stress that fluid is surveyed is zero, Gu the condition can be reduced to solid skeletal
Shear stress is zero;
Represent that pore media solid state shear stress, its expression formula are3) interface medium volume
Conservation;Represent solid skeletal displacement and pore-fluid displacement in porous media and it is equal to interface upper flow displacement
β be porous media porosity, μys2For pore media solid phase Normal Displacement, μyf2For pore media liquid phase Normal Displacement, μyL1
For fluid side Normal Displacement, its expression formula is respectively
4) interface fluid volume conservation;What is characterized is acoustic pressure conversion and the porous media that interface upper and lower fluid mutually moves product
The displacement relative equilibrium of middle Fluid-solid action;pL1-pf2=T β (μf2-μs2)
Represent that the relative flow rates of interface fluid are exchanged by two-way pressure to cause;Wherein β is situated between for saturation fluid porous
Matter porosity, T represent flow impedance;
According to two kinds of limit modes:T=0, i.e. pL1=pf2, correspond to saturation fluid porous media perforate situation, fluid in space
Can freedom swap with interface upper flow;T=∞, correspond to saturation fluid porous media closed pore situation, due to pL1、pf2
For infinite, then μf2=μs2;
Dielectric displacement, stress are tried to achieve using potential function relation, with reference to boundary condition, arrangement is obtained following equation:
Det (M)=0
Wherein M is:
4. it is a kind of based on the porosity of porous medium detection means for flowing solid Interface Wave, it is characterised in that to include:
Boundary wave experimental provision, for providing boundary wave experimental situation, completes multimetering or the array received of boundary wave, data
Collection;
Boundary wave excitation module, for exciting the boundary wave in boundary wave experimental provision;
Signal receiving module, for the directly reception interface ripple signal in water, while measuring acoustic pressure and quality vibration velocity in sound field
Each quadrature component, can completely reflect sound field information;
PC upper computer modules, for receiving signal that signal receiving module transmission comes and being analyzed and process;
Display module is calculated, the signal of PC upper computer modules feedback is calculated, and result is carried out with data or graphic form
Show;
The boundary wave excitation module and signal receiving module are arranged on boundary wave experimental provision, the boundary wave excitation module
It is connected with PC upper computer modules with signal receiving module, the PC upper computer modules connection calculates display module.
5. porosity of porous medium detection means according to claim 4, it is characterised in that the boundary wave experimental provision
Including experimental trough, the slide bar being arranged on experimental trough, slide bar angular adjustment frame and it is arranged on slide bar and can be on slide bar
The receiver module sliding bar of slip;Described slide bar one end is provided with for measuring the encoder of slide bar angle, the encoder connection
PC upper computer modules, described slide bar one end are movably arranged on experimental trough by encoder, its other end and the slide bar angle
Degree adjusting bracket is slidably connected, and solid dielectric is placed with the experimental trough, and described interface is arranged on the solid dielectric
Ripple excitation module;Described receiver module sliding bar one end is slidably connected on slide bar, and its other end is installed described signal and receives mould
Block.
6. porosity of porous medium detection means according to claim 5, it is characterised in that the slide bar angular adjustment frame
Including the cunning on elevating lever in base, the support bar arranged on base, the elevating lever that can be slided up and down in support bar and setting
Set, the sliding sleeve are slidably connected slide bar, and the elevating lever is slidably connected with support bar by slide block.
7. porosity of porous medium detection means according to claim 5, it is characterised in that the boundary wave excitation module
Transducer is used, the transducer is arranged on solid dielectric by wedge.
8. porosity of porous medium detection means according to claim 5, it is characterised in that the signal receiving module is adopted
It is vector hydrophone.
9. porosity of porous medium detection means according to claim 8, it is characterised in that the receiver module sliding bar
One end is installed on slide bar by sliding shoe, and the sliding shoe is provided with the motor of control slide block motion, the motor
Connection host computer, the slide bar are used screw thread slide bar, are provided with what is engaged with screw thread on screw thread slide bar on the sliding shoe
Gear, the gear connect motor.
10. porosity of porous medium detection means according to claim 5, it is characterised in that the calculating display module
The display module shown including the computing module calculated to signal and by result of calculation, the computing module are included
Data initialization module for initiation parameter and the algoritic module for calculating porosity of porous medium.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610924504.1A CN106546523B (en) | 2016-10-24 | 2016-10-24 | The detection device and method of porosity of porous medium based on stream solid Interface Wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610924504.1A CN106546523B (en) | 2016-10-24 | 2016-10-24 | The detection device and method of porosity of porous medium based on stream solid Interface Wave |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106546523A true CN106546523A (en) | 2017-03-29 |
CN106546523B CN106546523B (en) | 2019-08-30 |
Family
ID=58393061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610924504.1A Expired - Fee Related CN106546523B (en) | 2016-10-24 | 2016-10-24 | The detection device and method of porosity of porous medium based on stream solid Interface Wave |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106546523B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109871623A (en) * | 2019-02-25 | 2019-06-11 | 西安交通大学 | A kind of porous media model design method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005103766A2 (en) * | 2004-04-23 | 2005-11-03 | Schlumberger Canada Limited | Method and system for monitoring of fluid-filled domains in a medium based on interface waves propagating along their surfaces |
CN102914590A (en) * | 2011-08-02 | 2013-02-06 | 通用电气公司 | Method and apparatus for porosity measurement and defect detection |
CN104749082A (en) * | 2015-03-31 | 2015-07-01 | 南昌航空大学 | Ultrasonic multifunctional evaluation method and ultrasonic multifunctional evaluation device for void content |
-
2016
- 2016-10-24 CN CN201610924504.1A patent/CN106546523B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005103766A2 (en) * | 2004-04-23 | 2005-11-03 | Schlumberger Canada Limited | Method and system for monitoring of fluid-filled domains in a medium based on interface waves propagating along their surfaces |
CN102914590A (en) * | 2011-08-02 | 2013-02-06 | 通用电气公司 | Method and apparatus for porosity measurement and defect detection |
CN104749082A (en) * | 2015-03-31 | 2015-07-01 | 南昌航空大学 | Ultrasonic multifunctional evaluation method and ultrasonic multifunctional evaluation device for void content |
Non-Patent Citations (3)
Title |
---|
QINGBANG HAN ET AL.: "Propagation characteristics of interface waves between a porous medium and a sediment-containing two-phase fluid", 《ULTRASONICS》 * |
乔文孝 等: "多孔介质声学研究进展", 《物理学进展》 * |
刘仲一 等: "一种孔隙介质声波模型正演方法", 《地质评论》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109871623A (en) * | 2019-02-25 | 2019-06-11 | 西安交通大学 | A kind of porous media model design method |
CN109871623B (en) * | 2019-02-25 | 2020-10-30 | 西安交通大学 | Porous medium model design method |
Also Published As
Publication number | Publication date |
---|---|
CN106546523B (en) | 2019-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106226392A (en) | Water-oil phase flow containing rate measuring method based on ultrasonic attenuation mechanism model | |
Langlois et al. | Polynomial relations for quasi-static mechanical characterization of isotropic poroelastic materials | |
CN101696925B (en) | Device and method for testing performance of two-phase flow drag reducer | |
Reinsch et al. | A new technique for the simultaneous, real-time measurement of membrane compaction and performance during exposure to high-pressure gas | |
CN107636425A (en) | It is determined that assess method, computer program product and the device of the characteristic variable for the measurement arrangement for including clamp-on ultrasonic flowmeter and measurement pipe | |
CN106287240A (en) | A kind of pipeline leakage testing device based on acoustic emission and single-sensor localization method | |
CN102252816A (en) | Testing method of vortex-induced vibration of cylinder | |
CN102788738A (en) | Ultrasonic array detection method for multi-phase liquid density and concentration | |
Longo | Experiments on turbulence beneath a free surface in a stationary field generated by a Crump weir: free-surface characteristics and the relevant scales | |
CN103529237A (en) | Method and device for measuring sediment group settling velocity | |
CN101344505A (en) | On-line monitoring method and device for pollution of hollow fiber film | |
CN102606135A (en) | Detecting system for mass and flow of return solid particles in oil drilling shaft and detecting method | |
CN110530765A (en) | Underwater bubble group's size distribution parameter inversion method based on measuring non-linear parameters | |
CN105784707A (en) | Gas-liquid two-phase flow measurement method based on high-speed photography technology and device thereof | |
Li et al. | Cake-layer deposition, growth, and compressibility during microfiltration measured and modeled using a noninvasive ultrasonic technique | |
CN106546523A (en) | Detection means and method based on the porosity of porous medium of stream solid Interface Wave | |
CN103604728A (en) | Sand grain diameter detection device and method based on fluid/solid interface wave | |
CN114279940A (en) | Visual teaching system for shear seepage characteristics of fractured rock mass | |
CN110631958A (en) | Gas-liquid two-phase flow parameter measuring method | |
CN204514402U (en) | A kind of differential pressure mass flowmeter for vortex street | |
CN104568651B (en) | Calculating method of on-line double-pressure vibrating tube type drilling fluid density measurement instrument | |
CN110320341A (en) | Fluid phase state test device and experimental method in porous media | |
CN106597546A (en) | Method for detecting river bottom surface sediment physical properties | |
Jin et al. | A blind deconvolution method for attenuative materials based on asymmetrical Gaussian model | |
CA2868978A1 (en) | Speed of sound and/or density measurement using acoustic impedance |
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: 20190830 |