Three-phase distribution calculation method, three-phase distribution calculation device, and three-phase distribution calculation system for soil, and program

Abstract

To provide a technology for obtaining three-phase distribution of soil by using soil data without actually measuring the three-phase distribution.SOLUTION: Method (S1) for calculating three-phase distribution of soil includes dielectric constant acquisition step (S12) for acquiring a dielectric constant ε of the soil, heat transfer rate acquisition step (S13) for acquiring a heat transfer rate λ of the soil, and volume ratio calculation step (S14) for calculating volume ratios G, L, and S of the soil. The volume ratio calculation step uses a first function formula (1) in which the respective volume ratios of the gaseous phase, liquid phase and solid phase of the soil are G, L and S, and the dielectric constant ε of the soil is represented by dielectric constants εg, εl and εs of the gaseous phase, liquid phase and solid phase of the soil, respectively, a second function formula (2) in which the thermal conductivity λ of the soil is represented by thermal conductivities λg, λl and λs of the gaseous phase, liquid phase and solid phase of the soil, respectively, and a third function formula (3) in which the total of the respective volume ratios of the gaseous phase, liquid phase and solid phase of the soil is 1.SELECTED DRAWING: Figure 2

Description

[Technical Field]

[0001] The present invention relates to a method for calculating three-phase distribution in soil, a three-phase distribution calculation device, a three-phase distribution calculation system, and a program. [Background technology]

[0002] One type of soil data is the three-phase distribution, which is the volume ratio of the gas, liquid, and solid phases in the soil. The three-phase distribution is considered to be a factor that affects crop growth. Therefore, when growing crops in a field, information on the three-phase distribution of the soil is considered important, and methods have been developed to sample soil and actually measure the three-phase distribution. One such measurement method is known as the actual volume method. For example, Patent Document 1 discloses a measuring device that can mechanically measure the three-phase distribution of sampled soil.

[0003] Meanwhile, there are documents that examine the relationship between any of the volume fractions of the gas, liquid, or solid phases of soil and the thermal conductivity or dielectric constant of soil. For example, Non-Patent Document 1 examines the relationship between the thermal conductivity of soil and its porosity and water content. Furthermore, Non-Patent Document 2 describes that the dielectric constant of soil is related to its water content, soil texture, and density. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] Japanese Patent Application Laid-Open No. 2010-217156 [Non-patent literature]

[0005] [Non-Patent Document 1] Kobayashi, 1965, 3. Thermal conductivity of soil, Hokkaido University Geophysical Research Report 13:37-51. [Non-patent document 2] Shimobe, 2004, Relationship between soil dielectric constant and water content, 39th Geotechnical Engineering Research Conference, 307-308. Summary of the Invention [Problem to be solved by the invention]

[0006] However, the method described in Patent Document 1 involves sampling soil and measuring the three-phase distribution, which requires the effort of sampling soil. The above-mentioned Non-Patent Documents 1 and 2 state that the dielectric constant or thermal conductivity of soil is related to any one of the volume fractions of the gas, liquid, or solid phases, but do not disclose a method for determining the three-phase distribution using any soil data.

[0007] An object of one aspect of the present invention is to provide a technique for determining the three-phase distribution of soil using soil data without actually measuring it. [Means for solving the problem]

[0008] In order to solve the above-mentioned problems, a method for calculating a three-phase distribution of soil according to one aspect of the present invention includes a permittivity acquisition step of acquiring a permittivity ε of soil, a thermal conductivity acquisition step of acquiring a thermal conductivity λ of the soil, and a step of calculating a permittivity ε of the soil by dividing the permittivity ε of the soil by the volume fractions of the gas phase, liquid phase, and solid phase of the soil, where G, L, and S are the volume fractions of the gas phase, liquid phase, and solid phase of the soil. g , ε l , ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,ε l ,ε s ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by g , λ l , λ s The second functional formula (2) below is expressed as: λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third functional formula (3) is used, in which the total value of the volume ratios of the gas phase, liquid phase, and solid phase of the soil is 1: 1.0=G+L+S …(3) and a volume ratio calculation step of calculating volume ratios G, L, and S of the soil using the above formula.

