JP2015172494A - Environment load molecule source evaluation method, environment load molecule source evaluation system, and computer program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To correctly identify an environment load molecule source at an environment load molecule occurring at an environment load molecule detection part and affecting an environment load molecule as quickly as possible.SOLUTION: An environment load molecule generation source evaluation method is configured to: solve a simulation calculation model expression of a numeral fluid dynamics by letting concentration and environment condition of each molecule included in gas occurring in each environment load molecule generation source measured by a vacuum ultraviolet-single photon ionization-mass analysis method as a boundary condition to thereby derive a distribution of the concentration of each molecule in an analysis area on the basis of various boundary conditions; store the concentration of each molecule at each position obtained as a result of the derivation, and the concentration of each molecule and environment condition at each environment load molecule generation source when the result thereof is obtained as a simulation result; measure the concentration of each molecule at an environment load molecule detection part by the vacuum ultraviolet-single photon ionization-mass analysis method; and extract the concentration of each molecule measured at the environment load molecule detection part and the concentration of each molecule at each environment load molecule generation source corresponding to the environment condition from the simulation result.

Description

  The present invention relates to an environmental load molecule generation source evaluation method, an environmental load molecule generation source evaluation system, and a computer program, and particularly affects the environmental load molecule from the measurement result of gas molecules at a location where the environmental load molecule is detected. It is suitable for use in deriving an index for identifying the environmental load molecule generation source.

Environmentally damaging molecules may be generated in residential areas by gas generated from factory chimneys. In such a case, it is conceivable to measure the concentration of molecules in the atmosphere by a GC / MS (Gas Chromatography / Mass Spectrometry) analyzer at a location where environmental load molecules are generated such as in a residential area. In the following description, environmental impact molecules include NOx, SOx, benzene, and the like, as well as odors generated when various molecules are mixed in a complementary manner. In the following, odors are described mainly among the environmental load molecules. The “odor generation source such as a factory chimney” is referred to as “odor generation source” as necessary, and the “location where odor is detected” is referred to as “odor detection location” as necessary.
However, the GC / MS analyzer requires a measurement time of about 10 minutes. Odor changes from moment to moment depending on the wind direction, air volume, and the like. Therefore, it is difficult to apply the GC / MS analyzer to odor measurement.

Therefore, among the environmental load molecules, there are techniques described in Patent Documents 1 and 2 and Non-Patent Document 1 as techniques for measuring odor particularly.
In Patent Document 1, first, a reference odor vector S1 is created in a four-dimensional space (odor space) formed by four odor sensors. In addition, a measurement odor vector SX is created in the odor space from the measurement sample. Based on the angle θ between the reference odor vector S1 and the measurement odor vector SX, the similarity α between the coke odor of the reference sample and the odor of the measurement sample is calculated. Further, the orthographic projection of the measured odor vector SX with respect to the reference odor vector S1 is taken, and the intensity T of the coke odor contained in the odor of the measurement sample is calculated based on the length L of the orthographic projection.

In Patent Document 2, first, M sets of sample data strings are calculated based on N output values measured by N odor sensors for each of M measurement samples, and an odor group is associated with each sample data string. Further, the measurement data string is calculated based on N output values obtained by measuring the measurement object. Then, an odor group associated with one set of sample data strings having a good correlation with the measurement data string among the M sets of sample data strings is specified.
Non-Patent Document 1 discloses that when downwash occurs due to the positional relationship between an odor generation source and a building, the concentration of each molecule is derived by a diffusion model considering downwash.

JP 2009-156769 A JP 2007-218704 A JP 2011-252893 A

Industrial Environment Management Association website, "Ministry of Economy, Trade and Industry Low Rise Industrial Source Dispersion Model METI-LIS Model ver.3.02 Instruction Manual", March 2012 30th of month Yoshio Nagata, Norifumi Takeuchi, “Threshold Measurement Results of Odorous Substances Using the Three-Point Comparison Odor Bag Method”, Proceedings of the 29th Annual Meeting of the Air Pollution Society, p. 528, 1988 Yasuhisa Sakai, et al., “Odor threshold measurement of hydrocarbon components in automobile exhaust gas”, 28th Annual Meeting of the Air Pollution Society, p. 448, 1987

By the way, there are usually several candidates for the odor generation source described above, so if an odor is detected at the odor detection location, identify which odor generation source affects the odor detected at the odor detection location. However, it is necessary to take measures to suppress the generation of odors.
However, in the techniques described in Patent Documents 1 and 2, the odor is specified using the measurement result of the gas sensor as an index. Therefore, the technologies described in Patent Documents 1 and 2 cannot identify the odor-causing molecule, and as a result, identify the odor source that affects the odor generated at the odor detection location. I can't. In Non-Patent Document 1, although the odor-causing molecule can be identified, the concentration of each molecule is only derived from the diffusion equation, so it is difficult to accurately consider the influence of the wind. It is.

  The present invention has been made in view of the above problems, and it is possible to identify the environmental load molecule generation source that has an influence on the environmental load molecule detected at the environmental load molecule detection site as accurately as possible. Although it aims, it aims at pinpointing the odor generation source which has an influence on the odor detected especially in the odor detection location as much as possible. In the following description, odors are particularly described among the environmental load molecules, but the present invention is not limited to odors and is a technology applicable to all environmental load molecules.

  The environmental load molecule generation source evaluation method of the present invention is an environmental load molecule generation source evaluation method for deriving an index for identifying an environmental load molecule generation source that affects an environmental load molecule detected at an environmental load molecule detection location. By irradiating a vacuum ultraviolet laser beam, a plurality of molecules contained in the gas are ionized at once with a single photon, and the mass spectrum of each molecule is measured. Based on the first environmental load molecule measurement step and the mass and spectrum intensity of each molecule measured in the first environmental load molecule measurement step, the concentration of each molecule in the plurality of environmental load molecule generation sources A first molecular concentration deriving step for deriving the concentration, and a concentration and an environmental condition of each molecule in the plurality of environmental load molecule generation sources derived by the first molecular concentration deriving step. Deriving molecular concentration distribution, which is the distribution of the concentration of each molecule in the analysis region, which is a predetermined three-dimensional region, by solving the simulation calculation model formula of computational fluid dynamics as the boundary condition, with different boundary conditions A first molecular concentration distribution deriving step to be performed; a concentration of each molecule at each position in the analysis region obtained by the molecular concentration distribution derived by the first molecular concentration distribution deriving step; and the molecular concentration The concentration of each molecule in the plurality of environmental load molecule generation sources when the distribution is derived and at least a part of the environmental conditions when the molecular concentration distribution is derived are stored in a storage medium as a simulation result. Irradiating vacuum ultraviolet laser light and ionizing a plurality of molecules contained in the gas with a single photon. The environmental load molecule detection point is determined based on the second environmental load molecule measurement step performed at the environmental load molecule detection point and the mass of each molecule measured in the second environmental load molecule measurement step. A second molecular concentration deriving step for deriving the concentration of each molecule in the method, a concentration of each molecule at the environmental load molecule detection site derived by the second molecule concentration deriving step, and the environmental load molecule detection site; Extraction of the concentration of each molecule at the plurality of environmental load molecule generation sources corresponding to the environmental condition that can be regarded as an environmental condition when the environmental load molecule is detected at the detection site of the environmental load molecule from the result of the simulation And the environmental load detected at the environmental load molecule detection location based on the concentration of each molecule at the plurality of environmental load molecule generation sources extracted by the extraction step. An index deriving step for deriving an index for specifying an environmental load molecule generation source that affects the molecule, and an output step for outputting the index derived by the deriving step.

  The environmental load molecule generation source evaluation system of the present invention specifies a measurement device that measures the mass of molecules contained in a gas and an environmental load molecule generation source that affects the environmental load molecules detected at the environmental load molecule detection location. An environmental load molecule generation source evaluation system having an analysis device for deriving an index for the measurement, wherein the measurement device irradiates a vacuum ultraviolet laser beam so that a plurality of molecules contained in the gas are converted into one photon. Ionization at once and measuring the mass spectrum of each molecule was performed at each of a plurality of environmental load molecule generation sources and the environmental load molecule detection location, and the analysis device was measured by the measurement device The concentration of each molecule obtained based on the mass and spectral intensity of each molecule, wherein the concentration of each molecule in the plurality of environmental load source and the environmental condition are defined as boundary conditions. Then, by deriving a simulation calculation model formula of numerical fluid dynamics, deriving a molecular concentration distribution that is a concentration distribution of each molecule in an analysis region that is a predetermined three-dimensional region is performed with different boundary conditions. A concentration of each molecule at each position in the analysis region obtained from the molecular concentration distribution derived by the first molecular concentration distribution deriving means and the molecular concentration distribution derived from the first molecular concentration distribution deriving means; A storage for storing in the storage medium as a simulation result, the concentration of each molecule in the plurality of environmental load molecule generation sources when derived and at least a part of the environmental conditions when the molecular concentration distribution is derived Means, the concentration of each molecule obtained based on the mass of each molecule measured by the measuring device, and the concentration of each molecule at the environmental load molecule detection site The concentration of each molecule in the plurality of environmental load molecule generation sources corresponding to the environmental load molecule detection location and the environmental condition that can be regarded as the environmental condition when the environmental load molecule is detected at the environmental load molecule detection location. Extraction means for extracting from the result of the simulation, and environmental load molecules detected at the environmental load molecule detection location based on the concentration of each molecule in the plurality of environmental load molecule generation sources extracted by the extraction means And an index deriving unit for deriving an index for identifying an environmental load molecule generation source that affects the source, and an output unit for outputting the index derived by the index deriving unit.

  A computer program of the present invention is a computer program for causing a computer to derive an index for identifying an environmental load molecule generation source that affects an environmental load molecule detected at an environmental load molecule detection location. The mass and spectral intensity of each molecule measured with a measuring device that ionizes a plurality of molecules contained in the gas with a single photon by irradiating vacuum ultraviolet laser light and measures the mass spectrum of each molecule. The concentration of each molecule derived based on the above, the concentration of each molecule in the plurality of environmental load molecule generation source, and the environmental condition, by solving the simulation calculation model formula of numerical fluid dynamics as boundary conditions, Deriving a molecular concentration distribution that is a concentration distribution of each molecule in an analysis region that is a predetermined three-dimensional region, A first molecular concentration distribution deriving step performed with different boundary conditions, and each molecule at each position in the analysis region obtained by the molecular concentration distribution derived by the first molecular concentration distribution deriving step. A concentration, a concentration of each molecule in the plurality of environmental load molecule generation sources when the molecular concentration distribution is derived, and at least a part of an environmental condition when the molecular concentration distribution is derived; As a result, the storage step of storing in the storage medium, the concentration of each molecule derived based on the mass of each molecule measured by the measurement device, the concentration of each molecule at the environmental load molecule detection location, The environmental load molecule detection location and the environmental conditions that can be regarded as the environmental conditions when the environmental load molecule is detected at the environmental load molecule detection location, Based on the extraction step of extracting the concentration of the molecule from the result of the simulation and the concentration of each molecule at the plurality of environmental load molecule generation sources extracted by the extraction step, the concentration of the molecule is detected at the environmental load molecule detection location. Causing a computer to execute an index deriving step for deriving an index for identifying a source of an environmental load molecule that affects the environmental load molecule, and an output step for outputting the index derived by the index deriving step. It is characterized by.

  According to the present invention, the mass of each molecule in each of a plurality of environmental load molecule generation sources and environmental load molecule detection locations is measured by vacuum ultraviolet one-photon ionization mass spectrometry, and the concentration of each molecule is determined from the mass of each molecule. Is derived. From the results of solving the computational model equations of numerical fluid dynamics simulation using the concentration and environmental conditions of each molecule in multiple odor generating sources as boundary conditions, the concentration of each molecule in multiple environmental load molecular generating sources and each position at each position The concentration of molecules and environmental conditions are stored in association with each other. Extract the concentration of each molecule from a plurality of environmental load molecule generation sources corresponding to the concentration of each molecule at the environmental load molecule detection location, and based on the concentration of each molecule from the environmental load molecule generation source, An index for identifying a source of an environmental load molecule that affects the detected environmental load molecule is derived. Therefore, it is possible to identify the environmental load molecule generation source that affects the environmental load molecule detected at the environmental load molecule detection location as accurately as possible.

