ES2365653T3 - Procedimiento para determinar una dosis de radiación y procedimiento para determinar una curva de isodosis asociada. - Google Patents

Procedimiento para determinar una dosis de radiación y procedimiento para determinar una curva de isodosis asociada. Download PDF

Info

Publication number
ES2365653T3
ES2365653T3 ES08843228T ES08843228T ES2365653T3 ES 2365653 T3 ES2365653 T3 ES 2365653T3 ES 08843228 T ES08843228 T ES 08843228T ES 08843228 T ES08843228 T ES 08843228T ES 2365653 T3 ES2365653 T3 ES 2365653T3
Authority
ES
Spain
Prior art keywords
source
radiation
point
plane
dose
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
ES08843228T
Other languages
English (en)
Inventor
Véronique MASSE
Maurice Chiron
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.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
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 Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique CEA
Application granted granted Critical
Publication of ES2365653T3 publication Critical patent/ES2365653T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/169—Exploration, location of contaminated surface areas

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Radiation (AREA)
  • Radiation-Therapy Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

Procedimiento de determinación de dosis de radiación emitida, de manera real o simulada, por una fuente de material fisible en una instalación, que incluye las siguientes etapas: - determinación de un número de fisiones que se producen, en función del tiempo, en el material fisible de la fuente, - determinación de coeficientes de atenuación de materiales que constituyen las paredes verticales y los suelos de la instalación y, más generalmente, cualquier pantalla que pueda colocarse en la trayectoria de la radiación emitida, - determinación, a partir de datos geométricos descriptivos de la instalación, entre un plano fuente (PE) sensiblemente perpendicular a las paredes verticales de la instalación y que contiene un punto fuente (E) representativo de la fuente de radiación y un plano de visualización (PV) paralelo al plano fuente, de un conjunto de planos característicos (Pj) perpendiculares al plano fuente y que contienen, cada uno, el punto fuente (E) y al menos una arista de unión entre dos paredes verticales de la instalación, - barrido angular (θ) de los planos característicos alrededor de un eje (Zp) perpendicular al plano fuente y que pasa por el punto fuente (E) para definir al menos un plano de cálculo (Po), - determinación, para el plano de cálculo, de un conjunto de rectas características (Qj), pasando cada recta característica por el punto fuente (E) y por al menos un punto situado en la unión de dos aristas de unión, - en una recta de cálculo (D) situada en la intersección del plano de visualización y del plano de cálculo, determinación de posiciones de puntos de intersección entre la recta de cálculo y las rectas características, - entre los puntos de intersección presentes en la recta de cálculo, selección de los puntos de intersección Δj situados en zonas al aire libre de la instalación, - cálculo de la dosis de radiación d(Δj) presente en cada punto Δj, a partir del número de fisiones en función del tiempo, una distancia que separa el punto fuente (E) del punto Δj y los coeficientes de atenuación de los materiales constitutivos de las paredes verticales y/o de los suelos y/o de cualquier pantalla que separan el punto fuente (E) del punto Δj.

Description

