WO2016000666A1 - Method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere - Google Patents
Method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere Download PDFInfo
- Publication number
- WO2016000666A1 WO2016000666A1 PCT/CZ2015/000062 CZ2015000062W WO2016000666A1 WO 2016000666 A1 WO2016000666 A1 WO 2016000666A1 CZ 2015000062 W CZ2015000062 W CZ 2015000062W WO 2016000666 A1 WO2016000666 A1 WO 2016000666A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- radon
- thoron
- atmosphere
- bulk
- activity
- Prior art date
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims description 20
- 229910052704 radon Inorganic materials 0.000 claims abstract description 107
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 106
- 230000000694 effects Effects 0.000 claims abstract description 66
- XQVKLMRIZCRVPO-UHFFFAOYSA-N 4-[(2-arsonophenyl)diazenyl]-3-hydroxynaphthalene-2,7-disulfonic acid Chemical compound C12=CC=C(S(O)(=O)=O)C=C2C=C(S(O)(=O)=O)C(O)=C1N=NC1=CC=CC=C1[As](O)(O)=O XQVKLMRIZCRVPO-UHFFFAOYSA-N 0.000 claims abstract description 64
- 241000679125 Thoron Species 0.000 claims abstract description 64
- 238000005259 measurement Methods 0.000 claims abstract description 57
- 239000002689 soil Substances 0.000 claims abstract description 36
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 23
- 238000005070 sampling Methods 0.000 claims description 28
- 239000007789 gas Substances 0.000 claims description 23
- 238000011156 evaluation Methods 0.000 claims description 22
- 230000005855 radiation Effects 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 17
- 239000002680 soil gas Substances 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 11
- 239000003673 groundwater Substances 0.000 claims description 8
- 239000012528 membrane Substances 0.000 claims description 7
- 238000000926 separation method Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000008014 freezing Effects 0.000 claims description 5
- 238000007710 freezing Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000008901 benefit Effects 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 238000004611 spectroscopical analysis Methods 0.000 claims description 2
- 239000003643 water by type Substances 0.000 claims 4
- 230000008030 elimination Effects 0.000 claims 2
- 238000003379 elimination reaction Methods 0.000 claims 2
- 230000000750 progressive effect Effects 0.000 claims 2
- 238000012937 correction Methods 0.000 claims 1
- 238000012795 verification Methods 0.000 claims 1
- 239000003570 air Substances 0.000 description 35
- 238000010586 diagram Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000009792 diffusion process Methods 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 238000002156 mixing Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 206010044565 Tremor Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 229910052776 Thorium Inorganic materials 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000005494 condensation Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 230000005258 radioactive decay Effects 0.000 description 2
- 206010037844 rash Diseases 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000002344 surface layer Substances 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 229910052770 Uranium Inorganic materials 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000002274 desiccant Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 229930014626 natural product Natural products 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
Definitions
- the invention relates to the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere, namely for prediction of earthquakes.
- Beside predictions based on analysis of known historic earthquakes other predictions are based on changes of physical fields of the earth, especially electric, magnetic and gravity fields, monitoring of chemical changes of ground waters and gases, or changes of elevation of the earth's surface.
- Behaviour of animals is also studied for the prediction of earthquakes, as animals are sensitive to barometric and acoustic phenomena imperceptible to humans, which precede an imminent earthquake.
- the lithosphere is the solid envelope of the earth, consisting of the crust and the uppermost strata of the mantle.
- the lithosphere is not a compact envelope; it is split to large bodies known as the lithospheric plates, which "float” and collide with each other on top of the plastic layer of the earth's mantle, the asthenosphere.
- Monitoring of instantaneous bulk radon and thoron activity in soil gas is used as part of the monitoring of volcanic activity, e.g. in the Etna area (e.g. S.
- Document WO 2014/049408 describes the measuring equipment for the prediction and localization of earthquake epicentres by measuring anomalies of radon content in soil, which includes two interconnected measuring modules, of which one is installed in the soil, the other on the surface.
- the measured values of bulk radon activity are transferred to the server, where data from a number of independent measurement points are compared and processed using an algorithm that compares measured values from a network of ground and surface measurement points to the mean value of bulk radon activity, with the aim of predicting earthquakes and volcanic eruptions as well as calculating the epicentre and magnitude.
- the algorithm of data processing is not described in detail.
- the above measuring equipment only monitors and evaluates the time course of radon concentration in soil and on the surface, which does not seem adequate for reliable prediction of the time and place of a possible earthquake.
- the object of the invention is to find a new method of monitoring changes in the earth's lithosphere and atmosphere to enable an easier, more reliable and longer- term prediction of earthquakes than what today's state of technology allows.
- the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere are based on the understanding of the phenomena preceding the earthquake and its course in time.
- the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere are based on the understanding that changes preceding earthquakes are manifested in more natural products and more physical quantities than just anomalies of bulk radon activity in soil and on its surface.
- the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere as proposed by the invention are furthermore based on the fact that great stress occurs in the fault between lithospheric plates prior to an earthquake; the stress causes changes in the concentration of soil gases including radon and increased escape thereof; the changes thus provoked depend on the magnitude of the earthquake, its location and time. Radioactivity from soil gases, namely radon and its short-time decay products, ionizes air on a large scale and causes
- the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere as proposed by the invention are based on systematic continuous monitoring of important variables and analysis of their changes in time, including the velocity of their increase and decrease, systemic shift of time sequences, and further, more detailed analysis of the sequences using advanced methods of statistical data processing.
- Bulk radon activity in soil is understood as radon activity in a unit of volume (1m 3 ) of air collected from pores in soil.
- Areal velocity of radon exhalation from soil represents the activity of radon diffusing from a unit of area (soil in this case) into the atmosphere per unit of time.
- the two quantities may or may not be independent of one another and their ratio is another monitored parameter in the prediction of earthquake phenomena.
- Variation of other quantities in the atmosphere is related to the release of radon and/or thoron, as well as their quantification and identification of the trajectory of their movement in the atmosphere. Their interdependence in time is important for correct earthquake prediction.
- Simultaneous measurement of radon and/or thoron is conducted at different altitudes in the atmosphere and at different depths in soil, bedrock and water.
- monitoring point detectonically disturbed area, geological fault zone, groundwater resource, etc.
- Monitoring at different altitudes and depths makes it possible to track radon from its source in the bedrock to the ground surface, whereby significant differences of half-times of decay of radon and thoron can be helpful in general monitoring of the dynamism of soil gases.
- comphreensive monitoring of the above profiles can help to offset disturbances causing variation of radon and thoron concentration in time, which cannot be attributed to the phenomena of interest. (E.g. the effect of rainfall and snowfall, freezing of the surface layer of soil, temperature inversion.)
- EML equivalent mixing layer
- Continuous simultaneous measurement of atmospheric (mainly) radon at selected altitudes including its and thoron's DP, radon and/or thoron in water and soil air, areal velocity of radon and/or thoron exalation from soil, values of atmospheric dose inputs of gamma or beta radiation, solar radiation, gradients of temperature and humidity/moisture in air, soil and water, together with the measurement of precipitation, speed and trend of wind, constitute the base measurement profiles.
- the minima and maxima of bulk radon activity correspond to the minima and maxima of bulk activities of its DP.
- the measured values of BRA vary within well-measurable limits ca. 5-40 Bq/m 3
- equivalent bulk radon activity (EBRA) varies within well-measurable limits ca. 2-20 Bq/m 3 .
- BA of 218 Po varies in EBRA maxima on the order of tens of Bq/m 3
- BA of 214 Pb, 214 Bi varies on the order of single Bq/m .
- the well-measurable coefficient of imbalance varies from 0.5 to 0.7 in agreement with published data.
- Measured activities of atmospheric DP of thoron vary at the margin of measurability on the order of tenths of Bq/m 3 .
- FIG. 1 Block diagram of the equipment according to the invention Obr. 1
- Primary and filtered time sequences are processed using methods of mathematical statistics, whereby namely long-term trends are compared against the characteristic background of the site and remarkable changes in time (anomalies), which are the subject of interest in relation to changes in the earth's crust. Particularly significant are velocities of increase or decrease of the variables monitored, and the
- the modular monitoring system makes it possible at the same time to create calibration protocols for given site using long-term comparison of time sequences of the variables monitored and seismic activity on the site of interest.
- XI f(t, z) areal velocity of radon and/or thoron exhalation from soil;
- X2 f(t,z) bulk activity of radon and/or thoron from soil gas
- X3 f(t,z) bulk radon activity in the atmosphere
- X4 f(t, z) bulk activity of radon and/or thoron in water.
- the equipment based on the invention consists for example of
- CCU central control unit 1
- Fig. 1 shows an example of the configuration of the proposed instrument.
- Analog and digital outputs of measurement from modules A to F are stored in CCU (1) online in its external memory and, either all at a time or selected ones only, transferred in one-minute intervals, using built-in emitting modem (2) and protocol GPRS, to the user PC with a connected receiving modem (6).
- Transfer parameters including the selection of measured variables can be configured using the program Wincentral, installed in the user PC.
- Wincentral furthermore allows remote setting and diagnosing of selected parameters of key HW modules A to F from the user PC, using built-in commands.
- Wincentral makes it possible to control and read the entire network of stations, each of them including a CCU (1) with any combination of modules A to F.
- Visualization of files of data received in the user PC is provided by the program Visualis, which also allows to export them to common spreadsheets.
- CCU 1 is the heart of the entire system; it logs the relevant data from modules A to F directly in the required variables or converts raw measured data to the required variables by means of mathematical algorithms, using its built-in microprocessor.
- the overall concept of the station includes remote wireless access to all the CCU functions using the program Wincentral, which allows complete control and setting of all the measuring modules from a single PC.
- the program Wincentral makes it possible for the CCU to set up a network of measuring stations controlled from a single PC, of which each station must include a central control unit.
- Power source block can be fed from the grid or battery (4), or solar panels (5).
- measuring modules A to F are as follows:
- Module A offers the possibility of continuous spectrometric and non- spectrometric detection of atmospheric radon and/or thoron. Modul A predstavuje
- - Module D offers the possibility of measurement of photon dose equivalent in the range from natural background up to accident level, or of beta radiation using any kind of measuring instruument, e.g. one based on GM counter and/or coefficient of imbalance F between bulk radon activity and its equivalent bulk activity (EBRA) in the atmosphere.
- EBRA equivalent bulk activity
- - Module E includes a set of adequate voltage and amperage sensors, capable of measuring the necessary physical parameters of weather, soil and water (rain gauges, anemometers, thermometers, solar radiation, etc.).
- - Module F offers the possibility of simultaneous continuous measurement of bulk radon and throron activity in soil gas, which is sampled from a groundwater resource through a special separation unit with an immersed semi-permeable membrane.
- Figure 2 shows possible arrangement of measuring module A for the measurement of atmospheric radon and/or thoron, in the single-detector option for continuous collection of samples of atmospheric radon and/or thoron, based on non- spectrometric detection with a scintillation chamber with ZnS-based
- radon and/or thoron then causes scintillations within the sensitive volume of the chamber. These are registered by the photomultiplier of the scintillation chamber and taken by the counter as time pulses.
- the pulse count obtained per unit of time represents the rate of bulk activity of radon and/or thoron.
- the CCU uses published algorithms to calculate radon and/or thoron values. The results of the measurement are saved in CCU memory and at the same time transferred through the built-in modem online at one-minute intervals to the user PC by means of GPRS protocol.
- An advantage of the measuring module is that it offers the choice of manual or remote setting of the detector's own background using radon-free technical gas, contained in pressure cylinder (28), with the aid of gas management system (20, 22, 23, 24).
- Menu of the evaluation unit (11) can be used for manual setting, the evaluation unit consistsing of a two-track multichannel analyzer (42) (MCA, MC 2000).
- MCA two-track multichannel analyzer
- the remote option is provided by the CCU on command from the user PC.
- the two-track multichannel analyzer 42 furthermore allows another simultaneous detection of radon and/or thoron by connecting an additional, spectrometric or non-spectrometric detector of alpha, beta or gamma radiation.
- the stability of the detector in use is achieved by means of working standard 241 Am (29), which is inserted into the detection space of the detector when needed, and is controlled from the user PC by means of an electro-magnetic valve.
- module A To allow the use of module A at low temperature, its critical parts, i.e.
- electromagnetic valves (15, 20) are insulated and tempered using adjustable thermostats (21).
- Temperature sensor (13) helps to offset the effect of temperature variation, which affects namely the response of detector photomultiplier.
- Exchangeable filter (31) capturing air humidity is installed to eliminate the negative effect of condensation of air humidity inside the detector, and to improve the efficiency of detection and spectrometry of the measuring instruments; furthermore, optional heating of the adequately insulated detector is installed to react to outside air humidity, using the adjustable air humidity sensor of the hygrostat (30).
- Samples can be collected from a grommet of any size and shape, fixed in any suitable way and opening on the outside of the casing at a suitable collecting height above ground. Simultaneous measurement of the height profile of atmospheric radon and/or thoron is possible using an additional detector or setting up a network of measuring stations.
- Figure 3 shows a possible arrangement of measuring module B to measure samples of atmospheric bulk activities of DP of radon and thoron, captured on the endless filter (filtering tape).
- DP of radon and thoron in the atmosphere are fixed to aerosol. This is why they are detected applying the well-known principle of collection on filter during the suction of a defined volume of air sample at a known flow velocity through sensor (43) and filtering tape (50) using a powerful pump (44).
- the bulk activity of DP is then calculated from the known volume of air transiting the filter per unit of time and their activity determined from the filter using a spectroscopic semiconductor detector. Alpha radiation emanating from the DP is used for the detection, as it allows registering a sufficient count of signal pulses.
- the sample collected at filter is first subjected to the first detector (46), by whose means and using adequate timing of the measurement (choice of measurement times and intervals between them) the bulk activity of individual short-time radon DP ( 218 Po, 214 Pb, 214 Bi) is determined.
- the sample is then subjected to the second detector (47) for analogous determination, again with adequate timing, of the bulk activity of short-time thoron DP ( 212 Pb, 212 Bi).
- thoron DP measurement is finished, air samples are collected again from the "endless" filter (filter coil continuously winds off and then winds back to the beginning instantaneously) and the entire measurement cycle is repeated.
- the proposed algorithm provides for both spectrometric and non-spectrometric determination of the bulk activity of individual DP of radon and thoron, the necessary measurement of the background of the tape, simultaneous calibration of both detectors, and irradiation of the tape by means of the couple of double standards 41 Am (53, 54).
- a certain disadvantage of the use of detection of alpha radiation from the activity of the DP captured on the filter seems to be their self-absorption in the volume of the filter, which is a function of the change of its areal mass and subsequently causes deviations from the identified detection efficiency of the detectors using the double calibrating standard 241 Am (53).
- the module includes the second double standard 241 Am (54), which allows for the irradiation of the tape (always while using the two detectors (46, 47) and subsequent determination of the changes of areal mass of the tape (due to increasing contamination of the tape in its multiple use), based on the measured changes of amplitude spectra, obtained in irradiation.
- the calculation algorithm then allows correcting true detection efficiency to the changes found of the areal mass of the tape.
- the movement of the tape is controlled by the control unit (41) and the alpha spectra measured by means of the couples of detectors (46, 47) are evaluated using the two-track multichannel analyzer 42 (MC 2000).
- the proposed embodiment allows the control of the functions of tape movement, evaluation of the measured alpha spectra, measurement of the background, and the function of the 2 1 Am double standards either manually from the CCU keyboard, or remotely from the user PC.
- the collection duct for samples of active air (49), of appropriate size and length, can open in any suitable way from the cover box of the instrument. Simultaneous detection of radon and thoron DP is secured by the location of the network of measuring instruments at appropriate heights.
- Measuring module D provides for continuous measurement of the input of photon dose equivalent or beta radiation.
- Measurement of dose inputs from the module adequately shielded from the effect of radon and thoron deposits on ground surface, allows for quantitative assessment of atmospheric variability of radon and thoron DP.
- module E consists of amperage or voltage transducers of variables measured a) in the atmosphere, b) in soil and c) in water, which generally serve as auxiliary variables for the interpretation of diffusion of the profile of the key variables and the stability of the atmosphere, measured by modules A to C. All the transducers are connected directly into the analog and digital CCU inputs and are located either on a measuring mast, ca. 10 m high, at the measuring station, or in soil and water according to their use.
- Typical daily gradients ca. 0-3 °C in upward direction represent the lowest stability of the atmosphere with good mixing, explicitly accompanied by typical afternoon (12 to 3pm) daily minima of radon and/or thoron and their short-time DP.
- the gradient is of the opposite trend for the rest of the day, and the atmosphere becomes more stable with less mixing and consequently higher values of atmospheric radon and/or thoron and their DP.
- the choice of solar radiation sensor also matters because it allows direct quantification of the time of insolation.
- the speed and direction of the movement of air masses (or the monitored profile) at defined heights (10 m in this case, the height of the mast) should then be quantified by the transducer of wind speed and trend.
- Fig. 4 shows measuring module C.
- Module C enables simultaneous continuous measurement of bulk radon and thoron activity in soil gas, and the velocity of its exhalation from the surface layer of soil, using ionizing chambers through which air flows, coming from samples collected in ground probes or from exhalation vessel installed on the surface.
- the system is able to display data on the screen of the local control and evaluation unit ERM-4, which, following the concept of the measuring station, sends them on after processing through the central control unit 1 of the station E-log to the user PC.
- Block diagram of the module is shown in Fig. 4.
- each sampling track (67, 68, 69, 70, 71, 72), separated by a delay unit to eliminate thoron from the sampled gas (64, 65, 66), where the sampled gas is also the working gas of the detectors, and each sampling track includes a flow pump (76, 77, 78) with recording (73, 74, 75) of the flow of air along the sampling track;
- Fig. 5 shows block diagram of measuring module F.
- Module F enables
- Module F consists of
- evaluation and control unit (90) which, by means of measuring very low electric currents caused by ionization of air in the sensitive space of the detectors, allows for simultaneous evaluation of bulk radon and thoron activity from the sampling track for collection of soil gas, separated from groundwater resource (81).
- each samplmg track includes a flow pump (89) with recording of the flow of air along the sampling track (88);
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Radiation (AREA)
Abstract
Monitoring and evaluating changes the earth's lithosphere and atmosphere, namely with the intention to predict earthquakes, the following variables selected from group (A) are measured in independent measuring modules at defined time intervals: • - bulk radon and/or thoron activity in the atmosphere; • - bulk activity of decay products of radon and thoron in the atmosphere; • - bulk radon and thoron activity in soil or bedrock; • - areal exhalation of radon from soil, and • - bulk radon activity in water. The results of measurement are transferred to the central unit or computer, a curve of the measurement profile in the time interval is created from the values of the variables measured in the interval, the curves of measurement profiles are aligned into a successive time sequence, and the increments of planes under the curves of measurement profiles of the same time sequence are compared with each other or against the increment of the priority variable measured, selected from group (A).
Description
Method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere
Field of technology
The invention relates to the method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere, namely for prediction of earthquakes.
Present state of art
In spite of the progress made in the last decade in the research of the dynamism of the earth's body, prediction of earthquakes is still a very difficult matter. It is possible to determine the probability of tremors in certain area; it is also possible to determine the year in which expected tremors are likely to occur, thoug precise prediction of the hour or the day resists all efforts.
So far predictions of earthquakes have been based chiefly on analyses of data from past earthquakes, which allow compiling maps of seismic hazards, in which probable earthquake events are charted for defined areas. Three facts are taken into account in the determination of the degree of hazard, i.e. knowledge of historic earthquakes, knowledge of the geological constitution of the area, and knowledge of the situation in terms of geological engineering. The resulting data are plotted in maps, which show either the expected intensity of tremors in the earthquake zones, or probable acceleration of the earth's surface. In order to predict earthquakes, so- called foreshocks, i.e. minor shocks that occur prior the main movements, are also monitored.
Beside predictions based on analysis of known historic earthquakes, other predictions are based on changes of physical fields of the earth, especially electric, magnetic and gravity fields, monitoring of chemical changes of ground waters and gases, or changes of elevation of the earth's surface.
Behaviour of animals is also studied for the prediction of earthquakes, as animals are sensitive to barometric and acoustic phenomena imperceptible to humans, which precede an imminent earthquake.
It has also been found that unusual atmospheric phenomena precede great earthquakes. This is why monitoring stations are installed in earthquake zones, and
data describing events in the upper layer of the atmosphere and in the ionosphere during earthquakes are collected and transmitted via satellite.
Attention has focused recently on the monitoring of seismic waves before, during and after earthquakes. It has been found that the velocity of seismic waves differs by the degree of stress in consequence to the opening and closing of cracks in the earth's crust. The seismic signal transmitted from fault zones between the lithospheric plates precedes the earthquake by several tens of minutes. Mounting stress in the lithosphere, which heralds an imminent earthquake, is detected by measuring the changes of the velocity of seismic waves. The lithosphere is the solid envelope of the earth, consisting of the crust and the uppermost strata of the mantle. The lithosphere is not a compact envelope; it is split to large bodies known as the lithospheric plates, which "float" and collide with each other on top of the plastic layer of the earth's mantle, the asthenosphere.
Existing studies include those focusing on radon as a gas whose occurrence may precede earthquakes (T rique et al, 1999; Igarashi et al, 1995, 1993; Igarashi and Wakita, 1990; Wakita et. al, 1991; Liu et al, 1985; Hauksson and Goddard, 1981; Noguchi and Wakita, 1977; Shapiro et al, 1981; Teng, 1980; Wakita et al.,1980) However, no reliable earthquake prediction based on the content of radon in soil or in air has been published so far. An increase of radon concentration prior to the earthquake proper (typically days to weeks) has been detected in most cases, though in some cases, on the contrary, a drop of radon concentration was detected in some cases.
; Monitoring of instantaneous bulk radon and thoron activity in soil gas is used as part of the monitoring of volcanic activity, e.g. in the Etna area (e.g. S.
Giammanco, Measurements of220Rn and 222Rn and C02 emissions in soil and fumarole gases on Mt. Etna volcano (Italy): Implications for gas transport and shallow ground fracture, 2007).
)
The results of the above monitoring processes do not yield adequate data for precise and reliable prediction of earthquakes. No detailed physical description of relevant dynamic processes in the earth's crust, focusing on earthquake prediction, has been published so far.
General prediction of earthquakes is still hindered by the diversity of conditions preceding the earthquake. Each fault zone has a totally different geological
constitution and different history. The search is still going on for the set of conditions which evidently precede the earthquake.
Document WO 2014/049408 describes the measuring equipment for the prediction and localization of earthquake epicentres by measuring anomalies of radon content in soil, which includes two interconnected measuring modules, of which one is installed in the soil, the other on the surface. The measured values of bulk radon activity are transferred to the server, where data from a number of independent measurement points are compared and processed using an algorithm that compares measured values from a network of ground and surface measurement points to the mean value of bulk radon activity, with the aim of predicting earthquakes and volcanic eruptions as well as calculating the epicentre and magnitude. The algorithm of data processing is not described in detail. The above measuring equipment only monitors and evaluates the time course of radon concentration in soil and on the surface, which does not seem adequate for reliable prediction of the time and place of a possible earthquake.
The object of the invention is to find a new method of monitoring changes in the earth's lithosphere and atmosphere to enable an easier, more reliable and longer- term prediction of earthquakes than what today's state of technology allows.
Substance of the invention
The inadequacies of the present state of art are significantly eliminated, and the object of the invention is achieved by the method of monitoring changes in the earth's lithosphere and atmosphere according to Claim 1. Efficient executions of the method based on the invention are described in the dependent claims 2 to 4.
The inadequacies of the present state of art are significantly eliminated, and the object of the invention is achieved by the equipment formonitoring of changes in the earth's lithosphere and atmosphere according to independent Claim 5.
Efficient executions of the equipment based on the invention are described in dependent Claims 6 to 12.
The method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere are based on the understanding of the phenomena preceding the earthquake and its course in time.
The method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere are based on the understanding that changes preceding earthquakes are manifested in more natural products and more physical quantities than just anomalies of bulk radon activity in soil and on its surface.
The method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere as proposed by the invention are furthermore based on the fact that great stress occurs in the fault between lithospheric plates prior to an earthquake; the stress causes changes in the concentration of soil gases including radon and increased escape thereof; the changes thus provoked depend on the magnitude of the earthquake, its location and time. Radioactivity from soil gases, namely radon and its short-time decay products, ionizes air on a large scale and causes
condensation of water vapours. Temperature and infrared radiation also increase significantly due to this process.
The method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere as proposed by the invention are based on systematic continuous monitoring of important variables and analysis of their changes in time, including the velocity of their increase and decrease, systemic shift of time sequences, and further, more detailed analysis of the sequences using advanced methods of statistical data processing.
One of the seismic side effects before and during an earthquake is a change of the pressure of soil gases, which, geologically speaking, must inevitably contain radon and thoron due to their parent elements uranium and thorium, contained in the earth's crust. Another assumed side effect is the formation of cracks and fractures in the crust, including minor cracks. A change of pressure (increase) is manifested by an increased escape of radon and thoron at "favourable" points, and on the contrary by a reduced occurrence at other points, possibly in relation to the formation of cracks and fractures. In consequence to an increased discharge of radon and thoron, their short-time decay products (DP) also form over time, namely after diffusion of the gases in question to the atmosphere, because the DP formed in soil/rock are captured there.
An increased occurrence namely of radon, trailed in the atmosphere by an increase of their DP, is also manifested by increased count of ions at the site of discharge.
Due to the dilution and spontaneous radioactive decay in the atmosphere further away from the site of the discharge (eruption), the original bulk radon activity (BRA), the initial ion count, and bulk activities (BA) of individual DP gradually decrease, and, last but not least, the initial ratio (BA) of individual DP also changes significantly, as is easy to measure especially for radon, i.e. the ratio (BA) of radio nuclides 218Po:214Pb:214Bi. What changes principally is the most easily measurable ratio (BA) of DP of radon and BRA, known in literature as the coefficient of imbalance F (see Note 1 below). The trajectory of diffusion of atmospheric radon and DP of radon and thoron is influenced namely by the speed and trend of wind, which are to be measured in the station or network of stations together with radon and DP.
Bulk radon activity in soil is understood as radon activity in a unit of volume (1m3) of air collected from pores in soil. Areal velocity of radon exhalation from soil represents the activity of radon diffusing from a unit of area (soil in this case) into the atmosphere per unit of time. The two quantities may or may not be independent of one another and their ratio is another monitored parameter in the prediction of earthquake phenomena.
Whilst the physical half-time of decay of radon is ca. 3.8 days and the half-time of its longest DP 214Pb is ca. 27 minutes, the half-time of decay of thoron is short (ca. 1 min.), but the half-time of its longest DP 214Pb is as much as ca. 11 hours. This implies that beside the detection of radon and its DP, detection of DP of thoron could be quite important at points of expected discharges where geology suggests a content of thorium in soil. Simultaneous measurement of radon and thoron including their DP with totally different half-times therefore allows for sensitive monitoring of different velocities of changes of physical processes. Given the aforesaid longer physical half-time and with regard to the expected speed of air masses (metres to tens of metres per second) on the global scale, the detection of radon and its DP are more important.
The following factors are seen as key elements for the prediction of earthquakes using detection of radon and its DP according to the invention:
a) appropriate location of the station or network of stations in the area of identified geological faults;
b) long-term measurement of the "background", namely of radon and its DP, correlated with the stability of the atmosphere in view of its vertical mixing;
c) long-term measurement of atmospheric background of radon and its DP, as well as its variations in soil and water, correlated to selected soil and water parameters (water temperature, soil moisture and temperature at selected depth profiles) including exhalation of radon from soil.
The most frequent measurement is conducted in the top strata of soils (non- cohesive rock environment). Generally speaking, monitoring can be conducted also in the bedrock, confined or unconfined aquifers, etc.
Variation of other quantities in the atmosphere (input of photon dose equivalent, beta radiation, solar radiation, wind speed and trend, altitudinal temperature gradient, relative humidity, air pressure, and quantity of rainfall) is related to the release of radon and/or thoron, as well as their quantification and identification of the trajectory of their movement in the atmosphere. Their interdependence in time is important for correct earthquake prediction.
Simultaneous measurement of radon and/or thoron is conducted at different altitudes in the atmosphere and at different depths in soil, bedrock and water.
The selection of appropriate monitoring point (tectonically disturbed area, geological fault zone, groundwater resource, etc.). Monitoring at different altitudes and depths makes it possible to track radon from its source in the bedrock to the ground surface, whereby significant differences of half-times of decay of radon and thoron can be helpful in general monitoring of the dynamism of soil gases. Of utmost importance is the fact that comphreensive monitoring of the above profiles can help to offset disturbances causing variation of radon and thoron concentration in time, which cannot be attributed to the phenomena of interest. (E.g. the effect of rainfall and snowfall, freezing of the surface layer of soil, temperature inversion.)
Measurement of atmospheric radon at different altitudes is also important in terms of assessment of the altitude of equivalent mixing layer (EML) of the atmosphere, which is important for the evaluation of stability of the atmosphere with regard to the conditions of dispersion. Generally speaking, knowledge of the stability and variability of the atmosphere is necessary for correct interpretation of the changes of measured radon in relation to earthquake prediction. EML corresponds to the height of air column with assumed homogeneous distribution of radon throughout the volume of the column, and its bulk activity equal to that at the ground level.
Continuous simultaneous measurement of atmospheric (mainly) radon at selected altitudes including its and thoron's DP, radon and/or thoron in water and soil air, areal velocity of radon and/or thoron exalation from soil, values of atmospheric dose inputs of gamma or beta radiation, solar radiation, gradients of temperature and humidity/moisture in air, soil and water, together with the measurement of precipitation, speed and trend of wind, constitute the base measurement profiles.
General comparison of changes in time and dislocations of planes under "the curves of measurement profiles" is a totally new method of earthquake prediction. So far the existing state of technology has not offered any known instrument that would allow to create such a "base profile".
The following conclusions have been made from the measurement and evaluation of variables of group (A) and variables of group (B):
Experience of the measurement of atmospheric radon (module A) and its DP (module B), in agreement with published data, suggests typical minima at midday ca. from 12 to 4pm, and maxima in the early hours of the day (12 to 4am). The variability is of a typical shape though it is disturbed seasonally (spring, summer, autumn, winter), which is probably related to the aforesaid temperature gradient, i.e. atmospheric stability.
Generally speaking, the minima and maxima of bulk radon activity correspond to the minima and maxima of bulk activities of its DP. Whilst the measured values of BRA vary within well-measurable limits ca. 5-40 Bq/m3, equivalent bulk radon activity (EBRA) varies within well-measurable limits ca. 2-20 Bq/m3. Whilst BA of 218Po varies in EBRA maxima on the order of tens of Bq/m3, BA of 214Pb, 214Bi varies on the order of single Bq/m . The well-measurable coefficient of imbalance varies from 0.5 to 0.7 in agreement with published data. Measured activities of atmospheric DP of thoron vary at the margin of measurability on the order of tenths of Bq/m3.
Note 1. EBRA = 0.105 BA (218Po) + 0.516 BA ( 214Pb) + 0.379 BA (214Bi>, where BAs are directly measured bulk activities of radon DP from the filter (filtering tape).
List of illustrations in the drawings
The equipment for the monitoring of changes in the earth's lithosphere and atmosphere according to the invention is illustrated in drawings, which represent the following:
Fig. 1 Block diagram of the equipment according to the invention Obr. 1
Fig. 2 Block diagram of measurement module A
Fig 3 Block diagram of measurement module B
Fig. 4 Block diagram of measurement module C
Fig. 5 Block diagram of measurement module F
Examples of embodiments of the invention
Example 1 :
A modular station for the measurement of natural radioactivity and key disturbing factors allows on the location of interest to determine the characteristic
development in time of the monitored variables in relation to the changes of variables describing ambient effects. Such characteristic changes in time make it possible to eliminate disturbing environmental effects, thus creating a database of filters applicable to the primary time sequences of the variables of interest, which are subsequently analyzed in relation to earthquake prediction.
Primary and filtered time sequences are processed using methods of mathematical statistics, whereby namely long-term trends are compared against the characteristic background of the site and remarkable changes in time (anomalies), which are the subject of interest in relation to changes in the earth's crust. Particularly significant are velocities of increase or decrease of the variables monitored, and the
characteristic shift of time sequences. The modular monitoring system makes it possible at the same time to create calibration protocols for given site using long- term comparison of time sequences of the variables monitored and seismic activity on the site of interest.
Example 2
I. Optimal location of one or more stations in a network for on-line monitoring of the depth/altitude profile of the following continuously measured so-called key variables of interest:
XI = f(t, z) areal velocity of radon and/or thoron exhalation from soil;
X2 = f(t,z) bulk activity of radon and/or thoron from soil gas;
X3 = f(t,z) bulk radon activity in the atmosphere;
X4 = f(t, z) bulk activity of radon and/or thoron in water.
At the same time, the "equipment for the monitoring ..." or network thereof measures auxiliary (explanatory) variables from the following groups:
A) Soil:
- soil temperature profiles at two depths minimum Al(t,z)
- soil moisture profiles at two depths minimum A2 (t,z)
B) Air:
- wind trend Bl (t,z) and speed B2 (t.z) at two altitudes minimum
- absolute air humidity at two altitudes minimum B3(t,z)
- solar radiation B4(t,z)
- air temperature gradient, e.g. at 2 m/10 m B5 (t)
- coefficient of imbalance of radon and its DP in two altitude profiles minimum B6 (t,z)
- equivalent bulk thoron activity measured in two altitude profiles minimum B7(t,z)
- dose input of gamma radiation measured in two altitude profiles minimum B8(t,z)
- total beta radiation measured at two altitudes minimum B9(t,z)
C) Water:
- water temperature in two depth profiles minimum Cl(t,z) t
II. Furthermore, monitoring of the time course of key variables Xl(t,z) - X4(t,z) is conducted for long enough time (seasonally, annually, multiannually) against relevant influencing parameters, i.e. continuously measured variables of groups A, B, C, and their background values Xip (t,z) are fixed without disturbing seismic changes (measured with seismograph), using an adequate model, e.g. regression model with additive or multiplicative approach.
III. The following facts are found in relation to different earthquake phenomena: a) statistically significant deviations (shifts and variations) of the measured time sequences of the key variables Xi(t,z) from their relevant background values Xip (t,z), i.e. monitoring variation of time sequences Xi(t,z) vs Xip(t,z) for i = 1-4;
b) statistically significant correlations, shifts and variations between signal variables X#ij(t,z)pro i≠j, cleared of the background.
IV. A suitable model of earthquake prediction is designed subsequently, using both examples a) and b).
As Fig. 1 shows, the equipment based on the invention consists for example of
- central control unit 1 (CCU);
- emitting (2) and receiving (6) data modem;
- source block of output feeding voltage 12 V dc (3);
- measurement modules minimum A to F;
- transfer SW Wincentral and visualization SW Visualis;
- waterproof box to allow for field deployment of the instrument.
Fig. 1 shows an example of the configuration of the proposed instrument. Analog and digital outputs of measurement from modules A to F are stored in CCU (1) online in its external memory and, either all at a time or selected ones only, transferred in one-minute intervals, using built-in emitting modem (2) and protocol GPRS, to the user PC with a connected receiving modem (6). Transfer parameters including the selection of measured variables can be configured using the program Wincentral, installed in the user PC. Wincentral furthermore allows remote setting and diagnosing of selected parameters of key HW modules A to F from the user PC, using built-in commands. Thus the concept of the program Wincentral makes it possible to control and read the entire network of stations, each of them including a CCU (1) with any combination of modules A to F. Visualization of files of data received in the user PC is provided by the program Visualis, which also allows to export them to common spreadsheets.
CCU 1 is the heart of the entire system; it logs the relevant data from modules A to F directly in the required variables or converts raw measured data to the required variables by means of mathematical algorithms, using its built-in microprocessor. The overall concept of the station includes remote wireless access to all the CCU functions using the program Wincentral, which allows complete control and setting of all the measuring modules from a single PC. The program Wincentral makes it possible for the CCU to set up a network of measuring stations controlled from a single PC, of which each station must include a central control unit.
Power source block can be fed from the grid or battery (4), or solar panels (5).
The functions and connections of measuring modules A to F are as follows:
- Module A offers the possibility of continuous spectrometric and non- spectrometric detection of atmospheric radon and/or thoron. Modul A predstavuje
- Module B offers the possibility of continuous spectrometric measurement of bulk activity of air samples captured on the endless filter, using an appropriate combination of possible couple of detectors (alpha/alpha, alpha/beta,
alpha/gamma) with subsequent determination of bulk activities of individual DP of radon and thoron.
- Module C offers the possibility of simultaneous and continuous measurement of the bulk activity of radon and thoron in soil, and measurement of areal velocity of their exhalation.
- Module D offers the possibility of measurement of photon dose equivalent in the range from natural background up to accident level, or of beta radiation using any kind of measuring instruument, e.g. one based on GM counter and/or coefficient of imbalance F between bulk radon activity and its equivalent bulk activity (EBRA) in the atmosphere.
- Module E includes a set of adequate voltage and amperage sensors, capable of measuring the necessary physical parameters of weather, soil and water (rain gauges, anemometers, thermometers, solar radiation, etc.).
- Module F offers the possibility of simultaneous continuous measurement of bulk radon and throron activity in soil gas, which is sampled from a groundwater resource through a special separation unit with an immersed semi-permeable membrane.
Figure 2 shows possible arrangement of measuring module A for the measurement of atmospheric radon and/or thoron, in the single-detector option for continuous collection of samples of atmospheric radon and/or thoron, based on non- spectrometric detection with a scintillation chamber with ZnS-based
photomultiplier. The sample of active air, stripped of primary DP by filter (18), is pumped by pump (19) through drier (31) and detector (12) with arbitrarily fixed and known flow velocity, measured by airflow gauge (14). During its radioactive decay, radon and/or thoron then causes scintillations within the sensitive volume of the chamber. These are registered by the photomultiplier of the scintillation chamber and taken by the counter as time pulses. The pulse count obtained per unit of time represents the rate of bulk activity of radon and/or thoron. The CCU uses
published algorithms to calculate radon and/or thoron values. The results of the measurement are saved in CCU memory and at the same time transferred through the built-in modem online at one-minute intervals to the user PC by means of GPRS protocol.
An advantage of the measuring module is that it offers the choice of manual or remote setting of the detector's own background using radon-free technical gas, contained in pressure cylinder (28), with the aid of gas management system (20, 22, 23, 24). Menu of the evaluation unit (11) can be used for manual setting, the evaluation unit consistsing of a two-track multichannel analyzer (42) (MCA, MC 2000). The remote option is provided by the CCU on command from the user PC.
The two-track multichannel analyzer 42 (MCA) furthermore allows another simultaneous detection of radon and/or thoron by connecting an additional, spectrometric or non-spectrometric detector of alpha, beta or gamma radiation. The stability of the detector in use is achieved by means of working standard 241 Am (29), which is inserted into the detection space of the detector when needed, and is controlled from the user PC by means of an electro-magnetic valve.
To allow the use of module A at low temperature, its critical parts, i.e.
electromagnetic valves (15, 20) are insulated and tempered using adjustable thermostats (21).
Temperature sensor (13) helps to offset the effect of temperature variation, which affects namely the response of detector photomultiplier. Exchangeable filter (31) capturing air humidity is installed to eliminate the negative effect of condensation of air humidity inside the detector, and to improve the efficiency of detection and spectrometry of the measuring instruments; furthermore, optional heating of the adequately insulated detector is installed to react to outside air humidity, using the adjustable air humidity sensor of the hygrostat (30). Samples can be collected from a grommet of any size and shape, fixed in any suitable way and opening on the outside of the casing at a suitable collecting height above ground. Simultaneous measurement of the height profile of atmospheric radon and/or thoron is possible using an additional detector or setting up a network of measuring stations.
Figure 3 shows a possible arrangement of measuring module B to measure samples of atmospheric bulk activities of DP of radon and thoron, captured on the endless filter (filtering tape).
DP of radon and thoron in the atmosphere are fixed to aerosol. This is why they are detected applying the well-known principle of collection on filter during the suction of a defined volume of air sample at a known flow velocity through sensor (43) and filtering tape (50) using a powerful pump (44). The bulk activity of DP is then calculated from the known volume of air transiting the filter per unit of time and their activity determined from the filter using a spectroscopic semiconductor detector. Alpha radiation emanating from the DP is used for the detection, as it allows registering a sufficient count of signal pulses. In the proposed embodiment, the sample collected at filter is first subjected to the first detector (46), by whose means and using adequate timing of the measurement (choice of measurement times and intervals between them) the bulk activity of individual short-time radon DP (218 Po, 214 Pb, 214 Bi) is determined. The sample is then subjected to the second detector (47) for analogous determination, again with adequate timing, of the bulk activity of short-time thoron DP (212Pb, 212Bi). When thoron DP measurement is finished, air samples are collected again from the "endless" filter (filter coil continuously winds off and then winds back to the beginning instantaneously) and the entire measurement cycle is repeated. The proposed algorithm provides for both spectrometric and non-spectrometric determination of the bulk activity of individual DP of radon and thoron, the necessary measurement of the background of the tape, simultaneous calibration of both detectors, and irradiation of the tape by means of the couple of double standards 41Am (53, 54). A certain disadvantage of the use of detection of alpha radiation from the activity of the DP captured on the filter seems to be their self-absorption in the volume of the filter, which is a function of the change of its areal mass and subsequently causes deviations from the identified detection efficiency of the detectors using the double calibrating standard 241 Am (53). The module includes the second double standard 241Am (54), which allows for the irradiation of the tape (always while using the two detectors (46, 47) and subsequent determination of the changes of areal mass of the tape (due to increasing contamination of the tape in its multiple use), based on the measured changes of amplitude spectra, obtained in irradiation. The calculation algorithm then allows correcting true detection efficiency to the changes found of the areal mass of the tape.
The movement of the tape is controlled by the control unit (41) and the alpha spectra measured by means of the couples of detectors (46, 47) are evaluated using the two-track multichannel analyzer 42 (MC 2000). The proposed embodiment allows the control of the functions of tape movement, evaluation of the measured alpha spectra, measurement of the background, and the function of the 2 1 Am
double standards either manually from the CCU keyboard, or remotely from the user PC.
The collection duct for samples of active air (49), of appropriate size and length, can open in any suitable way from the cover box of the instrument. Simultaneous detection of radon and thoron DP is secured by the location of the network of measuring instruments at appropriate heights.
Module D
Measuring module D provides for continuous measurement of the input of photon dose equivalent or beta radiation.
Measurement of dose inputs from the module, adequately shielded from the effect of radon and thoron deposits on ground surface, allows for quantitative assessment of atmospheric variability of radon and thoron DP.
Module E
In principle, module E consists of amperage or voltage transducers of variables measured a) in the atmosphere, b) in soil and c) in water, which generally serve as auxiliary variables for the interpretation of diffusion of the profile of the key variables and the stability of the atmosphere, measured by modules A to C. All the transducers are connected directly into the analog and digital CCU inputs and are located either on a measuring mast, ca. 10 m high, at the measuring station, or in soil and water according to their use.
In the atmosphere. With regard to the interpretation of the diffusion of the monitored profile of the measured variables, the most important one is temperature difference (gradient = vector) between the ground collection layer of air, e.g. at 2 m, and that at 10 m. Its absolute value and trend provide for an important indicator of the stability of ground atmosphere. Typical daily gradients ca. 0-3 °C in upward direction represent the lowest stability of the atmosphere with good mixing, explicitly accompanied by typical afternoon (12 to 3pm) daily minima of radon and/or thoron and their short-time DP. The gradient is of the opposite trend for the rest of the day, and the atmosphere becomes more stable with less mixing and consequently higher values of atmospheric radon and/or thoron and their DP. As the gradient is primarily influenced by the intensity of insolation, the choice of solar radiation sensor also matters because it allows direct quantification of the time of insolation. The speed and direction of the movement of air masses (or the
monitored profile) at defined heights (10 m in this case, the height of the mast) should then be quantified by the transducer of wind speed and trend.
In water and soil.
Fig. 4 shows measuring module C. Module C enables simultaneous continuous measurement of bulk radon and thoron activity in soil gas, and the velocity of its exhalation from the surface layer of soil, using ionizing chambers through which air flows, coming from samples collected in ground probes or from exhalation vessel installed on the surface. The system is able to display data on the screen of the local control and evaluation unit ERM-4, which, following the concept of the measuring station, sends them on after processing through the central control unit 1 of the station E-log to the user PC.
Block diagram of the module is shown in Fig. 4.
Module C consists of
- evaluation and control unit (79), which, by means of measuring very low electric currents caused by ionization of air in the sensitive space of the detector, allows for simultaneous evaluation of bulk radon and thoron activity from individual sampling tracks of air gas collected at different depths (61, 62) and exhalation vessels installed on ground surface (63);
- two flow detectors, connected in series, for each sampling track (67, 68, 69, 70, 71, 72), separated by a delay unit to eliminate thoron from the sampled gas (64, 65, 66), where the sampled gas is also the working gas of the detectors, and each sampling track includes a flow pump (76, 77, 78) with recording (73, 74, 75) of the flow of air along the sampling track;
- closed sampling track with exhalation vessel, with specially designed collection hoses to prevent freezing of water vapour along the sampling track and
- open sampling track installed at defined depth in bedrock with specially designed i collection hoses to prevent freezing of water vapour along the sampling track.
Fig. 5 shows block diagram of measuring module F. Module F enables
simultaneous continuous measurement of bulk radon and thoron activity in soil gas, sampled from a groundwater resource through a special separation unit with ; semi-permeable membrane. The system is able to display data on the screen of its local control and evaluation unit ERM-4, which, following the concept of the
measuring station, sends them on after processing through the CCU E-log to the user PC.
Module F consists of
a) evaluation and control unit (90), which, by means of measuring very low electric currents caused by ionization of air in the sensitive space of the detectors, allows for simultaneous evaluation of bulk radon and thoron activity from the sampling track for collection of soil gas, separated from groundwater resource (81).
b) two flow detectors, connected in series, for each sampling track (85, 87), separated by a delay unit (86) to eliminate thoron from the sampled gas, where the sampled gas is also the working gas of the detectors, and each samplmg track includes a flow pump (89) with recording of the flow of air along the sampling track (88);
c) closed sampling track with special separation unit (84), providing for radon transfer from water through semi-permeable membrane (83) to the sampling track; d) separation vessel with special semi-permeable membrane of a defined length and diameter, enabling one-off or continuous collection of water from a groundwater resource.
List of reference characters
A measuring module A
B measuring module B
C measuring module C
D measuring module D
E measuring module E
F measuring module F
1 central control unit
2 emitting modem (GPRS)
3 source block
4 batteries
5 solar cable
6 PC 6 with receiving module
11 monitoring and evaluation unit
12 scintillation chamber 12 for radon detection
13 photomultiplier temperature sensor 13
14 air flow gauge
15 tempered electromagnetic valve
16 exit valve of scintillation chamber
17 entry valve of scintillation chamber
18 filter
19 pump
20 tempered electromagnetic valve 20 with a nozzle
21 Fluid thermos regulator 21
22 output closure 22 of reduction valve
23 reduction valve
24 pressure sensor
25 low pressure gauge
] 26 pressure cylinder manometer 26
27 pressure cylinder closure 27
28 pressure cylinder 28 (background gas)
29 working standard 29
30 hygrostat
; 31 drying agent
41 control unit
42 two-track multichannel analyser
air flow sensor
pump
filtering tape winding coil 45
, 47 detectors 46 (1) and 47 (2)
mechanical travel block 48 of filtering tape suction neck 49 for air sample collection filtering tape
spare coil 51 of filtering tape
preamplifiers 52 of detectors 1 ,2
couple of working standards 53 of stability of detectors 1, 2
couple of working standards 54 of tape contamination , 62 collection depths
exhalation vessels
, 65, 66 delay unit
-72 sampling track
, 74, 45 record 73, 74, 75 of air flow along sampling track, 77, 78 flow pump
evaluation and control unit
ground water source
record 82 of air flow along sampling track semi-permeable membrane
separation unit
, 87 sampling track
delay unit
record of air flow 88 along sampling track flow pump
evaluation and control unit
Claims
1. Method of monitoring and evaluating changes in the earth's lithosphere and atmosphere, particularly for earthquake forecast, according to which independent measurement of variables is conducted in independent measuring modules at selected time intervals:
- bulk activity of radon and/or thoron and/or their short live decay products in atmosphere and
- bulk activity of radon and thoron in lithosphere and/or
- bulk activity of radon in lithospheric waters,
measured results are transferred to a central unit or computer;
a curve of profile of the measurement in a given time interval is created from values of the selected variables, measured in the same interval;
curves of measurement profiles are arranged in a progressive time sequence;
increments of planes under the measurement profile curves of the same time sequence are compared with each other or against the increment of a priority variable, selected from the measured variables.
2. Method of monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 1,
characterized by that
the variables selected from following group are measured in the independent measuring modules at defined time intervals:
- input of photon dose equivalent and/or beta radiation in atmosphere;
- coefficient of imbalance F between bulk radon activity and equivalent bulk radon activity (EBRA) in atmosphere;
- areal velocity of radon exhalation from soil;
- solar radiation;
- wind speed and trend;
- temperature, relative humidity and ambient pressure.
3. Method of monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 1 ,
characterized by that
- bulk activity of products of radon and thoron in atmosphere is measured simultaneously at different altitudes;
- bulk activity of radon and thoron in the lithosphere is measured simultaneously at different depths;
- radon bulk activity in lithosphenc waters is measured simultaneously at different depths.
4. Method of monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claims 2 and 3,
characterized by that
- temperature and ambient pressure are measured simultaneously at different altitudes in the atmosphere and at different depths in lithosphenc waters, and in the lithosphere.
5. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to any of Claims 1 to 4,
characterized by that
it includes independent measuring modules (A, B, C, F) for measurement of variables of natural radioactivity
- module (A) for measurement of bulk activity of radon and/or thoron in the atmosphere;
- module (B) for measurement of bulk activity of the short live decay products of radon and thoron in the atmosphere;
- module (C) for measurement of bulk radon and thoron activity in lithosphere and rates of exhalation of radon from soil per area, and
- module (F) for measurement of bulk radon activity in lithospheric waters
- central unit (1) or computer with control and evaluation software to synchronize measuring intervals and evaluate the measured values of the selected variables.
6. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 5,
characterized by
including independent module (D) for measurement of the input of photon dose equivalent or beta radiation in the atmosphere.
7. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 5 or 6,
characterized by
including independent module (E) for measurement of at least one of the variables: - solar radiation - speed and trend of wind - temperature, relative humidity and ambient pressure.
8. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 5,
characterized by that
the module (A) is comprised of
I. evaluation unit (11), consisting of a two-track multichannel analyzer for simultaneous evaluation of bulk activity of atmospheric radon and/or thoron from one or two alpha and/or gamma, alpha and/or beta spectrometric or non- spectrometric measuring detectors;
II. detection measurement unit (12), consisting of one or two alpha/alpha or alpha gamma or alpha beta spectrometric or non-spectrometric detectors;,
III. collecting part for samples of atmospheric radon and/or thoron, consisting of a pump (19) and a system of tubes and valves for continuous collection of atmospheric samples;
IV. gas management system including a pressure cylinder (29) with inactive gas for the determination of the proper background of the detector or couple thereof;
V. a block (3) of power sources and central control, log and transfer unit (1.
9. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 5,
characterized by that
module (B) consists of
I. evaluation unit (41) which is a two-track multichannel analyzer (42) allowing simultaneous evaluation of bulk activity of short live decay products of radon and thoron from atmospheric samples, captured on endless filter with the advantage of a couple of measuring detectors;
II. detection measurement unit (46, 47), consisting of a couple of spectrometric semiconductor detectors with preamplifiers (52) for analysis using alpha/alpha and/or alpha/beta or alpha/gamma semiconductor spectrometry;
III. collecting part for samples of atmospheric short live decay products of radon and thoron on infinite filter, consisting of suction neck (49) and mechanical elements (48) for the control of filtering tape movement (50) with control encoder;
IV. block (3) of power sources and central control, log and transfer unit (1, with the advantage of being shared by all the modules of the station.
10. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 9,
characterized by that
the equipment includes a couple of double standards of alpha radiation, e.g. 241 Am, whereby one of the couple of double standards, modified for travel in the space between the filtering tape and the detector, is designed for simultaneous calibration and verification of the function of both detectors in use, and
the other double standard, modified for rotary motion behind the filtering tape above the two detectors, is designed to determine changes of mass per area of the filtering tape, caused by progressive contamination, and to conduct continuous correction of the efficiency of the two detectors in the monitoring of the change of mass per area.
11. Equipment for monitoring and evaluating changes in the earth's lithosphere and atmosphere according to Claim 5,
characterized by that
module (C) consists of
I. evaluation and measurement unit (79) for simultaneous evaluation of bulk activity of radon and thoron from individual sampling tracks for collection of soil gas at different collection depths and in exhalation vessels installed on ground surface;
II. two flow detectors (67, 68, 69, 70, 71, 72), connected in series, separated for each sampling track by a delay unit (64, 65, 66) for elimination of thoron in the gas sample, whereas the gas sample is also the working gas of the detectors, and a flow pump with flow control is provided for each sampling track;
III. closed sampling track with exhalation vessel (63) with collecting hoses that prevent water vapor freezing on the collection track;
IV. open sampling track with collecting probe (61, 62) for soil gas, arranged at a ; defined depth in the bedrock with collecting hoses that prevent freezing of water vapor in the collection track.
12. Equipment for monitoring and evaluating changes in the earth's litho sphere and atmosphere according to Claim 5,
characterized by that
module (F) consists of
I. evaluation and control unit (90) for simultaneous evaluation of bulk radon and thoron activity in the sampling track for the collection of air, separated from a ground water body;
II. two flow detectors, connected in series for each sampling track (85, 87) , separated by a delay unit (86) for elimination of thoron in the gas sample, whereby the gas sample is also working gas for the detectors, and each sampling track is furnished with a flow pump (89) with flow control (88);
III. a closed sampling track with a separation unit (84) for radon transfer from water through a semi-permeable membrane (83) to the sampling track and
IV. a separation vessel with special semi-permeable membrane (83) of a defined length and cross-section for one-off or continuous collection of water from a ground resource (81).
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CZPV2014-449 | 2014-06-29 | ||
CZ2014-449A CZ2014449A3 (en) | 2014-06-29 | 2014-06-29 | Method of monitoring changes in earth lithosphere and atmosphere and apparatus for making the same |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2016000666A1 true WO2016000666A1 (en) | 2016-01-07 |
Family
ID=53512919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CZ2015/000062 WO2016000666A1 (en) | 2014-06-29 | 2015-06-17 | Method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere |
Country Status (2)
Country | Link |
---|---|
CZ (2) | CZ28361U1 (en) |
WO (1) | WO2016000666A1 (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785428A (en) * | 2016-03-10 | 2016-07-20 | 吴清荣 | Underground radon gas monitoring method and system for seismic precursor monitoring |
WO2017120392A1 (en) * | 2016-01-08 | 2017-07-13 | The Regents Of The University Of Michigan | Ice crystals and volcanic ash detection system |
CN108646286A (en) * | 2018-06-01 | 2018-10-12 | 衡阳师范学院 | A kind of reference unit and method of continuous reliable regulation and control precipitation rate of radon |
CN109520892A (en) * | 2019-01-22 | 2019-03-26 | 南华大学 | A kind of electric field-temperature-wind pressure coupled field radon and its daughter diffusion And Movement device |
US10508952B1 (en) | 2018-10-31 | 2019-12-17 | The Regents Of The University Of Michigan | Optimum spectral bands for active vision systems |
US10621865B2 (en) | 2018-03-29 | 2020-04-14 | The Regents Of The University Of Michigan | Road condition monitoring system |
CN112882124A (en) * | 2021-01-13 | 2021-06-01 | 核工业北京地质研究院 | Method for delineating beneficial parts of three-dimensional uranium mineralization |
IT202000006313A1 (en) * | 2020-03-25 | 2021-09-25 | Univ Degli Studi Di Salerno | System for monitoring and predicting concentrations of radioactive gases and related method. |
CN113627023A (en) * | 2021-08-13 | 2021-11-09 | 东北石油大学 | Method and system for calculating self-absorption amount of different horizontal bedding shales |
CN113777660A (en) * | 2021-09-15 | 2021-12-10 | 核工业航测遥感中心 | Method for solving atmospheric radon correction parameter by energy spectrum ratio method |
CN114088313A (en) * | 2021-11-13 | 2022-02-25 | 衡阳师范学院 | Method for continuously measuring radon exhalation rate by changing leakage coefficient and anti-diffusion coefficient |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5654899A (en) * | 1993-10-04 | 1997-08-05 | Leonard; Bobby E. | Radiation instrument and method to measure radon and thoron concentration, air change rate and source emanation rate magnitudes |
US20140077074A1 (en) * | 2012-09-17 | 2014-03-20 | New York University | Seismic prediction with decay products |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08201523A (en) * | 1991-05-16 | 1996-08-09 | Kagaku Gijutsucho Hoshasen Igaku Sogo Kenkyusho | Installation-type radon and thoron discrimination and concentration measuring apparatus |
JPH11295435A (en) * | 1998-04-15 | 1999-10-29 | Fujita Corp | Measure of concentration of radon in underground water and measurement system |
ITRM20020638A1 (en) * | 2002-12-20 | 2004-06-21 | Gioacchino Giuliani | DEVICE AND METHOD OF DETECTION OF THE CONCENTRATION VARIATION OF GAS RADON IN THE ENVIRONMENT AND USE OF THE DEVICE IN THE FORECAST OF SEISMIC EVENTS. |
ES2425801B1 (en) * | 2012-03-16 | 2014-08-14 | Universitat Rovira I Virgili | IDENTIFICATION STATION AND REAL-TIME MEASUREMENT OF GAMMA ENVIRONMENTAL RADIOACTIVITY THROUGH SPECTROMETRY ON PAPER FILTER |
FR2996010B1 (en) * | 2012-09-27 | 2016-01-01 | Gregory Jean | DEVICE FOR PREDICTING SEISM AND VOLCANIC ERUPTIONS, LOCATION OF EPICENTERS AND MAGNITUDES BY MEASURING RADON IN THE SOIL |
-
2014
- 2014-06-29 CZ CZ2015-30994U patent/CZ28361U1/en not_active IP Right Cessation
- 2014-06-29 CZ CZ2014-449A patent/CZ2014449A3/en unknown
-
2015
- 2015-06-17 WO PCT/CZ2015/000062 patent/WO2016000666A1/en active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5654899A (en) * | 1993-10-04 | 1997-08-05 | Leonard; Bobby E. | Radiation instrument and method to measure radon and thoron concentration, air change rate and source emanation rate magnitudes |
US20140077074A1 (en) * | 2012-09-17 | 2014-03-20 | New York University | Seismic prediction with decay products |
Non-Patent Citations (3)
Title |
---|
HARI PRASAD JAISHI ET AL: "Radon and thoron anomalies along Mat fault in Mizoram, India", JOURNAL OF EARTH SYSTEM SCIENCE, vol. 122, no. 6, 1 December 2013 (2013-12-01), India, pages 1507 - 1513, XP055222767, ISSN: 0253-4126, DOI: 10.1007/s12040-013-0361-4 * |
MELVIN J D ET AL: "AN AUTOMATED RADON-THORON MONITOR FOR EARTHQUAKE PREDICTION RESEARCH", NUCLEAR INSTRUMENTS & METHODS,, vol. 153, no. 1, 1 July 1978 (1978-07-01), pages 239 - 251, XP001346040 * |
W. D. CROZIER ET AL: "Measurements of radon 220 (thoron) in the atmosphere below 50 centimeters", JOURNAL OF GEOPHYSICAL RESEARCH, vol. 71, no. 20, 15 October 1966 (1966-10-15), US, pages 4735 - 4741, XP055222771, ISSN: 0148-0227, DOI: 10.1029/JZ071i020p04735 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017120392A1 (en) * | 2016-01-08 | 2017-07-13 | The Regents Of The University Of Michigan | Ice crystals and volcanic ash detection system |
US10336465B2 (en) | 2016-01-08 | 2019-07-02 | The Regents Of The University Of Michigan | Ice crystals and volcanic ash detection system |
CN105785428A (en) * | 2016-03-10 | 2016-07-20 | 吴清荣 | Underground radon gas monitoring method and system for seismic precursor monitoring |
US10621865B2 (en) | 2018-03-29 | 2020-04-14 | The Regents Of The University Of Michigan | Road condition monitoring system |
CN108646286A (en) * | 2018-06-01 | 2018-10-12 | 衡阳师范学院 | A kind of reference unit and method of continuous reliable regulation and control precipitation rate of radon |
US10508952B1 (en) | 2018-10-31 | 2019-12-17 | The Regents Of The University Of Michigan | Optimum spectral bands for active vision systems |
CN109520892A (en) * | 2019-01-22 | 2019-03-26 | 南华大学 | A kind of electric field-temperature-wind pressure coupled field radon and its daughter diffusion And Movement device |
CN109520892B (en) * | 2019-01-22 | 2024-01-30 | 南华大学 | Electric field-temperature-wind pressure coupling field radon and daughter diffusion migration device thereof |
IT202000006313A1 (en) * | 2020-03-25 | 2021-09-25 | Univ Degli Studi Di Salerno | System for monitoring and predicting concentrations of radioactive gases and related method. |
CN112882124A (en) * | 2021-01-13 | 2021-06-01 | 核工业北京地质研究院 | Method for delineating beneficial parts of three-dimensional uranium mineralization |
CN112882124B (en) * | 2021-01-13 | 2024-02-20 | 核工业北京地质研究院 | Method for delineating beneficial parts of three-dimensional uranium ore formation |
CN113627023A (en) * | 2021-08-13 | 2021-11-09 | 东北石油大学 | Method and system for calculating self-absorption amount of different horizontal bedding shales |
CN113627023B (en) * | 2021-08-13 | 2023-08-22 | 东北石油大学 | Self-absorption amount calculation method and system for different horizontal layers of shale |
CN113777660A (en) * | 2021-09-15 | 2021-12-10 | 核工业航测遥感中心 | Method for solving atmospheric radon correction parameter by energy spectrum ratio method |
CN113777660B (en) * | 2021-09-15 | 2023-10-20 | 核工业航测遥感中心 | Method for solving atmospheric radon correction parameters by energy spectrum ratio method |
CN114088313A (en) * | 2021-11-13 | 2022-02-25 | 衡阳师范学院 | Method for continuously measuring radon exhalation rate by changing leakage coefficient and anti-diffusion coefficient |
CN114088313B (en) * | 2021-11-13 | 2024-05-07 | 衡阳师范学院 | Method for continuously measuring radon exhalation rate by changing leakage coefficient and back diffusion coefficient |
Also Published As
Publication number | Publication date |
---|---|
CZ28361U1 (en) | 2015-06-23 |
CZ305772B6 (en) | 2016-03-09 |
CZ2014449A3 (en) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2016000666A1 (en) | Method and equipment for the monitoring of changes in the earth's lithosphere and atmosphere | |
Zafrir et al. | Differentiation between the effect of temperature and pressure on radon within the subsurface geological media | |
Singh et al. | Radon monitoring in soil gas and ground water for earthquake prediction studies in North West Himalayas, India | |
Andreasen et al. | Modeling cosmic ray neutron field measurements | |
Barkat et al. | Time series analysis of soil radon in Northern Pakistan: Implications for earthquake forecasting | |
US6663012B2 (en) | Method for detecting moisture in soils using secondary cosmic radiation | |
US8712700B2 (en) | Method and apparatus for detection of the remote origin fraction of radon present in a measuring site | |
Cannelli et al. | Italian radon monitoring network (Iron): A permanent network for near real-time monitoring of soil radon emission in Italy | |
Piersanti et al. | Long term continuous radon monitoring in a seismically active area | |
Pérez et al. | Precursory subsurface 222 Rn and 220 Rn degassing signatures of the 2004 seismic crisis at Tenerife, Canary Islands | |
Zoran et al. | Radon levels assessment in relation with seismic events in Vrancea region | |
Ván et al. | Long term measurements from the Mátra Gravitational and Geophysical Laboratory | |
Mao et al. | The temporal variation of radon concentration at different depths of soil: A case study in Beijing | |
Amestoy et al. | Effects of environmental factors on the monitoring of environmental radioactivity by airborne gamma-ray spectrometry | |
Benkovitz et al. | The dynamics of Rn-222 cyclic flow within the shallow geological subsurface media as a daily temporal variated source for exhalation into the air | |
Galeriu et al. | Radon, as a tracer for mixing height dynamics-an overview and RADO perspectives | |
De Cicco et al. | Multi-parametric approach to the analysis of soil radon gas for its validation as geoindicator in two sites of the Phlegrean Fields caldera (Italy) | |
Baskaran et al. | Radon measurement techniques | |
Zafrir et al. | Response of Radon in a seismic calibration explosion, Israel | |
Chauhan et al. | Multiple linear regression analysis to estimate hydrological effects in soil Rn-222 at Ghuttu, Garhwal Himalaya, India: a prerequisite to identify earthquake precursors | |
Anisimov et al. | Background Volumetric Activity of Radon and Thoron in Soil and Surface Atmosphere Observed by a Radon-Recording Seismic Station | |
Yadav et al. | Imprint of Diurnal and Semidiurnal Cyclicity in Radon Time Series of MPGO, Ghuttu Garhwal Himalaya: Evidence Based on Singular Spectrum Analysis | |
Gavrilyuk et al. | Monitoring the 222 Rn concentration in the air of low-background laboratories by means of an ion-pulse ionization chamber | |
Anisimov et al. | Electricity of the Undisturbed Atmospheric Boundary Layer of Middle Latitudes | |
Gaware et al. | Indigenous development of online radon and thoron monitors for applications in Uranium mining and Thorium processing facilities |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 15744866 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 15744866 Country of ref document: EP Kind code of ref document: A1 |