It is a figure which shows notionally the outline | summary of the odor generation source evaluation method (odor generation source evaluation system). It is a figure which shows an example of the schematic structure of the whole vehicle-mounted odor measuring system. It is a figure which shows an example of schematic structure of an odor measuring apparatus. It is a figure which shows an example of schematic structure of a vacuum ultraviolet light generation part. It is a figure which shows an example of schematic structure of an ionization part and a time-of-flight mass spectrometry part. It is a figure which shows an example of schematic structure of a gas introduction part. It is a figure which shows an example of a functional structure of an odor generation source analyzer. It is a figure which shows an example of a mass spectrum. It is a figure which shows an example of molecular concentration distribution. It is a figure which shows an example of the data item of the database in which a value is given based on the result of simulation in a table format. It is a figure which shows an example of the olfactory sensitivity coefficient of each molecule | numerator. It is a flowchart explaining an example of operation | movement of the odor generation source analyzer at the time of memorize | storing the result of simulation. It is a flowchart explaining the 1st example of operation | movement of the odor generation source analyzer at the time of specifying the odor generation source which most affects the odor in an odor detection location. It is a flowchart explaining the 2nd example of operation | movement of the odor generation source analyzer at the time of specifying the odor generation source which most affects the odor in an odor detection location. It is a figure which shows the result of having simulated the diffusion of NO molecule and benzene molecule.

Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[First Embodiment]
First, the first embodiment will be described.
FIG. 1 is a diagram conceptually showing an outline of an odor generation source evaluation method (odor generation source evaluation system) in the present embodiment.
In FIG. 1, smoke (gas) containing molecular species b and c is generated from the chimney of the factory A 100a, and smoke (gas) containing molecular species a, c and e is generated from the chimney of the factory B 100b. It is assumed that smoke (gas) including molecular species a, b, and c is generated from the C factory 100c, and smoke (gas) including molecular species b and c is generated from the chimney of the D factory 100d. Further, it is assumed that odor is generated at the odor detection location 110. Here, since the flow of gas from each factory 100a to 100d to the odor detection location 110 and the amount of flow of each gas are unknown, in the present embodiment, these are derived as follows.

First, the density | concentration of each molecule | numerator in A factory 100a-D factory 100d (chimney) used as these odor generating sources is measured. Then, based on the measurement result, a numerical fluid dynamics simulation is performed to derive the concentration distribution of each molecule in the analysis region 120 that is a three-dimensional region (three-dimensional coordinate region) including the odor detection part 110. Are stored under various conditions and stored together with the environmental conditions (meteorological conditions) at that time. At this time, it is unknown that odor is generated at the odor detection location 110.
Thereafter, when it is found that an odor is generated at the odor detection location 110, the concentration of each molecule at the odor detection location 110 is measured. Then, the concentration of each molecule in each factory 100a to 100d (the chimney) corresponding to the measurement result and the environmental condition (meteorological condition) at that time is extracted from the simulation result described above, and the extracted result is used as the basis. In addition, an index for identifying a factory that affects the odor detected at the odor detection location 110 is derived.
Hereinafter, the odor generation source evaluation method (odor generation source evaluation system) in this embodiment will be described in detail. The odor generation source and the odor detection location shown in FIG. 1 are examples, and the odor generation source and the odor detection position of the present embodiment are not limited to those shown in FIG.

In the present embodiment, an in-vehicle odor measurement system is adopted so that an odor generation source affecting the odor generated at the odor detection location can be identified on the spot.
FIG. 2 is a diagram schematically illustrating an example of the overall schematic configuration of the in-vehicle odor measurement system. In the present embodiment, the vehicle-mounted odor measurement system is installed in the one box car 10. 2A is a view of the inside of the one box car 10 as viewed from above, and FIG. 2B is a view of the inside of the one box car 10 as viewed from the side.
As shown in FIG. 2, the in-vehicle odor measurement system includes a battery box 20, an odor measurement device 30, an odor generation source analysis device 40, a power cord 51 to 53, a gas introduction unit 60, and a power connector 70. Have.

<Battery box 20>
The battery box 20 is for supplying electric power to the odor measuring device 30, the odor source analysis device 40, and the like. The battery box 20 of the present embodiment includes a battery, an inverter, a battery charger, a breaker, and the like.
The DC power output from the battery is converted into AC power by the inverter. The AC power is output to the odor measuring device 30 and the odor generation source analyzing device 40 through the power cord 51.

<Odor measuring device 30>
The odor measuring apparatus 30 of this embodiment is a vacuum ultraviolet one-photon ionization mass spectrometer using a vacuum ultraviolet laser beam having an ionization potential of 10.5 eV. Each component of the odor measuring device 30 is fixed to the same frame attached inside the one box car 10.
FIG. 3 is a diagram schematically illustrating an example of a schematic configuration of the odor measuring device 30.
In FIG. 3, a laser generation unit 31, a mirror 32, a vacuum ultraviolet light generation unit 33, and a time-of-flight mass analysis unit 34 included in the odor measurement device 30 are shown.

FIG. 4 is a diagram (cross-sectional view) schematically illustrating an example of a schematic configuration of the vacuum ultraviolet light generation unit 33.
In FIG. 4, the vacuum ultraviolet light generator 33 includes a synthetic quartz lens 33a, a MgF ₂ lens 33b, a vacuum chamber 33c, and a Xe gas introduction pipe 33d.
In the vacuum ultraviolet light generation unit 33, an ultraviolet light (ultraviolet laser light) 35 generated by the laser generation unit 31 (oscillation) in a gas cell in which the Xe gas 36 is sealed at an appropriate pressure (the third harmonic of the Nd: YAG laser (355 nm). )) Is introduced through the mirror 32, and the ultraviolet light 35 is again converted into a third harmonic wave at the focal point to generate (oscillate) vacuum ultraviolet light (vacuum ultraviolet laser light) 37 having a wavelength of 118 nm. The larger the number of laser repetitions, the better. However, this number of repetitions depends on the number of repetitions of the ultraviolet laser that is the source of the generation of vacuum ultraviolet light. Note that the wavelength of the vacuum ultraviolet light generated from the vacuum ultraviolet light generator 33 is not limited to 118 nm. For example, the range of 100 nm to 150 nm can be used as the wavelength of vacuum ultraviolet light.

Here, when the ultraviolet light 35, which is the third harmonic of the Nd: YAG laser, is directly incident on the ionization portion, the measurement molecule is highly likely to be decomposed. As a countermeasure, the ultraviolet light 35 is separated by adjusting the angle and position of the MgF ₂ lens 33b, which is a condenser lens, using the difference in refractive index between the ultraviolet light 35 and the vacuum ultraviolet light 37. Thus, it is effective to design the optical system so that only the vacuum ultraviolet light 37 is incident on the ionization portion.
The relationship between the photon energy ε and the wavelength λ of the ultraviolet laser light is expressed by the following equation (1).
ε = hc / λ (1)
However, h is Planck's constant and c is the speed of light.
Therefore, when the ionization potential of a certain target molecule is IP, the photon energy ε is set so as to satisfy the following expressions (2) and (3).
IP <ε (2)
IP <hc / λ (3)
When these conditions are satisfied, it becomes possible to ionize target target molecules all at once with one photon. That is, all molecules having an ionization potential lower than the photon energy of the vacuum ultraviolet laser can be detected all at once. As described above, in this embodiment, since the ionization potential of the vacuum ultraviolet laser beam is 10.5 eV, molecules such as nitrogen, oxygen, argon, etc. that are present in the atmosphere at a constant rate are not ionized regardless of odor. Therefore, molecules that do not contribute to odor even if the amount of molecules present is large are not detected.