Ámbito técnico y técnica anterior
La invención se refiere a un procedimiento de determinación de dosis de radiación, así como a un procedimiento de determinación de curva de isodosis a partir de dosis determinadas por el procedimiento de la invención.
Las curvas de isodosis son, para una radiación dada, líneas o superficies ficticias a cuyo nivel una dosis de radiación es sensiblemente constante. En el resto de la descripción, el término “radiación” se entenderá como una radiación gamma, una emisión de neutrones o una emisión de partículas alfa.
Se conoce la publicación titulada “Use of simulation technology for prediction of radiation dose in nuclear power plant” (Institution of Electrical Engineers, Stevenage, BG, 2004).
Este documento describe un procedimiento de determinación de dosis de radiaciones emitidas por una fuente de material fisible en una instalación.
Este procedimiento requiere, implícitamente, una etapa de determinación del número de fisiones que se producen, en función del tiempo, en el material fisible, así como una etapa de determinación de los coeficientes de atenuación de los “materiales pantalla” de la instalación.
Las instalaciones que contienen material fisible presentan un riesgo de accidente de criticidad a pesar de todas las disposiciones aplicadas para limitar este riesgo. Las consecuencias de un accidente de criticidad son, entre otras, el riesgo de una exposición externa importante a radiaciones emitidas y el riesgo de una exposición interna por dispersión de productos radiactivos en la instalación y/o el entorno. La invención se aplica a la gestión del riesgo de exposición externa.
Cuando ocurre un accidente de criticidad en una instalación, es importante, por motivos de seguridad, determinar cómo se reparten las radiaciones nocivas resultantes del accidente. Según la técnica conocida, las dosis de radiación se calculan manualmente con la ayuda de ábacos. Estos cálculos por medio de ábacos se efectúan en planos horizontales de la instalación. Si los cálculos en el plano de la fuente (plano horizontal de la instalación que contiene la fuente emisora de las radiaciones nocivas) conducen a una correcta descripción del reparto de las dosis de radiación, los cálculos en los planos paralelos al plano de la fuente conducen a una descripción más compleja de este reparto. Además, debido al carácter manual de los cálculos, los riesgos de error no son despreciables y los tiempos de cálculo pueden ser largos.
La invención no presenta estos inconvenientes.
Descripción de la invención
En efecto, la invención se refiere a un procedimiento de determinación de dosis de radiación emitida por una fuente de material fisible en una instalación que incluye las siguientes etapas:
-
determinación de un número de fisiones que se producen, en función del tiempo, en el material fisible de la
fuente,
-
determinación de coeficientes de atenuación de materiales que constituyen las paredes verticales y los
suelos de la instalación y, más generalmente, cualquier pantalla que pueda colocarse en la trayectoria de la
radiación emitida,
-
determinación, a partir de datos geométricos descriptivos de la instalación, entre un plano fuente
sensiblemente perpendicular a
las paredes verticales de la instalación y que contiene un punto fuente
representativo de la fuente de radiación y un plano de visualización paralelo al plano fuente, de un conjunto
de planos característicos perpendiculares al plano fuente y que contienen, cada uno, el punto fuente y al
menos una arista de unión entre dos paredes verticales de la instalación,
-
barrido angular de los planos característicos alrededor de un eje perpendicular al plano fuente y que pasa
por el punto fuente para definir al menos un plano de cálculo,
-
determinación, para el plano de cálculo, de un conjunto de rectas características, pasando cada recta
característica por el punto fuente y por al menos un punto situado en la unión de dos aristas de unión,
-
en una recta de cálculo situada en la intersección del plano de visualización y del plano de cálculo,
determinación de posiciones de puntos de intersección entre la recta de cálculo y las rectas características,
-
entre los puntos de intersección presentes en la recta de cálculo, selección de los puntos de intersección Δj
situados en zonas al aire libre de la instalación,
-
cálculo de la dosis de radiación d(Δj) presente en cada punto Δj, a partir del número de fisiones en función
del tiempo, de una distancia que separa el punto fuente del punto Δj y de los coeficientes de atenuación de
los materiales constitutivos de las paredes verticales y/o de los suelos y/o de cualquier pantalla que separen el punto fuente del punto Δj.
La invención se refiere asimismo a un procedimiento de determinación de curvas de isodosis de una radiación emitida por una fuente de material fisible en una instalación, caracterizado porque incluye sucesivamente:
-un procedimiento de determinación de dosis según la invención, y -una comparación de las dosis d(Δj) determinadas por el procedimiento de la invención con intervalos de dosis predeterminados de tal manera que:
•
si dos dosis calculadas d(Δj) y d(Δj+1) para dos puntos de intersección seleccionados consecutivos Δjy Δj+1 pertenecen a un mismo intervalo de dosis, se asigna una misma zona de pertenencia entre los dos puntos calculados, y
•
de lo contrario, se busca, por dicotomía, uno o varios puntos Δk cuya dosis (Δk) es un límite de intervalo de dosis, asignándose una misma zona de pertenencia entre dos puntos consecutivos pertenecientes al mismo intervalo de dosis,
- una formación de la curva de isodosis a lo largo de la recta de cálculo, en función de las zonas de pertenencia asignadas a las dosis de radiación calculadas.
Al diseñar una instalación que debe contener material fisible, la aplicación del procedimiento de la invención permite ventajosamente efectuar un análisis de las consecuencias de un posible futuro accidente de criticidad, con objeto de:
-definir las vías de evacuación de la instalación, -determinar puntos de reunión en la instalación, -posicionar las sondas de detección de accidente de criticidad, -balizar las áreas con riesgo de exposición externa en las proximidades de la instalación, si existen.
Cuando ocurre un accidente de criticidad en una instalación, el procedimiento de la invención permite ventajosamente efectuar, en tiempo real:
-una evaluación de los potenciales riesgos de exposición externa de los trabajadores situados en la instalación o a proximidad de la misma, -la implantación de medios o procedimientos necesarios para la gestión de la crisis respecto de la instalación, complementando las actuaciones ligadas al plan de emergencia.
Los procedimientos de la invención se ponen en funcionamiento preferiblemente por medio de un ordenador.
Breve descripción de los dibujos
Otras características y ventajas de la invención aparecerán en una realización preferida, descrita con referencia a las figuras adjuntas, en las cuales:
- la figura 1 representa, de manera simbólica, un ejemplo de instalación que contiene material fisible en el
que es susceptible de producirse un accidente de criticidad, -la figura 2 representa un cuadro sinóptico del procedimiento de determinación de dosis de la invención, -la figura 3 representa un cuadro sinóptico del procedimiento de determinación de curva de isodosis de la
invención, -la figura 4 representa una vista en corte de un volumen de la instalación en el que se calculan las dosis según el procedimiento de la invención, -la figura 5 representa un conjunto de planos característicos útiles para el cálculo de dosis según el procedimiento de la invención, -la figura 6 representa un conjunto de rectas características útiles para el cálculo de dosis según el procedimiento de la invención, -la figura 7 representa un ejemplo de reparto de curvas de isodosis obtenidas en un plano de visualización, según el procedimiento de la invención.
En todas las figuras, las mismas referencias designan los mismos elementos.
Descripción detallada de modos de realización particulares
La figura 1 representa, de forma simbólica, un ejemplo de instalación en la que es susceptible de producirse un accidente de criticidad.
La instalación está constituida, por ejemplo, por un edificio de varios pisos, con varias estancias en cada piso. Distintos sensores de medición Cnm están repartidos en las distintas estancias de la instalación. Los sensores Cnm se destinan , en el caso de un accidente de criticidad real, a efectuar mediciones de las radiaciones que permiten identificar la fuente emisora de la radiación nociva. Los sensores Cnm son, por ejemplo, cámaras de ionización, sensores proporcionales, etc. La instalación se localiza en un triedro directo (x, y, z) de manera que el eje z es el eje vertical según el cual se define la altura de la instalación y el plano (x, y) es un plano horizontal para la instalación.
La figura 2 representa un cuadro sinóptico de las distintas etapas del procedimiento de determinación de dosis de la invención (etapas 1-9).
El procedimiento de la invención empieza con tres etapas de lectura de datos, es decir una etapa 1 de lectura de datos de geometría de la instalación, una etapa 2 de lectura de datos de fuente y una etapa 3 de lectura de datos de trazado. El orden en el que se efectúan las etapas de lectura 1 a 3 es indiferente, pudiendo efectuarse estas tres etapas simultáneamente.
Los datos geométricos de la instalación G leídos en la etapa 1 son representativos de la configuración volumétrica del edificio (las distintas estancias del edificio, la envuelta del edificio) y de la configuración geométrica de las pantallas de protección que están presentes en el edificio.
Los datos de fuente S leídos en la etapa 2 son datos relativos a la fuente que emite las radiaciones. Están constituidos por el número de fisiones que se producen, en función del tiempo, en el lugar del accidente, datos geométricos que describen la geometría de la cuba en la que se ha producido el accidente (fuente puntual o fuente volumétrica) y datos del medio que caracterizan el medio en el que se ha producido el accidente (medio líquido, polvo, metal). Estos datos se obtienen, por una parte, a partir de mediciones de radiación proporcionadas por los sensores que han detectado el accidente de criticidad y, por otra, a partir de información previamente grabada relativa a los distintos productos presentes en la instalación.
Los datos de trazado T incluyen:
-datos umbral que definen intervalos de dosis predeterminados en los que se repartirán las dosis calculadas, y -datos geométricos relativos a los puntos de la instalación donde se desea que se calculen las dosis (plano(s) de visualización, zona(s) particular(es) o punto(s) particular(es) de la instalación).
Las etapas de lectura 1, 2 y 3 van seguidas de una etapa 4 de evaluación de los coeficientes de atenuación K(Mk) (k=1, 2, …, n) de los distintos materiales (Mk (k=1, 2, …, n) que componen la instalación y de datos teóricos Do(P) que representan las dosis de radiación que estarían presentes, en ausencia de cualquier pared o pantalla, en distintos puntos P de la instalación. La etapa de cálculo 4 de los coeficientes K(Mk) y de los datos Do(P) se efectúa a partir de los datos G y S, de datos T y de datos internos I que incluyen un modelo matemático de coeficiente de atenuación para cada tipo de material. De manera preferente, un coeficiente de atenuación se presenta en forma de ecuación polinómica. A título de ejemplo no limitativo, un coeficiente de atenuación K(Mk) de un material Mk atravesado por una radiación se escribe con la forma:
K(Mk) = aX+ bY+ cXY+ dX2 + eY2 + fZ + gW
Los coeficientes a, b, c, d, e, f, g son parámetros conocidos de valor fijo que son característicos del material Mk del que se intenta evaluar el coeficiente de atenuación. Las magnitudes X, Y, Z son variables características de la fuente de radiación y la magnitud W es una variable que representa el grosor del material Mk atravesado (W se precisará posteriormente). De manera más precisa, la variable X depende del tipo de fuente (líquido, polvo, metal), la variable Y depende del volumen de la fuente y la variable Z depende del tiempo que ha transcurrido entre el accidente y el momento en que se determina el coeficiente. Los coeficientes a, b, c, d, e, f y g son datos que pertenecen al conjunto de datos I mencionado anteriormente. Los datos X, Y, Z son datos que pertenecen al conjunto de datos S y el dato W se calcula a partir de los datos geométricos G y de datos de trazado T.
Para un tipo de fuente dado, la cantidad aX + bY + cXY + dX2 + eY2 + fZ es un término constante Ko. De este modo, la magnitud K(Mk) se expresa en forma de función de la única variable W, es decir:
K(Mk) = g xW + Ko
De manera más general, los datos internos I incluyen, además de las ecuaciones matemáticas de los coeficientes de atenuación y los coeficientes a, b, c, d, e, f, g, los siguientes datos:
-las unidades en las que se desea que se calculen las dosis (dosis en el aire Gy o equivalente de dosis Sv), y -las condiciones del cálculo de los coeficientes de atenuación (coeficientes de corrección de distancia entre la fuente y los puntos de cálculo).
En paralelo a la etapa de cálculo 4 se efectúan cuatro etapas de cálculo elementales 5, 6, 7 y 8. La etapa 5 es una etapa de determinación de planos característicos útiles para el cálculo de dosis. A título de ejemplo no limitativo, la figura 5 representa un conjunto de planos característicos Pj. La figura 5 representa una vista en corte de la instalación según el plano horizontal PE que contiene el punto E al que está asimilada la fuente emisora de las radiaciones nocivas. Los planos característicos están construidos entre el plano PE y un plano de visualización PV. El plano de visualización PV es un plano paralelo al plano PE en el que se calculan las dosis (véase la figura 4). El plano de visualización PV es un dato fijado mediante los datos de trazado T. Cada plano característico Pj es un plano vertical, es decir un plano perpendicular a los planos horizontales PE y PV, que contiene el punto E al que está asimilada la fuente emisora de las radiaciones nocivas y al menos una arista de unión entre dos paredes verticales incluidas entre los planos PE y PV. El conjunto de todos los planos que es posible construir según la regla precisada anteriormente constituye los planos característicos de la invención. En consecuencia, están afectadas todas las aristas de todas las estancias incluidas entre PE y PV y que son perpendiculares a los planos PE y PV. El conjunto de planos característicos se selecciona a partir de los datos geométricos G.
En la etapa 6 del procedimiento de la invención, se efectúa a continuación un barrido entre los planos característicos Pj con el fin de determinar distintos planos de cálculo PC. Los planos de cálculo PC se obtienen entonces mediante rotación, con un paso angular θ, de los planos característicos Pj alrededor de un eje Zp perpendicular a los planos PE y PV y que pasan por el punto fuente E. Cada plano de cálculo PC es un plano en el que se efectúa un cálculo de dosis, según una dirección dada, como se describe a continuación, a título de ejemplo no limitativo, en un plano de cálculo particular, en referencia a la figura 6.
A la etapa 6 de determinación de los planos de cálculo le sigue una etapa 7 de determinación de rectas características Qj en cada plano de cálculo. Para un plano de cálculo dado, una recta característica Qj pasa por el punto fuente E y por al menos un punto situado en la unión de dos aristas situadas en el plano de cálculo. Todas las rectas que es posible construir según la regla precisada anteriormente constituyen el conjunto de rectas características Qj de la invención para el plano de cálculo considerado. Por cada construcción, un plano de cálculo PC se divide en dos semiplanos simétricos entre sí con relación al eje vertical Zp. El conjunto de rectas características relativas a un plano de cálculo se divide, por lo tanto, en dos semiconjuntos de rectas características. La figura 6 ilustra, a título de ejemplo no limitativo, un semiconjunto de rectas características para el plano de cálculo PC de la figura 5. El semiplano de cálculo corta el plano de visualización PV según una recta D de vector unitario ū. Se determina entonces un conjunto de puntos característicos Δj pertenecientes a la recta D (etapa 8 del procedimiento de la invención). Se obtiene un punto característico Δj mediante la intersección de una recta característica Qj y de la recta D. La figura 6 representa, a título de ejemplo, una sucesión de puntos característicos Δ0, Δ1, Δ2, …, Δn. Los puntos característicos Δj poseen una posición geométrica conocida en la instalación. La estructura de la instalación entre el punto fuente E y cada uno de los puntos Δj es asimismo conocida (véase la figura 6). De este modo, a partir de los datos calculados anteriormente D0 (P) y K(Mk) de la posición conocida de los puntos Δj con relación a la fuente emisora E y de la estructura conocida de la instalación entre la fuente E y los puntos Δj, se puede calcular la dosis de radiación d(Δj) presente en cada punto Δj (etapa 9 del procedimiento de la invención).
La recta de cálculo D está constituida por zonas al aire libre y zonas de pared o pantalla. El cálculo de las dosis sólo presenta un interés real en las zonas al aire libre. Por lo tanto, sólo se efectúa el cálculo de las dosis d(Δj) preferiblemente para los puntos Δj situados en las zonas al aire libre.
El cálculo de la dosis en un punto Δj se obtiene mediante la siguiente ecuación:
imagen1
-D0 (P) es la dosis calculada, en ausencia de paredes y pantallas, en un punto arbitrario predeterminado P situado, en la trayectoria de la radiación, a una distancia l0 del punto fuente E (en el caso de una fuente volumétrica, el punto E es el centro del volumen de la fuente),
-Cd es un coeficiente de corrección de distancia tal que:
imagen1
donde l0 es la distancia mencionada anteriormente e l la distancia desde el punto fuente E hasta el punto Δj, y
-K(Mk) es el coeficiente de atenuación del material Mk mencionado anteriormente. A continuación, se precisa el coeficiente de atenuación K(Mk). Como ya se ha mencionado anteriormente, el coeficiente de atenuación de un material Mk atravesado por la radiación se escribe:
K(Mk) = g xW + K0
donde la magnitud W representa la distancia recorrida por la radiación a través del material Mk. De manera preferente, se define la magnitud W en función del ángulo α que forma la dirección de la radiación que atraviesa la pared o pantalla de material Mk con la normal al plano de dicha pared o dicha pantalla:
-Para un ángulo α incluido entre 0º y un valor límite predeterminado αlim (0 < αlim < π/2), W es el grosor real de material atravesado, y -Para un ángulo α incluido entre el valor límite predeterminado αlim y π/2, W es el valor Wlim del grosor de la pared o pantalla que corresponde al ángulo αlim.
Se elige la cantidad αlim de manera a no subestimar la dosis d(Δj) para los ángulos importantes. Dicha cantidad αlim varía con el tipo de radiación.
La figura 3 representa un cuadro sinóptico del procedimiento de determinación de curvas de isodosis de la invención. El procedimiento de determinación de curvas de isodosis repite el conjunto de etapas 1-9 descritas anteriormente y una etapa 10 adicional de reparto de las dosis calculadas en intervalos de dosis de valores predeterminados.
A continuación, se describe un ejemplo de dicho reparto en el que las dosis d(Δj) se reparten en n intervalos de dosis |di, di+1| (i=1,2,…,n).
El reparto de las dosis calculadas en los distintos intervalos de dosis se efectúa así:
-si las dosis d(Δj) y d(Δj+1) calculadas para dos puntos característicos sucesivos Δjy Δj+1 de una misma zona al aire libre pertenecen a un mismo intervalo |di, di+1|, se asigna entonces una misma zona Z entre dichos puntos;
-de lo contrario, se calcula la dosis d ((δ D+δ D+ i) /2) en el punto medio (δ D+δ D+ i) /2 y se busca, por dicotomía, uno o varios puntos Δk cuya dosis d(Δk) es un límite de intervalo de dosis, asignándose una misma zona de pertenencia entre dos puntos consecutivos pertenecientes al mismo intervalo de dosis.
Es posible obtener entonces, para una misma recta de cálculo D, una curva de isodosis C(Zi) en función de distintas zonas Zi (etapa 10 del procedimiento). Obtenidas para el conjunto de rectas de cálculo, es decir para el conjunto de planos de cálculo, las curvas de isodosis C(Zi) constituyen una superficie de isodosis en el conjunto del plano de visualización PV. A título de ejemplo no limitativo, la figura 7 ilustra un reparto de las dosis calculadas en cinco zonas Z1-Z5.
En el caso particular en que el plano de visualización es el plano horizontal PE que contiene el punto fuente E, todas las paredes y pantallas son atravesadas perpendicularmente a su superficie (α = 0). Los valores de los coeficientes de atenuación son entonces valores constantes K. De este modo, los cálculos quedan ventajosamente muy simplificados.
Si las dosis d(Δj) y d(Δj+1) no pertenecen al mismo intervalo [di, di+1], la distancia l que separa el punto fuente E de un punto Δj donde la dosis de radiación d(Δj) corresponde a un límite de intervalo se expresa sencillamente mediante la ecuación:
imagen1

Claims (7)

  1. REIVINDICACIONES
    1. Procedimiento de determinación de dosis de radiación emitida, de manera real o simulada, por una fuente de material fisible en una instalación, que incluye las siguientes etapas:
    -determinación de un número de fisiones que se producen, en función del tiempo, en el material fisible de la fuente,
    -determinación de coeficientes de atenuación de materiales que constituyen las paredes verticales y los suelos de la instalación y, más generalmente, cualquier pantalla que pueda colocarse en la trayectoria de la radiación emitida,
    -determinación, a partir de datos geométricos descriptivos de la instalación, entre un plano fuente (PE) sensiblemente perpendicular a las paredes verticales de la instalación y que contiene un punto fuente
    (E) representativo de la fuente de radiación y un plano de visualización (PV) paralelo al plano fuente, de un conjunto de planos característicos (Pj) perpendiculares al plano fuente y que contienen, cada uno, el punto fuente (E) y al menos una arista de unión entre dos paredes verticales de la instalación,
    -barrido angular (θ) de los planos característicos alrededor de un eje (Zp) perpendicular al plano fuente y que pasa por el punto fuente (E) para definir al menos un plano de cálculo (Po),
    -determinación, para el plano de cálculo, de un conjunto de rectas características (Qj), pasando cada recta característica por el punto fuente (E) y por al menos un punto situado en la unión de dos aristas de unión,
    -en una recta de cálculo (D) situada en la intersección del plano de visualización y del plano de cálculo, determinación de posiciones de puntos de intersección entre la recta de cálculo y las rectas características,
    -entre los puntos de intersección presentes en la recta de cálculo, selección de los puntos de intersección Δj situados en zonas al aire libre de la instalación,
    - cálculo de la dosis de radiación d(Δj) presente en cada punto Δj, a partir del número de fisiones en función del tiempo, una distancia que separa el punto fuente (E) del punto Δj y los coeficientes de atenuación de los materiales constitutivos de las paredes verticales y/o de los suelos y/o de cualquier pantalla que separan el punto fuente (E) del punto Δj.
  2. 2. Procedimiento de determinación de curvas de isodosis de una radiación emitida, de manera real o simulada, por una fuente de material fisible en una instalación, caracterizado porque incluye sucesivamente: -un procedimiento de determinación de dosis según la reivindicación 1, y -una comparación de las dosis d(Δj) determinadas por el procedimiento de la reivindicación 1 con
    intervalos de dosis predeterminados de tal manera que:
    •
    si dos dosis calculadas d(Δj) y d(Δj+1) para dos puntos de intersección seleccionados consecutivos Δjy Δj+1 pertenecen a un mismo intervalo de dosis, se asigna una misma zona de pertenencia (Zi) entre los dos puntos calculados, y
    •
    de lo contrario, se calcula una dosis de radiación d( (Δj+Δj+1) /2) en el punto medio situado entre los dos puntos consecutivos Δjy Δj+1 y se busca, por dicotomía, uno o varios puntos Δk cuya dosis d(Δk) es un límite de intervalo de dosis, asignándose una misma zona de pertenencia entre dos puntos consecutivos pertenecientes al mismo intervalo de dosis,
    - una formación de la curva de isodosis a lo largo de la recta de cálculo, en función de las zonas de pertenencia asignadas a las dosis de radiación calculadas.
  3. 3.
    Procedimiento según la reivindicación 2, en el que el barrido angular se efectúa sobre 360 grados, de manera que se reúne un conjunto de curvas de isodosis establecidas a lo largo de un conjunto de rectas de cálculo, para constituir una representación de las isodosis en el conjunto del plano de visualización.
  4. 4.
    Procedimiento según una cualquiera de las reivindicaciones 1 a 3, en el que la dosis de radiación presente en el punto de intersección seleccionado viene dada por la ecuación:
    imagen1
    donde -D0 (P) es la dosis calculada, en ausencia de toda pantalla y de pantallas, en un punto arbitrario predeterminado (P) situado, en la trayectoria de la radiación que se propaga entre el punto fuente (E) y el punto Δj, a una distancia l0 del punto fuente (E), -Cd es un coeficiente de corrección de distancia tal que:
    imagen1
    donde l0 es la distancia mencionada anteriormente y l la distancia desde el punto fuente (E) hasta el punto Δj, y
    -K(Mk) es el coeficiente de atenuación calculado de un material Mk atravesado por una radiación que se propaga entre el punto fuente (E) y el punto Δj.
  5. 5. Procedimiento según la reivindicación 4, en el que el coeficiente de atenuación K(Mk) viene dado por la fórmula:
    K(Mk) = g xW + K0,donde
    5
    - W es una magnitud que representa el grosor atravesado del material Mk,
    - g es un coeficiente conocido característico del material Mk,
    - K0 es un término conocido que depende de la fuente de radiación y del material Mk.
  6. 6. Procedimiento según la reivindicación 5, en el que la magnitud W se define en función del ángulo α que 10 forma la dirección de la radiación con la normal a la pared vertical del material Mk, de tal manera que:
    -para un ángulo α incluido entre 0º y un valor límite predeterminado αlim (0 < αlim < π/2), W es el grosor real de material atravesado, y -para un ángulo α incluido entre el valor límite predeterminado αlim y π/2, W es el grosor de material 15 atravesado por una radiación cuya dirección con la normal a la pared vertical es el ángulo αlim.
  7. 7. Procedimiento según cualquiera de las reivindicaciones anteriores, en el que la radiación emitida es una radiación gamma o una emisión de neutrones o una emisión de partículas alfa.
    20
ES08843228T 2007-10-22 2008-10-20 Procedimiento para determinar una dosis de radiación y procedimiento para determinar una curva de isodosis asociada. Active ES2365653T3 (es)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0758465A FR2922656B1 (fr) 2007-10-22 2007-10-22 Procede de determination de dose de rayonnement et procede de determination de course isodose associe
FR0758465 2007-10-22

Publications (1)

Publication Number Publication Date
ES2365653T3 true ES2365653T3 (es) 2011-10-07

Family

ID=39645671

Family Applications (1)

Application Number Title Priority Date Filing Date
ES08843228T Active ES2365653T3 (es) 2007-10-22 2008-10-20 Procedimiento para determinar una dosis de radiación y procedimiento para determinar una curva de isodosis asociada.

Country Status (13)

Country Link
US (1) US8536515B2 (es)
EP (1) EP2201407B1 (es)
JP (1) JP5461413B2 (es)
CN (1) CN101836130B (es)
AT (1) ATE507493T1 (es)
DE (1) DE602008006584D1 (es)
DK (1) DK2201407T3 (es)
ES (1) ES2365653T3 (es)
FR (1) FR2922656B1 (es)
PL (1) PL2201407T3 (es)
PT (1) PT2201407E (es)
SI (1) SI2201407T1 (es)
WO (1) WO2009053324A1 (es)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2953041B1 (fr) * 2009-11-24 2011-12-30 Commissariat Energie Atomique Procede de determination d'une grandeur operationnelle permettant une gestion simplifiee et optimisee d'echantillons irradies dans des installations et garantissant le respect d'objectifs de surete
JP6215530B2 (ja) * 2012-12-26 2017-10-18 中部電力株式会社 放射線の表示方法、放射線表示装置
CN106526647B (zh) * 2015-09-09 2019-12-03 同方威视技术股份有限公司 放射源检测方法和系统
JP7399031B2 (ja) * 2020-06-09 2023-12-15 三菱重工業株式会社 線量推定装置、線量推定方法および線量推定プログラム

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2740899A (en) * 1956-04-03 cummings
JPH07202246A (ja) * 1994-01-10 1995-08-04 Fuji Electric Co Ltd 半導体放射線検出素子
JP2000221292A (ja) * 1999-02-01 2000-08-11 Hitachi Ltd 放射線照射施設内作業者被爆線量当量シュミレーション装置
US6795801B1 (en) * 1999-08-13 2004-09-21 Electric Power Research Institute, Inc. Apparatus and method for analyzing anisotropic particle scattering in three-dimensional geometries
US7233888B2 (en) * 2002-07-09 2007-06-19 General Electric Company Monte Carlo criticality-mode systems and methods for computing neutron and gamma fluence in a nuclear reactor
JP3741694B2 (ja) * 2003-06-04 2006-02-01 東京電力株式会社 放射線量算出システム、放射線量算出方法、およびプログラム
JP2005049148A (ja) * 2003-07-31 2005-02-24 Hitachi Ltd 放射線線量率分布の可視化装置
US7197404B2 (en) * 2004-03-01 2007-03-27 Richard Andrew Holland Computation of radiating particle and wave distributions using a generalized discrete field constructed from representative ray sets

Also Published As

Publication number Publication date
US20100312532A1 (en) 2010-12-09
FR2922656B1 (fr) 2009-12-11
US8536515B2 (en) 2013-09-17
JP5461413B2 (ja) 2014-04-02
DE602008006584D1 (de) 2011-06-09
FR2922656A1 (fr) 2009-04-24
SI2201407T1 (sl) 2011-08-31
DK2201407T3 (da) 2011-08-15
WO2009053324A1 (fr) 2009-04-30
CN101836130B (zh) 2012-11-21
ATE507493T1 (de) 2011-05-15
EP2201407B1 (fr) 2011-04-27
PT2201407E (pt) 2011-06-29
JP2011501171A (ja) 2011-01-06
EP2201407A1 (fr) 2010-06-30
CN101836130A (zh) 2010-09-15
PL2201407T3 (pl) 2011-11-30

Similar Documents

Publication Publication Date Title
Monsalve et al. Application of laser scanning for rock mass characterization and discrete fracture network generation in an underground limestone mine
Li et al. Impacts of wind fields on the distribution patterns of traffic emitted particles in urban residential areas
ES2365653T3 (es) Procedimiento para determinar una dosis de radiación y procedimiento para determinar una curva de isodosis asociada.
Xie et al. A study of diurnal and short-term variations of indoor radon concentrations at the University of Michigan, USA and their correlations with environmental factors
Yu et al. Long-term measurements of radon progeny concentrations with solid-state nuclear track detectors
JP2013113594A (ja) 空間線量評価装置および方法
JP2012251918A (ja) 放射能量分布情報生成装置
JP5608085B2 (ja) 危険を伴う施設内で発生する事故の時間的経過を求める方法
Gulan Analysis of Long-Term Monitoring of Radon Levels in a Low-Ventilated, Semi-Underground Laboratory—Dose Estimation and Exploration of Potential Earthquake Precursors
Röttger et al. The German thoron progeny chamber—concept and application
Nikezic et al. Theoretical basis for long-term measurements of equilibrium factors using LR 115 detectors
Burghele et al. Intercomparison between radon passive and active measurements and problems related to thoron measurements
CN107783174B (zh) 一种核设施液态流出物在线监测仪探测效率的校准方法
Sima Monte Carlo simulation of radon SSNT detectors
Otahal et al. Determination of the radon progeny activity size distribution at different types of workplaces
CN201732166U (zh) 基于虚拟源原理的航空放射性测量系统刻度装置
JP7406769B2 (ja) 赤外線測定画像解析装置及び赤外線測定画像解析方法
Rizo-Maestre et al. Radon gas as an indicator for air quality control in buried industrial architecture: Rehabilitation of the old Británica warehouses in alicante for a tourist site
Alexandrov et al. Principles of Data Processing of the Muographic Experiment in the Cave Complex of Pskovo-Pechersky Monastery
Navratilova Rovenska et al. Intercomparison of radon and radon progeny concentration measurements performed in the Historic Silver Mine in Tarnowskie Góry, Poland
CN115267872B (zh) 辐射区域内空间γ辐射剂量及表面β放射性的测量方法
Spring et al. An Experimental Evaluation of Roof Reduction Factors within a Multi-story Structure
CN114325795A (zh) 一种有害辐射场中城市内测量点分布和测量值获取方法
Brown Requist et al. Real-time contamination monitoring platform using fixed air sensing in underground mining environments
Afeni et al. Effect of Line-Of-Sight Obstruction on Distance Measurement Accuracy in Prism Monitoring Surveying