Carbon dioxide dynamic simulation device, computer program, and recording medium

Abstract

And the dynamic state of carbon dioxide in the agricultural greenhouse is mastered. The present invention is provided with: an absorption amount calculation unit (11) for determining the net photosynthesis amount of a plant cultivated in an agricultural greenhouse as the carbon dioxide absorption amount of the plant; and a utilization rate calculation unit (12) that determines, as the utilization rate of the carbon dioxide supplied to the agricultural greenhouse, the ratio of the absorption amount to the amount of carbon dioxide supplied from the carbon dioxide supply source to the interior of the agricultural greenhouse. The absorption amount and utilization rate of carbon dioxide supplied to an agricultural greenhouse can be obtained, and the dynamic state of carbon dioxide can be clearly grasped.

Description

Technical Field

[0001] The present invention relates to a technology for simulating the dynamics of carbon dioxide in an agricultural greenhouse. Background Art

[0002] Countermeasures to global warming are being studied extensively as a pressing issue. In particular, reducing greenhouse gas emissions such as carbon dioxide from power plants, substations, factories, and waste treatment plants is expected to have a significant impact on global warming countermeasures.

[0003] Patent Document 1 discloses a combined waste incineration system. This system consists of an incineration system equipped with an incinerator for incinerating wastes such as municipal garbage, sewage sludge, and industrial waste, and a methane fermentation system equipped with a methane fermentation tank for methane fermentation of wet waste. Methane gas and carbon dioxide are recovered and separated, and the separated carbon dioxide is reused within the system, for example, for cleaning ash generated in the incineration system, thereby conserving resources.

[0004] Patent Document 2 discloses a technology for preventing the release of hot exhaust gas from a thermal power plant containing a large amount of carbon dioxide into the atmosphere. Instead, the hot exhaust gas containing carbon dioxide is cooled to a given temperature by heat exchange using seawater, and then supplied to a plant factory, thereby suppressing the release of greenhouse gases into the atmosphere.

[0005] Patent document 3 discloses a technology for separating carbon dioxide from waste gas emitted from thermal power plants, waste incineration plants, etc. using a water membrane with carbon dioxide absorption capacity. It also records that the separated carbon dioxide is used for plant cultivation, thereby effectively utilizing carbon dioxide without releasing it into the atmosphere, thereby achieving a solution to global environmental problems.

[0006] Non-Patent Document 1 discloses, as one of the technologies for preventing global warming, a technology that supplies heat generated by combustion of gas or the like to agricultural greenhouses for heating, and extracts and supplies carbon dioxide from exhaust gas to promote plant growth.

[0007] Prior art literature

[0008] Patent Literature

[0009] Patent Document 1: Japanese Patent Application Laid-Open No. 2006-212524

[0010] Patent Document 2: Japanese Patent Application Laid-Open No. 3-236723

[0011] Patent Document 3: Japanese Patent Application Laid-Open No. 2021-133314

[0012] Non-patent literature

[0013] Non-Patent Literature 1: ENERGY FRONTLINE, April 13, 2021, Vol. 25, "Rapid growth of tomatoes using CO2! What's the latest technology that kills three birds with one stone?" URL: https: / / ene-fro.com / article / ef195_a1 / Summary of the Invention

[0014] -Problems to be solved by the invention-

[0015] The technologies disclosed in Patent Documents 1-3 and Non-Patent Document 1 fundamentally involve recovering carbon dioxide and heat emitted from various energy-generating facilities and utilizing them for other purposes, thereby suppressing the release of greenhouse gases into the atmosphere. Patent Documents 2-3 and Non-Patent Document 1 share a common approach of suppressing carbon dioxide emissions into the atmosphere by supplying carbon dioxide to plants.

[0016] Even if the carbon dioxide absorption capacity of plants is utilized, to effectively suppress carbon dioxide emissions into the atmosphere, it is necessary to construct agricultural greenhouses that can secure an appropriate cultivation area. However, the larger the cultivation area, the higher the cost of building an agricultural greenhouse. Therefore, if a greenhouse is to be constructed, it is desirable to secure a cultivation area large enough to maximize the effective utilization of carbon dioxide emitted from energy-emitting equipment.

[0017] However, the above-mentioned conventional technologies have only suggested the concept of consuming carbon dioxide through plant cultivation, and have not studied the extent to which carbon dioxide supplied to the agricultural greenhouse is absorbed by the plants.

[0018] Due to these circumstances, the method of preventing global warming by utilizing agricultural greenhouses is not yet widely used. In the case of building agricultural greenhouses as a measure to prevent global warming, if it is possible to predict in advance the extent of carbon dioxide reduction that can be achieved by using recyclable waste gas in the agricultural greenhouses, it will serve as an incentive for the construction of new agricultural greenhouses. In this case, if the effect of carbon dioxide reduction can be grasped in terms of monetary value, it can be expected that the effect will be further increased. Of course, if it is also possible to grasp the amount of carbon dioxide reduction and the economic effect when utilizing waste gas from nearby energy emission equipment in existing agricultural greenhouses, it is expected that the utilization of waste gas will be further promoted.

[0019] Meanwhile, in agricultural greenhouses, to promote photosynthesis, carbon dioxide generators are installed to supply carbon dioxide. Exhaust gas from heating systems is also being utilized. In other words, even without utilizing exhaust gas from other facilities like factories, CO2 is actively supplied using fossil fuel-powered systems to promote photosynthesis. However, it's difficult to numerically determine the extent to which plant photosynthesis is promoted by the supplied CO2, and this practice remains limited for practical application.

[0020] The present invention is made in view of the above situation, and its purpose is to provide a carbon dioxide dynamic simulation device, a computer program and a recording medium that can grasp the dynamics of carbon dioxide in agricultural greenhouses, and promote the recycling of carbon dioxide in the exhaust gas of power plants, factories, etc., which helps to reduce greenhouse gases.

[0021] -Methods for solving the problem-

[0022] In order to solve the above-mentioned problems, the carbon dioxide dynamic simulation device of the present invention simulates the dynamics of carbon dioxide in an agricultural greenhouse, and comprises: an absorption amount calculation unit, which calculates the net photosynthesis amount of plants cultivated in the agricultural greenhouse as the absorption amount of carbon dioxide by the plants; and a utilization rate calculation unit, which calculates the ratio of the absorption amount in a given period to the supply amount of carbon dioxide supplied from the carbon dioxide supply source to the agricultural greenhouse as the utilization rate of the carbon dioxide supplied to the agricultural greenhouse.

[0023] Preferably, the absorption amount calculation unit calculates the net photosynthesis amount using the total daily sunlight amount, the light transmittance of the agricultural greenhouse, the light receiving rate of the plant, and the indoor carbon dioxide concentration of the agricultural greenhouse.

[0024] Preferably, the indoor carbon dioxide concentration is calculated based on a ventilation rate of air between indoors and outdoors of the agricultural greenhouse, and the ventilation rate is calculated using a difference between outdoor specific enthalpy and indoor specific enthalpy of the agricultural greenhouse.

