JP2013024659A - Odor diffusion simulation method and odor diffusion simulation device - Google Patents

Abstract

PROBLEM TO BE SOLVED: To calculate odor concentration more faithful to reality in consideration that components and component ratio of an odorant change by a chemical reaction in a process that the odorant diffuses.SOLUTION: In a composite odor concentration DB 24, component ratio, air diffusion concentration, and odor concentration of component materials of various composite odors are accumulated by being associated with one another. An advective diffusion concentration calculation part 31 calculates the air diffusion concentration of discharged materials assuming that the discharged materials advect and diffuse without causing a chemical reaction; a chemical reaction concentration attenuation amount calculation part 32 calculates the air diffusion concentration of the discharged materials attenuated by a fact that the discharged materials cause the chemical reaction during advection; a composite odorous component calculation part 33 calculates the component ratio and the air diffusion concentration of the composite odor consisting of products by the chemical reaction and the remaining discharged materials; a composite odor concentration calculation part 34 performs multiple regression analysis for data having the same components as those of the composite odor from the composite odor concentration DB 24, calculates a prediction formula of the odor concentration of the composite odor, and calculates the odor concentration of the composite odor on the basis of the prediction formula.

Description

  The present invention relates to an odor diffusion simulation method and an odor diffusion simulation apparatus for simulating diffusion of odors discharged into the atmosphere.

  In order to take measures against the spread of odors (especially offensive odors) to the outside, or to grasp the odor diffusion state, it is necessary to carry out an odor diffusion simulation. Local governments have established standards for odor index at site boundaries from the viewpoint of preventing malodors, and facilities of a certain scale are obliged to meet the standards, and the impact on the environment must be evaluated.

  Substances that cause odor (bad odor) are different from ordinary air pollutants, and many substances whose concentration and components change due to chemical reaction. For this reason, by adding and synergizing various kinds of odorous substances, a bad odor (a bad odor due to a composite odor) that is strongly felt in human olfaction may occur. Therefore, in the odor diffusion simulation, it is particularly necessary to take into account such a bad odor due to the composite odor.

  Patent Document 1 discloses a technique for measuring an odor (bad odor) substance concentration and evaluating the measured result with an odor intensity that enables a sensuous capture of the measured result. However, in Patent Document 1, no consideration is given to the influence of secondary products and decomposition substances caused by the chemical reaction of the combined malodorous substances, or the change in the component ratio of the malodorous substances.

  In addition, simulation methods related to environmental assessment regarding bad odor are shown in “Odor Index No. 2 Regulation Standard” (for example, see Non-Patent Document 1) and “Guidelines for Investigation of Impact on Living Environment” (for example, see Non-Patent Document 2). Has been. In these non-patent documents 1 and 2, the plume / puff model, which is the same diffusion calculation method as that of general air pollutants, is used for the prediction of malodorous substance diffusion. In the plume / puff model, the concentration of air pollutants can be calculated by performing relatively simple numerical calculations.

  In addition, in the simulation methods shown in Non-Patent Documents 1 and 2, the parameters of the plume / puff model are adjusted in order to take into account the characteristics of odor (bad odor), and the composite odor is evaluated by human senses. Odor concentration is used. However, since the composite odor is evaluated as diffusion of the odor concentration obtained from the human sense, the concentration change or the component change due to the chemical change of the actual malodorous substance is not taken into consideration.

  As an atmospheric dispersion simulation model that takes into account chemical changes in substances, there is "ADMER (Atmospheric Dispersion Model for Exposure and Risk Assessment)" developed by the National Institute of Advanced Industrial Science and Technology (see, for example, Non-Patent Document 3). In “ADMER”, diffusion prediction is performed in consideration of concentration attenuation due to chemical change of a substance. Concentration decay is a concept based on the idea that a decomposition coefficient is set for a substance for which diffusion prediction is performed, and the rate at which the substance is decomposed by a chemical reaction increases with the progress of the diffusion prediction time of the substance. However, in “ADMER”, only the concentration decay of a single substance is calculated, and the influence of secondary generated substances and decomposed substances due to chemical changes is not considered.

Japanese Patent No. 4374723