[0009] A three-phase distribution calculation device according to one aspect of the present invention includes a data acquisition unit that acquires permittivity data and thermal conductivity data of soil, and a data acquisition unit that acquires permittivity data and thermal conductivity data of soil, and a data acquisition unit that acquires permittivity data and thermal conductivity data of soil, where G, L, and S are volume fractions of a gas phase, a liquid phase, and a solid phase of the soil, respectively. g , ε l , ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,ε l ,ε s ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by g , λ l , λ s The second functional formula (2) below is expressed as: λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third functional formula (3) is used, in which the total value of the volume ratios of the gas phase, liquid phase, and solid phase of the soil is 1: 1.0=G+L+S …(3) and a calculation unit that calculates the volume ratios G, L, and S of the soil using the above.

[0010] A three-phase distribution calculation system according to one aspect of the present invention includes the above-described three-phase distribution calculation device, a sensor that detects the dielectric constant and thermal conductivity of soil, and a data logger that can collect data measured by the sensor and transmit the data to the three-phase distribution calculation device.

[0011] The three-phase distribution calculation device or three-phase distribution calculation system according to each aspect of the present invention may be realized by a computer. In this case, the control program that causes the computer to operate as each part (software element) of the three-phase distribution calculation device or three-phase distribution calculation system to realize the three-phase distribution calculation device or three-phase distribution calculation system on the computer, and the computer-readable recording medium on which the control program is recorded, also fall within the scope of the present invention. [Effects of the Invention]

[0012] According to one aspect of the present invention, the three-phase distribution of soil can be determined using soil data without actually measuring it. [Brief explanation of the drawings]

[0013] [Figure 1] 1 is a block diagram showing the configuration of a three-phase distribution calculation system 100 and a three-phase distribution calculation device 1 according to an embodiment. [Figure 2] 10 is a flowchart showing the flow of a three-phase distribution calculation method S1 according to the embodiment. [Figure 3] 10 is a graph showing the correlation between the calculated value (estimated value) of the volume fraction of the solid phase in gray lowland soil obtained according to the three-phase distribution calculation method S1 of the embodiment and the actual measured value obtained using the actual volume method. [Figure 4] 10 is a graph showing the correlation between the calculated (estimated) volume fraction of the liquid phase in gray lowland soil obtained according to the three-phase distribution calculation method S1 of the embodiment and the measured value obtained using the actual volume method. [Figure 5] 10 is a graph showing the correlation between the calculated value (estimated value) of the volume fraction of the gas phase in gray lowland soil obtained according to the three-phase distribution calculation method S1 of the embodiment and the actual measured value obtained using the actual volume method. [Figure 6] 10 is a graph showing the correlation between the calculated value (estimated value) of the volume fraction of the solid phase in a field obtained according to the three-phase distribution calculation method S1 of the embodiment and the actual measured value obtained using the actual volume method. [Figure 7]10 is a graph showing the correlation between the calculated value (estimated value) of the volume fraction of the liquid phase in a field obtained according to the three-phase distribution calculation method S1 of the embodiment and the actual measured value obtained using the actual volume method. [Figure 8] 10 is a graph showing the correlation between the calculated value (estimated value) of the volume fraction of the gas phase in a field obtained according to the three-phase distribution calculation method S1 of the embodiment and the actual measured value obtained using the actual volume method. DETAILED DESCRIPTION OF THE INVENTION

[0014] [Embodiment 1] An embodiment of the present invention will be described in detail below with reference to the drawings. Fig. 1 is a block diagram showing the configuration of a three-phase distribution calculation system 100 and a three-phase distribution calculation device 1 according to this embodiment. The three-phase distribution calculation device 1 is a device that calculates a three-phase distribution from the measured dielectric constant and thermal conductivity of soil. The three-phase distribution calculation system 100 includes a sensor Bn and a data logger An in addition to the three-phase distribution calculation device 1. Although Fig. 1 shows the overall configuration of the three-phase distribution calculation system 100, it is also possible to use only the three-phase distribution calculation device 1 alone.

[0015] (Three-phase distribution calculation device 1) First, the three-phase distribution calculation device 1 will be described with reference to the drawings. As shown in Fig. 1, the three-phase distribution calculation device 1 includes a data acquisition unit 11, a calculation unit 12, and a control unit 20. The three-phase distribution calculation device 1 may also include an output unit 13, a communication unit 14, and a data table update unit 15.