[0025] Preferably, the absorption amount calculation unit includes: a greenhouse structure information acquisition unit, which acquires the structure information of the agricultural greenhouse including the light transmittance; a plant information acquisition unit, which acquires information of the plant of the cultivation object including the light reception rate; and a meteorological information acquisition unit, which acquires meteorological information of the area where the agricultural greenhouse is set up.

[0026] Preferably, the carbon dioxide dynamic simulation device comprises: a first greenhouse gas emission amount calculation unit, which calculates the first greenhouse gas emission amount when the supply amount of carbon dioxide supplied to the agricultural greenhouse is supplied by a device that generates carbon dioxide by burning fossil fuels as the carbon dioxide supply source; a second greenhouse gas emission amount calculation unit, which calculates the second greenhouse gas emission amount when the supply amount of carbon dioxide supplied to the agricultural greenhouse is supplied by carbon dioxide recovered from exhaust gas discharged from an exhaust gas generating source serving as the carbon dioxide supply source; and a greenhouse gas reduction amount calculation unit, which calculates the difference between the first greenhouse gas emission amount and the second greenhouse gas emission amount as the greenhouse gas reduction amount when the exhaust gas is used.

[0027] Preferably, the first greenhouse gas emission amount calculation unit obtains a difference between the carbon dioxide supply amount and the carbon dioxide absorption amount as the first greenhouse gas emission amount.

[0028] Preferably, the second greenhouse gas emission amount calculation unit calculates the second greenhouse gas emission amount by setting a value corresponding to the carbon dioxide absorption amount to a negative number.

[0029] Preferably, the carbon dioxide dynamics simulation device includes a cost calculation unit that compares the carbon dioxide supply amount when the fossil fuel is used and the carbon dioxide supply amount when the exhaust gas is used to calculate a reduction in carbon dioxide supply cost.

[0030] In addition, the present invention provides a computer program that enables a computer to function as a carbon dioxide dynamic simulation device in an agricultural greenhouse, and the computer program enables the computer to perform the following steps: calculating the net photosynthesis amount of plants cultivated in the agricultural greenhouse as the amount of carbon dioxide absorbed by the plants; and calculating the ratio of the absorption amount within a given period to the supply amount of carbon dioxide supplied from the carbon dioxide supply source to the agricultural greenhouse as the utilization rate of the carbon dioxide supplied to the agricultural greenhouse.

[0031] Preferably, in the step of determining the carbon dioxide absorption amount, the net photosynthesis amount is determined using the total sunlight amount, the light transmittance of the agricultural greenhouse, the light reception rate of the plant, and the indoor carbon dioxide concentration of the agricultural greenhouse.

[0032] Preferably, the indoor carbon dioxide concentration is calculated based on a ventilation rate of air between indoors and outdoors of the agricultural greenhouse, and the ventilation rate is calculated using a difference between outdoor specific enthalpy and indoor specific enthalpy of the agricultural greenhouse.

[0033] Preferably, in the step of calculating the amount of carbon dioxide absorbed, the computer program causes the computer to perform the following steps: obtaining structural information of the agricultural greenhouse including the light transmittance; obtaining information of the plant of the cultivated object including the light reception rate; and obtaining meteorological information of the area where the agricultural greenhouse is set up.

[0034] Preferably, the computer program further causes the computer to perform the following steps: calculating the first greenhouse gas emission amount when the carbon dioxide supply amount supplied to the agricultural greenhouse is supplied by a device that generates carbon dioxide by burning fossil fuels as the carbon dioxide supply source; calculating the second greenhouse gas emission amount when the carbon dioxide supply amount supplied to the agricultural greenhouse is supplied by carbon dioxide recovered from exhaust gas discharged from an exhaust gas generating source that serves as the carbon dioxide supply source; and calculating the difference between the first greenhouse gas emission amount and the second greenhouse gas emission amount as the greenhouse gas reduction amount when the exhaust gas is utilized.

[0035] Preferably, the computer program further causes the computer to execute the step of comparing the carbon dioxide supply amount when the fossil fuel is used and the carbon dioxide supply amount when the exhaust gas is used to supply the carbon dioxide, thereby calculating a reduction in carbon dioxide supply cost.

[0036] The present invention also provides a computer-readable recording medium having the computer program recorded thereon. The recording medium storing the computer program may be any non-transitory recording medium. Non-transitory recording media are not particularly limited, and examples thereof include floppy disks, hard disks, CD-ROMs, MO (magneto-optical disks), DVD-ROMs, and memory cards.

[0037] -Effects of the Invention-

[0038] According to the present invention, the amount of carbon dioxide absorbed by plants in an agricultural greenhouse can be measured. Therefore, the utilization rate of carbon dioxide supplied to the agricultural greenhouse can be calculated based on the supply and absorption of carbon dioxide within a given period. As a result, the dynamics of the supply and absorption of carbon dioxide in the agricultural greenhouse, as well as the emission of carbon dioxide outside the agricultural greenhouse during ventilation, can be quantitatively understood. In addition, by using waste gas as a source of carbon dioxide and having it absorbed by plants, it is possible to contribute to the reduction of greenhouse gases, and by calculating the supply of waste gas, the absorption of carbon dioxide by plants, and the utilization rate, it is possible to quantitatively and economically understand the reduction effect, and to contribute to the promotion of the installation of agricultural greenhouses that utilize waste gas from factories, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1This is a block diagram showing a schematic configuration of a carbon dioxide dynamics simulation device according to one embodiment of the present invention.

[0040] Figure 2 It is a diagram showing an example of an input screen and an output screen displayed on a display.

[0041] Figure 3 This is a flowchart for explaining the calculation process of the daytime heat dissipation amount in the absorption amount calculation unit.

[0042] Figure 4 This is a flowchart for explaining the calculation process of the indoor carbon dioxide concentration in the absorption amount calculation unit.

[0043] Figure 5 This is a flowchart for explaining the calculation process of the net photosynthesis amount and the carbon dioxide absorption amount of the plant in the absorption amount calculation unit.

[0044] Figure 6 This is a flowchart for explaining the calculation process of the carbon dioxide supply amount.

[0045] Figure 7 This is a block diagram showing a schematic configuration of a carbon dioxide dynamics simulation device according to another embodiment of the present invention. DETAILED DESCRIPTION

[0046] Hereinafter, embodiments of the present invention will be described in further detail based on the accompanying drawings. Figure 1 This figure shows the schematic configuration of a carbon dioxide dynamics simulation device 1 according to one embodiment of the present invention. As shown in the figure, the carbon dioxide dynamics simulation device 1 according to this embodiment is composed of a computer (the type of computer is not limited and includes personal computers, microcomputers, portable information terminals, etc.) equipped with a processor (CPU) 1a and a storage unit (the term "storage unit" in this specification includes any volatile or non-volatile recording medium such as a main storage device or memory, and is not limited to any one type) 1b.

