ES2952720A1 - Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo - Google Patents

Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo Download PDF

Info

Publication number
ES2952720A1
ES2952720A1 ES202330342A ES202330342A ES2952720A1 ES 2952720 A1 ES2952720 A1 ES 2952720A1 ES 202330342 A ES202330342 A ES 202330342A ES 202330342 A ES202330342 A ES 202330342A ES 2952720 A1 ES2952720 A1 ES 2952720A1
Authority
ES
Spain
Prior art keywords
medium
size distribution
refractive index
particle size
procedure
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
Application number
ES202330342A
Other languages
English (en)
Other versions
ES2952720B2 (es
Inventor
Chorro Miguel Angel Báez
Rodriguez Borja Vidal
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Universidad Politecnica de Valencia
Original Assignee
Universidad Politecnica de Valencia
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Universidad Politecnica de Valencia filed Critical Universidad Politecnica de Valencia
Priority to ES202330342A priority Critical patent/ES2952720B2/es
Publication of ES2952720A1 publication Critical patent/ES2952720A1/es
Priority to PCT/ES2024/070254 priority patent/WO2024223970A1/es
Application granted granted Critical
Publication of ES2952720B2 publication Critical patent/ES2952720B2/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02Investigating particle size or size distribution
    • G01N15/0205Investigating particle size or size distribution by optical means

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 Analysing Materials By Optical Means (AREA)

Abstract

Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo. Procedimiento y sistema para estimar la distribución del tamaño de las partículas combinando la información sobre módulo y fase del campo electromagnético esparcido. La detección de la señal esparcida hacia delante se realiza con medios capaces de obtener la información del módulo y la fase del campo eléctrico. A partir de estas medidas es posible derivar directamente el coeficiente de extinción y el índice de refracción. La distribución de tamaño se calcula entonces mediante un método iterativo de Twomey, donde el factor para cada iteración está modificado de manera que incluye información sobre el coeficiente de extinción y el índice de refracción.

Description

DESCRIPCIÓN
Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo
SECTOR DE LA TÉCNICA
La invención pertenece al sector del estudio de sistemas de partículas y en especial en agregados, aerosoles, emulsiones, polvos y suspensiones.
Más en particular, la presente invención se refiere a un procedimiento para la estimación de la distribución del tamaño de múltiples partículas a partir de la detección coherente de la luz esparcida hacia adelante.
ESTADO DE LA TÉCNICA
Existe la necesidad de comprender la naturaleza de los sistemas de partículas en numerosos campos, desde la farmacología a la ciencia de materiales. Con diferencia, la característica más importante de las muestras de partículas es la distribución del tamaño. Esta característica determina la calidad, la seguridad y el rendimiento del producto final.
El proceso de estimación de la distribución del tamaño de partículas (DTP) generalmente se organiza en los siguientes pasos:
(i) Preparación: procesamiento de muestras, configuración de instrumentos, calibración y alineamiento, instalación del sistema de entrega de muestras.
(ii) Medición: iluminar la muestra con una señal y medir la intensidad de la luz esparcida hacia adelante por la muestra para compararla con la intensidad sin muestra y de ahí derivar el coeficiente de extinción de la muestra.
(iii) Análisis directo: determinar el efecto del esparcimiento de la luz provocado por las partículas de la muestra bajo prueba sobre el coeficiente de extinción, asumiendo parámetros específicos para las partículas (forma, índice de refracción, concentración, etc.).
(iv) Análisis inverso: determinar las características desconocidas de la muestra utilizando los resultados del análisis directo y los datos de medición correspondientes. Un paso esencial de esta etapa es el análisis del contenido de información de las ecuaciones a invertir. Un análisis inverso generalmente se complementa con información a priori sobre la solución de interés.
(v) Estadísticas y análisis de datos: comprender los datos de cálculo, obtener la distribución del tamaño de las partículas utilizando una variedad de modelos (número/volumen/masa).
Para la modelización del coeficiente de extinción t ( v ) de la muestra en el paso de análisis directo (ii) es común utilizar la siguiente expresión:
Figure imgf000003_0001
0
Donde N0 es la cantidad total de partículas por unidad de volumen, D es el diámetro de las partículas, v es la frecuencia del campo electromagnético empleado en la medida, f (D) es la DTP normalizada del medio y Qext(v,D) es el factor de eficiencia de extinción para una sola partícula esférica, dado por la teoría de Mie.
Para obtener la DTP es necesario invertir la ecuación anterior que, tras la adecuada discretización de la integral, puede expresarse en forma de sistema de ecuaciones lineales:
Figure imgf000003_0002
Sin embargo, la obtención de la solución f mediante la inversión de la matriz K no es posible, ya que estos sistemas suelen presentar un mal condicionamiento, es decir, pequeños cambios en los valores de sus elementos generan soluciones radicalmente diferentes. Para la adecuada solución de este tipo de problema es necesario aplicar técnicas de inversión regularizada.
Un método conocido para la inversión de la distribución de tamaño de partículas a partir de medidas de extinción, (paso iv) es el método de Twomey. Éste corresponde a la categoría de métodos iterativos que obtienen una estimación de la solución tras la repetida multiplicación de la solución previa por un cierto factor partiendo desde una solución de prueba inicial. En concreto, la solución iterativa de la distribución para un diámetro Dj según el método de Twomey tiene la siguiente forma:
Figure imgf000004_0001
El factor multiplicativo en cada iteración depende, por un lado, del cociente entre la extinción medida (rmeas) y la calculada (t^ ) con el modelo de Mie (ecuación 1) empleando la solución de la iteración previa (fp). Por otro lado, este cociente está pesado por una función (Wtj ) calculada previamente a partir del modelo de Mie (ecuación 1) y cuya función es regular la contribución de cada tamaño de partícula proporcionalmente a su contribución en el valor de los parámetros ópticos. El factor iterativo finalmente se obtiene como el producto de las cantidades anteriores para cada frecuencia, v¿.
Una ventaja de este tipo de métodos iterativos es que aseguran que la solución tenga un valor positivo. Además, permiten trabajar con cantidades experimentales cuya descripción matemática no esté relacionada linealmente con la distribución de tamaños de partícula.
Puesto que conocer la distribución del tamaño de partículas es crítico en muchas aplicaciones, es deseable disponer de un método que consiga una mayor precisión en el cálculo de dicha distribución.
RESUMEN DE LA INVENCIÓN
La presente invención parte de la superación del prejuicio técnico de que la estimación de la distribución de tamaños de partícula de un medio heterogéneo a partir de la información sobre intensidad de la señal dispersada ofrece los mismos resultados que la estimación a partir de la información de fase y que la agregación de estas dos es redundante. Por el contrario, se ha encontrado que una modificación del método de Twomey para la resolución del problema inverso a partir de la luz esparcida hacia delante permite incorporar tanto la información de fase como la del módulo del campo para obtener la distribución de tamaños de partícula de forma más precisa.
La invención proporciona un procedimiento y sistema para estimar la distribución del tamaño de las partículas combinando pues la información del módulo y fase del campo electromagnético esparcido. Para ello, la detección de la señal esparcida hacia delante se debe realizar con medios capaces de obtener la información del módulo y fase del campo eléctrico. A partir de estas medidas es posible derivar directamente el coeficiente de extinción y el índice de refracción.
BREVE DESCRIPCIÓN DE LAS FIGURAS
Con objeto de ayudar a una mejor comprensión de las características de la invención y para complementar esta descripción, se acompañan como parte integrante de la misma las siguientes figuras, cuyo carácter es ilustrativo y no limitativo:
La Figura 1 muestra los elementos esenciales de la invención.
Las Figuras 2a-2d muestran resultados experimentales comparando la distribución obtenida mediante el procedimiento del estado de la técnica y aquella obtenida mediante el procedimiento de la invención.
DESCRIPCIÓN DETALLADA
El procedimiento para la estimación de la distribución del tamaño de múltiples partículas de la invención comprende los siguientes pasos:
- Iluminar la muestra con una fuente de radiación que incluya distintas frecuencias (Fig. 1) ;
- detectar una porción de radiación esparcida hacia adelante mediante un detector capaz de obtener el campo óptico complejo, es decir, capaz de extraer la información tanto del módulo como de la fase del campo electromagnético;
- medir tanto el efecto del medio heterogéneo en el campo complejo (campo muestra) como el campo en ausencia de dicho medio (campo referencia), es decir, medir el campo electromagnético de la luz esparcida hacia adelante por la muestra para compararla con el campo sin muestra
- calcular la transmisión, T, a través de la muestra como el valor absoluto del cociente de los módulos de los campos muestra y referencia. Calcular la diferencia de fase A0 como la resta de las fases de los campos muestra y referencia;
- derivar los parámetros ópticos del medio heterogéneo bajo prueba ( n exp, índice de refracción y rexp, extinción) a partir de la información experimental de transmisión y de la diferencia de fase calculadas anteriormente mediante las siguientes fórmulas:
Figure imgf000006_0003
donde d es el espesor de la muestra, v la frecuencia y c la velocidad de la luz en el vacío;
modelizar los parámetros ópticos de la muestra heterogénea, para lo que puede emplearse la teoría de Waterman-Truell para medios con grandes cantidades de partículas, que tiene en cuenta el esparcimiento múltiple de la luz, cuyo efecto es notable a partir de pequeñas fracciones volumétrica de partículas (>1%). Esta teoría permite calcular el índice de refracción y la extinción con la siguiente fórmula:
Figure imgf000006_0004
donde Re(-) y Im(-) hacen referencia a tomar la parte real e imaginaria del número complejo en el paréntesis. (S(0/n)) son las amplitudes de Mie para ángulos 0 y n para una sola partícula. kmed — - n med es la constante de propagación del medio matriz (es decir, sin tener en cuenta el efecto de las partículas agregadas) y nmed el índice de refracción de dicho medio. El símbolo {) indica un promedio de las amplitudes S pesadas con la DTP del medio;
obtener la distribución del tamaño de las partículas a partir de los parámetros ópticos mediante un nuevo algoritmo basado en la modificación del método de Twomey para incorporar al factor la contribución de las medidas de índice de refracción además de las del coeficiente de extinción:
Figure imgf000006_0001
donde nexp es el índice de refracción medido experimentalmente
Figure imgf000006_0002
es el índice de refracción calculado usando la ecuación 5 con la DTP de la iteración p. Wij son los pesos calculados mediante la ecuación. 3 y Wi] es el análogo empleando el índice de refracción
W n _ ncalcjv i>Dj ) ™ ij l t n calc(vi,DJ) ( ) - el parámetro a representa un equilibrio entre la importancia que tiene la extinción o el índice de refracción en la inversión. Este parámetro se elegirá como aquel que produzca una distribución de tamaños cuyos parámetros ópticos calculados se asemejan más a los experimentales. Para ello se elige a como aquel que produce un mínimo de la métrica ERM, donde ERM es
Figure imgf000007_0001
- esta métrica se basa en determinar el error relativo medio (ERM) que es la suma de la media de los residuos relativos de la extinción y el índice de refracción. Se ha determinado que los valores mínimos de esta métrica corresponden con valores de a que se encuentran próximos a los valores óptimos de forma que puede considerarse un buen indicador para determinar a sin más información que las medidas experimentales. En caso de que el valor de ERM no presente un mínimo muy pronunciado, es posible promediar las soluciones en un rango de valores de a alrededor del mínimo. Se ha observado que un rango 0.25 proporciona una estimación robusta de la verdadera DTP en la mayoría de los casos.
Resultados experimentales:
La técnica propuesta para la obtención del tamaño de partículas ha sido validada experimentalmente. En las figuras 2a y 2b se muestran parámetros ópticos de muestras granuladas hechas a partir de polvo de PTFE y microesferas de vidrio. La fracción de volumen acumulada de las esferas es de 5%. La figura 2c muestra los resultados de inversión de la DTP utilizando el método de Twomey original (línea punteada) y utilizando el método propuesto (línea continua) junto a la DTP dada por el fabricante. En la figura 2d se muestra el valor de la métrica utilizada para determinar el valor del parámetro a.
En las Figuras 2a y 2b se pueden observar los parámetros ópticos medidos (negro) así como los calculados a partir del modelo de esparcimiento de la luz (Tcatc, ncalc) usando como entrada la DTP invertida. Los resultados de la inversión empleando tanto el método de Twomey tradicional que sólo tiene en cuenta la extinción (línea punteada) como el método propuesto aunando información de extinción e índice de refracción (línea continua) se muestran en la inferior junto a la DTP dada por el fabricante. En la Fig. 2d se muestran los valores de ERM obtenidos con varios parámetros a y la posición del mínimo empleado en las inversiones. El método propuesto consigue disminuir la desviación respecto a la DTP tabulada consiguiendo un error menor que el método de Twomey. Concretamente, la mejora resulta de un 32.6%.
Esta validación experimental se ha efectuado en el infrarrojo lejano empleando un instrumento de Espectroscopía de Terahercios en Dominio Temporal que permite obtener el coeficiente de extinción y el índice de refracción de muestras en la banda de 0.2 a 1.6 THz. Un experto entenderá que sistemas que empleen antenas y detectores de RF o microondas, interferómetros espectrales de baja coherencia e interferómetros con fuentes coherentes sintonizables también son capaces de proporcionar la información sobre módulo y fase del campo electromagnético esparcido en sus correspondientes rangos espectrales de operación y por lo tanto son también aptos para su uso en el presente procedimiento.
A la vista de esta descripción y figuras, el experto en la materia podrá entender que la invención ha sido descrita según algunas realizaciones preferentes de la misma, pero que múltiples variaciones pueden ser introducidas en dichas realizaciones preferentes, sin exceder el objeto de la invención tal y como ha sido reivindicada.