[0016] The data acquiring unit 11 acquires soil permittivity data and thermal conductivity data. Hereinafter, the soil permittivity data and thermal conductivity data may be referred to as soil data or simply as data. The permittivity data and thermal conductivity data are actually measured soil data, and as will be described later, the data acquiring unit 11 acquires the soil data transmitted from a data logger An and received via the communication unit 14. Alternatively, the data acquiring unit 11 may acquire soil permittivity data and thermal conductivity data recorded in an external memory or database, etc., via the communication unit 14, which will be described later. Furthermore, when only the three-phase distribution calculation device 1 is used, the data acquiring unit 11 may acquire soil data previously recorded in the memory 22 by the user (i.e., soil data input by the user from an input device (not shown) of the three-phase distribution calculation device 1).

[0017] The calculation unit 12 calculates the three-phase distribution of the soil from the obtained dielectric constant data and thermal conductivity data. Specifically, when the volume fractions of the gas phase, liquid phase, and solid phase of the soil are expressed as G, L, and S, respectively, the calculation unit 12 solves the following simultaneous equations including G, L, and S.

[0018] The dielectric constant ε of the soil is calculated by dividing the dielectric constant ε of the gas, liquid, and solid phases of the soil by g , ε l , ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,ε l ,ε s ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by g , λ l , λ s The second functional formula (2) below is expressed as: λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third functional formula (3) is used, in which the total value of the volume ratios of the gas phase, liquid phase, and solid phase of the soil is 1: 1.0=G+L+S …(3) The volume ratios G, L, and S of the soil are calculated using the above formula.

[0019] The first and second functional formulas (1) and (2) are not particularly limited to specific formulas. For example, the following set of functional formulas (4), (5) and (3) may be used.

[0020]

number

[0021] Logarithm of The first function formula (4) expresses the dielectric constant ε of the soil as the dielectric constant ε of the gas, liquid, and solid phases of the soil. g , ε l , ε s The second function (5) expresses the thermal conductivity λ of the soil as a logarithm of the thermal conductivity λ of the gas, liquid, and solid phases of the soil. g , λ l , λ s The third functional formula (3) is a formula in which the total volume ratio of the gas, liquid, and solid phases of the soil is 1. Based on the prior art, the inventors considered that the permittivity of soil is the contribution of the permittivity of each phase in proportion to its own volume, and that the thermal conductivity of soil is the contribution of the thermal conductivity of each phase in proportion to its own volume. The above functional formulas (4), (5), and (3) are examples of formulas based on this consideration. Note that functional formulas (1) and (2) are not limited to the above functional formulas (4) and (5).

[0022] Dielectric constant ε of each soil phase g , ε l , ε s and the thermal conductivity of each phase of the soil, λ g , λ l , λ s The dielectric constant ε of each phase must be calculated in advance. This calculation method will be described later. g , ε l , ε s and the thermal conductivity of each phase λ g , λ l , λ s may be recorded in the memory 22 as a data table such as that shown in Table 1 below.

[0023] [Table 1]

[0024] The output unit 13 outputs information including the three-phase distribution calculated by the calculation unit 12 to an external display unit 30 or the like. The display unit 30 is a variety of displays, a printing device, or the like.

[0025] The communication unit 14 is an information communication interface that receives data from an external memory or a data logger An. The communication unit 14 may be, for example, Wi-Fi (registered trademark) or Bluetooth (registered trademark) for connecting to the Internet 50.

[0026] The data table update unit 15 updates the data table recorded in the memory 22. The data table stores data on the dielectric constant and thermal conductivity of each phase for each type of soil, as shown in Table 1. The data table update unit 15 may record new soil data in the data table, or may rewrite already recorded data. The new data may be data input by the user from the input device of the three-phase distribution calculation device 1, or may be various data recorded in an external database acquired via the communication unit.