FIG. 5 is a diagram schematically illustrating an example of a schematic configuration of the ionization unit 38 and the time-of-flight mass analysis unit 34.
The ionization unit 38 and the time-of-flight mass analysis unit 34 are changed from atmospheric pressure to a vacuum state using, for example, a rotary pump and two turbo molecular pumps (manufactured by Pfeiffer Vacuum).
FIG. 6 is a diagram (cross-sectional view) schematically illustrating an example of a schematic configuration of the gas introduction unit 60.
The gas suction port of the gas introduction part 60 is installed in the right rear side surface of the one box car 10, for example. The gas introduction unit 60 includes a pump that introduces the gas flowing in from the gas suction port into the ionization unit 38. Even when a small pump with an exhaust speed of about 70 L / s is used as this pump, the gas introduction unit 60 has a dust filter in order to introduce and measure the outside air as it is without being disturbed. From the collected gas to be measured, dust of 7 μm or more is removed with a dust filter, and the gas to be measured that has passed through the dust filter is ionized as a vacuum chamber from the gas outlet 60a (minimal pinhole). Introduced continuously into section 38.

The gas introduction unit 60 shown in FIG. 6 is an electrode that also serves as an ejection port. A direct injection type gas introduction unit that introduces ionized molecules into the time-of-flight mass analysis unit 34 coaxially with the ejection port without bending. Part.
As shown in FIG. 6, in the direct injection type, the nose electrode 60b (the measurement object gas injection port 60a) and the ion drawing electrode 60c (the ion drawing port 60d) are arranged in a straight line so as to face each other. Be placed. For this reason, ionization efficiency and efficiency of introducing ions into the time-of-flight mass spectrometer 34 can be increased.

  The measurement target gas G introduced into the ionization unit 38 from the measurement target gas ejection port 60a of the nose electrode 60b is photoionized at the ionization point P by being irradiated with the vacuum ultraviolet light 37, and becomes the ion I. The ions I enter the time-of-flight mass analyzer 34 while accelerating through the ion inlet 60d of the ion pulling electrode 60c due to the potential difference between the nasal electrode 60b and the ion pulling electrode 60c.

Returning to the description of FIG. 5, the time-of-flight mass analyzer 34 includes an MCP (multi-channel plate) detector 34 a and a reflectron 34 b.
The ion I (ionized molecule) is turned back by the reflectron 34b and detected by the MCP detector 34a, and the signal of the ion I is output as a current amount. A mass spectrum (mass number) is obtained by accelerating the ionized measurement target gas in a pulse manner and detecting a time difference until detection by the MCP detection unit 34a. In other words, if the amount of charge received by the ions I is constant, a mass spectrum is obtained by utilizing the fact that the larger the mass-to-charge ratio, the slower the flight time.

In the present embodiment, the above-described measurement of the ion I signal is individually performed in each of the plurality of odor generating sources.
For example, when the odor source is a chimney, the one box car 10 is moved to the vicinity of the chimney. Then, the tip of the gas inlet 60 on the gas suction port side is arranged inside the chimney, and the signal of the ion I of the molecule that affects the odor among the molecules contained in the gas taken in from the chimney is described above. Measure as you did. In this case, for example, a heat-resistant probe that introduces gas into the inside is disposed at the tip of the gas inlet 60 on the gas suction port side, and the heat-resistant probe is inserted into the chimney.

  Further, the above-described signal of the ion I is also measured at the odor detection location. As described above, the odor detection location is a location where odor is detected by the gas generated from the odor generation source. Residential areas become odor detection points. In this case, the one box car 10 is moved to the odor detection location. Then, the tip of the gas inlet 60 side of the gas introduction part 60 is directed to the atmosphere, and the signal of the ion I of the molecule that affects the odor among the molecules contained in the gas taken in from the atmosphere is measured.

<Odor source analysis device 40>
The odor source analysis apparatus 40 of the present embodiment is configured using a preamplifier, a picometer, and a commercially available notebook personal computer (notebook PC). Each component of the odor generation source analyzing device 40 is also fixed to the same frame attached to the inside of the one box car 10, similarly to the odor measuring device 30.
The ion I signal detected by the MCP detector 34a is amplified by a preamplifier (for example, VT120 manufactured by ortec) and input to a notebook PC via a picometer. FIG. 7 is a diagram illustrating an example of a functional configuration of the odor generation source analyzer 40. The functions of each block shown in FIG. 7 are realized by executing a dedicated program installed in a notebook PC, for example. Below, an example of the function which the odor generation source analyzer 40 has is demonstrated.

(Gas mass information acquisition unit 41)
The gas mass information acquisition unit 41 acquires signals of ions I continuously detected by the MCP detection unit 34a. In this embodiment, the gas mass information acquisition part 41 acquires the signal of the ion I continuously detected by the MCP detection part 34a every second.
In this embodiment, the gas mass information acquisition part 41 acquires separately the signal of the ion I continuously detected by the MCP detection part 34a in each of several odor generating sources. Moreover, the gas mass information acquisition part 41 also acquires the signal of the ion I continuously detected by the MCP detection part 34a in the odor detection location.

(Mass spectrum creation part 42)
The mass spectrum creation unit 42 creates a mass spectrum based on the ion I signal acquired by the gas mass information acquisition unit 41. In this embodiment, since the signal of the ion I is obtained every second, a mass spectrum every second can be obtained.
FIG. 8 is a diagram illustrating an example of a mass spectrum. In FIG. 8, the horizontal axis is the mass number (m / z), and the vertical axis is the signal amount (arb.unit) of the ion I. In the example shown in FIG. 8, peaks 801, 802, 803, and 804 exist where the mass numbers are 79, 92, 103, and 106. These peaks 801, 802, 803, and 804 correspond to benzene, toluene, styrene, and xylene, respectively.
In the present embodiment, the mass spectrum creation unit 42 individually creates a mass spectrum based on the ion I signal acquired by the gas mass information acquisition unit 41 in each of the plurality of odor generating sources. The mass spectrum creation unit 42 also creates a mass spectrum based on the ion I signal acquired by the gas mass information acquisition unit 41 at the odor detection location.

(Molecular concentration deriving unit 43)
The molecular concentration deriving unit 43 derives the concentration (ppm) of each molecule based on the mass spectrum (mass and spectral intensity of each molecule) created by the mass spectrum creating unit 42.
In the present embodiment, a gas having a known concentration of all molecules to be measured is introduced into the odor measuring device 30, a mass spectrum is created as described above, and the relationship between the signal amount and the concentration in the mass spectrum is determined. A calibration curve is created and stored in advance for each molecule. Then, the molecular concentration deriving unit 43 extracts a signal amount indicating a peak corresponding to each molecule from the mass spectrum created by the mass spectrum creating unit 42, and based on the extracted signal amount and a calibration curve for the molecule, The concentration of the molecule is derived. In the example shown in FIG. 8, the concentration of benzene was 8 ppm, the concentration of toluene was 2 ppm, the concentration of styrene was 0.4 ppm, and the concentration of xylene was 0.5 ppm.

As described above, in the present embodiment, the molecular concentration deriving unit 43 obtains the mass spectrum every second, so that the concentration of each molecule is obtained every second. The molecular concentration deriving unit 43 also derives an average value of the concentration of each molecule for a predetermined period (10 minutes in this embodiment). The molecular concentration deriving unit 43 also derives the maximum value and the minimum value for each predetermined period (in this embodiment, 1 minute) of the concentration of each molecule. In the following description, “average value of molecular concentration over a predetermined period” is referred to as “average value of molecular concentration” as necessary, and “maximum value / minimum value of molecular concentration per predetermined period” This is referred to as “maximum value / minimum value of molecular concentration” as necessary.
In the present embodiment, the molecular concentration deriving unit 43 calculates the instantaneous value (value per second) / average value of the concentration of each molecule in each of a plurality of odor generating sources from the mass spectrum created by the mass spectrum creating unit 42.・ Deriving maximum and minimum values. Similarly, the molecular concentration deriving unit 43 derives the instantaneous value / average value / maximum value / minimum value of the concentration of each molecule at the odor detection location from the mass spectrum created by the mass spectrum creating unit 42.

(First molecular concentration distribution deriving unit 44)
The first molecular concentration distribution deriving unit 44 derives the distribution of the concentration of each molecule in each predetermined three-dimensional region (the concentration of each molecule at each position) by solving the simulation calculation model formula of numerical fluid dynamics. In the following description, “concentration distribution of each molecule” is referred to as “molecular concentration distribution” as necessary, and “predetermined three-dimensional region” is referred to as “analysis region”. Of course, the analysis region is set so that the plurality of odor generation sources and odor detection locations described above are all included.

In the present embodiment, the first molecular concentration distribution deriving unit 44 derives a molecular concentration distribution by a finite volume method (SIMPLE method). An example of a simulation calculation model formula of numerical fluid dynamics is shown below.
First, the following formulas (1) and (2) are used as mass conservation formulas for gas in each cell.

Here, ρ is density, t is time, and v is a flow velocity vector. Also, p _op is pressure, R is gas constant, M _W is molecular weight, and T is temperature.
Next, the following equations (3) and (4) are used as gas momentum conservation equations in each cell.

Here, p _S is a static pressure, τ _eff (double upper line) is a stress tensor, ρ ₀ is an average density, and g is a gravitational acceleration vector. Μ _eff is a viscosity coefficient, I is a unit tensor, and k is turbulent energy.
Next, the following equations (5) to (9) are used as gas energy equations in each cell.

Here, h is the enthalpy, k _eff is the effective thermal conductivity, and J _j is the diffusion velocity of the molecule j. Y _j is the mass fraction of molecule j and h _j is the enthalpy of molecule j. k is the thermal conductivity, c _p is the constant pressure specific heat, μ _t is the turbulent viscosity coefficient, and Pr _t is the turbulent Prandtl number. c _p (upper line) is the average low pressure specific heat. c _{p, j} is the constant pressure specific heat of the molecule j.
Next, the following formulas (10) to (12) are used as molecular concentration conservation formulas in each cell.

Here, D _{j, m} is the diffusion coefficient of the molecule j in the mixed gas (m means mixed), Sc _t is the number of random Schmidt. D _t is a turbulent diffusion coefficient.
Next, the following formulas (13) and (14) are used as a turbulent flow model (standard k-ε model) in each cell.