Ministry of the Environment Water and Atmospheric Environment Bureau, Atmosphere Living Environment Office, “Understanding Odor Index Regulation No. 2 Standards-Regulation Standards of Offensive Odors at Gas Exhaust Ports of Chimneys, etc.”, [Search June 20, 2011], Internet <URL : Http://www.env.go.jp/air/akushu/kisei/index.html> Notification of Ministry of the Environment, “Guidelines for Survey on Impacts on Living Environment of Waste Management Facilities”, September 4, 2006, [Search June 20, 2011], Internet <URL: http: //www.env.go. jp / hourei / syousai.php? id = 11000364> National Institute of Advanced Industrial Science and Technology, “ADMER2: Atmospheric Dispersion Model for Exposure and Risk assessment”, [Search June 20, 2011], Internet <URL: http://www.aist-riss.jp/software /admer/en/index_en.html

  As described above, in the conventional odor diffusion simulation, the influence of a substance that is chemically changed in the process of malodorous substance diffusion and is secondarily generated is hardly taken into consideration. However, malodor is usually a complex odor generated from a plurality of odorous substances. The odor concentration of the composite odor is not always the sum of the odor concentrations of the odors generated from each odor substance. Since the odor has an interaction, the odor concentration of the composite odor becomes larger or smaller depending on the kind, combination, component ratio, etc. of the odor substance. Therefore, in the odor diffusion simulation, it is required to calculate the values of the components of the odor substance and the component ratio as close to reality as possible.

  Therefore, an object of the present invention is to consider an odor diffusion simulation method capable of calculating a more realistic odor concentration in consideration of the fact that components and component ratios change due to chemical changes in the process of odor substance diffusion. And providing an odor diffusion simulation apparatus.

  The odor diffusion simulation method and apparatus according to the present invention is a composite odor concentration data configured by associating the component ratio of each component substance of the composite odor, the atmospheric diffusion concentration, and the odor concentration with respect to each of the plurality of composite odors. An odor diffusion simulation method and apparatus for simulating the diffusion of odors discharged from an emission source using an arithmetic processing apparatus that can access a composite odor concentration database that stores therein, (1) the arithmetic processing apparatus includes: Based on a calculation model that inputs data and emission source data, and the emission substances that constitute the odor specified in the emission source data advect and diffuse in the atmosphere without a chemical reaction, at a calculation target point away from the emission source Advection diffusion concentration calculation processing for calculating the atmospheric diffusion concentration and advection time of the emission material, and (2) the calculated advection time The residual concentration, which is the atmospheric diffusion concentration of the emission material remaining when the emission material chemically reacts, is calculated, and the residual concentration is subtracted from the atmospheric diffusion concentration of the emission material calculated in the advection diffusion concentration calculation process. A chemical reaction concentration attenuation amount calculation process for calculating the difference amount as a concentration attenuation amount of the atmospheric diffusion concentration of the emission material due to the chemical reaction; and (3) calculating a component ratio of the products generated by the chemical reaction. The atmospheric diffusion concentration of the product is calculated based on the calculated component ratio and the concentration attenuation amount, and the emission material is calculated based on the calculated atmospheric diffusion concentration of the product and the residual concentration of the emission material. A composite odor component calculation process for calculating a component ratio and an atmospheric diffusion concentration of the composite odor generated by the advection diffusion and chemical reaction of (4) the composite odor concentration database; From the accumulated composite odor data, the composite odor data of the composite odor having the same components as the composite odor generated by the advection diffusion and chemical reaction of the emission is extracted, and the multiple regression analysis using the extracted composite odor data Substituting the component ratio and atmospheric diffusion concentration of the composite odor produced by the advection diffusion and chemical reaction of the exhaust into the prediction formula for predicting the odor concentration from the component ratio and atmospheric diffusion concentration of the composite odor obtained by And a composite odor concentration calculation process for calculating an odor concentration of the composite odor at the calculation target point.

  According to the present invention, it is possible to perform an odor diffusion simulation capable of calculating a more realistic odor concentration in consideration of changes in the components and component ratio due to chemical changes in the process of odor substance diffusion. .