[0027] The control unit 20 performs overall control of the three-phase distribution calculation device 1. The control unit 20 includes at least one processor 21 and at least one memory 22. The processor 21 can be configured using a general-purpose processor such as at least one MPU (Micro Processing Unit) or CPU (Central Processing Unit). The memory 22 may include multiple types of memory such as ROM (Read Only Memory) and RAM (Random Access Memory). As an example, the processor 21 realizes the functions of each unit by loading various control programs recorded in the ROM of the memory 22 into the RAM and executing them. Furthermore, the processor 21 may include a dedicated processor configured using an ASIC (Application Specific Integrated Circuit), FPGA (Field Programmable Gate Array), PLD (Programmable Logic Device), or the like.

[0028] (Three-phase distribution calculation system 100) Next, we will explain the three-phase distribution calculation system 100. As shown in Fig. 1, the three-phase distribution calculation system 100 includes a three-phase distribution calculation device 1, one or more sensors Bn that detect the dielectric constant and thermal conductivity of soil, and one or more data loggers An.

[0029] The sensor Bn may be a sensor that directly measures the permittivity and thermal conductivity of the soil, or one or more sensors that can calculate the permittivity and thermal conductivity, and the type may be arbitrary. For example, the sensor Bn may be a combination of a sensor that can measure the permittivity and a sensor that can measure the thermal conductivity. When two sensors, a permittivity measuring sensor and a thermal conductivity measuring sensor, are placed, they are placed at an appropriate distance so that the measured values ​​are not affected by each other. The sensor Bn may also be a combination of one or more sensors that acquire information that can be used to calculate the permittivity and thermal conductivity of the soil. The sensor Bn is embedded in the soil.

[0030] Data logger An is paired with sensor Bn and collects data measured by sensor Bn. Data logger An has a transmission function for transmitting data to the three-phase distribution calculation device 1 via the Internet 50, for example. Data logger An transmits data periodically or in response to an instruction from the three-phase distribution calculation device 1. The pair of sensor Bn and data logger An is placed in the soil to be measured. With this configuration, the dielectric constant and thermal conductivity of each soil can be obtained whenever necessary, and changes in these values ​​can be tracked.

[0031] (Dielectric constant and thermal conductivity of each phase) As mentioned above, the dielectric constant and thermal conductivity of each soil phase must be obtained in advance. This method is described below. First, the dielectric constant and thermal conductivity of the entire soil to be measured are measured at least once. Next, the soil is collected in a cylindrical container or similar, and the volume fraction of each phase is determined using conventional techniques. Any method for determining the volume fraction can be used. For example, when using the actual volume method, the volume and initial weight of the collected soil are measured. Next, the soil is dried by heating or other methods, and the weight after drying is measured. The moisture content can be determined by subtracting the dried weight from the initial weight. This moisture content is the liquid phase volume, and the volume is calculated as 1 g / cm. Next, the specific gravity and volume of the dried soil are determined. The volume of the gas phase is determined by subtracting the volume of water and the volume of soil from the volume of the collected soil. In this manner, the volume and volume fraction (G, L, S) of each phase can be determined. Alternatively, the volume fraction of each phase may be determined using the measuring device described in Patent Document 1.

[0032] Next, the measured dielectric constant, thermal conductivity, and volume fraction of the entire soil were substituted into the above function equations (1) and (2), and the dielectric constant ε of each phase of the soil was calculated by the optimization method. g , ε l , ε s and the thermal conductivity of each phase of the soil, λ g , λ l , λ sThe type of optimization method can be any, but the least squares method may be used, for example. Furthermore, when there is a large amount of actual measurement data for the volume fraction, dielectric constant, and thermal conductivity of each phase, an optimization program may be used to find the optimal (smallest error) dielectric constants εg, εl, and εs of each phase and the thermal conductivity λg, λl, and λs of each phase.

[0033] (Three-phase distribution calculation method S1) Next, a three-phase distribution calculation method S1 according to this embodiment will be described with reference to the drawings. Fig. 2 is a flowchart showing the flow of the soil three-phase distribution calculation method S1. As shown in the figure, the soil three-phase distribution calculation method S1 includes steps S11 to S14.

[0034] Step S11 is a preparatory step in which the dielectric constant and thermal conductivity are calculated in advance for each of the gas, liquid, and solid phases for each type of soil. These are obtained using the method described above. If data on the dielectric constant and thermal conductivity for each phase has already been obtained, step S11 may be omitted.