Claims (5)

REIVINDICACIONES
1. Procedimiento para el cálculo de la estimación de la distribución del tamaño de partículas en medios heterogéneos, donde el procedimiento comprende los pasos de: - a) iluminar el medio con una fuente de radiación a distintas frecuencias;
- b) medir la luz esparcida hacia adelante por el medio mediante un detector capaz de obtener el campo óptico complejo para compararla con la señal sin medio;
- c) calcular la transmisión, T, a través del medio como el valor absoluto del cociente de los módulos de los campos con medio y sin medio y calcular la diferencia de fase A<p como la resta de las fases de dichos campos muestra y referencia
- d) derivar el índice de refracción nexp y el coeficiente de extinción rexp, del medio heterogéneo a partir de la información experimental de transmisión y de la diferencia de fase calculadas anteriormente mediante las siguientes fórmulas:
Figure imgf000009_0001
donde d es el espesor de la muestra, v la frecuencia y c la velocidad de la luz en el vacío;
- e) modelizar los parámetros ópticos del medio heterogéneo empleando la teoría de Waterman-Truell para medios con grandes cantidades de partículas para obtener los valores del índice de refracción ncaic y el coeficiente de extinción r caíc;
- f) obtener la distribución del tamaño de las partículas mediante el método iterativo de Twomey, donde el factor para cada iteración está modificado según la siguiente expresión:
Figure imgf000009_0002
donde nexp es el índice de refracción medido experimentalmente y n f y es el índice de refracción calculado con la DTP de la iteración p.W-j son los pesos calculados mediante la ecuación. 3 y W-j es el análogo empleando el índice de refracción
Figure imgf000009_0003
y donde el parámetro a se elige como aquel que produce un mínimo de la métrica ERM, donde ERM es
Figure imgf000010_0001
2. Sistema para llevar a cabo el procedimiento de la reivindicación 1 que comprende: una fuente de radiación con distintas frecuencias, un detector capaz de obtener el campo óptico complejo y un procesador provisto de medios de programa para realizar los pasos c-f.
3. Sistema según la reivindicación 2 caracterizado por que el detector es un espectroscopio de terahercios en dominio temporal, un detector de RF o microondas, un interferómetro espectral de baja coherencia o un interferómetro con fuentes coherentes sintonizables.
4. Producto de programa de ordenador que comprende instrucciones para llevar a cabo los pasos c-f de la reivindicación 1.
5. Medio legible por ordenador que comprende almacenado el programa capaz de llevar a cabo los pasos c-f de la reivindicación 1.
ES202330342A 2023-04-28 2023-04-28 Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo Active ES2952720B2 (es)