Here, ε is the turbulent energy dissipation rate, σ _k is the turbulent Prandtl number of _k , and G _k represents the generation of turbulent energy due to the gradient of the average flow velocity. σε is the turbulent Prandtl number of ε, and G ₁ ε, G ₂ ε, C ₁ ε, and C ₃ ε are constants.

The first molecular concentration distribution deriving unit 44 solves the equations (1) to (14) by the finite volume method (SIMPLE method) to obtain the mass fraction Y _{j of the} molecule j (that is, the concentration of the molecule j). The steady solution at each position is derived as a molecular concentration distribution.
When solving the equations (1) to (14), the first molecular concentration distribution deriving unit 44 is a region that forms a boundary with the atmosphere included in the analysis region and affects the wind flow. Gas type, concentration, outflow, temperature, density at any position between the boundary surface (surface of chimney, building, bank, tree, sea, river, road, etc.) and the boundary surface of the analysis area And information on the wind speed, direction, and temperature of the atmosphere as boundary conditions. At this time, the initial values of the wind speed and direction of the atmosphere were fixed, and the wind direction and wind speed that changed thereby were used. Here, as the boundary condition of the concentration of each molecule at the position of the odor generating source, the concentration of each molecule in the odor generating source derived by the molecular concentration deriving unit 43 is given. As other boundary conditions, measured values and the like are appropriately given. In addition, known parameters (gas constant R and the like) in equations (1) to (14) are also given as appropriate. Since the analysis region is a three-dimensional region, the region where the boundary condition is set is also a three-dimensional region. FIG. 9 is a diagram schematically showing an example of the molecular concentration distribution of naphthalene in the analysis region (factory) obtained by solving the simulation calculation model formula of numerical fluid dynamics. In FIG. 9, naphthalene is present in a region 901 (it is difficult to see for the sake of notation, but the density varies depending on the density in the region 901). In addition, the present inventors have said odor obtained from the concentration of naphthalene at the odor detection location 902 obtained from the molecular concentration distribution in the analysis region shown in FIG. 9 and the mass spectrum measured by vacuum ultraviolet one-photon ionization mass spectrometry. Both of the naphthalene concentrations at the detection location 902 were about 100 ppb, and it was confirmed that the two coincided well.

  In the present embodiment, the first molecular concentration distribution deriving unit 44 derives the instantaneous value / average value / maximum value / concentration value of each molecule at each position of the plurality of odor generating sources derived by the molecular concentration deriving unit 43. Each of the minimum values is individually given as a boundary condition, and the above-described numerical simulation model of computational fluid dynamics is solved to derive the molecular concentration distribution of each molecule. Further, by providing various boundary conditions, the molecular concentration distribution of each molecule under various conditions is derived. Based on the instantaneous value, average value, maximum value, and minimum value of the concentration of each molecule derived from data observed at multiple odor sources, with a constant step size within the possible concentration range of each molecule. The concentration taken was given as a boundary condition. For the weather conditions (environmental conditions) of wind speed, wind direction, and temperature, the values taken at a constant step size within the possible range were given as boundary conditions. In this way, the simulation was performed under all possible conditions to derive the molecular concentration distribution of each molecule in the analysis region. At that time, the concentration of each molecule generated from the odor source may change depending on the situation, but in the same source, the relative ratio of the concentration of each molecule is not changed. Changing the concentration of each molecule while maintaining it may be a boundary condition.

(Data storage unit 45)
The data storage unit 45 stores, as a database, a simulation result based on the molecular concentration distribution of each molecule derived by the first molecular concentration distribution deriving unit 44 and the boundary condition when the molecular concentration distribution is derived. .
FIG. 10 is a diagram illustrating an example of a data item of a database (a database stored by the data storage unit 45) to which a value is given based on a simulation result in a table format.
In this embodiment, the position in the analysis region (position that is a candidate for the odor detection location), the concentration of each molecule at that position (molecular concentration), and each boundary condition set when the concentration is derived Information in which the concentration of the odor source (source concentration) and the wind direction, wind speed, and temperature of the atmosphere at a predetermined atmospheric observation position (area) are correlated with each other is stored as a simulation result by the data storage unit 45. The In FIG. 10, there are four odor generation sources “a”, “b”, “c”, “d”, and molecules “A”, “B”, “C”, “D”, “D”. The case where there are five of “E” is shown as an example. Further, the atmospheric observation position (region) is, for example, a predetermined position (region) at the boundary of the side surface of the analysis region. However, the atmospheric observation position (region) is not limited to this, and may be, for example, an arbitrary position in the analysis region (a position that is a candidate for an odor detection location).

  In the present embodiment, the value obtained by the first molecular concentration distribution deriving unit 44 is not stored as information on the simulation result, but the value is classified. Specifically, of the data items shown in FIG. 10, any one of 16 directions is set for “wind direction”. For other data items (“wind speed”, “temperature”, “source concentration”, “position”, and “molecular concentration”), any one of a plurality of preset ranges is set. As for “position”, a certain area where the same odor is considered to be generated is determined in advance as the same position, and information for specifying any of the areas is set. Also, representative values are set for each of “wind direction”, “wind speed”, “temperature”, “source concentration”, “position”, and “molecular concentration”. In the present embodiment, the representative value is used in subsequent calculations.

Further, as described above, in this embodiment, the instantaneous value, average value, maximum value, and minimum value of the concentration of each molecule at each position of a plurality of odor generating sources are given as boundary conditions. Therefore, even if “wind direction”, “wind speed”, “temperature”, “position”, and “molecular concentration” other than “source concentration” shown in FIG. 10 are the same, different information is set as “source concentration”. May be.
As described above, since the first molecular concentration distribution deriving unit 44 derives the molecular concentration distribution of each molecule under various conditions, the data storage unit 45 stores various simulation results. All combinations of information that can be set as “wind direction”, “wind speed”, “temperature”, “source concentration”, “position”, and “molecular concentration” are stored in the data storage unit 45. Most preferred.

(Data extraction unit 46)
The data extraction unit 46 starts operation when the data storage unit 45 stores the simulation result and the molecular concentration deriving unit 43 derives the concentration of each molecule at the odor detection location. In the following description, for convenience of description, “instantaneous value / average value / maximum value / minimum value of molecule concentration” is abbreviated as “molecule concentration” as necessary.

The data extraction unit 46 acquires information on the wind direction, wind speed, and temperature of the atmosphere at the above-described atmospheric observation position (region). As such information, for example, a result measured at the atmospheric observation position (region) at a time close to the time when the measurement at the odor detection location is performed by the odor measurement device 30 can be adopted. The data extraction unit 46, for example, corresponds to the concentration of each molecule at the odor detection location derived by the molecular concentration deriving unit 43 and the information on the wind direction, wind speed, and temperature of the atmosphere at the atmospheric observation position (region). The concentration (generation source concentration) of each molecule in the odor generation source is extracted from the simulation result stored in the data storage unit 45 (see FIG. 10).
The concentration of each molecule in each odor source (source concentration) includes an instantaneous value / average value / maximum value / minimum value, but the data extraction unit 46 determines the instantaneous value / average value / maximum value / Of the minimum values, the record of the one specified by the user is extracted.

(Odor source selection unit 47)
The odor generation source selection unit 47 selects all the combinations of odor generation sources except one of the plurality of odor generation sources one by one in order. Since the example shown in FIG. 10 shows the case where there are four odor generation sources that may affect the odor at the odor detection location, the odor generation source selection unit 47 is an unselected combination 3 One odor source is selected in turn.
(Second molecular concentration distribution deriving unit 48)
The second molecular concentration distribution deriving unit 48 is an odor among the conditions when a simulation result corresponding to the concentration of each molecule (generation source concentration) in each odor generation source extracted by the data extraction unit 46 is obtained. A molecular concentration distribution of each molecule is derived based on conditions excluding boundary conditions for one odor generation source not selected by the generation source selection unit 47, and the odor detection location is determined based on the derived molecular concentration distribution. Deriving the concentration of each molecule. That is, the second molecular concentration distribution deriving unit 48 assumes that no gas is discharged from one odor generation source that has not been selected by the odor generation source selection unit 47, and other conditions are extracted by the data extraction unit 46. The molecular concentration distribution of each molecule is derived using the same conditions as those obtained when the simulation result corresponding to the concentration (source concentration) of each molecule in each of the generated odor generation sources is obtained. Since the method for deriving the molecular concentration distribution is the same as that by the first molecular concentration distribution deriving unit 44, detailed description thereof is omitted here.

(Odor source identifying unit 49)
The odor generation source specifying unit 49 is stored in the data storage unit 45 in association with the concentration (generation source concentration) of each molecule in each odor generation source extracted by the data extraction unit 46 (at the odor detection location). ) Based on the concentration of each molecule (molecular concentration), the odor concentration at the odor detection location is derived. The odor concentration is expressed by the following equation (15).

Here, Y _j is the concentration of molecule j at the odor detection location. K _j is the olfactory sensitivity coefficient of molecule j. FIG. 11 is a diagram illustrating an example of the olfactory sensitivity coefficient of each molecule.
In FIG. 11, the olfactory sensitivity coefficient is obtained by quantifying the degree of influence of each molecule on the odor. The olfactory threshold represents the minimum value of the concentration that humans perceive as odor (see Non-Patent Documents 2 and 3). For example, when benzene reaches 2.7 ppm, humans perceive it as an odor. The olfactory threshold of each molecule when the olfactory threshold of benzene is used as a reference is the olfactory sensitivity coefficient. For example, the olfactory sensitivity coefficient of toluene is 2.7 ÷ 0.33 = 8.18≈8. Thus, the odor concentration is represented by the sum of values obtained by multiplying the concentration Y _j of each molecule at the odor detection location and the olfactory sensitivity coefficient K _j of each molecule for each identical molecule.