[0035] Step S12 is a step for acquiring the dielectric constant ε of the soil. Step S12 is executed by the data acquiring unit 11 of the three-phase distribution calculating device 1.

[0036] Step S13 is a step for acquiring the thermal conductivity λ of the soil. Step S13 is executed by the data acquisition unit 11 of the three-phase distribution calculation device 1. It is preferable that steps S12 and S13 are performed in parallel at the same timing. This is because if the measurements are taken at different times, the soil conditions may change. In the above-described three-phase distribution calculation system 100, the dielectric constant ε and the thermal conductivity λ of the soil can be acquired in parallel.

[0037] Step S14 calculates the dielectric constant ε of the soil by dividing the dielectric constants ε of the gas, liquid, and solid phases of the soil. g , ε l , ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,εl ,ε s ) …(1) The thermal conductivity of the soil λ is the thermal conductivity of the gas, liquid, and solid phases of the soil λ g , λ l , λ s The second function (2) below is expressed as λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third function formula (3) in which the total volume ratio of the gas phase, liquid phase, and solid phase of the soil is 1, 1.0=G+L+S …(3) Step S14 is a step of calculating the volume ratios G, L, and S of the soil using the above equations. Step S14 is executed by the calculation unit 12 of the three-phase distribution calculation device 1.

[0038] The specific formulas of the first functional formula (1) and the second functional formula (2) are not particularly limited. For example, the above-mentioned functional formulas (4) and (5) may be used instead of the above-mentioned first functional formula (1) and the second functional formula (2).

[0039] According to the three-phase distribution calculation device 1 and three-phase distribution calculation method S1 configured as described above, the three-phase distribution of soil can be calculated using soil data without actual measurement. By using the above-described three-phase distribution calculation device 1 or three-phase distribution calculation method S1, the three-phase distribution of soil can be calculated (estimated) in a simple manner. Conventionally, the three-phase distribution was calculated by sampling soil and actually measuring the volume fraction of each phase, but this work process can be greatly simplified. Therefore, it becomes easy to calculate the three-phase distribution of various soils at the required time. Furthermore, by continuously obtaining the three-phase distribution of soil, it can be useful for considering the type of crop suitable for each soil and for managing its growth. [Example]

[0040] An example in which the three-phase distribution was calculated using the above-mentioned method will be described. Figures 3 to 5 are graphs showing the correlation between the calculated (estimated) volume fractions of the solid, liquid, and gas phases in gray lowland soil, calculated according to the three-phase distribution calculation method S1, and the measured values ​​calculated using the actual volume method. The horizontal axis of the graph represents the measured values, and the vertical axis represents the estimated values. The correlation coefficient is shown in each graph. As shown in the figure, the correlation coefficient for the solid phase volume fraction is somewhat low, but the correlation coefficients for the liquid and gas phases are all 0.9 or higher, indicating that the estimations are highly accurate.

[0041] 6 to 8 are graphs showing the correlation between the calculated (estimated) values ​​of the volume fractions of the solid, liquid, and gas phases in a field obtained using a lysimeter according to the three-phase distribution calculation method S1 and the measured values ​​obtained using the actual volume method. As with FIGS. 3 to 5, the correlation coefficients between the liquid and gas phases are all 0.9 or higher, indicating that the estimations are highly accurate. As described above, it has been found that the three-phase distribution calculation method S1 of this embodiment can estimate a three-phase distribution equivalent to that obtained using the actual volume method.

[0042] [Software implementation example] The functions of the three-phase distribution calculation device 1 and the three-phase distribution calculation system 100 (hereinafter referred to as the "device") can be realized by a three-phase distribution calculation program, which is a program for making a computer function as the device, and which makes a computer function as each part executed by the control unit of the device.

[0043] In this case, the device includes a computer having at least one control device (e.g., processor 21) and at least one storage device (e.g., memory 22) as hardware for executing the three-phase distribution calculation program. The control device and storage device execute the three-phase distribution calculation program, thereby realizing the functions described in each of the above embodiments.

[0044] The three-phase distribution calculation program may be stored non-transitoryly in one or more computer-readable storage media. The storage media may or may not be included in the device. In the latter case, the program may be supplied to the device via any wired or wireless transmission medium.