The figure which showed the example of the structure of the odor diffusion simulation apparatus which concerns on embodiment of this invention. The figure which showed the example of the whole processing flow of the odor diffusion simulation which concerns on embodiment of this invention. The figure which showed the example of the processing flow of the chemical reaction concentration attenuation amount calculation process by the chemical reaction concentration attenuation amount calculation part. The figure which showed the example of the structure of chemical reaction density | concentration attenuation | damping rate DB. The figure which showed the mode of the density | concentration attenuation of the discharge | emission substance by a chemical reaction with the relationship graph of residual concentration and advection time. The figure which showed the example of the processing flow of the composite odor component calculation process by a composite odor component calculation part. The figure which shows the example of a structure of the product component data used by the composite odor component calculation process of FIG. 6, and memorize | stored in DB for work. The figure which showed the example of the structure of chemical reaction decomposition | disassembly rate and synthesis rate DB. The figure which showed the example of the detailed process flow of the chemical reaction product component ratio calculation process in the composite odor component calculation process of FIG. The figure which showed the example of the detailed process flow of step S672 in the chemical reaction product component ratio calculation process of FIG. The figure which showed the example of the detailed process flow of step S673 in the chemical reaction product component ratio calculation process of FIG. The figure which showed the example of the processing flow of the composite odor density | concentration calculation process by a composite odor density | concentration calculation part. The figure which showed the example of the structure of composite odor density DB.

  Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

  FIG. 1 is a diagram showing an example of the configuration of an odor diffusion simulation apparatus according to an embodiment of the present invention. As shown in FIG. 1, an odor diffusion simulation apparatus 100 includes an input device 10 composed of a keyboard and a mouse, a storage device 20 composed of a semiconductor memory, a hard disk device, etc., and an arithmetic processing device 30 composed of a microprocessor integrated circuit. And a display device 40 composed of a liquid crystal display or the like.

  On the storage device 20, a database such as a work DB (Database) 21, a chemical reaction concentration decay rate DB 22, a chemical reaction decomposition rate / synthesis rate DB 23, and a composite odor concentration DB 24 is configured. Above, functional blocks such as an advection diffusion concentration calculation unit 31, a chemical reaction concentration attenuation amount calculation unit 32, a composite odor component calculation unit 33, and a composite odor concentration calculation unit 34 are configured.

  Here, the work DB 21 configured on the storage device 20 stores various data used for calculation processing such as prediction conditions, weather data, and emission source conditions input via the input device 10, and data of calculation processing results. It is a database to memorize. The chemical reaction concentration decay rate DB 22 is a database in which the concentration decay rate, chemical reaction formula, chemical reaction formula composition rate, residual concentration calculation formula, and the like of exhaust substances obtained from experimental data and literature information are accumulated in advance. .

  The chemical reaction decomposition rate / synthesis rate DB 23 includes chemical reaction formulas of substances obtained from experimental data and literature information, substance decomposition rates, substance synthesis rates, substance determination concentrations, chemical reaction formula composition ratios, and the like. This is a database accumulated in advance. The composite odor concentration DB 24 is a database in which component ratios of composite odors obtained from experimental data, literature information, etc., atmospheric diffusion concentrations of composite odors, odor concentrations of composite odors, and the like are accumulated in advance.

  On the other hand, the advection diffusion concentration calculation unit 31 configured on the arithmetic processing unit 30 uses the input data (prediction conditions, weather data, emission source conditions, etc.) stored in the work DB 21 as the emission discharged from the emission source. Applied to a calculation model (for example, plume / puff model) in which the substance advects and diffuses in the atmosphere without chemical reaction, and calculates the advective diffusion of the emitted substance, and the atmospheric diffusion of the emitted substance at a predetermined calculation target point Calculate the concentration and advection time.

  In addition, the chemical reaction concentration attenuation amount calculation unit 32 calculates a residual concentration that is the atmospheric diffusion concentration of the emission material remaining after the emission material chemically reacts with the atmosphere or the like during the advection time calculated by the advection diffusion concentration calculation unit 31. To do. Further, the chemical reaction concentration attenuation amount calculation unit 32 calculates the difference amount obtained by subtracting the calculated residual concentration from the atmospheric diffusion concentration of the emission material calculated by the advection diffusion concentration calculation unit 31, and the atmospheric diffusion concentration of the emission material due to the chemical reaction. It is calculated as the density attenuation amount.

  In addition, the composite odor component calculation unit 33 calculates the component ratio of the product resulting from the chemical reaction of the discharged substances, using the data accumulated in the chemical reaction decomposition rate / synthesis rate DB 23. Then, based on the calculated component ratio and the concentration attenuation amount calculated by the chemical reaction concentration attenuation amount calculation unit 32, the atmospheric diffusion concentration of the product due to the chemical reaction is calculated. Furthermore, the composite odor component calculation unit 33 determines the component ratio and the component odor of the composite odor generated by the advection diffusion of the exhaust material and the chemical reaction based on the atmospheric diffusion concentration of the product by the chemical reaction and the residual concentration of the exhaust material. Calculate the atmospheric diffusion concentration.

  In addition, the composite odor concentration calculation unit 34 uses the composite odor data stored in the composite odor concentration DB 24 to provide composite odor data for the composite odor having the same components as the composite odor generated by advection diffusion and chemical reaction of emissions. Is extracted, and multiple regression analysis is performed using the extracted composite odor data, and a prediction formula for predicting the odor concentration is obtained from the component ratio of the composite odor and the atmospheric diffusion concentration. Then, the component ratio and atmospheric diffusion concentration of the composite odor generated by the advection diffusion and chemical reaction of the exhaust gas calculated by the composite odor component calculation unit 33 are substituted into the prediction formula, and the composite odor at the calculation target point is calculated. Calculate the odor concentration.

  Further, an odor concentration distribution creating unit (not shown) configured on the arithmetic processing device 30 is based on the odor concentration distribution at a number of calculation target points in a predetermined area calculated by using the composite odor concentration calculating unit 34. A diagram is created, and the created odor concentration distribution map is displayed on the display device 40 in a superimposed manner on a topographic map or a terrain map.

  Note that the functions of these functional blocks are realized by the arithmetic processing device 30 executing a predetermined program stored in the storage device 20 in advance. That is, the odor diffusion simulation apparatus 100 is realized by a so-called computer including the arithmetic processing device 30 and the storage device 20. In that case, the computer may be a plurality of computers connected to each other via a network or the like. Each database configured on the storage device 20 may be included in a separate computer.

  Further, the arithmetic processing unit 30 is not limited to a so-called CPU (Central Processing Unit) of a computer, and the advection diffusion concentration calculating unit 31, the chemical reaction concentration attenuation calculating unit 32, the complex odor component calculating unit 33, As long as the function of the complex odor concentration calculation unit 34 or the like is realized, a dedicated control circuit configured by an FPGA (Field Programmable Gate Array) or the like may be used.

  Further, the work DB 21, the chemical reaction concentration decay rate DB 22, the chemical reaction decomposition rate / synthesis rate DB 23, and the composite odor concentration DB 24 stored on the storage device 20 are respectively configured on the storage device 20 of different computers. It doesn't matter. Even if these databases are configured on the storage device 20 of the same computer as the arithmetic processing device 30 or are stored on the storage device 20 of another computer via a network, the processing unit 30 It only needs to be accessible.

  Subsequently, FIG. 2 is a diagram illustrating an example of the entire processing flow of the odor diffusion simulation according to the embodiment of the present invention. The odor diffusion simulation is realized by each functional block configured on the arithmetic processing device 30.

  As shown in FIG. 2, the arithmetic processing device 30 first uses the input device 10 to predict prediction condition data for odor diffusion simulation (latitude, longitude, X direction prediction range, Y direction prediction of the reference point of the prediction target area). The range, the X direction prediction mesh number, the Y direction prediction mesh number, etc.) are input (step S1) and stored in the work DB 21. Next, the arithmetic processing unit 30 inputs weather data (wind direction, wind speed, atmospheric stability, etc.) via the input device 10 (step S2) and stores it in the work DB 21. Next, the arithmetic processing unit 30 inputs the odor substance emission source data (e.g., the position of the emission source, the shape of the outlet, the discharge amount, the temperature, the discharge speed, and the discharge substance) via the input device 10 (step S3). And stored in the work DB 21.

  Subsequently, the arithmetic processing device 30 executes an advection diffusion concentration calculation process as a function of the advection diffusion concentration calculation unit 31 (step S4). That is, the arithmetic processing unit 30 uses the plume / puff model (see Non-Patent Document 1 or the like) based on the prediction condition data, weather data, and emission source data stored in the work DB 21 (hereinafter, referred to as odor concentration calculation target point). The atmospheric diffusion concentration and the advection time of the emission material at the target point are simply calculated, and the atmospheric diffusion concentration and the advection time are stored in the work DB 21. In the plume / puff model, the exhausted material is merely advected and diffused, and it is not assumed that the exhausted material chemically reacts with the atmosphere or the like.

  Next, the arithmetic processing unit 30 executes a chemical reaction concentration attenuation amount calculation process as a function of the chemical reaction concentration attenuation amount calculation unit 32 (step S5). That is, the arithmetic processing unit 30 reads the atmospheric diffusion concentration of the emission material calculated in the advection diffusion concentration calculation process of step S4 and the concentration of the emission material read from the chemical reaction concentration attenuation rate DB 22 for the emission material included in the emission source data. Attenuation rate data, chemical reaction formula, chemical reaction formula composition rate, and residual concentration calculation formula are used to calculate the concentration attenuation amount of the exhausted substance due to the chemical reaction at the target point and save it in the work DB 21. A more detailed processing flow of this chemical reaction concentration attenuation amount calculation processing will be separately described with reference to FIG.

  Next, the arithmetic processing unit 30 performs a composite odor component calculation process as a function of the composite odor component calculation unit 33 (step S6). That is, the arithmetic processing unit 30 calculates the concentration attenuation amount of the discharged substance calculated in the chemical reaction concentration attenuation amount calculation process in step S5, and the chemical reaction formula, decomposition rate, and synthesis of the substance read from the chemical reaction decomposition rate / synthesis rate DB23. Using the rate data, the component ratio and the atmospheric diffusion concentration of the composite odor generated by the advection and chemical reaction of the emission material are calculated and stored in the work DB 21. A more detailed processing flow of the composite odor component calculation process will be described separately with reference to FIG.

  Next, the arithmetic processing unit 30 executes a composite odor concentration calculation process as a function of the composite odor concentration calculator 34 (step S7). That is, the arithmetic processing unit 30 performs multiple regression analysis on the composite odor data including the component ratio of the composite odor accumulated in the composite odor concentration DB 24, the atmospheric diffusion concentration, and the composite odor concentration, and the obtained multiple regression analysis is performed. The composite odor component ratio and the atmospheric diffusion concentration calculated in the composite odor component calculation process in step S6 are substituted into the prediction formula to calculate the composite odor concentration, and the calculated composite odor concentration is stored in the work DB 21. A more detailed processing flow of the composite odor concentration calculation process will be described separately with reference to FIG.

  Next, details of the chemical reaction concentration attenuation amount calculation processing by the chemical reaction concentration attenuation amount calculation unit 32 will be described with reference to FIGS. Here, FIG. 3 is a diagram showing an example of the processing flow of the chemical reaction concentration attenuation amount calculation processing by the chemical reaction concentration attenuation amount calculation unit 32, and FIG. 4 shows an example of the configuration of the chemical reaction concentration attenuation rate DB22. FIG. 5 and FIG. 5 are graphs showing the state of decay of the concentration of exhausted substances due to a chemical reaction, in a relational graph between residual concentration and advection time.

  As shown in FIG. 3, in the chemical reaction concentration attenuation amount calculation process (step S5 in FIG. 2), the arithmetic processing unit 30 first uses the emission source data in step S3 (see FIG. 2) from the chemical reaction concentration attenuation rate DB22. Each data of the chemical reaction formula, chemical reaction formula composition rate, concentration decay rate, and residual concentration calculation formula associated with the designated emission material is read (step S51).

  Here, as shown in FIG. 4, in the chemical reaction concentration decay rate DB 22, a plurality of chemical reaction formulas related to the emission material are associated with the emission material, and further, for each chemical reaction equation, The concentration attenuation rate, chemical reaction composition rate, residual concentration calculation formula, and the like are associated with each other. In the example of FIG. 4, only NOx is exemplified as the emission material, and it is shown that three types of chemical reactions proceed with respect to the NOx.

  Next, in FIG. 3, the arithmetic processing unit 30 calculates the chemical reaction formula of the emission substance read out in step S51 and the atmospheric diffusion concentration of the emission substance calculated by the plume / puff model, the advection time (see step S4 in FIG. 2). Then, using the concentration attenuation rate and residual concentration calculation formula (see FIG. 4), the residual concentration after the chemical reaction of the emission material is calculated (step S52).