Further, the odor generation source specifying unit 49 derives the odor sensitivity at the odor detection location by the formula (15) based on the concentration of each molecule at the odor detection location derived by the second molecular concentration distribution deriving unit 48. To do.
In the following explanation, “the odor detection location derived based on the concentration of each molecule at the odor detection location corresponding to the concentration (generation source concentration) of each molecule in each odor generation source extracted by the data extraction unit 46. The “odor concentration at the bottom” is referred to as “odor concentration by all sources” as necessary. On the other hand, “the odor concentration at the odor detection location derived based on the concentration of each molecule at the odor detection location derived by the second molecular concentration distribution deriving unit 48” is set to “exclude specific sources” as necessary. This is referred to as “odor concentration by source”.

  As described above, the odor generation source selection unit 47 selects all the combinations of odor generation sources excluding one of the plurality of odor generation sources one by one in order, and for each of these combinations, the second The molecular concentration distribution deriving unit 48 derives the concentration of each molecule at the odor detection location. Therefore, the odor concentration by the source excluding the specific source is derived by the number of combinations selected by the odor source selection unit 47. On the other hand, only one odor concentration from all sources is derived.

  The odor source identifying unit 49 determines the odor concentration difference value, which is a value obtained by subtracting the odor concentration by the source excluding the specific source from the odor concentration by all the sources, and the odor concentration by the source excluding the specific source. Are derived individually. Next, the odor generation source specifying unit 49 selects an odor generation source that is not selected by the odor generation source selection unit 47 when the maximum odor concentration difference value is obtained from the derived odor concentration difference values. It is specified as an odor source that generates gas that most affects the odor.

(Output unit 50)
The output unit 50 outputs information including information on the odor generating source specified by the odor generating source specifying unit 49. For example, the output unit 50 displays information including information on the odor generation source specified by the odor generation source specifying unit 49 on the display device, or the storage medium in the odor generation source analysis device 40 or the odor generation source analysis device 40. The information including the information on the odor generating source specified by the odor generating source specifying unit 49 is output by storing the information in a portable storage medium connected to the odor generating device.

  In addition to the odor generation source information specified by the odor generation source specifying unit 49, the output unit 50 includes the concentration of each molecule in the odor generation source selected by the odor generation source selection unit 47 and the second molecular concentration distribution. Information specifying the molecular concentration distribution derived by the deriving unit 48 may be output. The output unit 50 also detects the odor detection location derived by the second molecular concentration distribution deriving unit 48 in addition to or instead of the information specifying the molecular concentration distribution derived by the second molecular concentration distribution deriving unit 48. The concentration of each molecule in may be output. Further, the output unit 50 may output at least one of the odor concentration by all the generation sources derived by the odor generation source specifying unit 49 and the odor concentration by the generation sources excluding the specific generation source. Furthermore, the output unit 50 may output other information derived by the odor generation source analyzer 40.

<Flowchart>
Next, an example of the operation of the odor source analysis apparatus 40 when storing the simulation result will be described with reference to the flowchart of FIG.
First, in step S1201, the gas mass information acquisition unit 41 acquires a signal (gas mass information) of ions I continuously detected by the MCP detection unit 34a using the odor generation source as a measurement location.
Next, in step S1202, the mass spectrum creation unit 42 creates a mass spectrum of the gas generated in the odor generating source based on the ion I signal acquired in step S1201 (see FIG. 8).

Next, in step S1203, the molecular concentration deriving unit 43 derives the concentration of each molecule of gas generated in the odor generating source based on the mass spectrum created in step S1202 and the calibration curve stored in advance.
Next, in step S1204, the first molecular concentration distribution deriving unit 44 determines whether or not the processing in steps S1201 to S1203 has been performed for all the preset odor generation sources. As a result of the determination, if the processes of steps S1201 to S1203 are not performed for all the odor generation sources, the process returns to step S1201 and the processes of steps S1201 to S1204 are repeated until the processes are performed for all the odor generation sources. .

And if the process of step S1201-S1203 is performed about all the odor generation sources, it will progress to step S1205.
In step S1205, the first molecular concentration distribution deriving unit 44 sets a plurality of predetermined boundary conditions including the “concentration of each molecule of gas generated from the odor generating source” derived in step S1203. ), One unset boundary condition (set) is set.
Next, in step S1206, the first molecular concentration distribution deriving unit 44 calculates the numerical calculation model equations of numerical fluid dynamics shown in equations (1) to (14) based on the boundary conditions set in step S1205. The molecular concentration distribution is derived by solving by the finite volume method (SIMPLE method) (see FIG. 9).

Next, in step S1207, the data storage unit 45 stores a simulation result based on the molecular concentration distribution derived in step S1206 and the boundary condition when the molecular concentration distribution is derived (see FIG. 10). ).
Next, in step S1208, the first molecular concentration distribution deriving unit 44 determines whether or not the processing in steps S1205 to S1207 has been performed for all predetermined boundary conditions (sets). As a result of this determination, if the processing in steps S1205 to S1207 has not been performed for all predetermined boundary conditions (sets), the process returns to step S1205, and processing is performed for all boundary conditions (sets). Until then, the processing of steps S1205 to S1207 is repeated.
Then, when the processes of steps S1205 to S1207 are performed for all (boundary) boundary conditions (sets), the process of the flowchart of FIG.

Next, an example of the operation of the odor generation source analysis device 40 when specifying the odor generation source generating the gas that most affects the odor at the odor detection location will be described with reference to the flowchart of FIG.
First, in step S1301, the gas mass information acquisition unit 41 acquires signals (gas mass information) of ions I continuously detected by the MCP detection unit 34a using the odor detection location as a measurement location.
Next, in step S1302, the mass spectrum creation unit 42 creates a mass spectrum of the gas at the odor detection location based on the ion I signal acquired in step S1301 (see FIG. 8).

Next, in step S1303, the molecular concentration deriving unit 43 derives the concentration of each molecule at the odor detection location based on the mass spectrum created in step S1302 and the calibration curve stored in advance.
Next, in step S1304, the data extraction unit 46 corresponds to the concentration of each molecule in the odor detection location derived in step S1303 and the information on the wind direction, wind speed, and temperature of the atmosphere at the atmospheric observation position (region). The concentration (generation source concentration) of each molecule in the odor generation source is extracted from the simulation result (see FIG. 10). In the present embodiment, the concentration (generation source concentration) of each molecule at each odor generation source includes an instantaneous value, an average value, a maximum value, and a minimum value. Only one specified by the user is selected.

Next, in step S1305, the odor source identifying unit 49 detects the odor stored as a simulation result in association with the concentration (source concentration) of each molecule in each odor source extracted in step S1304. Based on the concentration of each molecule at the location, the calculation of the equation (15) is performed to derive the odor concentration at the location where the odor is detected (odor concentration due to all sources).
Next, in step S1306, the odor generation source selection unit 47 selects an unselected combination from the combinations of odor generation sources excluding one of the plurality of odor generation sources.

Next, in step S1307, when the second molecular concentration distribution deriving unit 48 obtains a simulation result corresponding to the concentration (generation source concentration) of each molecule in each odor generation source extracted in step S1304. Among these conditions, the molecular concentration distribution of each molecule is derived based on the conditions excluding the boundary condition for one odor generation source not selected in step S1306, and the odor detection is performed based on the derived molecular concentration distribution. The concentration of each molecule at the location is derived.
Next, in step S1308, the odor generation source specifying unit 49 calculates the equation (15) based on the concentration of each molecule in the odor detection location derived in step S1307, and the odor concentration ( Derivation of odor concentration by sources excluding specific sources).

  Next, proceeding to step S1309, the odor generation source selection unit 47 determines whether or not all of the combinations of odor generation sources other than one of the plurality of odor generation sources have been selected. As a result of the determination, if not all the combinations of odor generation sources except for one of the plurality of odor generation sources have been selected, the process returns to step S1306. And the process of step S1306-S1309 is repeatedly performed until it selects all the combinations of the odor generation source except one of several odor generation sources.

When all the combinations of odor generation sources excluding one of the plurality of odor generation sources are selected, the process proceeds to step S1310. In step S1310, the odor generation source specifying unit 49 derives the odor concentration difference value, which is a value obtained by subtracting the odor concentration of the generation sources excluding the specific generation source from the odor concentration of all the generation sources in step S1308. Each odor concentration by the source excluding the specific source is derived individually, and the odor source that generates the gas that most affects the odor at the odor detection location is specified from the odor concentration difference value.
Next, in step S1311, the output unit 50 outputs information including information on the odor generation source specified in step S1310. And the process by the flowchart of FIG. 13 is complete | finished.

<Summary>
As described above, in this embodiment, the concentration of each molecule contained in the gas generated at each odor source is measured by vacuum ultraviolet one-photon ionization mass spectrometry using a vacuum ultraviolet laser beam having an ionization potential of 10.5 eV. To do. By solving the simulation calculation model of computational fluid dynamics using the concentration of each molecule contained in the gas generated from each odor source and the environmental conditions as boundary conditions, the analysis of each molecule in the analysis area, which is a predetermined three-dimensional area, is performed. Deriving the concentration distribution is performed under various boundary conditions. The concentration of each molecule at each position obtained as a result, the concentration of each molecule at each odor generation source when the result is obtained, and the environmental condition are stored as simulation results. Thereafter, the concentration of each molecule at the odor detection location is measured by vacuum ultraviolet one-photon ionization mass spectrometry using a vacuum ultraviolet laser beam having an ionization potential of 10.5 eV, and the concentration of each molecule at the measured odor detection location is measured. And the concentration of each molecule at each odor source corresponding to the environmental condition is extracted from the simulation result. Based on the concentration of each molecule at each extracted odor source, the odor concentration by all sources and the odor concentration by sources excluding specific sources are derived, and based on these, Derives and outputs information on the odor source that has the most influence on the odor that is present.

Therefore, first, the concentration of each molecule contained in the gas generated in each odor source is measured by vacuum ultraviolet one-photon ionization mass spectrometry using a vacuum ultraviolet laser beam having an ionization potential of 10.5 eV. Molecules that exist at a constant rate in the atmosphere, such as oxygen and argon, are not ionized regardless of odor. Therefore, molecules that do not contribute to odor even if the amount of molecules present is large are not detected.
In addition, by solving the computational calculation model of computational fluid dynamics, the distribution of the concentration of each molecule in the analysis area, which is a predetermined three-dimensional area, is derived. Reflecting this, the concentration of each molecule in the analysis region can be derived.
In addition, because the odor source that affects the odor generated at the odor detection location is identified based on the odor concentration, odor detection is performed by considering not only the amount of molecules but also the degree of influence on the odor. It is possible to identify the odor source that affects the odor generated at the location.