[0045] Furthermore, some or all of the functions of the control blocks can be realized by logic circuits. For example, integrated circuits in which logic circuits functioning as the control blocks are formed are also included in the scope of the present invention. In addition, the functions of the control blocks can also be realized by, for example, a quantum computer.

[0046] (Additional notes) A method for calculating a three-phase distribution of soil according to a first aspect of the present invention includes: a permittivity acquisition step of acquiring a permittivity ε of the soil; A thermal conductivity acquisition step of acquiring the thermal conductivity λ of the soil; The volume fractions of the gas, liquid, and solid phases of the soil are G, L, and S, respectively, and the dielectric constant ε of the soil is expressed as the dielectric constant ε of the gas, liquid, and solid phases of the soil. g , ε l , ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,ε l ,ε s ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by g , λ l , λ s The second functional formula (2) below is expressed as: λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third functional formula (3) is used, in which the total value of the volume ratios of the gas phase, liquid phase, and solid phase of the soil is 1: 1.0=G+L+S …(3) and a volume ratio calculation step of calculating volume ratios G, L, and S of the soil using the above formula.

[0047] A method for calculating a three-phase distribution in soil according to a second aspect of the present invention is the method according to the first aspect, wherein the first functional formula is the following formula (4): logε=G*logε g +L*logε l +S*logε s …(4) The second function formula is the following formula (5): logλ=G*logλ g +L*logλ l +S*logλ s …(5).

[0048] The method for calculating a three-phase distribution of soil according to a third aspect of the present invention is the method for calculating a three-phase distribution of soil according to the first or second aspect, wherein the dielectric constant ε g , ε l , ε s and the thermal conductivity λ g , λ l , λ s The method further includes the preliminary step of determining:

[0049] A method for calculating a three-phase distribution in soil according to a fourth aspect of the present invention is the same as in the third aspect, wherein the preparation step includes a step of deriving the volume fractions G, L, and S of the soil by an actual measurement method, and a step of measuring the dielectric constant ε and the thermal conductivity λ of the soil at least once, and calculating the dielectric constant ε that satisfies the first functional formula (1) by substituting the derived volume fractions G, L, and S. g , ε l , ε s and the thermal conductivity λ that satisfies the second functional formula (2) into which the volume ratios G, L, and S are substituted. g , λ l , λ s and determining by an optimization method.