Here, assuming that the atmospheric diffusion concentration of the emission material is C ₀ , the concentration decay rate is k _d , and the advection time is t, for example, the residual concentration of the emission material NO _{x at} the time t due to a chemical reaction of NO _x → NO ₂ + O _2. C _t is represented by the following formula (1). Note that C ₀ ≧ 0, C _t ≧ 0, and t ≧ 0.
C _t = C ₀ · exp (−k _d · t) (1)

The equation (1) the relationship between the residual concentration C _t and advection time t after the density attenuation by, for example, represented by the graph shown in FIG. That is, as the advection time t elapses, a chemical reaction of NO _x → NO ₂ + O ₂ proceeds, and the residual concentration of NO _x decreases. On the other hand, the concentrations of NO ₂ and O ₂ produced by this chemical reaction increase (not shown).

  Referring to FIG. 3 again, the arithmetic processing unit 30 determines the residual concentration after the chemical reaction according to each chemical reaction equation calculated in step S52 and the chemical reaction equation composition rate read from the chemical reaction concentration attenuation rate DB 22 in step S51. Is used to calculate the residual concentration of the emission material at the target point (step S53) and store it in the work DB 21.

  Next, the arithmetic processing unit 30 calculates the difference between the atmospheric diffusion concentration of the emission material at the target point calculated in step S3 and the residual concentration after the chemical reaction of the emission material at the target point calculated in step S53. The concentration attenuation amount of the discharged substance is calculated (step S54) and stored in the work DB 21.

  Subsequently, the details of the composite odor component calculation process by the composite odor component calculation unit 33 will be described with reference to FIGS. 6 is a diagram showing an example of the processing flow of the composite odor component calculation process by the composite odor component calculation unit 33, and FIG. 7 is used in the composite odor component calculation process of FIG. 6 and is stored in the work DB 21. FIG. 8 is a diagram showing an example of the configuration of the chemical reaction decomposition rate / synthesis rate DB 23.

  As shown in FIG. 6, in the composite odor component calculation process, the arithmetic processing unit 30 first sets an initial value of product component data composed of a product used in the process, its component ratio, and its atmospheric diffusion concentration. (See FIG. 7A) and stored in the work DB 21 (step S61). That is, in the product component data, the initial value of the product is the emission material itself, the initial value of the component ratio is represented by the component ratio of the emission material, and the atmospheric diffusion concentration is calculated by the chemical reaction concentration attenuation calculation process. It is represented by the calculated concentration attenuation amount of the emission material (see FIG. 3, step S54). Incidentally, in FIG. 7A, it is assumed that the emission material is only NOx, and the initial value of the component ratio is 1.00.

  Next, the arithmetic processing unit 30 reads product component data (product, component ratio, atmospheric diffusion concentration) from the work DB 21 (step S62). Further, the arithmetic processing unit 30 obtains, from the chemical reaction decomposition rate / synthesis rate DB 23, the product name included in the product component data read in step S62, the substance name of the substance related to the decomposition and synthesis of the product, the chemical reaction The equation and determination density are read (step S63).

  Here, as shown in FIG. 8, the chemical reaction decomposition rate / synthesis rate DB 23 includes the substance name of the substance related to the decomposition and synthesis of the substance, the determination concentration of the substance (minimum value of the concentration to be considered as odor), decomposition And a synthetic chemical reaction formula, a decomposition rate of a substance generated by decomposition, a synthesis rate of a substance used for a synthesis reaction, and the like.

  Next, in FIG. 6, the arithmetic processing unit 30 has an atmospheric diffusion concentration larger than the determination concentration in the data read from the chemical reaction decomposition rate / synthesis rate DB 23 in step S 63, and a decomposition or synthesis chemical reaction equation. It is determined whether or not there is a substance having (step S64).

  And when the substance applicable to the determination conditions exists (it is Yes at step S64), the arithmetic processing unit 30 selects the target substance which changes a component ratio by chemical reaction from the applicable substance ( In step S65), one chemical reaction formula for decomposition or synthesis of the target substance is selected from the data read from the chemical reaction decomposition rate / synthesis rate DB 23 in step S63 (step S66).

  Next, the arithmetic processing unit 30 calculates the component ratio of the product changed by the chemical reaction based on the selected chemical reaction formula (step S67). The component ratio calculation process will be described in detail separately with reference to FIGS. 9 to 11. In this process, product component data (product, component, etc.) stored in the work DB 21 is included. Ratio, atmospheric diffusion concentration) is updated as appropriate.

  Next, the arithmetic processing unit 30 determines whether another chemical reaction formula exists for the decomposition or synthesis of the target substance in the data read from the chemical reaction decomposition rate / synthesis rate DB 23 in step S63. (Step S68) If another chemical reaction formula exists (Yes in Step S68), the process returns to Step S66, and in Step S66, a chemical reaction formula that has not been selected in the previous processing is selected, and Step S67 is selected. The following processing is executed again. If another chemical reaction formula does not exist (No in step S68), the arithmetic processing unit 30 returns to step S62, and executes the processing in step S62 and subsequent steps again.

  On the other hand, in the determination in step S64, if no substance corresponding to the determination condition exists (No in step S64), the product component data (product, component ratio, atmospheric diffusion) stored in the work DB 21 at that time Concentration) and the residual concentration of the exhaust material calculated in the chemical reaction concentration attenuation calculation process (see step S53 in FIG. 3), including the advective and diffused exhaust material and the product generated by the chemical reaction of the exhaust material. The component ratio and atmospheric diffusion concentration of the odorous substance are calculated (step S69).

  In other words, the odorous substance mentioned here includes exhausted substances remaining after advection and diffusion up to this point, and products resulting from chemical reaction of the emitted substances with the atmosphere. Further, the product, the component ratio, and the atmospheric diffusion concentration of the product resulting from the chemical reaction of the emission material are represented by the product component data stored in the work DB 21 at this time (example in FIG. 7B). See).

  FIG. 9 is a diagram showing an example of a detailed process flow of the chemical reaction product component ratio calculation process (FIG. 6, step S67) in the composite odor component calculation process of FIG. In FIG. 9, the arithmetic processing unit 30 first determines whether or not the chemical reaction state of the chemical reaction formula selected in step S66 is “decomposition” (step S671), and the chemical reaction state is “decomposition”. In this case (Yes in step S671), the component ratio of the product after the decomposition reaction of the target substance is calculated (step S672). On the other hand, when the chemical reaction state is not “decomposition”, that is, “synthesis”. (No in step S671), the component ratio of the product after the synthesis reaction of the target substance is calculated (step S673). An example of a more detailed processing flow of steps S672 and S673 will be separately described in detail with reference to FIGS.

  Next, the processor 30 reads the atmospheric diffusion concentration and the component ratio of all the products from the work DB 21 that stores the product component data reflecting the calculation results in steps S672 and S673 (step S674). It is determined whether or not the sum of the atmospheric diffusion concentrations of all the read products is equal to the concentration attenuation amount calculated in step S5 (step S675).

  As a result of the determination in step S675, when the sum of the atmospheric diffusion concentrations of all products is not equal to the concentration attenuation amount (No in step S675), the arithmetic processing unit 30 determines the component ratio of the products read in step S674. Is used to correct the atmospheric diffusion concentration of the product so that the sum of the atmospheric diffusion concentration of the product and the concentration attenuation calculated in step S5 are equal (step S676), and the work DB 21 (product component data) And the process of FIG. 9 is terminated. If the result of determination in step S675 is that the sum of atmospheric diffusion concentrations of all products is equal to the concentration attenuation amount (Yes in step S675), the processing in FIG.

  FIG. 10 is a diagram showing an example of a detailed processing flow of step S672 in the chemical reaction product component ratio calculation processing of FIG. In FIG. 10, the arithmetic processing unit 30 first reads out the decomposition rate from the chemical reaction equation decomposition rate / synthesis rate DB 23 for the chemical reaction equation selected in step S66 (see FIG. 6) (step S6721). Next, the arithmetic processing unit 30 calculates the component ratio and the atmospheric diffusion concentration of the product after the decomposition reaction from the decomposition rate of the read chemical reaction formula and the atmospheric diffusion concentration of the target substance (step S6722).

  Next, the processor 30 determines whether or not the product after the decomposition reaction is included in the product component data of the work DB 21 (step S6723), and if included (in step S6723). Yes), the component ratio and atmospheric diffusion concentration of the product in the product component data of the work DB 21 are updated using the component ratio and atmospheric diffusion concentration of the product after the decomposition reaction calculated in step S6722 (step S6724). ).

  On the other hand, when it is determined in step S6723 that the product after the decomposition reaction is not included in the product component data of the work DB 21 (No in step S6723), the arithmetic processing unit 30 uses the product for work. Register in the product component data of DB21 (step S6725). Next, the arithmetic processing unit 30 uses the component ratio and atmospheric diffusion concentration of the product after the decomposition reaction calculated in step S6722 to use the component ratio and atmospheric diffusion concentration of the product in the product component data of the work DB 21. Is updated (step S6726).

  FIG. 11 is a diagram showing an example of a detailed processing flow of step S673 in the chemical reaction product component ratio calculation processing of FIG. In FIG. 11, the arithmetic processing unit 30 first reads the synthesis rate from the chemical reaction rate decomposition rate / synthesis rate DB 23 for the chemical reaction formula selected in step S66 (see FIG. 6) (step S6731). From the work DB 21 (product component data), the component ratio and atmospheric diffusion concentration of the substances used for the chemical reaction of the synthesis by the selected chemical reaction formula are read (step S6732). Next, the arithmetic processing unit 30 synthesizes the chemical reaction formula read out in step S6731, the atmospheric diffusion concentration of the target substance, the component ratio of the substance used for the synthesis read out in step S6732, and the atmospheric diffusion concentration. The component ratio and atmospheric diffusion concentration of the product after the reaction are calculated (step S6733).

  Next, the arithmetic processing unit 30 determines whether or not the product after the synthesis reaction is included in the product component data of the work DB 21 (step S6734), and if included (in step S6734). Yes), using the component ratio and atmospheric diffusion concentration of the product after the synthesis reaction calculated in step S6733, the component ratio and atmospheric diffusion concentration of the product in the product component data of the work DB 21 are updated (step S6735). ).

  On the other hand, if it is determined in step S6734 that the product after the synthesis reaction is not included in the product component data in the work DB 21 (No in step S6734), the arithmetic processing unit 30 uses the product for work. Register in the product component data of DB21 (step S6736). Next, the arithmetic processing unit 30 uses the component ratio and atmospheric diffusion concentration of the product after the synthesis reaction calculated in step S6733, and the component ratio and atmospheric diffusion concentration of the product in the product component data of the work DB 21. Is updated (step S6737).

  Next, with reference to FIG. 12 and FIG. 13, the details of the composite odor concentration calculation process (see FIG. 2, step S7) by the composite odor concentration calculator 34 will be described. Here, FIG. 12 is a diagram showing an example of the processing flow of the composite odor concentration calculation process by the composite odor concentration calculator 34, and FIG. 13 is a diagram showing an example of the configuration of the composite odor concentration DB 24.

  First, the arithmetic processing unit 30 extracts composite odor concentration data of a composite gas composed of the same components as the composite odor obtained in the composite odor component calculation process (FIG. 2, step S6) from the composite odor concentration DB 24 ( Step S71).

  Here, as shown in FIG. 13, the composite odor concentration DB 24 accumulates composite odor concentration data in which various composite gases (composite odors) are associated with the component ratio, atmospheric diffusion concentration, and odor concentration of the composite gas. Database. Note that the composite odor concentration data stored in the composite odor concentration DB 24 may be experimental data obtained through experiments, or may be literature data described in published literature. Furthermore, calculation result data calculated by the odor diffusion simulation apparatus 100 according to the present embodiment may be used. In addition, the odor concentration of the composite gas referred to here is usually a numerical value determined by a person such as an odor determination person.

  In the example of FIG. 13, only one example of the composite gas is described for each of the experimental data, the literature data, and the calculation result data. However, there are many types of the composite gas, and the composite gas has the same component. However, a large number of components having different component ratios and atmospheric diffusion concentrations are accumulated.

  Therefore, in the process of step S71 of FIG. 12, the composite odor in the present embodiment (that is, a gas containing a discharge substance discharged from the discharge source and a product in which the discharge substance is chemically changed) from the composite odor concentration DB 24. Compound odor concentration data for a complex gas (composite odor) having the same components as is extracted.

  Therefore, the arithmetic processing unit 30 applies multiple regression analysis to the composite odor concentration data (consisting of the component ratio of the composite gas (composite odor), the atmospheric diffusion concentration, and the odor concentration) extracted in step S71, A multiple regression prediction formula for obtaining the odor concentration from the component ratio and the atmospheric diffusion concentration is derived (step S72).