[0047] Specifically, in the CO2 dynamics simulation device 1 of this embodiment, a computer program that executes the steps for causing the computer of the CO2 dynamics simulation device 1 to function as the absorption amount calculation unit 11 and the utilization rate calculation unit 12 is stored in the storage unit 1b. The computer program is typically stored on a non-volatile recording medium such as a hard disk or SSD, internal or external to the computer (CO2 dynamics simulation device 1), and is read and executed by the processor 1a. Furthermore, various data may be stored in a storage unit connected via a communication line, in addition to the internal or external storage unit of the CO2 dynamics simulation device 1.

[0048] The absorption amount calculation unit 11 calculates the net photosynthesis amount of the plants grown in the agricultural greenhouse as the amount of carbon dioxide (CO2) absorbed by the plants. The net photosynthesis amount P is the amount of carbon dioxide actually absorbed by the plants and can be calculated based on the total daily sunlight (MJ / m 2 ) (Io), light transmittance of the covering material of the agricultural greenhouse (τ), light reception rate of the plant (φ), indoor carbon dioxide concentration of the agricultural greenhouse (indoor CO2 concentration: C in ), conversion coefficient to photosynthesis rate (m 3 / J)(μ), the factor involved in the reverse reaction (δ). Total sunlight amount per day (MJ / m 2 )(I o ) × light transmittance of the covering material of the agricultural greenhouse (τ) becomes the incident solar heat of the agricultural greenhouse (MJ / m 2 ).

[0049] Specifically, it can be obtained using the following simultaneous equations (1) to (4).

[0050] (1) Gains and losses of carbon dioxide (indoor and outdoor carbon dioxide exchange rate (E))

[0051] E=V(C in -C out )