[0050] A three-phase distribution calculation device according to a fifth aspect of the present invention includes: a data acquisition unit that acquires dielectric constant data and thermal conductivity data of soil; The volume fractions of the gas, liquid, and solid phases of the soil are G, L, and S, respectively, and the dielectric constant ε of the soil is expressed as follows: g , ε l, ε s The first function formula (1) below is expressed as ε=f(G,L,S,ε g ,ε l ,ε s ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by g , λ l , λ s The second functional formula (2) below is expressed as: λ=g(G,L,S,λ g ,λ l ,λ s ) …(2) The following third functional formula (3) is used, in which the total value of the volume ratios of the gas phase, liquid phase, and solid phase of the soil is 1: 1.0=G+L+S …(3) and a calculation unit that calculates the volume ratios G, L, and S of the soil using the above.

[0051] A three-phase distribution calculation system according to a sixth aspect of the present invention includes the three-phase distribution calculation device according to the fifth aspect, a sensor for detecting the dielectric constant and thermal conductivity of soil, and a data logger capable of collecting data measured by the sensor and transmitting the data to the three-phase distribution calculation device.

[0052] A three-phase distribution calculation program according to aspect 7 of the present invention is a three-phase distribution calculation program for causing a computer to function as the three-phase distribution calculation device according to aspect 5, and is a program for causing a computer to function as the data acquisition unit and the calculation unit.

[0053] A non-transitory recording medium according to an eighth aspect of the present invention is a computer-readable non-transitory recording medium having the three-phase distribution calculation program according to the seventh aspect recorded thereon.

[0054] The present invention is not limited to the above-described embodiments, and various modifications are possible within the scope of the claims. Embodiments obtained by appropriately combining the technical means disclosed in different embodiments are also included in the technical scope of the present invention. [Explanation of symbols]

[0055] 1... Estimation device 11…Acquisition part 12...Generation section 13…Estimation part 14,30...Control unit 2...Data generation device 21…Acquisition part 22...Conversion unit 23...Generation section 24...Label assignment section 25...Size conversion section 26...Data Augmentation Section 27…Input / output section 60...Mechanical model 100…Display device

Claims

1. a permittivity acquisition step of acquiring the permittivity ε of the soil; A thermal conductivity acquisition step of acquiring the thermal conductivity λ of the soil; The volume fractions of the gas phase, liquid phase, and solid phase of the soil are G, L, and S, respectively, and the dielectric constant ε of the soil is expressed as the dielectric constant ε of the gas phase, liquid phase, and solid phase of the soil. ｇ , ε ｌ , ε ｓ The following first functional formula (1) is expressed as: ε=f(G,L,S,ε) ｇ ,he ｌ ,he ｓ ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by ｇ , λ ｌ , λ ｓ The following second function formula (2) is expressed as: λ=ﺇ(G, L, S, λ ｇ ,l ｌ ,l ｓ ) …(2) The following third functional formula (3) is defined as a sum of the volume ratios of the gas phase, liquid phase, and solid phase of the soil, which is 1: 1.0=G+L+S…(3) a volume ratio calculation step of calculating volume ratios G, L, and S of the soil using the above formula; A method for calculating three-phase distribution in soil, including:

2. The first functional formula is the following formula (4): logε=G*logε ｇ +L*logε ｌ +S*logε ｓ …(4) The second function formula is the following formula (5): logλ=G*logλ ｇ +L*logλ ｌ +S*logλ ｓ …(5) The method for calculating three-phase distribution in soil according to claim 1 .

3. For each type of soil, the dielectric constant ε ｇ , ε ｌ , ε ｓ and the thermal conductivity λ ｇ , λ ｌ , λ ｓ The method for calculating a three-phase distribution in soil according to claim 1 or 2, further comprising a preparation step of determining:

4. The preparation step includes a step of deriving the volume ratios G, L, and S of the soil by an actual measurement method; The dielectric constant ε and the thermal conductivity λ of the soil are each measured at least once, and the dielectric constant ε that satisfies the first functional formula (1) is calculated by substituting the derived volume fractions G, L, and S. ｇ , ε ｌ , ε ｓ and the thermal conductivity λ that satisfies the second function formula (2) into which the volume ratios G, L, and S are substituted. ｇ , λ ｌ , λ ｓ and calculating the three-phase distribution of soil according to claim 3 by an optimization method.

5. a data acquisition unit that acquires dielectric constant data and thermal conductivity data of soil; The volume fractions of the gas phase, liquid phase, and solid phase of the soil are G, L, and S, respectively, and the dielectric constant ε of the soil is expressed as follows: ｇ , ε ｌ , ε ｓ The following first functional formula (1) is expressed as: ε=f(G,L,S,ε) ｇ ,he ｌ ,he ｓ ) …(1) The thermal conductivity λ of the soil is calculated by dividing the thermal conductivity λ of the gas phase, liquid phase, and solid phase of the soil by ｇ , λ ｌ , λ ｓ The following second function formula (2) is expressed as: λ=ﺇ(G, L, S, λ ｇ ,l ｌ ,l ｓ ) …(2) The following third functional formula (3) is defined as a sum of the volume ratios of the gas phase, liquid phase, and solid phase of the soil, which is 1: 1.0=G+L+S…(3) A calculation unit that calculates the volume ratios G, L, and S of the soil using the above formula. A three-phase distribution calculation device comprising:

6. The three-phase distribution calculation device according to claim 5 ; a sensor for detecting the dielectric constant and thermal conductivity of the soil; a data logger capable of collecting data measured by the sensor and transmitting the data to the three-phase distribution calculation device; A three-phase distribution calculation system comprising:

7. A three-phase distribution calculation program for causing a computer to function as the three-phase distribution calculation device according to claim 5, the three-phase distribution calculation program causing a computer to function as the data acquisition unit and the calculation unit.

8. A computer-readable non-transitory recording medium having the three-phase distribution calculation program according to claim 7 recorded thereon.

Images