Here, the multiple regression analysis is a technique for predicting the value of another variable (objective variable) from the distribution values of a plurality of variables (explanatory variables), where the explanatory variables are X _{1 to} X _n and the objective variable is Y. At this time, the value of the objective variable Y is predicted by Expression (2), which is a basic expression of multiple regression analysis (hereinafter referred to as multiple regression prediction expression).
Y = b ₀ + b ₁ · X ₁ + b ₂ · X ₂ +... + B _n · X _n (2)

In the multiple regression prediction formula (2), b ₀ , b ₁ ,..., B _n are called partial regression coefficients and are obtained by a known general multiple regression analysis process. Therefore, in the case of the present embodiment, the component ratio and the atmospheric diffusion concentration of the composite odor correspond to the explanatory variables X _{1 to} X _n , and the odor concentration corresponds to the objective variable Y. Therefore, if the component ratio and the atmospheric diffusion concentration of the composite odor are known, the odor concentration of the composite odor can be calculated by the multiple regression prediction equation (2).

  Therefore, the arithmetic processing unit 30 substitutes the composite odor component ratio and the atmospheric diffusion concentration calculated in the composite odor component calculation process (FIG. 2, step S6) into the multiple regression prediction formula (2), thereby The odor concentration is calculated (step S73). Next, the arithmetic processing unit 30 registers the composite odor data including the component ratio of the composite odor, the atmospheric diffusion concentration and the odor concentration calculated by the above processing in the composite odor concentration DB 24 (step S74).

  As described above, in the present embodiment, since the composite odor data including the component ratio, the atmospheric diffusion concentration, and the odor concentration of the composite odor calculated by the odor diffusion simulation apparatus 100 is registered in the composite odor concentration DB 24, Data that can be used in the regression analysis will increase sequentially, and accordingly, the accuracy of the odor concentration of the complex odor predicted by the multiple regression prediction equation (2) is improved.

  The arithmetic processing unit 30 performs processing for displaying the display device 40 by combining the odor concentration at the target point obtained as described above with the topographic map or the topographic map. Omitted.

  As described above, according to the present embodiment, the emission substance (odor substance) discharged from the emission source attenuates the atmospheric diffusion concentration due to the chemical reaction with the atmosphere while being advected and diffused, and the chemical reaction The product component ratio and atmospheric diffusion concentration are calculated, and the odor concentration of the resulting composite odor is calculated. Therefore, in the present embodiment, it is possible to obtain a highly accurate odor concentration that is closer to reality than in a conventional plume / puff model.