[0052] (V: air ventilation volume, C out : atmospheric carbon dioxide concentration (can be a constant)

[0053] (2) Carbon dioxide absorption (P)

[0054] P=I o ·τ·φ·μ·C in / (C out +δ)

[0055] (3) Carbon dioxide supply (S)

[0056] S=P+E

[0057] (4) Heat budget equation

[0058] I o τ = k cover ·1 / β·(T in -T out )+V·(H in -H out )·ρ

[0059] (k cover : Heat transfer coefficient of the covering material, β: Insulation ratio (ratio of the ground area to the covering area), T in: Indoor temperature, T out : Outdoor temperature, H in : Indoor specific enthalpy, H out : outdoor specific enthalpy, ρ: air density)

[0060] The absorption amount calculation unit 11 includes a greenhouse information acquisition unit 111 , a plant information acquisition unit 112 , and a weather information acquisition unit 113 .

[0061] The greenhouse information acquisition unit 111 acquires structural information of the agricultural greenhouse to be simulated. The structural information includes the type of roof material of the agricultural greenhouse, the average heat dissipation coefficient associated with the light transmittance and the type of roof material, the designed ground area, the designed covered area, and the like. In this embodiment, the dimensions of the agricultural greenhouse, such as the designed ground area and the designed covered area, can be fixed dimensions for simulation, or can be roughly conceivable dimensions specified based on the area of ​​land on which construction is planned. In addition, it is also possible to simulate the designed ground area and the designed covered area of ​​an agricultural greenhouse of a size that can supply the heat required for heating a given plant for a given period of time, based on the premise of utilizing waste heat from factories, etc. in the area where construction is planned, and use this data. The structural information is stored in the structural information database 31, and the greenhouse information acquisition unit 111 accesses the structural information database 31 to obtain the required information.

[0062] The plant information acquisition unit 112 acquires information about the plant to be cultivated. The plant information includes the light receiving rate and the conversion coefficient (m 3 / J), and the daytime ventilation temperature (daytime ventilation temperature) suitable for the cultivation of the plant. This information is stored by plant type in the plant information database 32. Furthermore, the "ventilation temperature" is the daytime management temperature suitable for the plant to be cultivated and is used to calculate the predicted value of the daytime indoor temperature described later.

[0063] The weather information acquisition unit 113 acquires weather information for the area where the agricultural greenhouse is located. This weather information includes the area's daily average daytime temperature, daily maximum temperature, daily minimum temperature, daily total sunshine, and daily day length. This information is stored in the weather information database 33, corresponding to each area. The weather information database 33 can be created by importing necessary information from, for example, the area's regional weather observation system (AMeDAS). In this embodiment, the weather information acquisition unit 113 accesses the weather information database 33 to read this information. However, a mechanism that directly accesses the area's regional weather observation system and reads the information can also be used.

[0064] In addition, place names, plant types, types of agricultural greenhouse roof materials, and carbon dioxide (CO2) supply amounts are input using appropriate input devices and displayed on the Figure 2 The input screen of such a display is shown ( Figure 2 left half of the screen).

[0065] The utilization rate calculation unit 12 calculates the ratio of the carbon dioxide absorption amount calculated by the absorption amount calculation unit 11 to the carbon dioxide supply amount supplied from the carbon dioxide supply source to the agricultural greenhouse during a given period as the utilization rate. The carbon dioxide supply amount is set by inputting the supply amount assumed from the carbon dioxide supply source into the input screen as described above. The carbon dioxide absorption amount, supply amount, and utilization rate are output to the corresponding display. Figure 2 The right half of the screen is output.

[0066] The carbon dioxide supply source includes devices that generate carbon dioxide by burning fossil fuels, such as carbon dioxide generators, and, when utilizing exhaust gas, exhaust gas generation sources such as power plants, substations, factories, and garbage disposal plants.

[0067] Next, an example of simulation of carbon dioxide dynamics using the carbon dioxide dynamics simulation apparatus 1 of this embodiment will be described.

[0068] [Calculation of carbon dioxide absorption]

[0069] In this embodiment, the carbon dioxide absorption amount is the net photosynthesis amount, which can be calculated using the above formulas (1) to (4). An example of the calculation procedure of the carbon dioxide absorption amount using the above formulas (1) to (4) is described below.

[0070] (Acquisition of information)

[0071] The simulation operator inputs Figure 2 Enter the required information on the input screen shown ( Figure 3 As input items, the place name (e.g., Sapporo), the type of plant (e.g., tomato), and the type of roofing material (e.g., FRA) are input. When this information is input, the greenhouse information acquisition unit 111 accesses the structural information database 31 (S302), reads the average heat dissipation coefficient corresponding to the type of roofing material (e.g., 5), and obtains the designed ground area of ​​a given size agricultural greenhouse used for simulation (e.g., 1000m 2 ), designed coverage area (e.g., 1890.58m 2) (S303). In addition, the "designed ground area" is the design value of the ground area of ​​the planned agricultural greenhouse to be built, and in order to reduce the calculation load, it is a value calculated assuming a square ground shape, for example. The "designed covered area" is the design value of the area covered by the covering material after combining the side surface, end surface, top surface forming the roof, and tangent surface of the agricultural greenhouse. The "average heat dissipation coefficient" is a value obtained by dividing the heating heat of one night by the number of night temperature hours and the covering area for multiple types of films (roofing materials, lining covering materials (curtains)) during the experimental period, and averaging them according to the type category of the film (roofing materials, lining covering materials (curtains)). In this embodiment, three greenhouses A to C (surface area: 262m2) built at the Oyama Factory Farm of Seiwa Co., Ltd. from November 9, 1980 to February 3, 1981 were used. 2 , ground area: 120m 2 This value was obtained through an experiment in which four growing beds were set up in a greenhouse (with the following parameters: Outer film: Agricultural vinyl film; Heating method: Warm air heater; Curtain type: Greenhouse A... Single layer of agricultural vinyl film; Greenhouses B and C... Biaxial double layer; Curtain position (lower layer): Greenhouses A and B... 1.8m (from the ground); Greenhouse C... 1.85m (from the ground); Interlayer distance for double layer: 20cm). In the experiment, each growing bed was covered with black film only, and nozawa, lettuce, green peas, and lentils were planted in each growing bed. The average heating load factor for each greenhouse was calculated for each film type. When the heating load factor of Greenhouse A was set to 100%, the percentage of the heating load factor for Greenhouses B and C, depending on the film combination, was calculated. This value was defined as the "average heat dissipation factor."

[0072] The plant information acquisition unit 112 accesses the plant information database 32 ( S304 ), and acquires the daytime ventilation temperature (23° C.) corresponding to the input cultivation target plant “tomato” ( S305 ).

[0073] The weather information acquisition unit 113 accesses the weather information database 33 (or from public weather data, etc.) (S306) to obtain the daily average daytime temperature of "Sapporo" (the average outside temperature calculated on a daily basis using the daily maximum temperature and daily minimum temperature obtained from the weather data of the past year), and the daily day length (s) = 32400s (a fixed value is used in this embodiment, but a value taking into account the latitude and season can also be used) (S307).

[0074] After acquiring these information, the absorption amount calculation unit 11 further executes the following steps.

[0075] (Calculation of daytime heat dissipation)

[0076] First, the daily daytime indoor temperature is calculated using the daily average daytime temperature and the daytime ventilation temperature (23°C) (S308). The daily daytime indoor temperature can be calculated using a machine learning model that uses an approximate equation derived through regression analysis using the daily daytime indoor temperature as the target variable and the daily average daytime temperature and the daytime ventilation temperature (23°C) as explanatory variables. In this embodiment, based on the derived approximate equation, if the daily daytime average temperature is less than 5°C, the daily daytime indoor temperature is calculated as (daytime ventilation temperature - 2). If the daily daytime average temperature is 5°C or higher, the daily daytime indoor temperature is calculated as (-0.8224 × daily daytime average temperature + 1.299 × daytime ventilation temperature + 0.03229 × daily daytime average temperature squared - 18.94 / daily daytime average temperature). In this example, the average daytime temperature is -0.7291°C, and thus "21°C" is calculated using the former formula (S308).

[0077] Next, the heat preservation ratio is calculated ( S309 ). The heat preservation ratio is calculated as the design floor area / design covering area, and in the above example, it is calculated as 1000 / 1890.58=0.5289.

[0078] Next, calculate the daily heat dissipation (MJ / m 2 )(S310).

[0079] Daily daytime heat dissipation (MJ / m 2 ) = (Daily daytime indoor temperature (°C) - Daily average daytime temperature (°C)) / Insulation ratio × Average heat dissipation coefficient of roof material × Daily day length (s) / 3600 × 4.186 / 1000 = 7.7384 (MJ / m 2 )

[0080] (Calculation of daily solar heat deduction)

[0081] Next, the total daily sunshine amount (MJ / m 2 ) (calculated based on the meteorological data for one year) is multiplied by the light transmittance of the roofing material obtained from the structural information database 31 to obtain the daily incident solar heat (MJ / m 2 )( Figure 4 S401). Using the daily incident solar heat (MJ / m 2 ) and by Figure 3 The daily heat dissipation (MJ / m 2 ), calculate the daily solar heat deduction (MJ / m2 )(S402).

[0082] Daily daytime heat dissipation (MJ / m 2 )>Daily incident solar heat(MJ / m 2 ), it is calculated as "0".

[0083] Daily daytime heat dissipation (MJ / m 2 )≤ Daily incident solar heat (MJ / m 2 ) is calculated by “daily daytime heat dissipation (MJ / m 2 ) - Daily incident solar heat (MJ / m 2 )" to find out.

[0084] (Calculation of outdoor specific enthalpy)

[0085] Next, the daily maximum and minimum temperatures of the region acquired from the weather information database 33 or public weather data are used to determine the daily outdoor specific enthalpy (J / kg) ( S403 ). The outdoor specific enthalpy is determined using the following learned model.

[0086] First, using the average daytime temperature and absolute humidity obtained from meteorological authorities at multiple observation locations, the daytime specific enthalpy was calculated using the following formula (where Td is the temperature and χ is the absolute humidity).

[0087] H=1.006·Td+(1.86)·Td+2501)·χ

[0088] The obtained specific enthalpy was used as the target variable, and the maximum and minimum temperatures were obtained daily and used as explanatory variables. This process was repeated for each of the aforementioned observation locations, obtaining approximately two years of data from official meteorological authorities. A learned model was generated using linear regression. As an example, the learned model used in this embodiment is represented as follows.

[0089] Outdoor specific enthalpy (J / kg): H = 4211 + 1389·T h +382.4·T l -16.26·T h 2 +55.90T l 2 (T h : Daily maximum temperature, T l : Daily minimum temperature) (Coefficient of determination = 0.9683)

[0090] Thus, the outdoor specific enthalpy is obtained using the daily maximum temperature and daily minimum temperature of the region.

[0091] (Calculation of indoor specific enthalpy)

[0092] Used by Figure 3 The indoor specific enthalpy (J / kg) is calculated using the daily daytime indoor air temperature obtained in S308 (S404). A learned model based on linear regression is generated using data from actual agricultural greenhouses to determine the indoor air temperature for determining the indoor specific enthalpy. Specifically, in an agricultural greenhouse operated by Tomato Park Co., Ltd., an associated company of the present applicant, the daytime indoor air temperature (target variable), the daytime average external air temperature (explanatory variable), and the ventilation temperature (explanatory variable) are measured for several dozen weeks, and a learned model is generated using linear regression. Furthermore, the relative humidity used to determine χ is set to a fixed value, and the indoor air temperature (Td) is divided into a constant range and subjected to regression analysis to generate a learned model with the indoor specific enthalpy as the target variable. As an example, the learned model used in this embodiment is represented as follows.

[0093] Indoor specific enthalpy (J / kg): H = 63.242Td 2 +181.18Td+16504 (In the range of indoor temperature Td = 10 to 40°C, determination coefficient = 0.9996)

[0094] (Calculation of ventilation rate)

[0095] Next, calculate the daily ventilation rate (m 3 / m 2 ) (S405). In addition to the daily solar heat deduction (MJ / m 2 ), outdoor specific enthalpy (J / kg) and indoor specific enthalpy (J / kg), and the hourly interstitial ventilation rate (m2) specified according to the type of roofing material obtained from the structural information database 31 by the greenhouse information acquisition unit 111 is also used. 3 / hm 2 ) (Example: 1.8m 3 / hm 2 ) and the above-mentioned daily length of day (s) (32400s) associated with "Sapporo" (S406, S407) to calculate the daily ventilation rate.