<Modification>
In the present embodiment, in the process of identifying the odor generation source generating the gas that most affects the odor at the odor detection location, the odor concentration by all the generation sources and the odor concentration by the generation source excluding the specific generation source (See steps S1305 and S1306 to S1309). However, it is not always necessary to derive these in the process of identifying the odor generation source generating the gas that most affects the odor at the odor detection location. That is, for all of the simulation results (all of the records in FIG. 10), the odor concentration by all the sources and the odor concentration by the sources other than the specific source may be derived and stored. . In this case, the data storage amount increases, but the processing time can be shortened.

  Further, in this embodiment, as an index for identifying the odor generation source that affects the odor generated at the odor detection location, the odor generation source that most affects the odor generated at the odor detection location. Was derived. However, the index for identifying the odor generation source that affects the odor generated at the odor detection location is not limited to this. For example, it occurs at odor detection locations for multiple odor sources according to the odor concentration difference value, which is the value obtained by subtracting the odor concentrations from the sources excluding specific sources from the odor concentrations from all sources. A rank may be given in order from the one that has a great influence on the odor, and the information on the rank may be used as an index for identifying the odor generation source that affects the odor generated at the odor detection location.

Further, in the present embodiment, a steady solution of a computational calculation model formula of numerical fluid dynamics is derived. However, an unsteady solution may be derived. In this case, the boundary condition is changed with the passage of time. In such a case, the first molecular concentration distribution deriving unit 44 can derive the instantaneous value, average value, maximum value, and minimum value of the concentration of each molecule at the odor detection location.
If the notebook PC is also placed on the one-box car 10 as in the present embodiment, an indicator for identifying the odor generation source that affects the odor generated at the odor detection location is used as the odor detection. Can be derived at points. However, it is not always necessary to place the notebook PC on the one-box car 10.
In the present embodiment, the concentration of each molecule is derived by the odor generation source analyzer 40 (molecular concentration deriving unit 43), but it is not always necessary to do so. For example, the odor measuring device 30 may be provided with a function for deriving the concentration of each molecule, and the concentration of each molecule derived by the odor measuring device 30 may be input to the odor generation source analyzing device 40.
The odor measurement method and odor measurement system according to the first embodiment have been described above as an example of the environmental load molecule measurement method and the environmental load molecule measurement system. Odor is due to molecules. The first embodiment can be applied to general molecules regardless of molecules that emit odor. The first embodiment of the present invention can also be applied when the measurement molecule is an environmental load molecule.

[Second Embodiment]
Next, a second embodiment of the present invention will be described.
In the first embodiment, the odor generation source that is not selected when the odor concentration difference value, which is a value obtained by subtracting the odor concentration by the generation source excluding the specific generation source from the odor concentration by all the generation sources, is the largest. It was derived as the odor source that has the most influence on the odor generated at the odor detection location. On the other hand, in this embodiment, one odor generation source is selected from a plurality of odor generation sources, and each molecule is assumed to have gas discharged from only the selected odor generation source among the plurality of odor generation sources. The molecular concentration distribution is derived to derive the odor concentration. Then, the odor generation source selected when the highest odor concentration is obtained among the derived odor concentrations is derived as the odor generation source that most affects the odor generated at the odor detection location. As described above, the present embodiment and the first embodiment are mainly different in part of the processing when specifying the odor generation source that most affects the odor generated at the odor detection location. Therefore, in the description of the present embodiment, the same parts as those in the first embodiment are denoted by the same reference numerals as those in FIGS.

FIG. 14 is a flowchart for explaining an example of the operation of the odor generation source analyzer 40 when the odor generation source generating the gas that most affects the odor at the odor detection location is specified. In the present embodiment, the odor source analysis apparatus 40 executes the process according to the flowchart of FIG. 14 instead of the flowchart of FIG.
Steps S1401 to S1404 are the same as steps S1301 to S1304 in FIG.
That is, in step S1401, the gas mass information acquisition unit 41 acquires the ion I signal (gas mass information) at the odor detection location, and in step S1402, the mass spectrum creation unit 42 determines the gas mass spectrum at the odor detection location. In step S1403, the molecular concentration deriving unit 43 derives the concentration of each molecule at the odor detection location. In step S1404, the data extraction unit 46 extracts the concentration (generation source concentration) of each molecule at each odor generation source from the simulation result.

In step S1405, the odor generation source selection unit 47 selects one unselected odor generation source from the plurality of odor generation sources.
Next, in step S1406, the second molecular concentration distribution deriving unit 48 obtains a simulation result corresponding to the concentration (generation source concentration) of each molecule in each odor generation source extracted in step S1404. Among these conditions, the molecular concentration distribution of each molecule is derived based on the conditions excluding the boundary condition for the odor source not selected in step S1405, and the odor detection location is determined based on the derived molecular concentration distribution. The concentration of each molecule is derived.

Next, in step S1407, the odor generation source specifying unit 49 calculates the equation (15) based on the concentration of each molecule at the odor detection location derived in step S1406, and calculates the odor concentration at the odor detection location. To derive. In the following description, this odor concentration is referred to as “odor concentration by a specific source” as necessary.
Next, in step S1408, the odor generation source selection unit 47 determines whether or not all of the plurality of odor generation sources have been selected. As a result of the determination, if all of the plurality of odor generating sources are not selected, the process returns to step S1405. And the process of step S1405-S1408 is repeated until all the some odor generating sources are selected.

When all of the plurality of odor generating sources are selected, the process proceeds to step S1409. In step S1409, the odor generation source selected in step S1405 when the odor concentration by the largest specific generation source among the odor concentrations by the specific generation source derived in step S1407 is obtained is determined at the odor detection location. Identified as the source of odor generating gas that most affects odor.
Next, in step S1410, the output unit 50 outputs information including information on the odor generation source specified in step S1409. And the process by the flowchart of FIG. 14 is complete | finished.
Even if it does as mentioned above, there exists the same effect as what was demonstrated in 1st Embodiment.

<Modification>
As in the first embodiment, this embodiment has the most influence on the odor generated at the odor detection location as an index for identifying the odor generation source that affects the odor generated at the odor detection location. The source of odor generation is derived. However, even in this embodiment, a plurality of odor generation sources are ranked in order from the one that has a great influence on the odor generated at the odor detection location, and information on the ranking is generated at the odor detection location. It is good also as a parameter | index for pinpointing the odor generation source which affects the odor which it has. In this case, the ranking is performed in descending order of the odor concentration by the specific generation source. In addition, the ratio of each odor source that contributes to the odor generated at the odor detection location shall be derived as an index for identifying the odor generation source that affects the odor generated at the odor detection location. Can do. For example, a value obtained by dividing the odor concentration by the specific generation source by the odor concentration by all the generation sources can be derived as the ratio of the odor generation source contributing to the odor generated at the odor detection location.
In addition, also in the present embodiment, various modifications described in the first embodiment can be employed.
The odor measurement method and odor measurement system according to the second embodiment have been described above as an example of the environmental load molecule measurement method and the environmental load molecule measurement system. Odor is due to molecules. The second embodiment can be applied to general molecules regardless of molecules that emit odor. The second embodiment of the present invention can also be applied when the measurement molecule is an environmental load molecule.

[Relationship with Claims]
<Claims 1, 7, 13>
The first environmental load molecule measuring step is realized, for example, when the odor measuring device 30 detects a signal of ions I in each odor generating source.
The first molecular concentration deriving step is realized, for example, when the molecular concentration deriving unit 43 executes the process of step S1203 in FIG.
The first molecular concentration distribution deriving step (means) is realized, for example, when the first molecular concentration distribution deriving unit 44 executes the processing of steps S1205 and S1206 in FIG. At least a part of the environmental conditions corresponds to, for example, the wind direction, wind speed, and temperature of the atmosphere at a predetermined atmosphere observation position (region).
The storage step (means) is realized, for example, by the data storage unit 45 executing the process of step S1207 in FIG.
The second environmental load molecule measurement step is realized, for example, when the odor measuring device 30 detects the ion I signal at the odor detection location.
The second molecular concentration deriving step is realized, for example, when the molecular concentration deriving unit 43 executes the process of step S1303 in FIG. 13 or the process of step S1403 in FIG.
The extraction process (means) is realized, for example, by the data extraction unit 46 executing the process of step S1304 in FIG. 13 or the process of step S1404 in FIG.
In the index deriving step (means), for example, the second molecular concentration distribution deriving unit 48 executes the process of step S1307 in FIG. 13 or the process of step S1407 in FIG. 14, and the odor source identifying unit 49 in FIG. This is realized by executing the processing of steps S1305, S1306, S1308 to S1310 of FIG. 14 or steps S1405 and S1407 to S1409 of FIG.
The output process (means) is realized, for example, when the output unit 50 executes the process in step S1311 in FIG. 13 or the process in step S1410 in FIG.
The odor generation source evaluation system is realized by, for example, an in-vehicle odor measurement system.
<Claims 2 and 8>
The second molecular concentration distribution deriving step (means) is realized, for example, when the second molecular concentration distribution deriving unit 48 executes the process of step S1307 in FIG.
The first odor concentration deriving step (means) is realized, for example, when the odor generation source specifying unit 49 executes the process of step S1305 in FIG.
The second odor concentration deriving step (means) is realized, for example, when the odor generation source specifying unit 49 executes the process of step S1308 in FIG.
The odor concentration difference deriving step (means) is realized, for example, when the odor generation source specifying unit 49 executes the process of step S1310 of FIG.
<Claims 3 and 9>
The second molecular concentration distribution deriving step (means) is realized, for example, when the second molecular concentration distribution deriving unit 48 executes the process of step S1406 in FIG.
The first odor concentration deriving step (means) is realized, for example, when the odor generation source specifying unit 49 executes the process of step S1407 in FIG.
<Claims 6 and 12>
Information indicating the degree of contribution to the environmental load molecule detected at the environmental load molecule detection location is attached to a plurality of odor generation sources in order, for example, from the one having a great influence on the odor generated at the odor detection location. It is realized by the information of the ranking. In addition, the information indicating the degree of contribution to the environmental load molecule detected at the environmental load molecule detection location is also realized by, for example, information indicating the ratio of the odor generation source contributing to the odor generated at the odor detection location. Is done.

