NL1024984C2 - Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. - Google Patents
Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. Download PDFInfo
- Publication number
- NL1024984C2 NL1024984C2 NL1024984A NL1024984A NL1024984C2 NL 1024984 C2 NL1024984 C2 NL 1024984C2 NL 1024984 A NL1024984 A NL 1024984A NL 1024984 A NL1024984 A NL 1024984A NL 1024984 C2 NL1024984 C2 NL 1024984C2
- Authority
- NL
- Netherlands
- Prior art keywords
- reflection
- probability
- measured
- particles
- factor
- Prior art date
Links
- 239000002245 particle Substances 0.000 title claims description 126
- 238000009826 distribution Methods 0.000 title claims description 66
- 238000000034 method Methods 0.000 title claims description 34
- 239000007788 liquid Substances 0.000 title claims description 23
- 238000005259 measurement Methods 0.000 claims description 37
- 230000014509 gene expression Effects 0.000 claims description 26
- 238000007476 Maximum Likelihood Methods 0.000 claims description 7
- 230000000873 masking effect Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 239000002699 waste 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/02—Investigating particle size or size distribution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02416—Solids in liquids
Landscapes
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Separation Of Solids By Using Liquids Or Pneumatic Power (AREA)
Description
»
Titel: Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof
De uitvinding betreft een werkwijze en inrichting waarin, met behulp van weerkaatsing van een ultrasone pulsen door macroscopische deeltjes zoals van oliedruppels of zandkorrels in water, informatie over de grootteverdeling van de deeltjes in een vloeistof bepaald wordt.
5
Europese octrooiaanvrage No. 0801305 beschrijft een inrichting voor het karakteriseren van een suspensie. De inrichting wekt ultrasone pulsen op in een bundel in een vloeistof. Reflectie van een puls aan een individueel deeltje in de vloeistof resulteert in een echo. Als het deeltje zich 10 op een gegeven positie bevindt is de amplitude van de reflectie één op één aan de grootte van het deeltje gekoppeld. Bij metingen aan reflecties van achtereenvolgende pulsen op deeltjes van verschillende grootte worden zodoende verschillende amplitudes gemeten.
De bekende inrichting maakt een histogram van tellingen van het 15 aantal keren dat reflecties met verschillende amplitude voorkomen. Het histogram bevat informatie over de deeltjesgrootteverdeling, dat wil zeggen de concentratie van deeltjes als functie van de grootte van de deeltjes. Met informatie over de deeltjesgrootteverdeling op zijn beurt kan bijvoorbeeld de totale concentratie van deeltjes in de vloeistof bepaald worden.
20 Europese octrooiaanvrage No. 0801305 beschrijft een techniek om informatie over de grootte uit de verdeling van de amplitudes te halen. Er wordt een uitdrukking gebruikt voor de waarschijnlijkheid dat reflecties met verschillende amplitudes optreden. Deze uitdrukking relateert de waarschijnlijkheid aan de deeltjesgrootteverdeling. Als de reflecterende 25 deeltjes zich altijd op dezelfde plaats in de bundel zouden bevinden dan zou het aantal reflecties met een bepaalde amplitude evenredig zijn met de fractie aan deeltjes met de grootte die tot die amplitude leidt.
1024984 2 «
De informatie over de deeltjesgrootte verdeling zit echter versmeerd over de amplitudes, omdat de reflectieamplitude van een deeltje behalve van de grootte van het deeltjes ook afhankelijk is van de positie van het deeltje in de bundel. Hoe verder het deeltje zich bij reflectie uit het 5 centrum van de bundel bevindt, hoe kleiner de amplitude. De gebruikte uitdrukking voor de waarschijnlijkheid van een reflectiemeting met een bepaalde amplitude versmeert de deeltjesgrootteverdeling over de amplitudeverdeling om dit effect tot uitdrukking te brengen.
Met behulp van een maximum likelihood techniek worden 10 parameters van de deeltjesgrootte verdeling geschat die de zo uitgedrukte waarschijnlijkheid van de daadwerkelijk gemeten aantallen deeltjes maximaliseren. Zodoende wordt de versmering ongedaan gemaakt. In concreto worden de parameters zo gekozen dat een som geminimaliseerd wordt van de kwadraten van de verschillen tussen gemeten aantallen 15 reflecties met amplitudes in verschillende amplitudebereiken en voorspelde aantallen. Een dergelijke som is een indicatie voor de waarschijnlijkheid van de combinatie van gemeten reflecties, maar uiteraard kunnen ook andere indicaties van de waarschijnlijkheid gebruikt worden, zoals het product van de waarschijnlijkheden van de verschillende gemeten 20 amplitudes.
Een noodzakelijk voorwaarde voor de gebruikte techniek ie dat de reflecties van individuele deeltjes onderscheiden kunnen worden. Als reflecties van meerdere deeltjes door elkaar gemeten worden, dan raakt een deel van de reflecties gemaskeerd. Zo kan reflectie van een groter deeltje 25 reflectie van een kleiner deeltje onzichtbaar maken, maar ook kan reflectie van een kleiner deeltje reflectie van een groter deeltje maskeren als het grotere deeltje zich voldoende veel verder uit het centrum van de bundel bevindt dan het kleinere deeltje.
Zodoende geeft de techniek alleen bij voldoende lage concentraties 30 betrouwbare resultaten. De maximaal bruikbare concentratie kan verhoogd
102 498/J
♦ 3 worden door het minimaliseren van het volume waarin deeltjes ononderscheidbaar gemeten worden, bijvoorbeeld door het gebruik van een gefocusseerde bundel en kleine tijdswindows waarin echo's geaccepteerd worden. Er zijn echter grenzen aan de toepasbaarheid van dergelijke 5 technieken. Daarom beperkt de eis van gescheiden waarneming van individuele deeltjes de toepasbaarheid van de techniek.
Het is een doel van de uitvinding om het gebruik van ultrasone metingen voor de karakterisatie van deeltjesgrootteverdeling in vloeistoffen 10 bruikbaar te maken voor hogere concentraties.
De huidige uitvinding berust op de realisatie.dat het mogelijk is om met een statistisch model het effect van het missen van deeltjes te corrigeren en dat daardoor de bekende meettechniek ook voor grotere concentraties gebruikt kan worden.
15 De uitvinding voorziet in een werkwijze volgens de conclusie 1. In de werkwijze worden de parameters van de deeltjesgrootte verdeling geschat met een maximum likelihood techniek, waarbij gebruik gemaakt wordt van een uitdrukking voor de waarschijnlijkheid van meting van reflectieamplitudes in termen van twee factoren. De eerste factor drukt de 20 waarschijnlijkheid P0(A) uit van een reflectiemeting waarvan een reflectie met de waarde A onderdeel van uitmaakt, onafhankelijk van de vraag of deze waarde gemaskeerd wordt doordat tevens een reflectie met een grotere waarde optreedt. Deze eerste factor kom in wezen overeen met de hele in de stand der techniek gebruikte uitdrukking voor de waarschijnlijkheid.
25 De tweede factor in de volgens de uitvinding gebruikte uitdrukking omvat de waarschijnlijkheid dat niet tevens een reflectie met een dominerende waarde van de eigenschap, die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting. Zodoende wordt het effect van maskering bij de schatting van de parameters meegenomen. Bij de 1024984 4 ι maximum ïikelihood schatting wordt zodoende zowel het effect van versmering als het effect van maskering samen meegenomen.
Deze en andere doelen en voordelige aspecten van de uitvinding zullen δ worden beschreven aan de hand van de volgende figuren
Figuur 1 toont een inrichting voor het karakteriseren van een vloeistof Figuur 2 toont aantallen reflecties met verschillende amplitudes Figuur 3 toont een gereconstrueerde deeltjesgrootteverdeling 10 Figuur 4 toont een berekende concentratie
Figuur 5 toont een flow chart van een werkwijze om parameters van een deeltjesgrootteverdeling te schatten »
Figuur 1 toont een inrichting voor het karakteriseren van een vloeistof, De 15 inrichting omvat een vloeistofkanaal 11, controle eenheid 10, een signaalgenerator 12, een ultrasone transducer 14, een detector 16 en een rekeneenheid 18. Controle eenheid 10 is aan signaalgenerator 12, detector 16 en rekeneenheid 18 gekoppeld. Signaalgenerator 12 hééft een uitgang die aan transducer 14 gekoppeld is. Detector 16 heeft een ingang die aan 20 transducer 14 gekoppeld is. Detector 16 heeft een uitgang die aan rekeneenheid 18 gekoppeld is. Transducer 14 is aan kanaal 11 gekoppeld voor het opwekken van een ultrasone bundel.
In bedrijf stroomt een vloeistof met daarin zwevende deeltjes door kanaal 11. Controle eenheid 10 triggert achtereenvolgende reflectiemetingen in elk 25 waarvan signaalgenerator 12 een puls genereert die door transducer 14 in een ultrasone golf in kanaal 11 wordt omgezet. Transducer 14 vangt reflecties van de puls op en voert resulterende signalen aan detector 16 toe. Detector 11 selecteert signalen die in een bepaald tijdswindow ten opzichte van de opgewekt puls optreden. Als gevolg van de geluidssnelheid in 30 vloeistof definieert het tijdswindow een zone vanwaaruit deeltjes aanleiding 1024984 5 kunnen geven tot reflecties. Detector 16 meet de amplitude van het reflectiesignaal (indien aanwezig) in het tijdswindow en voert de gemeten waarde van de amplitude toe aan rekeneenheid 18.
Rekeneenheid 18 verzamelt histogram informatie van de amplitudes.
5 Hiertoe gebruikt rekeneenheid een aantal amplituderanges en accumuleert rekeneenheid tellingen van het aantal reflectiemetingen die reflectieamplitudes in de respectieve ranges hebben opgeleverd. Nadat zodoende een groot aantal reflectiemetingen ie uitgevoerd berekent rekeneenheid 18 uit de histogram informatie verdere informatie over de 10 deeltjesgrootte verdeling.
Figuur 2 toont een gesimuleerd voorbeeld van histogrammen voor vloeistoffen waarin verschillende concentraties van deeltjes aanwezig zijn, terwijl die deeltjesgrootteverdeling fb(D) telkens hetzelfde is (dat wil zeggen dat de fractie van de deeltjes met een grootte tussen D en D+dD telkens 15 fo(D)dD is voor elke waarde van D). Zichtbaar is dat, behalve een toename van de aantallen reflecties ook de vorm van de amplitude verdeling verandert bij toenemende concentratie.
Figuur 3 toont met gestippelde lijnen schattingen van de deeltjesgrootte verdeling die uit de amplitude histogrammen kan worden 20 berekend op basis van de bekende techniek, voor een aantal verschillende concentraties. Tevens wordt de werkelijke verdeling getoond. De bekende techniek houdt geen rekening met het effect van de concentratie op de vorm van de amplitudeverdeling. Het resultaat is dat de geschatte deeltjesgrootteverdelingen onderling verschillen en bij hogere concentraties 25 afwijken van de echte deeltjesgrootteverdeling.
De. uitvinding voorziet in een verbetering van de bekende techniek om de deeltjesgrootteverdeling te schatten, zodat ook voor het effect van de concentratie van de deeltjes gecorrigeerd kan worden. Deze verbetering berust op de realisatie dat de afwijking van de vorm van de histogrammen 30 van figuur 2 een gevolg is van maskering van reflectieamplitudes die bij de 1024984 ♦ 6 een reflectiemeting wel optreden, maar in de telling gemist worden omdat bij dezelfde reflectiemeting één of meer reflecties met een grotere amplitude optreden.
Daartoe wordt een uitdrukking gebruik voor de 5 waarschijnlijkheidsverdeling P(A)dA dat een reflectiemeting een gemeten amplitude oplevert die in een infmitesimaal interval tussen A en A+dA ligt. Deze uitdrukking beeft de vorm van P(A)=P0(A) Q (A) 10
Hierin correspondeert P0(A) met de ongemaskeerde, waarschijnlijkheidsverdeling, dat wil zeggen waarschijnlijkheid van aanwezigheid van een reflectie, los van de vraag of deze gemaskeerd wordt door een andere reflectie met een grotere amplitude. De factor Q(A) drukt de 15 waarschijnlijkheid uit dat er in de dezelfde reflectiemeting geen deeltjes reflecties met grotere of gelijke amplitude dan A veroorzaken.
De factor Po(A) komt overeen met de waarschijnlijkheidsverdeling die in de stand der techniek gebruikt wordt voor de gehele waarschijnlijkheid verdeling P(A) dat een amplitude A in het interval A tot 20 A+dA ook daadwerkelijk gemeten wordt (dwz onder verwaarlozing van de mogelijkheid dat een reflectie gemaskeerd wordt). Deze factor P0(A) kan op op zich bekende wijze in de deeltjesgrootteverdeling fb(D) worden uitgedrukt. Een uitdrukking voor deze factor is bijvoorbeeld 25 P0(A)=CjdD&(D)G(A|D)Vmea8(D)
Hierin is C de concentratie van de deeltjes (het gemiddeld aantal deeltjes per volume-eenheid). G(A| D) is de conditionele waarschijnlijkheid dat er een reflectie met amplitude A detecteerbaar is, als een deeltje met een 30 grootte D een detecteerbare reflectie oplevert. De conditionele 1024984 7..
waarschijnlijkheid G(A | D) drukt het effect uit dat een deeltje met een grootte D op verschillende afstanden van het centrum van de bundel aanleiding geeft tot verschillende reflectie amplitudes. Vmeas(D) is het totale vloeistofvolume vanwaaruit een deeltje met grootte D een detecteerbare 5 reflectie kan opleveren. Dit volume wordt in de diepte gedefinieerd door het tijdswindow vanwaaruit tijdens de meting reflecties geaccepteerd worden en in de breedte door de afval van de bundelintensiteit tot een punt waarop de amplitude van de opgewekte reflecties niet meer boven het ruisniveau An uitsteekt.
10 In het veel voorkomende geval van een Gaussisch verdeelde bundel kan het volume Vmea8(D) worden uitgedrukt als
Vmeas(D)=K Az log(Ao(D)/ Αη)/2κ 15 Hierin is Az het diepteinterval waaruit reflecties geaccepteerd worden, κ is de standaard deviatie die bepaald wordt door de snelheid waarmee de intensiteit van de bundel als functie van de afstand tot het centrum van de bundel afvalt. Ao(D) is de maximale amplitude die een reflectie van een deeltje met grootte D zou opleveren als het deeltje zich in het centrum van 20 de bundel zou bevinden.
Voor een Gaussische bundel kan ook een uitdrukking voor de conditionele waarschijnlijkheid G(A | D) worden gegeven: G(A | D)= l/{ A*log (Ao(D)/An)} 25
Deze uitdrukking geldt mits Ao(D)>A>An. Voor amplitudes groter dan de maximaal met een deeltjesgrootte D haalbare amplitude Ao(D) of kleiner dan het ruisniveau An is de conditionele kans op detectie G(A | D) gelijk nul.
1024984 / 8 ι
Het zal overigens duidelijk zijn dat voor bundels met een andere dan Gaussische verdeling overeenkomstig andere uitdrukkingen voor Vmeas(D) en G(A | D) gebruikt zullen worden.
De nieuw geïntroduceerde factor Q(A), die het effect van maskering 5 van reflecties door andere reflecties met een grotere of gelijke amplitude dan de amplitude A uitdrukt, kan eveneens worden uitgedrukt in termen van de concentratie G en de deeltjesgrootteverdeling fb(D). Maatgevend hiervoor is de integraaluitdrukking 10 W (A) = IdA’idD fb(D) G(A' | D)
Hierin wordt de integraal over de amplitude genomen van het ruisniveau tot aan de amplitude A, waarvoor de kans op detectie berekend wordt. De factor Q kan worden benaderd door 15 Q(A)= exp[-C(1-W(A))]
Of bij benadering door een factor Q (A) = [ 1-p +p (1-W(A)/Vmeas) ]*> 20
Waarin de som over k genomen wordt van 1 tot η; n is een geheel getal dat niet groter is dan het maximaal aantal deeltjes dat fysisch past in het volume vanwaaruit meetbare reflecties optreden, maar ruimschoots groter dan het gemiddeld aantal deeltjes dat meetbare reflecties oplevert en 25 p=CVmeas/n, waarbij Vmeas het volume is van waaruit de grootste deeltjes nog echos boven het ruisniveau kunnen opwekken.
Een alternatieve uitdrukking voor de P 0 (A) factor is bijvoorbeeld
Po (A)= C jdD fb(D) dV(A,D)/ dA 30 1024984 9 *
Hierin ie C de concentratie van de deeltjes (het gemiddeld aantal deeltjes per volume-eenheid). V(A,D) is het volume waaruit een deeltje met grootte D een echo met een amplitude groter of gelijk A oplevert. Dit volume wordt in de diepte gedefinieerd door het tijdswindow vanwaaruit tijdens de meting 5 reflecties geaccepteerd worden en in de breedte door de afval van de bundelintensiteit met de afstand tot het centrum van de bundel.
In het veel voorkomende geval van een Gaussisch verdeelde bundel kan het volume V(A,D) worden uitgedrukt als 10 V(A, D)=rc Δζ log(Ao(D)/ A)/2k
Hierin is Δζ het diepteinterval waaruit reflecties geaccepteerd worden, κ is een parameter die de snelheid karakteriseert waarmee de intensiteit van de bundel als functie van de afstand tot het centrum van de bundel afvalt.
15 Ao(D) is de maximale amplitude die een reflectie van een deeltje met grootte D zou opleveren als het deeltje zich in het centrum van de bundel zou bevinden. Het zal overigens duidelijk zijn dat voor bundels met een andere dan Gauesische verdeling overeenkomstig andere uitdrukkingen voor V(A, D) gebruikt zal worden.
20 De nieuw geïntroduceerde factor Q(A), drukt de waarschijnlijkheid uit dat er naast de echo met amplitude A geen echo's met grotere amplitude optreden. Dit beschrijft het effect van maskering van reflecties door andere reflecties met een grotere of gelijke amplitude dan de amplitude A. De factor Q(A), kan eveneens worden uitgedrukt in termen van de concentratie C en 25 de deeltjesgrootteverdeling fü(D). Maatgevend hiervoor is de integraaluitdrukking W (A) = JdD fb(D) V(A,D) 30 De factor Q kan worden benaderd door 1024984 10 Q(A)= exp[ -CW(A) ]
Aan de hand van de uitdrukking P(A)= P0(A) Q voor de 5 waarschijnlijkheidsverdeling van metingen van reflecties met amplitude A die niet gemaskeerd worden door reflecties met grotere amplitude wordt uit het stelsel van reflectie amplitudes dat in een groot aantal metingen gemeten is de deeltjes grootte verdeling fb(D). Dit gebeurt bijvoorbeeld door een schatting fü(D) te zoeken die de som E minimaliseert van de kwadraten 10 tussen de voorspelde aantallen metingen Npred van amplitudes in een reeks amplitude intervallen IAi en de daadwerkelijk gemeten aantallen Nmeae E — Σΐ (Nmeas(IAi)· Npred(IAi))2 15 Deze som is in feite een indicatie voor de waarschijnlijkheid van het gemeten stelsel van amplitudes, zodat de deeltjes grootte verdeling fo(D) gekozen wordt die dit stelsel het meest waarschijnlijk maakt. Uiteraard kan in plaats van de som van de kwadraten ook met behulp van andere indicaties van de waarschijnlijkheid gemaximaliseerd worden, bijvoorbeeld 20 met behulp van een product Πι Po(A0 Q (Ai) van de waarschijnlijkheden van de verschillende metingen Ai.
25 Bij deze schatting wordt bijvoorkeur een parametrisatie van de deeltjesgrootteverdehng fb(D) gebruikt, bijvoorbeeld in de vorm van een multimodale verdeling fb(D) = li Ci exp(-(D-Di)2/2oi2 30 1 024984 ‘_ 11 (i=l, 2..M). In dit geval worden de parameters Ci Di en ai2 zo geschat dat de volgens de gebruikte indicatie de waarschijnlijkheid van de het gemeten stelsel reflectie amplitudes gemaximaliseerd wordt. Technieken voor dergelijke maximalisaties zijn in hun algemeenheid bekend, en door deze 5 toe te passen op een uitdrukking P0(A) Q (A) voor de waarschijnlijkheid van gemeten (en niet gemaskeerde) reflectieamplitudes Verkrijgt men een schatting voor de parameters.
Figuur 3 toont tevens verschillende deeltjesgrootteverdelingen die zodoende uit de aantallen reflecties in de voorbeelden van figuur 2 berekend 10 worden. Deze verdelingen vallen binnen de nauwkeurigheid van de figuur nagenoeg samen met de werkelijke verdeling. Het zal duidelijk zijn dat zodoende een betere schatting van de deeltjesgrootteverdeling verkregen wordt.
Figuur 4 toont de gevolgen van de schatting van de deeltjesgrootte 15 verdeling voor de concentratie. De concentratie is de fractie van de vloeistof die door de deeltjes wordt gevormd. Horizontaal staat de werkelijke concentratie en verticaal de uit de (gesimuleerde) metingen berekende concentratie. De getrokken lijn geeft aan waar berekende en werkelijke concentratie gelijk zijn. De cirkels geven de met de nieuwe techniek 20 berekende concentratie en de kruisjes de met de stand der techniek berekende concentraties. Uit de figuur blijkt dat met name voor hogere concentraties duidelijk betere resultaten verkregen worden.
Het zal duidelijk zijn dat de uitvinding niet gebonden is aan deze specifieke wijze van schatten van de deeltjesgrootteverdeling. Zo kan 25 bijvoorbeeld een andere vorm van parametriseren van de deeltjesgrootte verdeling gebruikt worden, bijvoorbeeld aan de hand van de waarden (en/of afgeleiden) van de deeltjesgrootteverdeling voor een aantal deeltjesgroottes en een interpolatie van de deeltjesgrootteverdeling tussen deze deeltjesgroottes, of bijvoorbeeld in de vorm van een quotiënt van twee 30 polynomen, waarvan de coëfficiënten de parameters vormen. Ook kunnen 1024984_ 12 ι benaderingen voor de verschillende termen in de wiskundige formule voor de kans P(IA) gebruikt worden.
Verder kan bijvoorbeeld eerst een schatting gemaakt worden van de factor Q(A) als functie van de amplitude Δ (bijvoorbeeld door provisorisch 5 de deeltjesgrootte verdeling te schatten en vervolgens de factor met die provisorisch geschatte verdeling uit te rekenen) en vervolgens kan men de gemeten aantallen reflecties Nmeas(A) met verschillende amplitudes corrigeren met deze factor Ncorr(A)= Nmeas(A)/Q(A). Uit de gecorrigeerde aantallen kan vervolgens op de uit de stand der techniek bekende wijze de 10 deeltjesgrootte verdeling geschat worden.
Hoewel de uitvinding beschreven is in termen van de expliciete deeltjesgrootteverdeling fb(D) zal het duidelijk zijn dat het voor bepaalde toepassingen niet nodig is om deze deeltjesgrootte verdeling ook expliciet uit te rekenen. Wanneer alleen een grootheid zoals bijvoorbeeld de 15 massadichtheid aan deeltjes nodig is, die volgt uit de integraal
JdD fo(D) m(D) dan kan deze grootheid, die een parameter van de deeltjesgrootte verdeling 20 is ook impliciet berekend worden bij het schatten van de parameters van de deeltjesgrootteverdeling die het stelsel van gemeten reflectie amplitudes het meest waarschijnlijk maakt.
De gebruikte techniek is niet beperkt tot metingen aan amplitudes, maar kan uitgebreid worden op elke vórm van metingen van eigenschappen 25 van de reflecties waarbij, bij aanwezigheid van meerdere deeltjes die reflecties veroorzaken, het effect van een dominerend deeltje het effect van de overige deeltjes maskeert.
De techniek is met name ook toepasbaar op vloeistoffen die verschillende soorten deeltjes bevatten, wanneer onderscheid gemaakt kan 30 worden van welk soort deeltje tot een reflectiemeting aanleiding geeft. Dit 1024934 13 onderscheid kan bijvoorbeeld aan het reflectiesignaal zelf zichtbaar zijn, i bijvoorbeeld door het teken van de eerste piek in een ontvangen ultrasoon reflectieeignaal in response op een ultrasone puls met een bepaald teken. In het geval van water waarin oliedruppels en zandkorrels zijn op gelost, 5 bijvoorbeeld, kunnen reflecties aan zandkorrels en oliedruppels onderscheiden worden aan het teken van de reflectie.
De gecorrigeerde waarschijnlijkheidsverdeling Pj(A)dA voor een meting van een reflectie van een deeltje van type j, met een amplitude in een bereik tussen A en A+dA kan in dit geval beschreven worden met 10
Pj(A)=Pjo(A) Πί Qi (A)
Hierin beschrijven verschillende factoren Qi (A) voor elk verschillend soort deeltjes i de waarschijnlijkheid dat geen deeltjes van het soort i tot detectie 15 leiden met grotere amplitude dan A. De factoren Q i (A) worden elk op de bovenbeschreven manier uitgedrukt in de deeltjesgrootteverdeling fbi(D) van een respectievelijke soort deeltjes i, evenals de waarschijnlijkheden Pj0(A). Met deze uitdrukkingen voor Pj(A) voor verschillende soorten deeltjes en een stelsel metingen van reflectieamplitudes die geklassificeerd zijn naar het 20 betrokken soort deeltjes kunnen op de bovenbeschreven wijze parameters van de deeltjesgrootteverdelingen fbi(D) geschat worden. Zo kunnen bijvoorbeeld deeltjesgrootteverdelingen geschat worden die een samengestelde som minimaliseren van kwadraten van afwijkingen 25 E = Xm (Nimeas(IAm)- NlpredflAm))^ van de aantallen Nimeas van gemeten reflecties van deeltjes van soort m in verschillende amplitudebereiken IAm en de voorspelde aantallen Nipred gemeten reflecties in deze bereiken.
1024984 t 14
Figuur 5 toont samenvattend de werkwijze voor het schatten van de deeltjesgrootteverdeling. In een eerst stap 51 wordt een puls opgewekt en in een tweede stap 52 wordt gedetecteerd of er in een bepaald tijdswindow na de puls een reflectie boven het ruisniveau terugkomt. In een derde stap 5 53 wordt bepaald in welk van een aantal intervallen IA de amplitude valt en wordt in een geheugen een telling voor het betrokken interval opgehoogd. Vanuit een vierde stap 54 worden de eerdere stappen herhaald totdat een bepaald aantal N reflectie experimenten is uitgevoerd. Vervolgens wordt een vijfde stap 55 uitgevoerd waarin parameters van een initiële 10 deeltjesgrootteverdeling gekozen wordt. Hiermee wordt in een zesde stap 56 met behulp van de formule P(A)=P0(A) Q (A) een voorspelling van de gemeten aantallen uitgerekend en wordt een verschil met de gemeten aantallen bepaald. Als het verschil voldoende klein is stopt de werkwijze in een zevende stap 57. Anders worden in een achtste stap 58 de parameters 15 aangepast, zo dat het verschil naar verwachting zal verminderen en wordt de werkwijze vanaf de zesde stap 56 herhaald. Uiteraard wordt schatting in de werkwijze bijvoorkeur met een computer uitgevoerd.
Het zal duidelijk zijn dat de uitvinding niet is beperkt tot de beschreven werkwijze, zo hoeft bijvoorbeeld in de zesde stap 56 de het verschil niet 20 expliciet uitgerekend te worden, maar kan in plaats daarvan informatie uitgerekend worden die nodig is om de aanpassing van de parameters te kiezen en kan de werkwijze stoppen in de zevende stap 57 als de aanpassing kleiner dan een drempel is.
1024984
Claims (8)
1. Werkwijze voor het meten van parameters die een deeltjesgrootteverdeling van deeltjes in een vloeistof beschrijven, welke werkwijze de stappen omvat van - het uitvoeren van een reeks reflectiemetingen, in elk waarvan een 5 signaalbundel in de vloeistof wordt opgewekt en een waarde (A) van een eigenschap van een reflectie aan een deeltje in de signaalbundel wordt gemeten; - het uitvoeren van een maximum likelihood schatting van de parameters gezien een combinatie van de gemeten waardes (A), op basis van een 10 uitdrukking voor een waarschijnlijkheid van de gemeten waardes als functie van de gemeten waardes, welke uitdrukking een eerste factor (P) bevat voor de waarschijnlijkheid van een reflectiemeting waarvan een reflectie met de gemeten waarde onderdeel van uitmaakt, gecorrigeerd met een tweede factor (Q) voor de waarschijnlijkheid dat niet tevens een reflectie met een 15 dominerende waarde van de eigenschap, die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting.
2. Werkwijze volgens conclusie 1, waarin de eerste factor (P) de deeltjesgrootte verdeling omvat, versmeerd met een kansverdeling dat een deeltje met een bepaalde grootte leidt tot een reflectiemeting waarvan een 20 reflectie met de gemeten waarde onderdeel van uitmaakt.
3. Werkwijze volgens conclusie 1 of 2, waarin de tweede factor (Q) de waarschijnlijkheid dat een reflectie met een andere dan de gemeten waarde onderdeel van uitmaakt van een reflectiemeting, geïntegreerd over een bereik aan andere dan de gemeten waardes.
4. Werkwijze volgens conclusie 3, waarin de tweede factor nagenoeg overeenkomt met exp(- C JdA'JdD fb(D) G(A' | D)), waarin D een deeltjesgrootte is, C een concentratie van de deeltjes, fb(D) een dichtheid 1024984____ » van deeltjes 'met deeltjesgrootte D en G(A | D) is de conditionele waarschijnlijkheid dat er een reflectie met amplitude A detecteerbaar is, als een deeltje met een grootte D een detecteerbare reflectie oplevert.
5. Werkwijze volgens één der voorafgaande conclusies, waarin de 5 maximum likelihood schatting het maken van tellingen omvat van aantallen reflectiemetingen waarin de gemeten waarde in respectieve waardeintervallen vallen, en de schatting zo gekozen wordt dat een samenstel van afwijkingen tussen de tellingen in de verschillende intervallen volgens de waarschijnlijkheid als functie van de gemeten 10 waardes voorspelde tellingen geminimaliseerd wordt.
6. Werkwijze volgens één der voorafgaande conclusies, waarin bij de reflectiemetingen onderscheid gemaakt wordt tussen verschillende soorten deeltjes die de reflecties veroorzaken en bij het uitvoeren van de maximum i likelihood schatting de uitdrukking gecorrigeerd wordt met een product van 15 respectieve tweede factoren voor de waarschijnlijkheid dat niet tevens een reflectie met een dominerende waarde van de eigenschap door respectievelijke soorten deeltjes, die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting.
7. Inrichting voor het meten van parameters die een 20 deeltjesgrootteverdeling van deeltjes in een vloeistof beschrijven, welke inrichting omvat - een vloeistofkanaal; - middelen om een ultrasone bundel op te wekken in het vloeistof kanaal; - middelen om een eigenschap van een reflectie van de bundel aan een 25 deeltje in het vloeistofkanaal te meten; - een verwerkingseenheid ingericht voor het uitvoeren van een maximum likelihood schatting van de parameters, gezien een combinatie van de gemeten waardes, op basis van een uitdrukking voor een waarschijnlijkheid van de gemeten waardes als functie van de gemeten waardes, welke 30 uitdrukking een eerste factor bevat voor de waarschijnlijkheid van een 1024984 reflectiemeting waarvan een reflectie met de gemeten waarde onderdeel van uitmaakt, gecorrigeerd met een tweede factor voor de waarschijnlijkheid dat niet tevens een reflectie met een dominerende waarde van de eigenschap, die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting.
8. Computer programma product met instructies voor het meten van parameters die een deeltjesgrootteverdeling van deeltjes in een vloeistof beschrijven, aan de hand van van een reeks reflectiemetingen, in elk waarvan een signaalbundel in de vloeistof wordt opgewekt en een waarde (A) van een eigenschap van een reflectie aan een deeltje in de signaalbundel 10 wordt gemeten; welke instructies zijn ingericht voor het uitvoeren van een maximum likelihood schatting van de parameters gezien een combinatie van de gemeten waardes (A), op basis van een uitdrukking voor een waarschijnlijkheid van de gemeten waardes als functie van de gemeten waardes, welke uitdrukking een eerste factor (P) bevat voor de 15 waarschijnlijkheid van een reflectiemeting waarvan een reflectie met de gemeten waarde onderdeel van uitmaakt, gecorrigeerd met een tweede factor (Q) voor de waarschijnlijkheid dat niet tevens een reflectie met een dominerende waarde van de eigenschap, die de gemeten waarde zou maskeren, deel uitmaakt van de reflectiemeting. 20 1024934
Priority Applications (10)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NL1024984A NL1024984C2 (nl) | 2003-12-10 | 2003-12-10 | Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. |
PCT/NL2004/000857 WO2005057183A1 (en) | 2003-12-10 | 2004-12-09 | Method and apparatus for obtaining information about the size distribution of macroscopic particles in a liquid |
JP2006543751A JP4823918B2 (ja) | 2003-12-10 | 2004-12-09 | 液体中の肉眼で見える粒子の粒度分布に関する情報を得るための方法および装置 |
CA2548984A CA2548984C (en) | 2003-12-10 | 2004-12-09 | Method and apparatus for obtaining information about the size distribution of macroscopic particles in a liquid |
CN2004800415502A CN1914499B (zh) | 2003-12-10 | 2004-12-09 | 用于获得关于液体中肉眼可见的颗粒大小分布的信息的方法和装置 |
EP04808773.8A EP1697721B1 (en) | 2003-12-10 | 2004-12-09 | Method and apparatus for obtaining information about the size distribution of macroscopic particles in a liquid |
RU2006124556/28A RU2376581C2 (ru) | 2003-12-10 | 2004-12-09 | Способ и устройство для получения информации о распределении макроскопических частиц в жидкости по размеру |
AU2004297478A AU2004297478A1 (en) | 2003-12-10 | 2004-12-09 | Method and apparatus for obtaining information about the size distribution of macroscopic particles in a liquid |
US10/582,282 US7610822B2 (en) | 2003-12-10 | 2004-12-09 | Method and apparatus for obtaining information about the size distribution of macroscopic particles in a liquid |
NO20062920A NO338984B1 (no) | 2003-12-10 | 2006-06-22 | Fremgangsmåte og anordning for å oppnå informasjon om størrelsesdistribusjon av makroskopiske partikler i en væske |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NL1024984A NL1024984C2 (nl) | 2003-12-10 | 2003-12-10 | Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. |
NL1024984 | 2003-12-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
NL1024984C2 true NL1024984C2 (nl) | 2005-06-13 |
Family
ID=34676034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NL1024984A NL1024984C2 (nl) | 2003-12-10 | 2003-12-10 | Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. |
Country Status (10)
Country | Link |
---|---|
US (1) | US7610822B2 (nl) |
EP (1) | EP1697721B1 (nl) |
JP (1) | JP4823918B2 (nl) |
CN (1) | CN1914499B (nl) |
AU (1) | AU2004297478A1 (nl) |
CA (1) | CA2548984C (nl) |
NL (1) | NL1024984C2 (nl) |
NO (1) | NO338984B1 (nl) |
RU (1) | RU2376581C2 (nl) |
WO (1) | WO2005057183A1 (nl) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070168354A1 (en) * | 2005-11-01 | 2007-07-19 | Jorey Ramer | Combined algorithmic and editorial-reviewed mobile content search results |
US9714893B2 (en) | 2009-11-13 | 2017-07-25 | Stable Solutions Llc | Particle size distribution profiles and use thereof to adjust a dispersion characteristic |
DE102010031129A1 (de) * | 2010-07-08 | 2012-01-12 | Endress + Hauser Flowtec Ag | Ultraschall-Partikelmesssystem |
DE102010031128A1 (de) * | 2010-07-08 | 2012-01-12 | Endress + Hauser Flowtec Ag | Ultraschall-Partikelmesssystem |
JP5915758B2 (ja) * | 2012-09-12 | 2016-05-11 | 株式会社島津製作所 | 粒度分布データ処理装置及びこれを備えた粒度分布測定装置、並びに、粒度分布データ処理方法及び粒度分布データ処理プログラム |
DE102013213362A1 (de) * | 2013-07-08 | 2015-01-08 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren zur Identifizierung und Quantifizierung von emittierenden Teilchen in Systemen |
CN109983334B (zh) * | 2016-06-14 | 2022-07-19 | 尼玛克股份有限公司 | 定量测定非均质材料混合物的组分的浓度或粒径的方法、设备及该设备的用途 |
RU2694448C1 (ru) * | 2018-12-13 | 2019-07-15 | Общество с ограниченной ответственностью "Нониус Инжиниринг" (ООО "Нониус Инжиниринг") | Способ измерения плотности пульпы |
US11204317B2 (en) * | 2019-02-26 | 2021-12-21 | Mitsubishi Electric Research Laboratories, Inc. | Tomographic imaging system |
CN112525785B (zh) * | 2020-11-25 | 2023-03-24 | 长江水利委员会长江科学院 | 一种基于新型多频超声探头的悬移质含沙量测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3779070A (en) * | 1971-11-23 | 1973-12-18 | Autometrics Co | Particle size and percent solids monitor |
EP0989397A1 (en) * | 1998-08-31 | 2000-03-29 | Lucent Technologies Inc. | Process for determining characteristics of suspended particles |
WO2003102550A1 (en) * | 2002-05-30 | 2003-12-11 | Nederlandse Organisatie Voor Toegepast- Natuurwetenschappelijk Onderzoek Tno | Ultrasonic characterization of high concentration particles |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5121629A (en) * | 1989-11-13 | 1992-06-16 | E. I. Du Pont De Nemours And Company | Method and apparatus for determining particle size distribution and concentration in a suspension using ultrasonics |
NL1003595C2 (nl) * | 1996-04-10 | 1997-10-14 | Tno | Werkwijze en inrichting voor het karakteriseren van suspensies. |
KR20040020869A (ko) * | 2000-12-18 | 2004-03-09 | 이 아이 듀폰 디 네모아 앤드 캄파니 | 서스펜션 입자의 초음파 사이징을 위한 방법 및 장치 |
CN1208613C (zh) * | 2002-11-14 | 2005-06-29 | 上海交通大学 | 一种可同时测定颗粒粒径和浓度的光学传感器 |
-
2003
- 2003-12-10 NL NL1024984A patent/NL1024984C2/nl not_active IP Right Cessation
-
2004
- 2004-12-09 RU RU2006124556/28A patent/RU2376581C2/ru active
- 2004-12-09 US US10/582,282 patent/US7610822B2/en active Active
- 2004-12-09 CN CN2004800415502A patent/CN1914499B/zh active Active
- 2004-12-09 CA CA2548984A patent/CA2548984C/en active Active
- 2004-12-09 WO PCT/NL2004/000857 patent/WO2005057183A1/en active Application Filing
- 2004-12-09 AU AU2004297478A patent/AU2004297478A1/en not_active Abandoned
- 2004-12-09 JP JP2006543751A patent/JP4823918B2/ja active Active
- 2004-12-09 EP EP04808773.8A patent/EP1697721B1/en active Active
-
2006
- 2006-06-22 NO NO20062920A patent/NO338984B1/no unknown
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3779070A (en) * | 1971-11-23 | 1973-12-18 | Autometrics Co | Particle size and percent solids monitor |
EP0989397A1 (en) * | 1998-08-31 | 2000-03-29 | Lucent Technologies Inc. | Process for determining characteristics of suspended particles |
WO2003102550A1 (en) * | 2002-05-30 | 2003-12-11 | Nederlandse Organisatie Voor Toegepast- Natuurwetenschappelijk Onderzoek Tno | Ultrasonic characterization of high concentration particles |
Also Published As
Publication number | Publication date |
---|---|
US7610822B2 (en) | 2009-11-03 |
US20070091301A1 (en) | 2007-04-26 |
RU2376581C2 (ru) | 2009-12-20 |
NO338984B1 (no) | 2016-11-07 |
CA2548984C (en) | 2013-10-22 |
EP1697721B1 (en) | 2016-03-16 |
WO2005057183A1 (en) | 2005-06-23 |
EP1697721A1 (en) | 2006-09-06 |
CA2548984A1 (en) | 2005-06-23 |
RU2006124556A (ru) | 2008-01-20 |
JP4823918B2 (ja) | 2011-11-24 |
JP2007514937A (ja) | 2007-06-07 |
CN1914499B (zh) | 2011-01-19 |
CN1914499A (zh) | 2007-02-14 |
NO20062920L (no) | 2006-09-11 |
AU2004297478A1 (en) | 2005-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NL1003595C2 (nl) | Werkwijze en inrichting voor het karakteriseren van suspensies. | |
NL1024984C2 (nl) | Werkwijze en inrichting voor het verkrijgen van informatie over de grootteverdeling van macroscopische deeltjes in een vloeistof. | |
David et al. | Modern spectral analysis techniques for blood flow velocity and spectral measurements with pulsed Doppler ultrasound | |
JP3207426B2 (ja) | パルスエコー方式で動作する間隔センサの受信信号における有効信号と障害信号の識別および分離方法 | |
US7894874B2 (en) | Method and apparatus for enhancing the detecting and tracking of moving objects using ultrasound | |
KR102156396B1 (ko) | 초음파반사율 및 수심을 이용한 다중회귀분석을 통해 부유사농도를 추정하는 방법 | |
Larrarte et al. | Acoustic Doppler flow-meters: a proposal to characterize their technical parameters | |
CN112327293A (zh) | 一种基于分形特征智能学习的海面目标检测方法 | |
CN110927735A (zh) | 基于多通道全波形激光雷达数据的多目标距离测量方法 | |
NL1020714C2 (nl) | Ultrasone karakterisatie van hoge concentratie deeltjes. | |
Pintado-Patiño et al. | Hydrodynamics and sediment transport under a dam-break-driven swash: An experimental study | |
Park et al. | A standard criterion for measuring turbulence quantities using the four-receiver acoustic Doppler velocimetry | |
Birjandi et al. | Improvement of Acoustic Doppler Velocimetry in bubbly flow measurements as applied to river characterization for kinetic turbines | |
CN109613505A (zh) | 一种抑制双重杂波的装置及其抑制方法 | |
Jones et al. | Measures of nonlinearity for single target tracking problems | |
Słota | Decomposition techniques for full-waveform airborne laser scanning data | |
RU2351947C2 (ru) | Способ измерения начальной скорости снаряда | |
CN114660558A (zh) | 基于改进keystone的等离子鞘套包覆目标相参积累方法、***、设备及存储介质 | |
Zedel | Noise in coherent Doppler sonar velocity measurements: Where does it come from and what can you do about it? | |
Wang et al. | Robust Multi-Target Tracking Algorithm Based on Automotive Millimeter-Wave Radar | |
Vogt et al. | Measurement and multivariate analysis of ultrasound parameters for quantitative characterization of liquids | |
Hurther et al. | Multistatic acoustic Doppler profilers for fine-scale studies of velocity and particle flux processes | |
Wu et al. | An Eulerian segment-matching method for bloodflow velocity estimation using pulsed ultrasound RF echoes | |
Vanovsky | Doppler radar: Retrieval of information about turbulence in rain | |
Fisher et al. | Spectral reconstruction method for liquid velocity measurement beyond the Nyquist limit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PD2B | A search report has been drawn up | ||
MM | Lapsed because of non-payment of the annual fee |
Effective date: 20180101 |