[0096] Specifically, first, compare the following two values:

[0097] Divide by the difference between the outdoor specific enthalpy (J / kg) and the indoor specific enthalpy (J / kg) (EA)

[0098] =100000×daily solar heat deduction (MJ / m 2 ) / (outdoor specific enthalpy (J / kg) - indoor specific enthalpy (J / kg)) / 1.2(air density kg / m3 ),

[0099] The gap air exchange rate (m 3 / hm 2 ) and the value obtained by multiplying the daily length (s) (EB)

[0100] = Daily air exchange rate (m 3 / m 2 ) = interstitial air exchange rate (m 3 / hm 2 ) × daily day length (s) / 3600.

[0101] When the comparison result is EA<EB, the daily air exchange rate (m 3 / m 2 ) = interstitial air exchange rate (m 3 / hm 2 ) × daily day length (s) / 3600.

[0102] In the case of EA≥EB, 100000× daily solar heat deduction (MJ / m 2 ) / (outdoor specific enthalpy (J / kg) - indoor specific enthalpy (J / kg)) / 1.2(air density kg / m 3 ).

[0103] (Calculation of indoor carbon dioxide concentration)

[0104] The absorption amount calculation unit 11 further calculates the carbon dioxide supply amount per unit area based on the input value of the carbon dioxide supply amount and the designed floor area (S408), and multiplies this value by the daily length (s) = 32400s to calculate the carbon dioxide supply amount per unit area per day (S409). In addition, when the carbon dioxide supply amount is in kg, it is divided by the carbon dioxide 1m3 at 1 atmosphere and 25°C. 3 The mass, converted to m 3 The daily carbon dioxide supply per unit area is calculated using the unit ( S410 ). In addition, the weather information acquisition unit 113 acquires the outdoor carbon dioxide concentration (410 ppm) of “Sapporo” from the weather information database 33 ( S411 , S412 ).

[0105] The plant information acquisition unit 112 of the absorption calculation unit 11 accesses the plant information database 32 and obtains the factor (δ), light reception rate (φ), and conversion coefficient involved in the reverse reaction for "tomato" (S413-S415). The conversion coefficient is used to convert the light reception (MJ) into the photosynthesis amount (m 3 ).

[0106] Upon obtaining the above information, the absorption amount calculation unit 11 calculates the daily indoor carbon dioxide concentration using the following formula ( S417 ).