Priority Applications (2)

Application Number Priority Date Filing Date Title
ES202330342A ES2952720B2 (es) 2023-04-28 2023-04-28 Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo
PCT/ES2024/070254 WO2024223970A1 (es) 2023-04-28 2024-04-23 Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
ES202330342A ES2952720B2 (es) 2023-04-28 2023-04-28 Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo

Publications (2)

Publication Number Publication Date
ES2952720A1 true ES2952720A1 (es) 2023-11-03
ES2952720B2 ES2952720B2 (es) 2024-11-11

Family

ID=88558708

Family Applications (1)

Application Number Title Priority Date Filing Date
ES202330342A Active ES2952720B2 (es) 2023-04-28 2023-04-28 Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo

Country Status (2)

Country Link
ES (1) ES2952720B2 (es)
WO (1) WO2024223970A1 (es)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0464337A2 (en) * 1990-07-02 1992-01-08 International Business Machines Corporation Measurement of size and refractive index of particles using the complex forward-scattered electromagnetic field

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0464337A2 (en) * 1990-07-02 1992-01-08 International Business Machines Corporation Measurement of size and refractive index of particles using the complex forward-scattered electromagnetic field

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SAWYER N B E ET AL. Amplitude and phase microscopy for sizing of spherical particles. Applied Optics 1 Aug. 2003 Opt. Soc. America USA. , 01/08/2003, Vol. 42, Páginas 4488 - 4498 ISSN 0003-6935 (print), *
SWIRNIAK GRZEGORZ ET AL. Forward and inverse analysis for particle size distribution measurements of disperse samples: A review. MEASUREMENT, 20211019 INSTITUTE OF MEASUREMENT AND CONTROL. LONDON, GB. , 19/10/2021, Vol. 187, ISSN 0263-2241, (DOI: doi:10.1016/j.measurement.2021.110256) Apartados 1, 2.3, 4.3.2 *
TAUBENBLATT M A ET AL. Measurement of the size and refractive index of a small particle using the complex forward-scattered electromagnetic field. Applied Optics 20 Nov. 1991 USA. , 20/11/1991, Vol. 30, Páginas 4972 - 4979 ISSN 0003-6935 (print), *
ZHOU J ET AL. Particle size and refractive index retrieval from the backscattering spectrum of white light using the Twomey iterative method: Simulation and experiment. Applied Optics 20060910 OSA - The Optical Society usa. , 10/09/2006, Vol. 45, Páginas 6876 - 6885 ISSN 1559-128X (print) ISSN 1539-4522 (electronic), (DOI: doi:10.1364/AO.45.006876) Apartado 2 *