The embodiment of the present invention described above can be realized by a computer executing a program. Further, a computer-readable recording medium in which the program is recorded and a computer program product such as the program can also be applied as an embodiment of the present invention. As the recording medium, for example, a flexible disk, a hard disk, an optical disk, a magneto-optical disk, a CD-ROM, a magnetic tape, a nonvolatile memory card, a ROM, or the like can be used.
In addition, the embodiments of the present invention described above are merely examples of implementation in carrying out the present invention, and the technical scope of the present invention should not be construed as being limited thereto. Is. That is, the present invention can be implemented in various forms without departing from the technical idea or the main features thereof.

Next, examples of the present invention will be described.
NOx and benzene are molecules that can employ a modification of the odor measurement method and odor measurement system described in the first embodiment.
Among NOx having a high environmental load, NO molecules and benzene, which is required to have an annual average of 3 μg / m ³ or less by the basic environmental method, are effective molecules to which the present invention can be applied. Normally, NOx that is evaluated only by the total amount can be evaluated for each molecule by applying the present invention. Moreover, the odor measuring system in each embodiment has high sensitivity, and can measure benzene having a concentration of 1 μg / m ³ or less with an average of 10 seconds. Therefore, by using the odor measurement system in the first embodiment, it is possible to apply the odor measurement system as atmospheric environmental load monitoring of NO and benzene. An example is shown.

Here, the factory chimney A and the factory chimney B at the location indicated by a star (☆) at the location indicated by a circle (◯) in FIG. 15 are used as NO, benzene, and toluene emission sources, and NO and The diffusion of benzene was simulated. From the chimney A, NO molecules were controlled and discharged in a range of 1 ppm to 100 ppm, benzene molecules from 100 ppb to 10 ppm, and toluene molecules from 100 ppb to 10 ppm. The discharge ratio of NO, benzene, and toluene from the chimney A was 1: 2.5: 3. On the other hand, from the chimney B, NO molecules are controlled to be discharged in a range of 1 ppm to 100 ppm, benzene molecules are controlled in a range of 100 ppb to 10 ppm, and toluene molecules are controlled in a range of 100 ppb to 10 ppm. The discharge ratio of NO, benzene, and toluene from the chimney B was 4: 1: 1. The amount of exhaust gas from the chimneys A and B hardly changed, the amount of exhaust gas from the chimney A was 1500 m ³ / min, and the amount of exhaust gas from the chimney B was 2300 m ³ / min.

FIG. 15 shows that the height of the chimney A is 13 m, the height of the chimney B is 20 m, the wind direction is north wind, the wind speed is 3 m / s, and the height is 2 m according to the actual factory layout conditions and weather conditions. It is a figure which shows the result of the simulation which shows the spreading | diffusion of NO molecule | numerator in the position. Similar simulations were performed with benzene and toluene, and similar results were obtained with only the difference in concentration.
On the other hand, NO molecules, benzene, and toluene were measured at an actual site height of 2 m located at the square mark (□) in FIG. The wind direction at the time of measurement was north wind, the wind speed was 3 m / s, the measurement result of NO molecules was 4 ppb, the measurement result of benzene molecules was 5.5 ppb, and the measurement result of toluene molecules was 6.5 ppb.

  In the simulation, NO molecules are 4ppb, benzene molecules are 5.5ppb, and toluene molecules are 6.5ppb at the square marks (□). The emission effects of chimney A (○) and chimney B (☆) The degree was found to be a 5: 1 ratio case. Moreover, when the generation | occurrence | production density | concentration was calculated from this result by simulation, the emission density | concentration of each molecule | numerator from the chimney A was NO 1ppm, benzene 2.5ppm, and toluene 3ppm. On the other hand, the discharge concentration of each molecule from the chimney B was 8 ppm for NO, 2 ppm for benzene, and 2 ppm for toluene.

  Actually, the emission concentration of each molecule from the chimneys A and B at the time of operation is consistent with the simulation value, and the amount from the source of NO, benzene, and toluene is predicted by this system that combines actual measurement and simulation. It was clarified that it was appropriate. However, the example described here is an example, and is not limited to this example.

DESCRIPTION OF SYMBOLS 10 One box car 20 Battery box 30 Odor measurement apparatus 40 Odor generation | occurrence | production source analysis apparatus 41 Gas mass information acquisition part 42 Mass spectrum preparation part 43 Molecular concentration derivation part 44 First molecular concentration distribution derivation part 45 Data storage part 46 Data extraction part 47 Odor generating source selection unit 48 Second molecular concentration distribution deriving unit 49 Odor generating source specifying unit 50 Output unit 51 to 53 Power cord 60 Gas introducing unit 70 Power connector 100a to 100d Odor generating source 110 Odor detection point

Claims (17)