[0107] Indoor carbon dioxide concentration = (-(light reception rate × daily incident solar heat × conversion coefficient + daily ventilation rate × δ - daily carbon dioxide supply (per unit area) (m 3 / m 2 )-daily ventilation rate × indoor carbon dioxide concentration) + ((light reception rate × daily incident solar heat × conversion coefficient + daily ventilation rate × δ-daily carbon dioxide supply (per unit area) (m 3 / m 2 ) - daily air change rate × outdoor carbon dioxide concentration) 2 -4× daily air exchange rate×(-δ× daily carbon dioxide supply (per unit area) (m 3 / m 2 )-δ×daily air change rate×square root of outdoor carbon dioxide concentration)) / 2 / daily air change rate (m 3 / m 2 )

[0108] Figure 4 The figure shows an example of calculation results obtained by specifically applying numerical values ​​to obtain the indoor carbon dioxide concentration by the absorption amount calculation unit 11. In this way, for example, the indoor carbon dioxide concentration is obtained as 0.004432509.

[0109] (Calculation of daily net photosynthesis)

[0110] The plant information acquisition unit 112 accesses the plant information database 32 and acquires the factor (δ), light reception rate (φ), and conversion coefficient ( Figure 5 Use the conversion coefficient to convert the amount of light received (MJ) into the amount of photosynthesis (m 3 )(S504).

[0111] The greenhouse information acquisition unit 111 accesses the structural information database 31 to acquire the designed floor area ( S505 ).

[0112] The daily incident solar heat (MJ / m 2 ) (S401), daily indoor carbon dioxide concentration (S417), and daily net photosynthesis amount (S506).

[0113] In this embodiment, it is obtained by the following formula.

[0114] Daily net photosynthesis (m 3 ) = Daily incident solar heat (MJ / m 2)×light receiving rate×conversion coefficient (m 3 / J) × daily indoor carbon dioxide concentration × designed floor area / (δ + daily indoor carbon dioxide concentration)

[0115] (Calculation of carbon dioxide absorption)

[0116] The net photosynthetic rate per day (m 3 ) is converted into the daily carbon dioxide absorption amount, which is multiplied by the conversion coefficient to convert to kg (S507), thereby calculating the daily carbon dioxide absorption amount in kg (S508).

[0117] By summing the daily CO2 absorption during the cultivation period (for example, cultivation starts on September 1st and ends on June 1st of the following year), the CO2 absorption for the entire cultivation period is calculated. As a simulation, the CO2 absorption is calculated on a monthly basis (S509), and the annual CO2 absorption is calculated by accumulating the data for one year (S510). By calculating the CO2 absorption for a specific period, such as the actual cultivation period, month, or year, the utilization rate and reduction of CO2 can be compared in the utilization rate calculation unit 12, the first greenhouse gas emission calculation unit 13, the second greenhouse gas emission calculation unit 14, the greenhouse gas reduction calculation unit 15, and the cost calculation unit 16, which will be described later.

[0118] Figure 5 An example of calculation results of the carbon dioxide absorption amount calculated by the absorption amount calculation unit 11 is shown. In this example, the net photosynthetic amount per day is calculated as 9.877358432 m 3 = 17.76936782 kg, daily carbon dioxide absorption = 17.76936782 kg, carbon dioxide absorption during one month (the cumulative example from January 1 to January 31 in the above cultivation period) = 0.65249752 t, and annual carbon dioxide absorption = 10.10189154 t. The absorption calculation unit 11 outputs "10.1 t" to Figure 2 The “CO 2 absorption amount” column of the output screen is displayed ( S511 ).

[0119] [Calculation of utilization rate]

[0120] The utilization rate calculation unit 12 calculates the utilization rate using the carbon dioxide absorption amount obtained by the absorption amount calculation unit 11 .

[0121] The utilization rate is obtained by comparing it with the carbon dioxide supply amount. The carbon dioxide supply amount is obtained as follows.

[0122] First, suppose that 15 kg / h is input as the carbon dioxide supply rate (refer to Figure 2 Divide the carbon dioxide supply by the above-mentioned designed floor area of ​​1000m 2 Calculate the hourly supply per unit area to be 0.015 kg / hm 2 ( Figure 6 Next, the day length of "Sapporo" (32400s) is read from the weather information database 33 (S602), and the hourly supply per unit area per day is calculated (S603), which is further multiplied by the design ground area of ​​1000m 2 , the daily carbon dioxide supply is calculated as 135 kg (S604).

[0123] Once the daily CO2 supply is calculated, the daily CO2 supply is accumulated over an appropriate period, and the period CO2 supply is calculated using the actual cultivation period, month, year, etc., similarly to the CO2 absorption described above (S605-S607). Thus, for example, the period CO2 supply for January is calculated as 4.185 t, and the period CO2 supply for the year is calculated as 40.77 t. The utilization rate calculation unit 12 outputs "40.8 t" to the Figure 2 The “CO 2 supply amount” column of the output screen is displayed ( S607 ).

[0124] For example, the utilization rate calculation unit 12 calculates the utilization rate as 10.10 / 40.77=24.8% when comparing between years, and calculates the utilization rate as 0.652 / 4.185=15.6% when comparing between January.

[0125] exist Figure 2 In the example of the output screen shown, "CO2 absorption: 10.1t" is displayed as the annual carbon dioxide absorption amount, and "CO2 supply: 40.8t" is displayed as the annual carbon dioxide supply amount, but the utilization rate calculation unit 12 further outputs the annual utilization rate as "25%" in the "CO2 utilization rate" column and displays it.

[0126] The carbon dioxide dynamics simulation device of this embodiment can determine the absorption and utilization rate of carbon dioxide supplied to an agricultural greenhouse. This allows for a clear understanding of the dynamics of carbon dioxide within the agricultural greenhouse. Furthermore, since the absorption and utilization rates can be determined, it is possible to suppress excessive carbon dioxide supply that is inconsistent with the utilization rate. This helps optimize the operating frequency and output of a device that supplies carbon dioxide through the combustion of fossil fuels. This, in turn, helps suppress greenhouse gas emissions caused by operating the device that supplies carbon dioxide more than necessary.

[0127] Next, a carbon dioxide dynamic simulation device 1 involved in another embodiment of the present invention is described. In addition to calculating the absorption amount and utilization rate in order to understand the dynamics of carbon dioxide, the carbon dioxide dynamic simulation device 1 involved in the other embodiment also compares the case where a device that produces carbon dioxide by using fossil fuels is used as a carbon dioxide supply source and the case where exhaust gas from factories, etc. is used, and can simulate the contribution to the reduction of greenhouse gases when exhaust gas is used.

[0128] In this embodiment, if Figure 7 As shown, similarly to the above-mentioned embodiment, there is an absorption amount calculation unit 11 and a utilization rate calculation unit 12, and a computer program for executing steps of functioning as a first greenhouse gas emission amount calculation unit 13, a second greenhouse gas emission amount calculation unit 14, a greenhouse gas reduction amount calculation unit 15 and a cost calculation unit 16 is stored in the storage unit 1b.

[0129] The first greenhouse gas emission calculation unit 13 has the function of calculating the first greenhouse gas emission amount when the carbon dioxide supplied to the agricultural greenhouse is supplied using a device that generates carbon dioxide through the combustion of fossil fuels. The second greenhouse gas emission calculation unit 14 has the function of calculating the second greenhouse gas emission amount when the supply amount is supplied using carbon dioxide separated and recovered from exhaust gas emitted from the exhaust gas generation source. The greenhouse gas reduction calculation unit 15 calculates the difference between the first greenhouse gas emission amount and the second greenhouse gas emission amount as the greenhouse gas reduction amount when exhaust gas supplied from the exhaust gas generation source is used.

[0130] The first greenhouse gas emission amount is equivalent to the amount obtained by subtracting the amount absorbed by plants from the amount supplied to the agricultural greenhouse. Regarding the second greenhouse gas emission amount, since exhaust gas is used, the amount emitted to the outside during ventilation does not constitute new carbon dioxide emissions. As a result, greenhouse gas emissions are reduced by the amount absorbed by plants. When the total amount supplied to the agricultural greenhouse is supplied by carbon dioxide separated / recovered from the exhaust gas, the second greenhouse gas emission amount calculation unit 14 calculates the amount equivalent to the absorption amount calculated by the absorption amount calculation unit 11 as a negative number as the second greenhouse gas emission amount. Therefore, the greenhouse gas reduction amount calculation unit 15 outputs the value obtained by adding the absolute value of the second greenhouse gas emission amount to the value of the first greenhouse gas emission amount as the reduction amount. As a result, when the total amount of carbon dioxide supplied is supplied by exhaust gas, the carbon dioxide supply amount becomes the difference between the two, i.e., the reduction amount, compared to an apparatus using fossil fuels.

[0131] The first greenhouse gas emission amount calculation unit 13, the second greenhouse gas emission amount calculation unit 14 and the greenhouse gas reduction amount calculation unit 15 calculate the first greenhouse gas emission amount, the second greenhouse gas emission amount and the greenhouse gas reduction amount, and output the values ​​to the respective Figure 2 The output screen shown is displayed.

[0132] The cost calculation unit 16 calculates the cost required for supplying carbon dioxide based on the carbon dioxide supply amount, compares the carbon dioxide supply cost when using fossil fuels with the carbon dioxide supply cost when using exhaust gas, and calculates the difference as the cost of reducing the carbon dioxide supply cost.

[0133] The carbon dioxide supply cost when using a device that burns fossil fuels is calculated using the unit price of the fuel, the calorific value of the fuel, and the carbon emission coefficient.

[0134] The unit price of carbon dioxide (CO2) when utilizing waste gas can be determined through agreements between the companies generating the waste gas, the builders of agricultural greenhouses, and the government. Alternatively, the market price of CO2 emissions trading can be referenced. Overall, the CO2 cost should be significantly lower than when using fossil fuels.

[0135] Next, an example of simulating the reduction status of greenhouse gases using the carbon dioxide dynamics simulation device 1 of this embodiment will be described.

[0136] [Calculation of greenhouse gas reduction amount]

[0137] The difference between the greenhouse gas emissions when a device generating carbon dioxide by burning fossil fuels is used as a carbon dioxide supply source and the greenhouse gas emissions when carbon dioxide from exhaust gas is used is calculated as the greenhouse gas reduction amount.

[0138] (When using a device that generates carbon dioxide by using fossil fuels)

[0139] The first greenhouse gas emission amount calculation unit 13 calculates the greenhouse gas emission amount using the following equation.

[0140] Greenhouse gas emissions = Carbon dioxide supply - Carbon dioxide absorption

[0141] Using the above example to compare on an annual basis, the first greenhouse gas emissions are calculated as:

[0142] 40.77t-10.10t=30.67t.

[0143] (When the total amount of carbon dioxide supplied is calculated as exhaust gas)

[0144] The second greenhouse gas emission amount calculation unit 14 calculates the greenhouse gas emission amount using the following equation.

[0145] Second greenhouse gas emissions = - Carbon dioxide absorption

[0146] Using the above example to compare between years, we can find:

[0147] Second greenhouse gas emissions = -10.10t.

[0148] (Greenhouse gas reduction)

[0149] The greenhouse gas reduction amount calculation unit 15 calculates the amount using the following equation.

[0150] Greenhouse gas reduction = primary greenhouse gas emissions - secondary greenhouse gas emissions

[0151] In the above example, the greenhouse gas reduction amount is calculated as:

[0152] 30.67t-(-10.10t)=40.77t.

[0153] In this example, the total amount of carbon dioxide supplied using exhaust gas is assumed, so the difference between the first greenhouse gas emissions and the second greenhouse gas emissions, i.e., the greenhouse gas reduction, is equal to the carbon dioxide supply. However, if, for example, exhaust gas from a factory or other facility is insufficient during a period, carbon dioxide generated by burning fossil fuels may be used to compensate for the deficiency. In this case, assuming that fossil fuels are used for one month out of a year during the cultivation period and exhaust gas is used for the remaining periods, the greenhouse gas reduction for the period is calculated using "carbon dioxide supply minus carbon dioxide absorption" similarly to the first greenhouse gas emissions calculation unit 13 for the one-month period. For the remaining periods, the negative of the carbon dioxide absorption is calculated as the greenhouse gas emissions for the period, and the sum of these two values ​​is calculated as the second greenhouse gas emissions.

[0154] For example, assuming the annual CO2 supply is 40.77 tons and the annual CO2 absorption is 10.10 tons, and assuming the monthly CO2 supply is 3.5 tons and the monthly CO2 absorption is 2.1 tons, the greenhouse gas emissions from the fossil fuel supply for one month are 3.5 - 2.1 = 1.4 tons. The greenhouse gas emissions (CO2 absorption) from the waste gas supply for the remaining 11 months are 10.10 - 2.1 = -8 tons. Therefore, the sum of these two amounts, -6.6 tons, represents the second greenhouse gas emissions, and the greenhouse gas reduction is 30.67 - (-6.6) = 37.27 tons.

[0155] For this reason, Figure 2 As shown in FIG. 1 , it is preferable to input the ratio of exhaust gas utilization in the carbon dioxide supply amount as the circulation rate as the input information.

[0156] In this case, in the absorption amount calculation unit 11, the daily carbon dioxide absorption amount is accumulated excluding the period corresponding to the circulation rate, and in the second greenhouse gas emission amount calculation unit 14, as described above, the second greenhouse gas emission amount according to the circulation rate can be calculated by adding up the period of supply using fossil fuel and the period of supply using exhaust gas.

[0157] exist Figure 2 In the output screen shown, in the column displayed as "CO2 absorption amount", the first greenhouse gas emission amount calculation unit 13, the second greenhouse gas emission amount calculation unit 14, and the greenhouse gas reduction amount calculation unit 15 respectively output "30.7t" as the annual greenhouse gas emission amount when a device that generates carbon dioxide through the use of fossil fuels is used (without resource circulation), output "-10.1t", which is a negative number equivalent to the value of the annual carbon dioxide absorption amount, as the annual greenhouse gas emission amount when exhaust gas is used at a circulation rate of 100% (with resource circulation), and output "40.8t" as the greenhouse gas reduction amount, and display them.

[0158] [Calculation of Carbon Dioxide Supply Cost Reduction Expenses]

[0159] The cost calculation unit 16 compares the carbon dioxide supply cost when the carbon dioxide supply amount is supplied by a device that generates carbon dioxide by burning fossil fuels and the carbon dioxide supply cost when the carbon dioxide supply amount is supplied using exhaust gas.

[0160] (When using a device that generates carbon dioxide by using fossil fuels)

[0161] For example, in the case of a device using LPG as fuel, if the unit price of LPG is set at ¥256 / kg, then based on a calorific value of 50.08MJ / kg and a carbon emission coefficient of 0.01637kgC / MJ, the unit price of carbon dioxide by mass is:

[0162] 256 / (0.01637×50.08×44 / 12)=¥85.

[0163] According to the above example, assuming that the carbon dioxide supply amount during January is 4.185 tons and the carbon dioxide supply amount during the year is 40.77 tons, the carbon dioxide supply amount is calculated as follows.

[0164] Monthly CO2 supply cost for January = ¥85 x 4185 kg = ¥355,725

[0165] Annual CO2 supply cost = ¥85 x 40,770 kg = ¥3,465,450

[0166] (When the total amount of carbon dioxide supplied is calculated as exhaust gas)

[0167] Based on the transaction price of exhaust gas, the unit price of carbon dioxide is assumed to be ¥38 and the calculation is performed in the same manner as above.

[0168] Monthly CO2 supply cost for January = ¥38 x 4185kg = ¥159,030

[0169] Annual CO2 supply cost = ¥38 x 40,770 kg = ¥1,549,260

[0170] (CO2 supply cost reduction costs)

[0171] Monthly CO2 reduction costs for January = ¥355,725 - ¥159,030 = ¥196,695

[0172] Annual CO2 reduction costs = ¥3,465,450 - ¥1,549,260 = ¥1,916,190

[0173] Furthermore, when the circulation rate is set, the carbon dioxide supply cost and the carbon dioxide supply cost reduction fee are calculated based on the circulation rate, similarly to the case of determining the greenhouse gas emission amount and greenhouse gas reduction amount described above.

[0174] exist Figure 2In the output screen shown, in the column displayed as "CO2 cost", the cost calculation unit 16 outputs "¥3,465,450" as the annual carbon dioxide supply cost when a device that generates carbon dioxide through the use of fossil fuels is used (without resource circulation), outputs "¥1,549,260" as the annual carbon dioxide supply cost when exhaust gas is used at a circulation rate of 100% (with resource circulation), and outputs "¥1,916,190" as the reduction cost, and displays them.

[0175] Based on the above, the carbon dioxide dynamic simulation device 1 of this embodiment can quantitatively and clearly show to what extent the amount of carbon dioxide supply can be reduced compared with the case of using a device that uses fossil fuels to supply carbon dioxide when the carbon dioxide supply is supplied by exhaust gas. In addition, it can quantitatively and clearly show to what extent the supply cost of carbon dioxide can be reduced.

[0176] Therefore, in this embodiment, the advantages of building agricultural greenhouses that recycle waste gas can be clarified, which serves as an incentive to promote the construction of agricultural greenhouses that can contribute to the reduction of greenhouse gases. In addition, greenhouse gas reduction can also be achieved in factories that emit waste gas, which can be conducive to the activation of carbon dioxide emission trading. That is, from the perspective of factories and the like that are the waste gas supply side, "¥1,549,260" as an example of the annual carbon dioxide supply cost in the case of "resource recycling" mentioned above becomes a new sales amount based on the waste gas that was previously only released into the atmosphere. Therefore, in Figure 2 In the output screen, it is preferable to set an item that displays this amount as the total sales of the exhaust gas supplier. Being able to grasp such a value before the construction of the agricultural greenhouse will serve as an incentive for the exhaust gas supplier to achieve cooperation with the agricultural greenhouse.

[0177] -Description of Reference Numerals-

[0178] 1 Carbon dioxide dynamic simulation device

[0179] 11 Absorption amount calculation unit

[0180] 111 Greenhouse Information Acquisition Department

[0181] 112 Plant Information Acquisition Department

[0182] 113 Meteorological Information Acquisition Department

[0183] 12 Utilization Rate Calculation

[0184] 13. 1st Greenhouse Gas Emissions Calculation Unit

[0185] 14. Second greenhouse gas emission calculation unit

[0186] 15 Greenhouse Gas Reduction Calculation Department

[0187] 16 Cost Calculation Department.

Claims

1. A carbon dioxide dynamic simulation device to simulate the dynamics of carbon dioxide in agricultural greenhouses. The carbon dioxide dynamic simulation device is characterized by having: an absorption amount calculation unit that obtains a net photosynthesis amount of plants cultivated in the agricultural greenhouse as an absorption amount of carbon dioxide by the plants; and The utilization rate calculation unit calculates a ratio of the absorption amount to the supply amount of carbon dioxide supplied from a carbon dioxide supply source into the agricultural greenhouse during a predetermined period as a utilization rate of the carbon dioxide supplied to the agricultural greenhouse.

2. The carbon dioxide dynamic simulation device according to claim 1, wherein: The absorption amount calculation unit calculates the net photosynthesis amount using the total daily sunlight amount, the light transmittance of the agricultural greenhouse, the light receiving rate of the plant, and the indoor carbon dioxide concentration of the agricultural greenhouse.

3. The carbon dioxide dynamic simulation device according to claim 2, wherein: The indoor carbon dioxide concentration is calculated based on the indoor and outdoor ventilation rate of the air in the agricultural greenhouse. The ventilation rate is calculated using the difference between the outdoor specific enthalpy and the indoor specific enthalpy of the agricultural greenhouse.

4. The carbon dioxide dynamic simulation device according to claim 2, wherein: The absorption amount calculation unit includes: a greenhouse structure information acquisition unit, which acquires the structural information of the agricultural greenhouse including the light transmittance; a plant information acquisition unit configured to acquire information of the plant to be cultivated including the light reception rate; and The weather information acquisition unit acquires weather information of an area where the agricultural greenhouse is installed.

5. The carbon dioxide dynamic simulation device according to claim 1, wherein: The carbon dioxide dynamic simulation device has: a first greenhouse gas emission amount calculation unit for calculating a first greenhouse gas emission amount when a device that generates carbon dioxide by burning fossil fuels is used as the carbon dioxide supply source to supply the carbon dioxide to the agricultural greenhouse; a second greenhouse gas emission amount calculation unit for calculating a second greenhouse gas emission amount when the amount of carbon dioxide supplied to the agricultural greenhouse is supplied by carbon dioxide recovered from exhaust gas emitted from an exhaust gas generation source serving as the carbon dioxide supply source; and The greenhouse gas reduction amount calculation unit calculates a difference between the first greenhouse gas emission amount and the second greenhouse gas emission amount as a greenhouse gas reduction amount when the exhaust gas is used.

6. The carbon dioxide dynamic simulation device according to claim 5, wherein: The first greenhouse gas emission amount calculation unit calculates a difference between the carbon dioxide supply amount and the carbon dioxide absorption amount as the first greenhouse gas emission amount.

7. The carbon dioxide dynamic simulation device according to claim 5, wherein: The second greenhouse gas emission amount calculation unit calculates the second greenhouse gas emission amount by setting a value corresponding to the carbon dioxide absorption amount to a negative number.

8. The carbon dioxide dynamic simulation device according to claim 5, wherein: The carbon dioxide dynamics simulation device includes a cost calculation unit that compares the carbon dioxide supply amount when the fossil fuel is used and the carbon dioxide supply amount when the exhaust gas is used to calculate a reduction cost for carbon dioxide supply.

9. A computer program for causing a computer to function as a device for simulating the dynamics of carbon dioxide in an agricultural greenhouse. The computer program is characterized in that it causes the computer to execute the following steps: determining the net photosynthesis of the plants cultivated in the agricultural greenhouse as the carbon dioxide absorption amount of the plants; and The ratio of the absorption amount to the supply amount of carbon dioxide supplied from the carbon dioxide supply source to the agricultural greenhouse during a predetermined period is determined as the utilization rate of the carbon dioxide supplied to the agricultural greenhouse.

10. The computer program according to claim 9, wherein In the step of determining the carbon dioxide absorption amount, the net photosynthesis amount is determined using the total daily sunlight amount, the light transmittance of the agricultural greenhouse, the light reception rate of the plant, and the indoor carbon dioxide concentration of the agricultural greenhouse.

11. The computer program according to claim 10, wherein The indoor carbon dioxide concentration is calculated based on the indoor and outdoor ventilation rate of the air in the agricultural greenhouse. The ventilation rate is calculated using the difference between the outdoor specific enthalpy and the indoor specific enthalpy of the agricultural greenhouse.

12. The computer program according to claim 10, wherein In the step of determining the amount of carbon dioxide absorbed, the computer program causes the computer to execute the following steps: acquiring structural information of the agricultural greenhouse including the light transmittance; acquiring information of the plant to be cultivated including the light reception rate; and Acquire weather information of the area where the agricultural greenhouse is installed.

13. The computer program according to claim 9, wherein The computer program further causes the computer to perform the following steps: determining a first greenhouse gas emission amount when a device that generates carbon dioxide by burning fossil fuels is used as the carbon dioxide supply source to supply the carbon dioxide to the agricultural greenhouse; determining a second greenhouse gas emission when the amount of carbon dioxide supplied to the agricultural greenhouse is supplied by carbon dioxide recovered from exhaust gas emitted from an exhaust gas generation source serving as the carbon dioxide supply source; and A difference between the first greenhouse gas emission amount and the second greenhouse gas emission amount is calculated as a greenhouse gas reduction amount when the exhaust gas is utilized.

14. The computer program according to claim 13, wherein The computer program further causes the computer to perform the following steps: The amount of carbon dioxide supplied when the fossil fuel is used is compared with the amount of carbon dioxide supplied when the exhaust gas is used to calculate a reduction in carbon dioxide supply cost.

15. A computer-readable recording medium, characterized in that The computer program according to any one of claims 9 to 14 is recorded.

Images