Also Published As

Publication number Publication date
ES2952720B2 (es) 2024-11-11
WO2024223970A1 (es) 2024-10-31

Similar Documents

Publication Publication Date Title
Boruah et al. Cosmic flows in the nearby Universe: new peculiar velocities from SNe and cosmological constraints
Wax et al. Cellular organization and substructure measured using angle-resolved low-coherence interferometry
Basse et al. Dark energy properties from large future galaxy surveys
Hill Review of optical scintillation methods of measuring the refractive-index spectrum, inner scale and surface fluxes
Moskalensky et al. Accurate measurement of volume and shape of resting and activated blood platelets from light scattering
Vega et al. Latex particle size distribution by dynamic light scattering: novel data processing for multiangle measurements
Riviere et al. Multispectral polarized BRDF: design of a highly resolved reflectometer and development of a data inversion method.
Gutiérrez-Reyes et al. Extension of Fresnel’s formulas for turbid colloidal suspensions: a rigorous treatment
CN110887800A (zh) 一种用于光谱法水质在线监测系统的数据校准方法
Kuz’min et al. Photon diffusion in random media and anisotropy of scattering in the Henyey–Greenstein and Rayleigh–Gans models
Gao et al. GMOPNet: A GAN-MLP two-stage network for optical properties measurement of kiwifruit and peaches with spatial frequency domain imaging
ES2952720B2 (es) Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejo
Ma et al. The measurement and modeling investigation on the spectral polarized BRDF of brass
Chen et al. Optical aerosol sizing method without prior refractive index
Stergar et al. Effects of phantom microstructure on their optical properties
CN106290094B (zh) 应用于扬尘颗粒物在线监测的mie散射快速计算方法
Sviridenkov et al. Parameterization of asymmetry factor of atmospheric aerosol scattering phase function
Ustinov On inverse reconstruction problems of the erythrocyte size distribution in laser diffractometry
Buccini et al. Weakly constrained Lucy–Richardson with applications to inversion of light scattering data
CN108645817B (zh) 一种多类型混合颗粒物质量浓度在线测量方法
Goulden Spectroturbidimetry of emulsions
Mroczka Method of moments in light scattering data inversion in the particle size distribution function
JP3645760B2 (ja) 粒子径分布解析方法
Dunn Angularly-Resolved Light Scattering for Sizing Organelle Distributions in Single Cells
RU206033U1 (ru) Устройство для определения количества частиц и распределения их по скоростям в жидких биологических средах

Legal Events

Date Code Title Description
BA2A Patent application published

Ref document number: 2952720

Country of ref document: ES

Kind code of ref document: A1

Effective date: 20231103

FG2A Definitive protection

Ref document number: 2952720

Country of ref document: ES

Kind code of ref document: B2

Effective date: 20241111