An environmental load molecule generation source evaluation method for deriving an index for identifying an environmental load molecule generation source that affects an environmental load molecule detected at an environmental load molecule detection location,
By irradiating with vacuum ultraviolet laser light, a plurality of molecules contained in the gas are ionized at once with one photon, and the mass spectrum of each molecule is measured at each of the plurality of environmental load molecule generation sources. A first environmental impact molecule measurement step;
A first molecular concentration derivation step for deriving the concentration of each molecule in the plurality of environmental load molecule generation sources based on the mass and spectral intensity of each molecule measured in the first environmental load molecule measurement step;
By solving the simulation calculation model formula of numerical fluid dynamics using the concentration and environmental conditions of each molecule in the plurality of environmental load molecule generation sources derived by the first molecular concentration deriving step as boundary conditions, a predetermined 3 A first molecular concentration distribution deriving step in which deriving a molecular concentration distribution that is a concentration distribution of each molecule in an analysis region that is a dimension region is performed with different boundary conditions;
The concentration of each molecule at each position in the analysis region obtained by the molecular concentration distribution derived by the first molecular concentration distribution deriving step, and the plurality of environments when the molecular concentration distribution is derived A storage step of storing the concentration of each molecule in the load molecule generation source and at least a part of the environmental conditions when the molecular concentration distribution is derived in a storage medium as a result of simulation;
A second environment in which a plurality of molecules contained in a gas are ionized at once with one photon by irradiation with vacuum ultraviolet laser light, and the mass of each molecule is measured at the environmental load molecule detection site. Loading molecule measurement process;
A second molecular concentration derivation step for deriving the concentration of each molecule at the environmental load molecule detection site based on the mass of each molecule measured in the second environmental load molecule measurement step;
Derived by the second molecular concentration derivation step, the concentration of each molecule at the environmental load molecule detection location, the environmental load molecule detection location, and the environmental load molecule detected at the environmental load molecule detection location An extraction process for extracting the concentration of each molecule in the plurality of environmental load molecule generation sources corresponding to the environmental condition that can be regarded as an environmental condition from the result of the simulation;
Based on the concentration of each molecule in the plurality of environmental load molecule generation sources extracted in the extraction step, an environmental load molecule generation source that affects the environmental load molecules detected at the environmental load molecule detection location is identified. An index deriving step for deriving an index for
An output step of outputting the index derived by the derivation step;
A method for evaluating an environmental load molecule generation source, comprising: The first molecular concentration distribution deriving step includes:
The relative ratio of the concentration of each molecule in the environmental load molecule generation source is assumed to be equal to the relative ratio of the concentration of each molecule in the environmental load molecule generation source obtained in the first molecular concentration derivation step. The environmental load molecule generation source evaluation method according to claim 1, wherein the molecular concentration distribution is derived with different boundary conditions.   The indicator for specifying the environmental load molecule generation source that affects the environmental load molecule detected at the environmental load molecule detection location is for specifying the odor generation source that affects the odor detected at the odor detection location. 3. The environmental load molecule generation source evaluation method according to claim 1 or 2, which is an index. The index derivation step includes
The environmental conditions stored as a result of the simulation corresponding to the concentration of each molecule in the plurality of odor generation sources extracted in the extraction step, and the plurality of odor generation sources extracted in the extraction step By solving the simulation calculation model of numerical fluid dynamics using the concentration of each molecule at the odor generation source excluding one odor generation source among the concentration of each molecule of A second molecular concentration distribution deriving step for deriving a molecular concentration distribution which is a concentration distribution by sequentially selecting a combination of odor generating sources excluding one of the plurality of odor generating sources;
A concentration of each molecule at the odor detection location obtained from the molecular concentration distribution derived by the second molecular concentration distribution deriving step, and an olfactory sensitivity coefficient obtained by quantifying the degree of influence on the odor of each molecule. A first odor concentration derivation step for individually deriving an odor concentration, which is a sum of values multiplied for each same molecule, for each combination of odor generation sources excluding one of the plurality of odor generation sources;
The concentration of each molecule at the odor detection location stored as a result of the simulation corresponding to the concentration of each molecule at the plurality of odor generation sources extracted by the extraction step, and the odor of each molecule A second odor concentration deriving step for deriving an odor concentration that is a sum of values obtained by multiplying the olfactory sensitivity coefficient obtained by quantifying the degree of the influence on the same molecule for each molecule;
The difference between the odor concentration derived by the second odor concentration derivation step and the odor concentration derived by the first odor concentration derivation step is calculated as an odor generation source excluding one of the plurality of odor generation sources. An odor concentration difference deriving step for individually deriving each of the combinations;
The index for identifying the odor generation source that affects the odor detected at the odor detection location is derived based on the difference derived by the odor concentration difference deriving step. Of environmental impact molecule generation source evaluation method. The index derivation step includes
The environmental conditions stored as a result of the simulation corresponding to the concentration of each molecule in the plurality of odor generation sources extracted in the extraction step, and the plurality of odor generation sources extracted in the extraction step By solving the simulation calculation model expression of numerical fluid dynamics using the concentration of each molecule at one odor generation source of the concentration of each molecule in the boundary condition as a boundary condition, the distribution of the concentration of each molecule in the analysis region A second molecular concentration distribution deriving step of deriving a molecular concentration distribution by sequentially selecting one of the plurality of odor generating sources;
A concentration of each molecule at the odor detection location obtained from the molecular concentration distribution derived by the second molecular concentration distribution deriving step, and an olfactory sensitivity coefficient obtained by quantifying the degree of influence on the odor of each molecule. A first odor concentration derivation step for individually deriving an odor concentration that is a sum of values multiplied for each same molecule for each of the plurality of odor generation sources,
An index for deriving an odor generation source that affects an odor detected at the odor detection location is derived based on the odor concentration derived in the first odor concentration deriving step. 4. The method for evaluating an environmental load molecule generation source according to 3. The first molecular concentration derivation step and the second molecular concentration derivation step include:
6. The method according to claim 1, wherein at least one of an average value of the concentration of each molecule at a predetermined time and a maximum value and a minimum value of the concentration of each molecule at a predetermined time are derived. The environmental load generation source evaluation method according to claim 1.   The environmental load molecule generation source evaluation method according to any one of claims 1 to 6, wherein the ionization potential of the vacuum ultraviolet laser light is 10.5 eV.   An index for identifying an environmental load molecule generation source that affects the environmental load molecule detected at the environmental load molecule detection location is detected at the environmental load molecule detection location among the plurality of environmental load molecule generation sources. Information indicating the environmental load molecule generation source that most influences the environmental load molecule and the degree of contribution of each of the plurality of environmental load molecule generation sources to the environmental load molecule detected at the environmental load molecule detection location The environmental load molecule generation source evaluation method according to any one of claims 1 to 7, comprising at least one of information and information. A measuring device for measuring the mass of molecules contained in the gas;
An analysis device for deriving an index for identifying a source of an environmental load molecule that affects the environmental load molecule detected at the environmental load molecule detection point;
An environmental impact molecular source evaluation system comprising:
The measuring device is
By irradiating vacuum ultraviolet laser light, ionizing a plurality of molecules contained in a gas at once with one photon and measuring the mass spectrum of each molecule, Performing with the environmental load molecule detection point,
The analysis device includes:
The concentration of each molecule obtained based on the mass and spectral intensity of each molecule measured by the measuring device, wherein the concentration of each molecule in the plurality of environmental load molecule generation sources and the environmental condition are boundary conditions. As described above, the molecular concentration distribution that is the distribution of the concentration of each molecule in the analysis region that is the predetermined three-dimensional region is derived by solving the simulation calculation model expression of the numerical fluid dynamics as follows. 1 molecular concentration distribution deriving means;
The concentration of each molecule at each position in the analysis region obtained by the molecular concentration distribution derived by the first molecular concentration distribution deriving means, and the plurality of environments when the molecular concentration distribution is derived Storage means for storing the concentration of each molecule in the load molecule generation source and at least a part of the environmental conditions when the molecular concentration distribution is derived in a storage medium as a result of simulation;
The concentration of each molecule obtained based on the mass of each molecule measured by the measuring device, the concentration of each molecule at the environmental load molecule detection site, the environmental load molecule detection site, and the environmental load molecule An environmental condition that can be regarded as an environmental condition when an environmental load molecule is detected at a detection location; and an extraction means for extracting the concentration of each molecule in the plurality of environmental load molecule generation sources corresponding to the environmental condition from the simulation result;
Based on the concentration of each molecule at the plurality of environmental load molecule generation sources extracted by the extraction means, the environmental load molecule generation source that affects the environmental load molecules detected at the environmental load molecule detection location is specified. Index deriving means for deriving an index for
Output means for outputting the index derived by the index deriving means;
An environmental load molecule generation source evaluation system characterized by comprising: The first molecular concentration distribution deriving means includes:
The boundary condition is varied assuming that the relative ratio of the concentration of each molecule in the environmental load molecule generation source is equal to and constant with the relative ratio of the concentration of each molecule in the environmental load molecule generation source obtained in the measurement apparatus. The molecular load distribution evaluation system according to claim 9, wherein the molecular concentration distribution is derived.   The indicator for specifying the environmental load molecule generation source that affects the environmental load molecule detected at the environmental load molecule detection location is for specifying the odor generation source that affects the odor detected at the odor detection location. The environmental load molecule generation source evaluation system according to claim 9 or 10, wherein the system is an index. The index derivation means includes
The environmental conditions stored as a result of the simulation corresponding to the concentration of each molecule in the plurality of odor generation sources extracted by the extraction means, and the plurality of odor generation sources extracted by the extraction means By solving the simulation calculation model of numerical fluid dynamics using the concentration of each molecule at the odor generation source excluding one odor generation source among the concentration of each molecule of Second molecular concentration distribution deriving means for deriving a molecular concentration distribution which is a concentration distribution by sequentially selecting a combination of odor generating sources excluding one of the plurality of odor generating sources;
A concentration of each molecule at the odor detection location obtained from the molecular concentration distribution derived by the second molecular concentration distribution deriving means, and an olfactory sensitivity coefficient obtained by quantifying the degree of influence of each molecule on the odor First odor concentration deriving means for individually deriving an odor concentration that is a sum of values multiplied for each same molecule for each combination of odor generation sources excluding one of the plurality of odor generation sources;
The concentration of each molecule at the odor detection location stored as a result of the simulation corresponding to the concentration of each molecule at the plurality of odor generating sources extracted by the extraction means, and the odor of each molecule Second odor concentration deriving means for deriving an odor concentration that is a sum of values obtained by multiplying the olfactory sensitivity coefficient obtained by quantifying the degree of the influence on the same molecule for each same molecule;
The difference between the odor concentration derived by the second odor concentration deriving means and the odor concentration derived by the first odor concentration deriving means is calculated as an odor generation source excluding one of the plurality of odor generation sources. Odor concentration difference deriving means for individually deriving each combination, and
12. An index for identifying an odor generation source that affects an odor detected at the odor detection location is derived based on the difference derived by the odor concentration difference deriving unit. Environmental impact molecular source evaluation system. The index derivation means includes
The environmental conditions stored as a result of the simulation corresponding to the concentration of each molecule in the plurality of odor generating sources extracted by the extracting means, and the plurality of odor generating sources extracted by the extracting means By solving the simulation calculation model expression of numerical fluid dynamics using the concentration of each molecule at one odor generation source of the concentration of each molecule in the boundary condition as a boundary condition, the distribution of the concentration of each molecule in the analysis region Second molecular concentration distribution deriving means for deriving a molecular concentration distribution by sequentially selecting one of the plurality of odor generating sources;
A concentration of each molecule at the odor detection location obtained from the molecular concentration distribution derived by the second molecular concentration distribution deriving means, and an olfactory sensitivity coefficient obtained by quantifying the degree of influence of each molecule on the odor First odor concentration deriving means for individually deriving an odor concentration that is a sum of values multiplied for each same molecule for each of the plurality of odor generating sources,
An index for identifying an odor generation source that affects an odor detected at the odor detection location is derived based on the odor concentration derived by the first odor concentration deriving means. The environmental load molecule generation source evaluation system according to 11. First molecular concentration deriving means for deriving the concentration of each molecule in the plurality of environmental load molecule generation sources based on the mass of each molecule measured by the measuring device;
Second molecular concentration deriving means for deriving the concentration of each molecule at the environmental load molecule detection site based on the mass of each molecule measured by the measuring device;
The first molecular concentration deriving means and the second molecular concentration deriving means are:
14. The method according to claim 9, wherein at least one of an average value of the concentration of each molecule at a predetermined time and a maximum value and a minimum value at a predetermined time of the concentration of each molecule are derived. The environmental load molecule generation source evaluation system according to claim 1.   The environmental load molecule generation source evaluation system according to any one of claims 9 to 14, wherein an ionization potential of the vacuum ultraviolet laser light is 10.5 eV.   An index for identifying an environmental load molecule generation source that affects the environmental load molecule detected at the environmental load molecule detection location is detected at the environmental load molecule detection location among the plurality of environmental load molecule generation sources. Information indicating the environmental load molecule generation source that most influences the environmental load molecule and the degree of contribution of each of the plurality of environmental load molecule generation sources to the environmental load molecule detected at the environmental load molecule detection location The environmental load molecule generation source evaluation system according to any one of claims 9 to 15, comprising at least one of information and information. A computer program for causing a computer to derive an index for identifying an environmental load molecule generation source that affects an environmental load molecule detected at an environmental load molecule detection location,
By irradiating vacuum ultraviolet laser light, a plurality of molecules contained in the gas are ionized at once with a single photon, and the mass and spectrum intensity of each molecule measured by a measuring device that measures the mass spectrum of each molecule A concentration of each molecule derived on the basis of each of the plurality of environmental load molecule generation sources and the environmental condition as a boundary condition. Deriving a molecular concentration distribution that is a distribution of the concentration of each molecule in an analysis region that is a three-dimensional region of
The concentration of each molecule at each position in the analysis region obtained by the molecular concentration distribution derived by the first molecular concentration distribution deriving step, and the plurality of environments when the molecular concentration distribution is derived A storage step of storing the concentration of each molecule in the load molecule generation source and at least a part of the environmental conditions when the molecular concentration distribution is derived in a storage medium as a result of simulation;
The concentration of each molecule derived based on the mass of each molecule measured by the measuring device, the concentration of each molecule at the environmental load molecule detection site, the environmental load molecule detection site, and the environmental load molecule An extraction process for extracting the concentration of each molecule in the plurality of environmental load molecule generation sources corresponding to the environmental condition that can be regarded as an environmental condition when an environmental load molecule is detected at a detection location;
Based on the concentration of each molecule in the plurality of environmental load molecule generation sources extracted in the extraction step, an environmental load molecule generation source that affects the environmental load molecules detected at the environmental load molecule detection location is identified. An index deriving step for deriving an index for
An output step of outputting the indicator derived by the indicator derivation step;
A computer program for causing a computer to execute.