10 Input device 20 Storage device 21 Work DB
22 Chemical reaction concentration decay rate DB
23 Chemical reaction decomposition rate and synthesis rate DB
24 Compound odor concentration DB
30 Arithmetic processor 31 Advection diffusion concentration calculation part (Advection diffusion concentration calculation means)
32 Chemical reaction concentration attenuation calculation section (Chemical reaction concentration attenuation calculation means)
33 Compound odor component calculation part (Compound odor component calculation means)
34 Complex odor concentration calculation unit (Compound odor concentration calculation means)
40 display device 100 odor diffusion simulation device

Claims (6)

Computation that allows access to a complex odor concentration database that accumulates complex odor concentration data that is configured by associating the component ratio of the complex odor components, atmospheric diffusion concentration, and odor concentration with each of multiple complex odors An odor diffusion simulation method for simulating diffusion of odor discharged from an emission source using a processing device,
The arithmetic processing unit includes:
A calculation target point that is distant from the emission source based on a calculation model in which meteorological data and emission source data are input, and emission substances constituting the odor specified by the emission source data are advected and diffused in the atmosphere without a chemical reaction Advection diffusion concentration calculation processing for calculating the atmospheric diffusion concentration and advection time of the emission material in
A residual concentration, which is an atmospheric diffusion concentration of the exhaust material remaining when the exhaust material chemically reacts during the calculated advection time, is calculated, and the atmospheric diffusion concentration of the exhaust material calculated in the advection diffusion concentration calculation process is used to calculate the residual concentration. A chemical reaction concentration attenuation amount calculation process for calculating a difference amount obtained by subtracting the residual concentration as a concentration attenuation amount of the atmospheric diffusion concentration of the emission material due to the chemical reaction;
The component ratio of the product generated by the chemical reaction is calculated, the atmospheric diffusion concentration of the product is calculated based on the calculated component ratio and the concentration attenuation amount, and the calculated atmosphere of the product is calculated. Based on the diffusion concentration and the residual concentration of the emission material, a composite odor component calculation process for calculating the component ratio and the atmospheric diffusion concentration of the composite odor generated by advection diffusion and chemical reaction of the emission material,
From the composite odor data accumulated in the composite odor concentration database, the composite odor data of the composite odor having the same components as the composite odor generated by the advection diffusion and chemical reaction of the emission is extracted, and the extracted composite odor data The prediction formula for predicting the odor concentration from the component ratio of the composite odor and the atmospheric diffusion concentration obtained by the multiple regression analysis using the odor, the component ratio of the composite odor generated by the advection diffusion and chemical reaction of the emission and the atmosphere Complex odor concentration calculation processing for substituting diffusion concentration and calculating the odor concentration of the complex odor at the calculation target point;
An odor diffusion simulation method characterized in that The arithmetic processing unit includes:
For each of a plurality of substances, the chemical reaction formula in which the substance reacts in the atmosphere, the residual concentration calculation formula for calculating the atmospheric diffusion concentration of the substance remaining by the chemical reaction, and the residual concentration calculation formula are used. The chemical reaction concentration decay rate database that stores the chemical reaction concentration decay rate data that is configured by associating the concentration decay rate with the chemical reaction composition rate of the chemical reaction formula for the substance is further configured to be accessible. And
In the chemical reaction concentration attenuation amount calculation process, when calculating the residual concentration of the emission material, the chemical reaction concentration attenuation rate data for the emission material is read from the chemical reaction concentration attenuation rate database, and the chemical reaction concentration attenuation rate is read. 2. The odor diffusion simulation method according to claim 1, wherein the residual concentration of the emission material is calculated using a chemical reaction formula, a residual concentration calculation formula, a concentration decay rate, and a chemical reaction composition rate included in the data. . The arithmetic processing unit includes:
For each of a plurality of substances, decomposition rate data including a chemical reaction formula of the decomposition reaction of the substance and a decomposition rate of the substance generated by the decomposition reaction, and synthesis by the substance generated by the decomposition reaction Chemical reaction decomposition rate, which is a combination of chemical reaction decomposition rate and synthesis rate data composed by associating the chemical reaction formula of the reaction and the synthesis rate data including the synthesis rate of the product generated by the synthesis reaction The composition rate database is further configured to be accessible,
In the composite odor component calculation process, when calculating the component ratio of the product, the chemical reaction decomposition rate / synthesis rate data associated with the emission is read from the chemical reaction decomposition rate / synthesis rate database, and the chemical The decomposition rate of the product generated by the decomposition reaction of the emission material included in the reaction decomposition rate / synthesis rate data, and the generation rate of the product generated by the synthesis reaction by the product generated by the decomposition reaction The odor diffusion simulation method according to claim 2, wherein a component ratio of the product is calculated. For each of a plurality of composite odors, a composite odor concentration database that stores composite odor concentration data configured by associating the component ratio of the component substances of the composite odor, the atmospheric diffusion concentration, and the odor concentration is provided. An odor diffusion simulation device for simulating the diffusion of odor discharged from
A calculation target point that is distant from the emission source based on a calculation model in which meteorological data and emission source data are input, and emission substances constituting the odor specified by the emission source data are advected and diffused in the atmosphere without a chemical reaction Advection diffusion concentration calculating means for calculating the atmospheric diffusion concentration and advection time of the emission material in
The residual concentration, which is the atmospheric diffusion concentration of the exhaust material remaining when the exhaust material chemically reacts during the calculated advection time, is calculated, and the atmospheric diffusion concentration of the exhaust material calculated by the advection diffusion concentration calculating means is used to calculate the residual concentration. A chemical reaction concentration attenuation amount calculating means for calculating a difference amount obtained by subtracting the residual concentration as a concentration attenuation amount of the atmospheric diffusion concentration of the emission material due to the chemical reaction;
The component ratio of the product generated by the chemical reaction is calculated, the atmospheric diffusion concentration of the product is calculated based on the calculated component ratio and the concentration attenuation amount, and the calculated atmosphere of the product is calculated. A composite odor component calculation means for calculating a component ratio and an atmospheric diffusion concentration of a composite odor generated by advection diffusion and chemical reaction of the exhaust material based on a diffusion concentration and a residual concentration of the exhaust material;
From the composite odor data accumulated in the composite odor concentration database, the composite odor data of the composite odor having the same components as the composite odor generated by the advection diffusion and chemical reaction of the emission is extracted, and the extracted composite odor data The prediction formula for predicting the odor concentration from the component ratio of the composite odor and the atmospheric diffusion concentration obtained by the multiple regression analysis using the odor, the component ratio of the composite odor generated by the advection diffusion and chemical reaction of the emission and the atmosphere A composite odor concentration calculating means for substituting a diffusion concentration and calculating an odor concentration of the composite odor at the calculation target point;
An odor diffusion simulation apparatus characterized by comprising: For each of a plurality of substances, the chemical reaction formula in which the substance reacts in the atmosphere, the residual concentration calculation formula for calculating the atmospheric diffusion concentration of the substance remaining by the chemical reaction, and the residual concentration calculation formula are used. A chemical reaction concentration attenuation rate database that accumulates chemical reaction concentration attenuation rate data configured to correspond to the concentration attenuation rate and the chemical reaction composition rate of the chemical reaction formula for the substance;
The chemical reaction concentration decay amount calculating means is:
When calculating the residual concentration of the emission material, the chemical reaction concentration attenuation rate data for the emission material is read from the chemical reaction concentration attenuation rate database, and the chemical reaction formula and the residual concentration included in the chemical reaction concentration attenuation rate data are read. The odor diffusion simulation apparatus according to claim 4, wherein the residual concentration of the emission substance is calculated using a calculation formula, a concentration decay rate, and a chemical reaction composition rate. For each of a plurality of substances, decomposition rate data including a chemical reaction formula of the decomposition reaction of the substance and a decomposition rate of the substance generated by the decomposition reaction, and synthesis by the substance generated by the decomposition reaction Chemical reaction decomposition rate, which is a combination of chemical reaction decomposition rate and synthesis rate data composed by associating the chemical reaction formula of the reaction and the synthesis rate data including the synthesis rate of the product generated by the synthesis reaction A synthesis rate database is further provided.
The complex odor component calculating means is:
When calculating the component ratio of the product, the chemical reaction decomposition rate / synthesis rate data associated with the emission is read from the chemical reaction decomposition rate / synthesis rate database, and the chemical reaction decomposition rate / synthesis rate data is read. A component of the product using a decomposition rate of a product generated by a decomposition reaction of the contained emission substance and a generation rate of a product generated by a synthesis reaction by a product generated by the decomposition reaction The odor diffusion simulation apparatus according to claim 5, wherein the ratio is calculated.

Images