WO2023063044A1 - 測定値温度係数補正装置及び測定値温度係数補正プログラム - Google Patents

測定値温度係数補正装置及び測定値温度係数補正プログラム Download PDF

Info

Publication number
WO2023063044A1
WO2023063044A1 PCT/JP2022/035234 JP2022035234W WO2023063044A1 WO 2023063044 A1 WO2023063044 A1 WO 2023063044A1 JP 2022035234 W JP2022035234 W JP 2022035234W WO 2023063044 A1 WO2023063044 A1 WO 2023063044A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
temperature
temperature coefficient
conductivity
brix
Prior art date
Application number
PCT/JP2022/035234
Other languages
English (en)
French (fr)
Inventor
敏夫 英
Original Assignee
ディー・クルー・テクノロジーズ株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ディー・クルー・テクノロジーズ株式会社 filed Critical ディー・クルー・テクノロジーズ株式会社
Priority to JP2023555064A priority Critical patent/JPWO2023063044A1/ja
Publication of WO2023063044A1 publication Critical patent/WO2023063044A1/ja

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q11/00—Accessories fitted to machine tools for keeping tools or parts of the machine in good working condition or for cooling work; Safety devices specially combined with or arranged in, or specially adapted for use in connection with, machine tools
    • B23Q11/10—Arrangements for cooling or lubricating tools or work
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18—Water

Definitions

  • the present invention relates to a measured value temperature coefficient correction device and a measured value temperature coefficient correction program, for example, a measured value temperature coefficient correction device for correcting temperature-dependent measured values among the characteristic values of a solution according to liquid temperature. and a measured value temperature coefficient correction program.
  • coolant liquid (cutting oil) is used as a coolant or lubricating oil in machining processes such as cutting and polishing processes.
  • water-soluble coolants and water-insoluble coolants are used after being diluted about 5 to 20 times.
  • the diluted coolant is mixed with water in a state where the undiluted coolant is emulsified.
  • the coolant diluted to an appropriate concentration is set to be alkaline with a pH of about 8 to 9, which is effective in preventing rust and bacteria.
  • a circulation system is normally built in which the coolant liquid is stored in a coolant tank, sent to the processing chamber through a pipe by a supply pump when used, and returned to the coolant tank through a recovery pipe after use.
  • the factory can reuse the coolant liquid over a long period of time.
  • the concentration of the water-soluble coolant liquid increases as the water content evaporates due to the heat generated during machining, and conversely, the coolant component adheres to the cutting waste and is removed together with the removal filter, resulting in a lower concentration.
  • Patent Document 1 discloses a technique for measuring the concentration of a coolant liquid using a conductivity meter.
  • the coolant liquid management device described in Patent Document 1 is obtained from an electric field generated in the secondary side toroidal transformer by the action of a full-circle current loop generated in the coolant liquid when an electric field is applied to the primary side toroidal transformer immersed in the coolant liquid.
  • An electromagnetic induction type conductivity meter that measures the amount of conductivity of the coolant liquid based on the voltage difference, a calculation unit that calculates the BRIX value corresponding to the amount of conductivity by applying the amount of conductivity to a conversion formula given in advance, and BRIX and an external communication interface for outputting the value to the outside.
  • One aspect of the measured value temperature coefficient correction device is a temperature measurement value obtained from a temperature sensor that measures the temperature of a solution, and a characteristic value of the solution measured by an electrical conduction method, the value of which varies depending on the temperature.
  • a measured value memory for storing in chronological order in association with measured solution values having varying temperature dependence, and a temperature derived from statistical processing of the measured solution values and the measured temperature values stored in the measured value memory.
  • a temperature coefficient calculator for calculating a coefficient; and a measured value temperature corrector for outputting the temperature-corrected solution measured value obtained by correcting the solution measured value using the temperature coefficient.
  • One aspect of the measured value temperature coefficient correction program is a measured value temperature coefficient correction program that corrects the measured value using the temperature coefficient by being executed by the calculation unit, and measures the temperature of the solution.
  • a temperature measurement value obtained from a temperature sensor is associated with a solution measurement value that is a characteristic value of the solution measured by the electrical conduction method and has a temperature dependence in which the value changes depending on the temperature, and the measured value memory is stored in chronological order.
  • a temperature coefficient calculation process for calculating a temperature coefficient derived from statistical processing for the solution measurement value and the temperature measurement value stored in the measurement value memory; and using the temperature coefficient and a measurement value temperature correction process of outputting the temperature-corrected solution measurement value obtained by correcting the solution measurement value.
  • the measured value temperature coefficient correction device and the measured value temperature coefficient correction program according to the present invention perform statistical processing on the solution measured value and the temperature measured value to calculate the temperature coefficient, thereby correcting the deviation of the solution measured value due to temperature. to correct.
  • the device for correcting the temperature coefficient of measured values and the program for correcting the temperature coefficient of measured values according to the present invention can correct the deviation of the measured conductivity caused by the liquid temperature.
  • FIG. 1 is a block diagram of a coolant liquid management device according to a first embodiment
  • FIG. 3 is a block diagram of a conductivity meter of the coolant management device according to the first embodiment
  • FIG. 4 is a block diagram of a correlation BRIX calculator according to the first embodiment
  • FIG. 4 is a graph for explaining the relationship between electrical conductivity and BRIX using a conversion formula used in the correlation BRIX calculator according to the first embodiment
  • 7 is a graph for explaining a temporal change function used in the correlation BRIX calculator according to the first embodiment
  • FIG. 11 is a block diagram of a correlation BRIX calculation unit according to the second embodiment
  • FIG. 9 is a graph for explaining the relationship between conductivity and BRIX using a conversion formula used in the correlation BRIX calculator according to the second embodiment
  • 9 is a graph for explaining a time-dependent change function used in a correlation BRIX calculator according to the second embodiment
  • FIG. 11 is a block diagram of a correlation BRIX calculation unit according to a third embodiment
  • FIG. FIG. 11 is a block diagram of a coolant management device according to a fourth embodiment
  • FIG. FIG. 12 is a block diagram of a correlation BRIX calculation unit according to a fourth embodiment
  • FIG. 11 is a block diagram of a measured value temperature coefficient correction unit according to a fourth embodiment
  • FIG. 13 is a diagram for explaining a format of data stored in a measured value memory of a measured value temperature coefficient correction unit according to the fourth embodiment
  • FIG. 11 is a diagram for explaining a format of data stored in a matrix memory of a measured value temperature coefficient correction unit according to the fourth embodiment
  • FIG. 13 is a diagram for explaining a format of data stored in a matrix memory of a measured value temperature coefficient correction unit according to the fourth embodiment
  • the programs described above include instructions (or software code) that, when read into a computer, cause the computer to perform one or more functions described in the embodiments.
  • the program may be stored in a non-transitory computer-readable medium or tangible storage medium.
  • computer readable media or tangible storage media may include random-access memory (RAM), read-only memory (ROM), flash memory, solid-state drives (SSD) or other memory technology, CDs - ROM, digital versatile disc (DVD), Blu-ray disc or other optical disc storage, magnetic cassette, magnetic tape, magnetic disc storage or other magnetic storage device.
  • the program may be transmitted on a transitory computer-readable medium or communication medium.
  • transitory computer readable media or communication media include electrical, optical, acoustic, or other forms of propagated signals.
  • FIG. 1 shows a schematic diagram for explaining the outline of a coolant liquid management device 1 according to the first embodiment.
  • two boxes ie, a control box 10 and a sensor box 20 are connected by a suspension wire Wi and a connecting jig.
  • the control box 10 is provided with a display section IND.
  • the coolant liquid management device 1 is connected to the weight Wt by a suspension wire Wi hung on a pulley. By moving the suspension wire Wi up and down, the coolant liquid management device 1 moves up and down in the coolant tank.
  • the sensor box 20 is equipped with various sensors for measuring the electrical conductivity and the like of the coolant while it is immersed in the coolant.
  • the sensor box 20 is equipped with a conductivity meter.
  • the conductivity meter is called a conductivity meter.
  • the control box 10 controls various sensors mounted in the sensor box 20, and is equipped with a communication interface and the like for communicating with the outside. Further, the control box 10 is provided at a position where the sensor box 20 is not immersed in the coolant liquid even when the sensor box 20 is immersed in the coolant liquid.
  • FIG. 2 shows a block diagram of the coolant management device 1 according to the first embodiment.
  • the control box 10 has a power source section 11, a computing section 12, a communication interface 13, an external communication interface 14, and a display section IND.
  • the sensor box 20 also has a power supply unit 21 , a correlation BRIX calculator 22 , a control signal generator 23 , a conductivity meter 24 , an analog-to-digital conversion circuit 25 and a communication interface 26 .
  • the power supply unit 11 distributes power supplied from an external AC adapter to each processing block of the control box 10 (not shown).
  • the power supply unit 11 also supplies power to the power supply unit 21 of the sensor box 20 .
  • the power supply unit 21 distributes the power supplied to each processing block of the sensor box 20 (not shown).
  • the computing unit 12 is, for example, a computing unit capable of executing a program such as a microcontroller.
  • the calculation unit 12 controls, for example, the measurement sequence of the coolant liquid BRIX measurement process executed in the coolant liquid management device 1 . Further, the calculation unit 12 performs control related to communication processing via the communication interface 13 and the external communication interface 14, and information display processing on the display unit IND.
  • the communication interface 13 is an interface circuit that performs specific communication processing for communication between the calculation unit 12 and the correlation BRIX calculation unit 22 via the communication interface 26 on the sensor box 20 side.
  • the external communication interface 14 is an interface circuit for performing wired or wireless communication with a host system (local server, cloud server, etc.).
  • the display unit IND is for displaying the operating status of the coolant management device 1, measurement information, and the coolant status.
  • the correlation BRIX calculation unit 22 is, for example, a calculation unit such as a microcontroller that can execute a program.
  • the correlation BRIX calculation unit 22 performs, for example, a BRIX calculation process for calculating BRIX based on the conductivity value obtained from the conductivity meter 24, a conductivity for converting the signal obtained from the conductivity meter 24 into a specific numerical value, Calculation processing, etc. are performed. The details of this BRIX calculation process will be described later. Further, the correlation BRIX calculator 22 instructs the control signal generator 23 whether or not to operate the conductivity meter 24 .
  • the control signal generator 23 gives specific instructions for operation to the conductivity meter 24 based on the instructions from the correlation BRIX calculator 22 .
  • the conductivity meter 24 measures the conductivity of the coolant while being immersed in the coolant.
  • the analog-to-digital conversion circuit 25 converts the output signal having an analog value output from the conductivity meter 24 into a signal having a digital value.
  • the communication interface 26 is an interface circuit that performs specific processing for the correlation BRIX calculator 22 to communicate with the calculator 12 .
  • FIG. 3 shows a block diagram of the conductivity meter of the coolant management device according to the first embodiment.
  • the conductivity meter 24 has a sine wave generator 30 , an amplifier 31 , a primary side toroidal transformer 32 , a secondary side toroidal transformer 33 , a low noise amplifier 34 and a signal amplitude detector 35 .
  • the sine wave generator 30 receives an instruction from the control signal generator 23 and outputs a sine wave signal for operating the primary side toroidal transformer 32 .
  • the amplifier 31 converts the sine wave signals supplied from the amplifier 31 into differential signals having phases different from each other by 180 degrees to drive the primary side toroidal transformer 32 .
  • the primary side toroidal transformer 32 generates an electric field based on the differential signal given from the amplifier 31 .
  • an electric field generated in the primary side toroidal transformer 32 causes a current loop that circulates through the inner circumferences of the primary side toroidal transformer 32 and the secondary side toroidal transformer 33 in the coolant liquid.
  • the secondary-side toroidal transformer 33 outputs to the low-noise amplifier 34 an output signal corresponding to the electric field generated due to the current loop generated by the primary-side toroidal transformer 32 .
  • the low noise amplifier 34 amplifies the output signal of the secondary side toroidal transformer 33 and supplies it to the signal amplitude detector 35 .
  • the signal amplitude detector 35 detects the amplitude of the signal given from the low noise amplifier 34 and outputs the detected signal to the analog-to-digital conversion circuit 25 as a measurement result.
  • the BRIX value is calculated using the conductivity obtained by the conductivity meter 24 . Therefore, the calculation method of BRIX in Embodiment 1 will be described in detail.
  • BRIX is a physical quantity used as a sugar content mainly in the food industry to measure the content of so-called sugars such as sucrose, fructose, invert sugar and glucose.
  • the BRIX value is defined as a value corresponding to the mass percentage of a sucrose solution at 20°C. When 100 g of an aqueous solution containing only 1 g of sucrose as a solute is measured with a BRIX refractometer, the reading BRIX value is 1%.
  • the conductivity is completely different from the value measured by the refractometer. Calculate the corresponding correlated BRIX value. Therefore, the correlation BRIX calculator 22 will be described in detail below. In the following description, among the functions of the correlation BRIX calculation unit 22, the part related to the correlation BRIX value will be described in detail, and the description of the part related to the control of the conductivity meter 24 will be omitted.
  • the correlation BRIX calculation unit 22 measures the conductivity and the BRIX value of the unused coolant immediately after the coolant is changed (total replacement of the coolant), and shows the relationship between the conductivity and the BRIX.
  • a regression function is created, and a BRIX approximation coefficient applied to a conversion formula (for example, a BRIX approximation function) used to calculate the correlation BRIX value Bd is calculated from this regression function.
  • a correlated BRIX value Bd is calculated from the conductivity using a BRIX approximation function to which the calculated BRIX approximation coefficient is applied.
  • the correlated BRIX calculator 22 collects parameter information to be used for measurement as a preparation before starting the measurement process for outputting the correlated BRIX value Bd of the coolant.
  • This parameter information includes the date and time D0 when the coolant liquid was added to the coolant tank, the reference conductivity which is the conductivity of the coolant liquid in the coolant tank measured by the conductivity meter 24 at a specific reference date and time such as before the start of measurement.
  • the correlation BRIX calculator 22 calculates the temporal change coefficient k using these parameter information before starting the measurement.
  • the correlation BRIX calculator 22 acquires the conductivity measurement value Cd of the coolant liquid with the conductivity meter 24, and corrects the acquired conductivity measurement value Cd by applying the coefficient of change over time k. to generate a corrected conductivity Cdn, and a correlated BRIX value Bd calculated by applying this corrected conductivity to the BRIX approximation function is output as the BRIX value of the coolant.
  • FIG. 4 shows a block diagram of the correlation BRIX calculator 22 according to the first embodiment.
  • the correlation BRIX calculation unit 22 includes a data storage unit 41, a parameter storage unit 42, an internal clock 43, an elapsed time counter 44, a temporal change coefficient calculation unit 45, a conductivity conversion unit 46, a correlation BRIX conversion unit. 47.
  • the correlation BRIX calculator 22 is, for example, a calculator that can execute a program like a microcomputer.
  • the data storage unit 41 and the parameter storage unit 42 can use a hard disk in the calculation unit, a non-volatile storage device such as an SSD (Solid-State Drive), or a volatile memory such as a DRAM.
  • other configuration blocks can be implemented by programs (for example, BRIX measurement program).
  • the data storage unit 41 is a storage unit that stores information on the coolant liquid as a library.
  • the data storage unit 41 holds information in which the BRIX approximation coefficient of the coolant liquid investigated in advance is associated with the manufacturer name and model number (or part number) of the coolant liquid.
  • the BRIX approximation coefficient is a coefficient included in the correlation BRIX approximation function created by using a separately prepared conductivity meter and BRIX meter through prior verification.
  • the correlation BRIX calculation unit 22 uses the quadratic function represented by Equation (1) as the correlation BRIX approximation function.
  • FIG. 5 shows a graph for explaining the relationship between conductivity and BRIX using a conversion formula used in the correlation BRIX calculator according to the first embodiment.
  • the BRIX approximation coefficient is determined so that equation (1) follows the curve of the BRIX approximation function obtained from a plurality of measurement points measured during measurement. The starting point of the curve shown in FIG.
  • the correlation BRIX calculation unit 22 applies the temporal change coefficient k to maintain a constant relationship between the electrical conductivity obtained from the equation (1) and the correlation BRIX value Bd.
  • the parameter storage unit 42 stores parameters used by the elapsed time counter 44 and the conductivity conversion unit 46. These parameters store coolant information, fluid change date and time D0, reference measurement date and time Dr1, reference conductivity Cr, dilution water conductivity Cw, reference BRIX value Br, and aging constants D1 and D2. Even if these parameters automatically collect information obtained from the conductivity meter 24, the separately prepared BRIX meter, the coolant information stored in the data storage unit 41, and the internal clock 43, the control box 10 side may be given from an input device (not shown) such as the display unit IND. Note that the fluid change date and time D0 and the aging constants D1 and D2 are information input from the display unit IND.
  • the reference measurement date Dr1 is the date and time when the reference conductivity Cr and the reference BRIX value Br are measured.
  • the reference conductivity Cr is the conductivity of the coolant injected into the coolant tank measured by the conductivity meter 24 on the reference measurement date and time Dr1.
  • the reference BRIX value Br is a BRIX value measured by an external BRIX meter on the reference measurement date and time Dr1 for the coolant liquid injected into the coolant tank.
  • the dilution water conductivity Cw is the conductivity measured by the conductivity meter 24 of the dilution water injected when the coolant injected into the coolant tank is diluted.
  • the time-dependent change constants D1 and D2 are values set when using the system. For example, the time-dependent change constant D1 is a value of about 50 days, and the time-dependent change constant D2 is a value of about 300 days. It is preferable to specify periods that vary in degree.
  • the built-in clock 43 counts the date and time and outputs current time information Dn. A part of this current time information Dn can also be used as the reference measurement date and time Dr1.
  • a time-dependent change coefficient calculation unit 45 calculates a time-dependent change coefficient k that corrects the amount of change in the conductivity from the reference date to the present using a time-dependent change function created in advance from the measurement results of the change in the conductivity of the coolant over time. do.
  • the temporal change coefficient calculator 45 applies the reference conductivity Cr, the reference BRIX value Br, the dilution water conductivity Cw, and the reference measurement date and time Dr1 to the temporal change function to calculate the temporal change coefficient k.
  • FIG. 6 shows a graph for explaining the temporal change function used in the correlation BRIX calculator according to the first embodiment.
  • the temporal change function ⁇ (D) is a function in which D is the number of days that have passed since the liquid change date and time D0.
  • the temporal change function ⁇ (D) has a large rate of change in a period in which the number of elapsed days D is small, and decreases as the number of elapsed days D increases.
  • Equation (2) includes the corrected conductivity Cdn output by the conductivity converter 46 .
  • This corrected conductivity Cdn is a value corresponding to the conductivity of the coolant liquid measured when formula (1) was derived.
  • the conductivity measurement value Cd is the conductivity of the coolant liquid measured at the elapsed time D. Equation (3) can be approximated to Equation (4) using the temporal change coefficient k and the temporal change constants D1 and D2.
  • Crn in equation (5) is the reference initial conductivity, which can be calculated by applying the reference BRIX value Br to equation (1) and from the solution formula for equation (1).
  • the reference initial conductivity Crn can be calculated by the formula (6).
  • the time-dependent change coefficient calculator 45 calculates the reference conductivity Cr, the diluted water conductivity Cw, the time-dependent change constants D1 and D2, and the reference BRIX value Br read from the parameter storage unit 42 by the equations (5) and (6). Apply the equation to calculate the coefficient of change over time k.
  • the conductivity conversion unit 46 acquires the conductivity measurement value Cd from the conductivity meter 24 that measures the conductivity of the coolant liquid, applies the coefficient of change over time to the conductivity measurement value Cd, and converts the conductivity measurement value Cd , the corrected conductivity Cdn is calculated from which the influence of the concentration change of the electrolytic component over time is removed. Specifically, the conductivity converter 46 calculates the corrected conductivity Cdn by performing the following calculation.
  • the corrected conductivity Cdn can be expressed by the formula (7) by modifying the formula (2) and using the formula (5).
  • the correlated BRIX converter 47 applies the corrected conductivity Cdn to a predetermined conversion formula and outputs a correlated BRIX value Bd corresponding to the corrected conductivity Cdn.
  • the predetermined conversion formula is obtained by applying the pre-calculated BRIX approximation coefficient to formula (1), and is represented by formula (8).
  • the correlation BRIX calculation unit 22 derives from the change over time of the coolant liquid, and uses the change coefficient k with time, whose value changes with time, to determine the electrolytic component from the conductivity measurement value Cd. Eliminates the effects of concentration changes. As a result, the correlation BRIX calculator 22 can output a highly accurate correlation BRIX value Bd with little deviation from the actual BRIX of the coolant liquid without changing the correlation BRIX approximation function prepared in advance.
  • the temporal change coefficient k calculated by the correlation BRIX calculation unit 22 is a function of time, and can be automatically changed with time by counting the measurement time. That is, by using the correlation BRIX calculator 22 according to the first embodiment, it is possible to continuously measure the BRIX of the coolant liquid without adjusting parameters.
  • Embodiment 2 describes a correlation BRIX calculation unit 22a, which is another form of the correlation BRIX calculation unit 22 described in Embodiment 1.
  • FIG. 2 the same components as those described in Embodiment 1 are denoted by the same reference numerals as in Embodiment 1, and descriptions thereof are omitted.
  • unused coolant is used in calculating the BRIX approximation function. Created based on the coolant liquid.
  • the reference date is a day when the coolant is already in use, such as the installation date of the coolant management device 1, and the electrical conductivity and BRIX of the coolant are measured. day.
  • FIG. 7 shows a block diagram of the correlation BRIX calculator 22a according to the second embodiment.
  • the correlation BRIX calculation unit 22a according to the second embodiment includes a parameter storage unit 52, a built-in clock 43, an elapsed time counter 54, a temporal change coefficient calculation unit 55, a conductivity conversion unit 56, a correlation BRIX conversion It has a portion 57 .
  • the parameter storage unit 52 stores parameters that the temporal change coefficient calculation unit 55 uses for calculation.
  • the parameters stored in the parameter storage unit 52 include the first reference measurement date and time Dr1, the first reference conductivity Cr1, the post-dilution conductivity Cr1′, the first reference BRIX value Br1, the post-dilution BRIX value Br1′, the second reference Measurement date and time Dr2, second reference conductivity Cr2, second reference BRIX value Br2, and dilution water conductivity Cw are included.
  • the reference date is set to the date when the coolant is already in use, such as the date when the coolant management device 1 is installed, and the date when the electrical conductivity and BRIX of the coolant are measured. do.
  • the date and time of this measurement date is defined as a first reference measurement date and time Dr1.
  • the post-dilution conductivity Cr1 ' is the conductivity of a diluted coolant obtained by further diluting the coolant for which the first reference conductivity Cr1 was measured to a preset concentration (for example, a concentration of 1/2).
  • the first reference BRIX value Br1 is the BRIX value of the coolant measured on the reference date.
  • the post-dilution BRIX value Br1' is the BRIX value of the coolant for which the post-dilution conductivity Cr1' was measured.
  • the correlation BRIX calculation unit 22a on a day after the first reference measurement date and time Dr1 and on a day when it is determined that the difference between the correlation BRIX value Bd and the actual BRIX value has increased, Carry out various measurements.
  • This remeasurement date is the reference measurement date Dr2.
  • the second reference conductivity Cr2 is the conductivity measurement value Cd measured on the reference measurement date and time Dr2.
  • the second reference BRIX value Br2 is the BRIX measured at the reference measurement date and time Dr2.
  • the diluent water conductivity Cw is the conductivity of the diluent used when diluting the coolant.
  • the elapsed time counter 54 outputs the length of time that has elapsed since the first reference measurement date and time Dr1 as the elapsed time D.
  • the temporal change coefficient calculator 55 converts the second reference conductivity Cr2, the second reference BRIX value Br2, the dilution water conductivity Cw, the reference measurement date Dr1, and the reference measurement date Dr2 into a temporal change function ⁇ (D ) to calculate the coefficient of change over time k.
  • the temporal change coefficient is set to 0, that is, the temporal change function ⁇ (D) is set to 0.
  • the temporal change coefficient calculator 55 also calculates BRIX approximation coefficients ar, br, and cr that are used in the correlated BRIX converter 57 .
  • the correlated BRIX calculator 22a requires the BRIX approximation coefficients ar, br, and cr in order for the correlated BRIX converter 57 to calculate the correlated BRIX value Bd at the start of operation. Therefore, first, a method for calculating the BRIX approximation coefficients ar, br, and cr will be described.
  • a BRIX approximation function used in the correlated BRIX transform unit 57 is represented by equation (9).
  • Br in the formula (9) is the BRIX of the coolant at the first reference measurement date and time Dr1
  • Cr is the conductivity of the coolant at the first reference measurement date and time Dr1.
  • FIG. 8 shows a graph for explaining the relationship between conductivity and BRIX using a conversion formula used in the correlation BRIX calculator 22a according to the second embodiment.
  • the BRIX approximation function in Embodiment 2 has a dilution water conductivity Cw as a starting point, a first point determined by a first reference conductivity Cr1 and a first reference BRIX value Br1, and a post-dilution conductivity It is set along a curve connecting Cr1' and a second point determined by the post-dilution BRIX value Br1'.
  • the correlation BRIX calculation unit 22a converts the conductivity measurement value Cd into the corrected conductivity Cdr using the temporal change coefficient k, so that the relationship between the conductivity measurement value Cd and the correlated BRIX value Bd at the reference measurement date Dr1 is Correction processing is performed so that the relationship of the BRIX approximation function is obtained.
  • k 0 in a period in which the elapsed time D is shorter than the elapsed time up to the second reference measurement date Dr2. Then, when the elapsed time D is greater than or equal to the elapsed time up to the second reference measurement date and time Dr2, a non-zero temporal change coefficient k is used. In the following description, a method for calculating k that is not 0 will be described.
  • Equations (13) and (14) include the corrected conductivity Cdr output by the conductivity converter 56 .
  • This corrected conductivity Cdr is a value corresponding to the conductivity of the coolant measured when deriving the formula (9) (for example, the first reference conductivity Cr1).
  • the conductivity measurement value Cd is the conductivity of the coolant liquid measured at the elapsed time D.
  • equation (14) can be approximated to equation (15) using the temporal change coefficient k.
  • Cr2n in equation (16) is the second reference date initial conductivity, which can be calculated by applying the second reference BRIX value Br2 to equation (9) and from the solution formula for equation (9).
  • the reference initial conductivity Cr2n can be calculated by the formula (17).
  • the temporal change coefficient calculation unit 55 applies the parameters read from the parameter storage unit 52 to the equations (16) and (17) to calculate the temporal change coefficient k.
  • FIG. 9 shows a graph for explaining the temporal change function used in the correlation BRIX calculator 22a according to the second embodiment.
  • the temporal change coefficient ⁇ (D) according to the second embodiment is in the period from the first reference measurement date and time Dr1 to the second reference measurement date and time Dr2. will be 0 if
  • the elapsed time D exceeds the second reference measurement date and time Dr2
  • the slope rises linearly with the passage of time with the temporal change coefficient k.
  • the conductivity conversion unit 56 acquires a conductivity measurement value Cd from a conductivity meter that measures the conductivity of the coolant liquid, applies the coefficient of change over time to the conductivity measurement value Cd, and converts the conductivity measurement value Cd into A corrected conductivity Cdr is calculated from which the influence of the concentration change of electrolytic components over time is removed. Specifically, the conductivity converter 56 calculates the corrected conductivity Cdr by performing the following calculation.
  • the corrected conductivity Cdr can be expressed by the formula (18) by modifying the formula (13) and using the formula (16).
  • the temporal change coefficient k is 0 when the elapsed time is D ⁇ Dr2.
  • the correlated BRIX converter 57 applies the corrected conductivity Cdr to a predetermined conversion formula and outputs a correlated BRIX value Bd corresponding to the corrected conductivity Cdr.
  • the predetermined conversion formula is obtained by applying the BRIX approximation coefficient calculated in advance to the formula (9), and is represented by the formula (19).
  • the correlation BRIX calculation unit 22a derives from the change over time of the coolant liquid, and uses the change coefficient k with time, whose value changes over time, to determine the electrolytic component from the measured conductivity value Cd. Eliminates the effects of concentration changes. As a result, the correlation BRIX calculation unit 22a can output the correlation BRIX value Bd with little deviation from the actual BRIX value of the coolant liquid and high accuracy without changing the correlation BRIX approximation function created on the day of installation or the like. can.
  • the temporal change function ⁇ (D) according to the second embodiment calculates the temporal change coefficient k from the linearly approximated temporal change characteristic
  • the correlation calculated by the correlation BRIX calculation unit 22a according to the second embodiment The BRIX value Bd is less accurate than the correlation BRIX calculator 22 according to the first embodiment.
  • the correlation BRIX calculation unit 22a calculates the BRIX approximation coefficient and the coefficient of change over time k using the conductivity and BRIX value that can be obtained from the coolant in use.
  • the coolant management device using the correlation BRIX calculation unit 22a does not need to prepare in advance the BRIX approximation characteristics of the coolant in use when not in use, unlike the first embodiment.
  • the time and effort required for installation can be greatly simplified compared to the coolant liquid management device according to
  • Embodiment 3 describes a correlation BRIX calculation unit 22b, which is another form of the correlation BRIX calculation unit 22 described in Embodiment 1.
  • FIG. in the description of the third embodiment the same components as those described in the first embodiment are assigned the same reference numerals as in the first embodiment, and the description thereof is omitted.
  • FIG. 10 shows a block diagram of the correlation BRIX calculator 22b according to the third embodiment.
  • the correlation BRIX calculation unit 22b according to the third embodiment is obtained by adding a dilution water correction unit 68 to the correlation BRIX calculation unit 22 according to the first embodiment. Further, in the correlation BRIX calculation unit 22b, with the addition of the dilution water correction unit 68, the parameter storage unit 42, the time-dependent change coefficient calculation unit 45, and the correlation BRIX conversion unit 47 are replaced with the parameter storage unit 62 and the time-dependent change coefficient calculation unit 65. , has a correlation BRIX transform unit 67 .
  • Coolant information, liquid change date and time D0, reference measurement date and time Dr1, reference conductivity Cr, dilution water conductivity Cw, reference BRIX value Br, temporal change constants D1 and D2, specified parameters are stored in the parameter storage unit 62.
  • Conductivity C1 is included.
  • the temporal change coefficient calculator 65 calculates the reference conductivity C0 in addition to calculating the temporal change coefficient k.
  • the dilution water correction unit 68 replaces the dilution water conductivity, which indicates the conductivity of the dilution water obtained by diluting the coolant liquid, with the reference conductivity C0 when the conversion formula (e.g., formula (8)) of the correlated BRIX conversion unit is determined.
  • the conductivity of the diluent water obtained by diluting the coolant actually poured into the coolant tank is defined as the diluent water conductivity Cw, and the difference is corrected.
  • the dilution water correction unit 68 outputs a dilution water corrected conductivity Cdnw as the corrected conductivity.
  • the correlated BRIX converter 67 calculates the reference BRIX value Bd using the diluted water corrected conductivity Cdnw as the corrected conductivity Cdn of the BRIX approximation function shown in equation (8).
  • the reference conductivity C0 is the conductivity of diluted water obtained by diluting the coolant measured when creating the formula (8).
  • the designated conductivity C1 is a value input from an input unit such as the display unit IND as a fixed value.
  • This designated conductivity C1 is a parameter for adjusting the slope of the BRIX approximation function, and is determined by verifying the BRIX approximation function after correction by the dilution water correction unit 68 .
  • the reference conductivity C0 is derived from the formula for solving the equation (8) using the BRIX approximation coefficients an, bn, and cn, and is calculated based on the equation (20).
  • This reference conductivity C0 is the conductivity at which the correlation BRIX value Bd calculated by the equation (8) is zero.
  • the dilution water correction unit 68 performs correction to eliminate the difference between the reference conductivity C0 and the dilution water conductivity Cw. Specifically, the dilution water correction unit 68 performs correction processing based on the equation (21) and outputs the dilution water corrected conductivity Cdnw. On the other hand, the calculation of the temporal change coefficient k also requires correction. (5) with reference to formula, Here, Crnw in the formula (22) corresponds to the reference initial conductivity Crn in the formula (8). Make corrections. Here, Crn in equation (23) can be obtained by equation (6).
  • the conductivity between the time of creating the conversion formula used in the dilution water correction unit 68 and the dilution water of the coolant liquid in the coolant tank Even if there is a difference in , the accuracy of the BRIX value calculated from the conductivity can be maintained at a high level.
  • Embodiment 4 a measurement value temperature coefficient correction unit that cancels the temperature dependence of the measurement value will be described.
  • This measured value temperature coefficient correction unit corrects the temperature coefficient of the measured value of the solution, which is a characteristic value of the solution measured by the electrical conductivity method and has a temperature dependency in which the value changes depending on the temperature. It can be applied to various measurements, not limited to conductivity measurements obtained from 24. Measuring instruments that output such measured values include, for example, a conductivity meter, a salinity meter, and a TDS (total dissolved solids) meter.
  • the conductivity measurement values handled in Embodiments 1 to 3 are used as examples of measurement values to be corrected by the measurement value temperature coefficient correction unit.
  • FIG. 11 shows a block diagram of a coolant liquid management device 4 according to the fourth embodiment.
  • the coolant management device 4 according to the fourth embodiment has a temperature sensor TS added to the coolant management device 1 according to the first embodiment. This temperature sensor TS measures the liquid temperature of the coolant liquid in the coolant tank. Further, the coolant liquid management device 4 according to the fourth embodiment replaces the correlation BRIX calculation unit 22 of the coolant liquid management device 1 according to the first embodiment with a correlation BRIX calculation unit 22c.
  • FIG. 12 shows a block diagram of the correlation BRIX calculator 22c according to the fourth embodiment.
  • the correlation BRIX calculation unit 22c is obtained by adding a measured value temperature coefficient correction unit 71 to the correlation BRIX calculation unit 22 according to the first embodiment.
  • the measured value temperature coefficient correction unit 71 corrects the conductivity measurement value Cd according to the temperature of the coolant liquid, and outputs the temperature-corrected conductivity measurement value Cd to the conductivity conversion unit 46. do. That is, the measured value temperature coefficient corrector 71 uses the measured value temperature coefficient corrector 71 to perform temperature correction processing for removing the temperature fluctuation component included in the conductivity measured value Cd.
  • FIG. 13 shows a block diagram of the measured value temperature coefficient correction unit 71 according to the fourth embodiment. Details of the measured value temperature coefficient correction unit 71 will be described with reference to FIG. As shown in FIG. 13, the measured value temperature coefficient correction unit 71 includes a measured value memory 81, a temperature coefficient calculation parameter storage unit 82, an internal clock 83, a temperature coefficient calculation unit 84, a matrix memory 85, a matrix memory 86, a conductivity temperature It has a corrector 87 .
  • the code of the conductivity measurement value before temperature correction is C
  • the code of the conductivity measurement value after temperature correction is Cd.
  • the measured value memory 81 associates the measured temperature value T and the measured conductivity value C acquired from the temperature sensor that measures the temperature of the coolant and stores them in time series. In addition to the measured temperature value T and the measured conductivity value C, the measured value memory 81 stores data related to the amount of change in conductivity. Therefore, data stored in the measured value memory 81 will be described.
  • FIG. 14 shows a diagram for explaining the format of data stored in the measured value memory of the measured value temperature coefficient correction unit according to the fourth embodiment.
  • the measurement time, the temperature measurement value T, and the conductivity measurement value C before temperature correction are associated and stored as one entry in the measurement value memory 81 .
  • a difference in conductivity measurement value between the previous entry and the current entry is stored as a conductivity difference.
  • the measured value memory 81 stores 1 when the conductivity difference is within the range specified by the variation width A specified by the temperature coefficient calculator 84, and when the variation rate difference is outside the range.
  • a fluctuation width flag F of 0 is stored.
  • the fluctuation rate flag shows only the fluctuation width A of ⁇ 1.0%, but the fluctuation width A may be set to other values.
  • the measured value memory 81 stores the value of the stable period counter W in which the count value counted up while the fluctuation width flag F is continuously set to 1 is recorded. This count value is calculated by the temperature coefficient calculator 84 . Note that the example shown in FIG. 14 shows an example in which each entry is saved at intervals of 10 minutes.
  • the temperature coefficient calculation parameter storage unit 82 stores the temperature coefficient calculation parameters input by the operator from the outside.
  • the temperature coefficient calculation parameters include the solution reference temperature Tref, the initial temperature coefficient ⁇ 0, the temperature coefficient range ⁇ min/ ⁇ max, and the calculation period Dmin/Dmax.
  • the built-in clock 83 outputs the measurement time of the measured value memory 81 and the current time information Dn used in the period calculation in the temperature coefficient calculator 84 .
  • the temperature coefficient calculation unit 84 calculates the temperature coefficient ⁇ derived from statistical processing performed on the conductivity measurement value C and the temperature measurement value T stored in the measurement value memory.
  • a matrix memory 85 and a matrix memory 86 are used in this statistical processing. Therefore, this statistical processing will be described in detail.
  • the temperature coefficient calculator 84 extracts the data stored in the measured value memory 81 while switching the combination of the variation width A and the stable period count value W of the stable period Tc.
  • the variation width A indicates the rate of change between the conductivity measurement value C of the previous entry and the conductivity measurement value C of the current entry, and is, for example, 0.1%, 0.2%, 0.1%, 0.2%, or 0.2%.
  • a size of about 3%, 0.5%, 0.7%, 1%, 2%, and 3% can be set.
  • the stable period count value W can be set to values of 2, 4, 8, 12, 24, 48, 72, and 96 (hours), for example.
  • the temperature coefficient calculation unit 84 refers to the measured value memory 81, and the variation width A is 0.2% or less.
  • the conductivity measurement value C and the temperature measurement value T in the part where the count value W continues to count up for 12 hours or more are extracted as data of the stabilization period Tc.
  • the temperature coefficient calculator 84 performs statistical processing on the extracted data and writes the result of the statistical calculation into the matrix memory 85 . Note that the temperature coefficient calculator 84 generates a data table for each combination of the variation width A and the stable period count value W.
  • FIG. 15 shows a diagram for explaining the format of data stored in the matrix memory 85 of the measured value temperature coefficient corrector 71 according to the fourth embodiment.
  • the matrix memory 85 stores data in a matrix.
  • the matrix memory 85 stores data in the row direction for each stable period Tc.
  • the results of statistical processing for each stable period Tc are stored in the column direction.
  • the number of the stable period Tc is stored in the CL1 column.
  • the second column CL stores the entry number of the measured value memory 81 at which the stable period Tc starts, and the third column stores the entry number of the measured value memory 81 at which the stable period Tc ends.
  • the length w of the stable period Tc is stored in the fourth column CL.
  • the fifth column CL stores the start time of each stable period Tc.
  • the sixth row of CL stores the average temperature Xav of each stable period Tc.
  • the seventh column of CL stores the average conductivity Yav of each stable period Tc.
  • the eighth column of CL stores the temperature deviation Sx for each stable period Tc.
  • the ninth column of CL stores the conductivity deviation Sy of each stable period Tc.
  • the 10th column CL stores the covariance Sxy of each stable period Tc.
  • the CL11 column stores the correlation coefficient R2 of each stable period Tc.
  • the 12th column of CL stores the primary coefficient a.
  • the 0th-order coefficient b is stored in the 13th column of CL.
  • the zero temperature t0 is stored in the 14th column CL.
  • the temperature coefficient candidate value ⁇ tc is stored in the 15th column CL.
  • the temperature coefficient limit Z is stored in the 16th column of CL.
  • the temperature coefficient limit Z is a flag that becomes 1 if the temperature coefficient candidate value ⁇ tc is within the temperature coefficient range ⁇ min/ ⁇ max, and becomes 0 if the temperature coefficient candidate value ⁇ tc is outside the temperature coefficient range ⁇ min/ ⁇ max. be. Entries with a temperature coefficient limit Z of 0 are excluded from processing using a matrix memory 86, which will be described later.
  • the 17th column CL stores the correlation coefficient average R 2 av.
  • the 18th column of CL stores the covariance average Sxyav/A.
  • a weighting factor K is stored in the 19th column CL.
  • the CL 20th column stores the temperature coefficient weighted average ⁇ av.
  • the correlation coefficient average R 2 av, the covariance average Sxyav/A, and the temperature coefficient weighted average ⁇ av stored in the lower rows are the temperature coefficient limit values of the entries in the upper rows.
  • the temperature coefficient calculator 84 generates the table data stored in the matrix memory 85 with respect to the variation range A and the stable period count value W, which are subject to statistical processing, and then aggregates the statistical data in the matrix memory 86. After that, the data stored in the matrix memory 85 is cleared and used for the next combination of the variation width A and the stable period count value W.
  • the temperature coefficient calculator 84 calculates the correlation coefficient average R 2 av and the covariance average Sxyav/ A, the temperature coefficient weighted average ⁇ av, the number of valid entries (the number of entries whose flag value of temperature coefficient limit Z is 1, the number of valid Tc), and the total number of entries (total number of Tc) are output.
  • FIG. 16 shows a diagram for explaining the format of data stored in the matrix memory 86 of the measured value temperature coefficient corrector 71 according to the fourth embodiment.
  • the temperature coefficient calculation unit 84 stores the aggregated correlation coefficient average R 2 av, covariance average Sxyav /A, the temperature coefficient weighted average ⁇ av, the number of valid entries (the number of entries whose flag value of the temperature coefficient limit Z is 1, the number of valid Tc), and the total number of entries (total number of Tc) are stored.
  • the temperature coefficient calculation unit 84 performs statistical processing and tabulation processing for all combinations of the candidates for the variation width A and the candidates for the stable period count value W, and fills the entries in the matrix memory 86 .
  • the temperature coefficient calculator 84 searches for an entry that satisfies a preset effective coefficient condition.
  • Valid coefficient conditions are, for example, entries that satisfy the following five conditions. 1st condition: Effective Tc number is 2 or more 2nd condition: Correlation average is 0.5 or more 3rd condition: Covariance average is 0.5 or more 4th condition: Entry with maximum covariance average 5th condition: 3rd condition If there are multiple entries that satisfy the first to fourth conditions, the entry with the largest correlation average is selected.
  • the temperature coefficient calculator 84 verifies the above five conditions for each entry, determines the temperature coefficient weighted average ⁇ av included in the entries satisfying all the conditions as the temperature coefficient ⁇ , and outputs it to the conductivity temperature corrector 87 . On the other hand, if there is no entry that satisfies the above five conditions, the temperature coefficient calculation unit 84 outputs the initial temperature coefficient ⁇ 0 stored in the temperature coefficient calculation parameter storage unit 82 as the temperature coefficient ⁇ . Further, the temperature coefficient calculator 84 reads out the solution reference temperature Tref from the temperature coefficient calculation parameter storage unit 82 together with the temperature coefficient ⁇ , and gives it to the conductivity temperature corrector 87 .
  • the conductivity temperature correction unit 87 corrects the conductivity measurement value C using the temperature coefficient ⁇ and the solution reference temperature Tref, and outputs the conductivity measurement value Cd.
  • the conductivity temperature correction unit 87 performs temperature correction based on, for example, formula (24).
  • the measurement value temperature coefficient correction unit 71 performs correction to cancel the temperature fluctuation component of the conductivity measurement value Cd.
  • the correlation BRIX value Bd output by the correlation BRIX calculator 22c according to the fourth embodiment does not include an error due to temperature.
  • the accuracy of the correlation BRIX value Bd can be improved by using the correlation BRIX calculator 22c according to the fourth embodiment.
  • a time-dependent change coefficient calculation unit that calculates a time-dependent change coefficient for correcting the amount of change in the conductivity from a reference date to the present using a time-dependent change function created in advance from the measurement result of the electrical conductivity of the coolant liquid over time;
  • a conductivity measurement value is obtained from a conductivity meter that measures the conductivity of the coolant liquid, the time-dependent change coefficient is applied to the conductivity measurement value, and the change in electrolytic component concentration over time is calculated from the conductivity measurement value.
  • a conductivity converter that calculates a corrected conductivity with the effect removed; a correlated BRIX converter that applies the corrected conductivity to a predetermined conversion formula to calculate a correlated BRIX value corresponding to the corrected conductivity;
  • BRIX measurement device with
  • a dilution water conductivity that is the conductivity of the diluent mixed with the coolant, a measurement date after the reference date, and the conductivity of the coolant measured on the measurement date further comprising a parameter storage unit that stores a reference conductivity and a second reference BRIX value that is the BRIX of the coolant measured on the measurement date;
  • the time-varying function is a function that linearly approximates the change in the conductivity measurement value and whose slope changes in a predetermined period
  • the temporal change coefficient calculation unit applies the second reference conductivity, the second reference BRIX value, the dilution water conductivity, the reference date, and the measurement date to the temporal change function, BRIX measurement device according to appendix 1, which calculates a coefficient of variation.
  • (Appendix 5) Dilution water conductivity indicating conductivity of dilution water obtained by diluting the coolant when the conversion formula of the correlated BRIX conversion unit is determined, and conductivity of dilution water obtained by diluting the coolant injected into the coolant tank 5.
  • the BRIX measurement device according to any one of Appendices 1 to 4, further comprising a dilution water correction unit that corrects the difference in the .
  • the measured value temperature coefficient correction unit a measurement value memory that associates and stores temperature measurement values obtained from a temperature sensor that measures the temperature of the coolant and the conductivity measurement values in chronological order; a temperature coefficient calculation unit that calculates a temperature coefficient derived from statistical processing performed on the conductivity measurement value and the temperature measurement value stored in the measurement value memory; and a conductivity temperature correction unit that corrects the conductivity measurement value using the temperature coefficient and outputs the temperature-corrected conductivity measurement value to the conductivity conversion unit.
  • BRIX measuring device according to item 1.
  • a solution reference temperature indicating a reference value of the temperature of the coolant liquid, an initial temperature coefficient indicating an initial value of the temperature coefficient, a temperature coefficient range indicating an upper limit and a lower limit range of the temperature coefficient, and a temperature measurement to be calculated.
  • the temperature coefficient calculation unit With reference to the measured value memory, the length of the stable period in which the difference between the two conductivity measured values that are consecutive in time series is equal to or less than a preset fluctuation width is calculated for each different fluctuation width, calculating a temperature coefficient weighted average using the average value, deviation, covariance, and correlation coefficient of the conductivity measurement value and the temperature measurement value included in the stable period for each stable period; Among the combinations of the length of the stable period and the magnitude of the variation range in which both the average value of the covariance and the average value of the correlation coefficient are equal to or greater than a predetermined value, the average value of the covariance is calculated for the combination with the largest value.
  • the BRIX measuring device wherein the temperature coefficient weighted average is calculated as the temperature coefficient.
  • the temperature coefficient calculation unit The temperature coefficient candidate value used to calculate the temperature coefficient weighted average invalidates the data of the stable period in which the temperature coefficient candidate value is outside the temperature coefficient range, and the temperature coefficient candidate value is within the temperature coefficient range. validating the data of the stable period including the temperature coefficient candidate value; 7.
  • the BRIX measuring apparatus according to appendix 7, wherein the temperature coefficient weighted average calculated from the valid data of the stable period is used as the temperature coefficient calculation target.
  • the temperature coefficient calculation unit Output the initial temperature coefficient as the temperature coefficient when there is no combination of the length of the stable period and the magnitude of the variation range in which both the average value of the covariance and the average value of the correlation coefficient are equal to or greater than a predetermined value.
  • the BRIX measurement device according to appendix 7 or 8.
  • a BRIX measurement program for calculating a BRIX value of the coolant liquid from a conductivity measurement value obtained from a conductivity meter that measures the conductivity of the coolant liquid by being executed by a calculation unit, A time-dependent change coefficient calculation process for calculating a time-dependent change coefficient for correcting the amount of change in the electrical conductivity from a reference date to the present using a time-dependent change function created in advance from the measurement results of the electrical conductivity of the coolant liquid over time; A conductivity conversion process of applying the coefficient of change over time to the measured conductivity value to calculate a corrected conductivity obtained by removing the influence of changes in electrolytic component concentration over time from the measured conductivity value by correcting; A correlated BRIX conversion process of applying the corrected conductivity to a predetermined conversion formula to calculate a correlated BRIX value corresponding to the corrected conductivity; BRIX measurement program that performs
  • Reference Signs List 1 4 coolant management device 10 control box 11 power supply unit 12 calculation unit 13 communication interface 14 external communication interface 15 control signal generation unit 20 sensor box 21 power supply unit 22 correlation BRIX calculation unit 23 control signal generation unit 24 conductivity meter 25 analog Digital conversion circuit 26 communication interface 30 sine wave generator 31 amplifier 32 primary side toroidal transformer 33 secondary side toroidal transformer 34 low noise amplifier 35 signal amplitude detector 41 data storage section 42, 52, 62 parameter storage section 43 built-in clock 44, 54 Elapsed time counters 45, 55, 65 Temporal change coefficient calculators 46, 56 Conductivity converters 47, 57, 67 Correlation BRIX converters 68 Dilution water corrector 71 Measured value temperature coefficient corrector 81 Measured value memory 82 Temperature coefficient calculation parameter Storage unit 83 Built-in clock 84 Temperature coefficient calculation unit 85 Matrix memory 86 Matrix memory 87 Conductivity temperature correction unit TS Temperature sensor IND Display unit

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Mechanical Engineering (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

従来のクーラント液管理装置では、溶液の測定値の液温に起因する測定値のズレが生じる問題がある。本発明にかかる測定値温度係数補正装置の一態様は、溶液の温度を計測する温度センサから取得した温度測定値（Ｔ）と、電気伝導方式によって計測される溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値（Ｃ）と、を関連付けて時系列に記憶する測定値メモリ（８１）と、測定値メモリ（８１）に記憶された溶液測定値（Ｃ）と温度測定値（Ｔ）とについて統計処理から導き出した温度係数（α）を算出する温度係数算出部（８４）と、温度係数（α）を用いて溶液測定値（Ｃ）を補正した温度補正済みの溶液測定値（Ｃｄ）を出力する測定値温度補正部（８７）と、を有する。

Description

測定値温度係数補正装置及び測定値温度係数補正プログラム
　本発明は測定値温度係数補正装置及び測定値温度係数補正プログラムに関し、例えば、溶液の特性値のうち温度依存性を有する測定値に対して液温に応じた補正を行う測定値温度係数補正装置及び測定値温度係数補正プログラムに関する。
　工場では、切削工程や研磨工程などの機械加工工程において冷却剤又は潤滑油としてクーラント液（切削油）が用いられる。このクーラント液には、水溶性のものと不水溶性のものがあり、水溶性クーラントは５倍～２０倍程度に希釈して使用する。希釈したクーラントは、水分中にクーラント原液が乳化した状態で混ざり合う状態となる。また、適切な濃度に希釈したクーラントはｐＨが８～９程度のアルカリ性になるよう設定されており、錆びやバクテリアの発生を防ぐ効果を発揮する。
　工場では通常、クーラント液はクーラント槽に貯蔵され、使用する際に供給ポンプによって加工室に配管を通して送られ、使用後は回収用配管を通してクーラント槽へ戻されるという循環システムが構築されている。これにより、工場では、クーラント液を長期間にわたって繰り返し利用することが可能となっている。しかしながら、クーラントは長期間使用しているうちに、フィルタでは除去できない切削屑や摺動油が混入し汚れが増えてくる。また、水溶性クーラント液は、加工時の発熱により水分が蒸発し濃度が濃くなったり、逆に切削屑などにクーラント成分が付着し除去フィルタで一緒に除去され濃度が薄くなったりする。このようなクーラント液の濃度変化を放置すると、本来の加工性能(冷却性、潤滑性)が得られなかったり、分離やｐＨ低下による錆やバクテリアの発生に起因したクーラント液の変色、濁り、腐敗、ヘドロの発生が生じる。このようなクーラント液の変化は、クーラント液の劣化と言われ、加工性能低下だけでなく工場内環境面の悪化も招いてしまう。
　このようなことから、加工性能を維持してクーラント液を利用するために、工場では、日常的にクーラント液の濃度等を管理している。そこで、クーラント液を管理する方法として屈折率計（例えば、ＢＲＩＸ計）を用いたクーラント液の濃度管理がある。このＢＲＩＸ計を用いた濃度管理は、週１回程度クーラント液を汲み上げ、濃度測定を行い、その結果に基づき規程の濃度となるようにクーラント液の濃度を調整するものである。また、濁りや臭いなどの他の項目については、視覚や嗅覚により確認することが行われている。そのため、クーラント液の濃度管理には多くの手間と経験を要する問題がある。そこで、クーラント液の濃度を導電量計を用いて計測する技術が特許文献１に開示されている。
　特許文献１に記載のクーラント液管理装置は、クーラント液に浸した１次側トロイダルトランスに電界を与えた際にクーラント液に生じる一巡電流ループの作用により二次側トロイダルトランスに生じる電界から得られる電圧差に基づきクーラント液の導電量を計測する電磁誘導型の導電量計と、導電量を予め与えられた変換式に適用して、導電量に対応するＢＲＩＸ値を算出する演算部と、ＢＲＩＸ値を外部に出力する外部通信インタフェースと、を有する。
特開２０２１－４１４８７号公報
　しかしながら、溶液の１つであるクーラント液は、導電率計を用いて濃度変化を検出しようとした場合、液温によって特性値にズレが生じる。許文献１に記載の技術では、このような溶液から測定される導電率の温度変化について考慮されておらず、計測される導電率に液温に起因したにズレが生じる問題がある。
　本発明にかかる測定値温度係数補正装置の一態様は、溶液の温度を計測する温度センサから取得した温度測定値と、電気伝導方式によって計測される前記溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値と、を関連付けて時系列に記憶する測定値メモリと、前記測定値メモリに記憶された前記溶液測定値と前記温度測定値とについて統計処理から導き出した温度係数を算出する温度係数算出部と、前記温度係数を用いて前記溶液測定値を補正した温度補正済みの前記溶液測定値を出力する測定値温度補正部と、を有する。
　本発明にかかる測定値温度係数補正プログラムの一態様は、演算部で実行されることで、温度係数を用いて測定値を補正する測定値温度係数補正プログラムであって、溶液の温度を計測する温度センサから取得した温度測定値と、電気伝導方式によって計測される前記溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値と、を関連付けて時系列に測定値メモリに記憶する測定値記憶処理と、前記測定値メモリに記憶された前記溶液測定値と前記温度測定値とについて統計処理から導き出した温度係数を算出する温度係数算出処理と、前記温度係数を用いて前記溶液測定値を補正した温度補正済みの前記溶液測定値を出力する測定値温度補正処理と、を前記演算部に実行させる。
　本発明にかかる測定値温度係数補正装置及び測定値温度係数補正プログラムは、溶液測定値と温度測定値とに対して統計処理を施して温度係数を算出することで温度による溶液測定値のズレを補正する。
　本発明にかかる測定値温度係数補正装置及び測定値温度係数補正プログラムは、計測される導電率に液温に起因したにズレを補正することができる。
実施の形態１にかかるクーラント液管理装置の外観を説明する図である。 実施の形態１にかかるクーラント液管理装置のブロック図である。 実施の形態１にかかるクーラント液管理装置の導電量計のブロック図である。 実施の形態１にかかる相関ＢＲＩＸ算出部のブロック図である。 実施の形態１にかかる相関ＢＲＩＸ算出部で利用される変換式を用いた導電率とＢＲＩＸとの関係を説明するグラフである。 実施の形態１にかかる相関ＢＲＩＸ算出部で利用される経時変化関数を説明するグラフである。 実施の形態２にかかる相関ＢＲＩＸ算出部のブロック図である。 実施の形態２にかかる相関ＢＲＩＸ算出部で利用される変換式を用いた導電率とＢＲＩＸとの関係を説明するグラフである。 実施の形態２にかかる相関ＢＲＩＸ算出部で利用される経時変化関数を説明するグラフである。 実施の形態３にかかる相関ＢＲＩＸ算出部のブロック図である。 実施の形態４にかかるクーラント液管理装置のブロック図である。 実施の形態４にかかる相関ＢＲＩＸ算出部のブロック図である。 実施の形態４にかかる測定値温度係数補正部のブロック図である。 実施の形態４にかかる測定値温度係数補正部の測定値メモリに蓄積されるデータのフォーマットを説明する図である。 実施の形態４にかかる測定値温度係数補正部の行列メモリに蓄積されるデータのフォーマットを説明する図である。 実施の形態４にかかる測定値温度係数補正部のマトリクスメモリに蓄積されるデータのフォーマットを説明する図である。
　説明の明確化のため、以下の記載及び図面は、適宜、省略、及び簡略化がなされている。また、様々な処理を行う機能ブロックとして図面に記載される各要素は、ハードウェア的には、ＣＰＵ（Central Processing Unit）、メモリ、その他の回路で構成することができ、ソフトウェア的には、メモリにロードされたプログラムなどによって実現される。したがって、これらの機能ブロックがハードウェアのみ、ソフトウェアのみ、またはそれらの組合せによっていろいろな形で実現できることは当業者には理解されるところであり、いずれかに限定されるものではない。なお、各図面において、同一の要素には同一の符号が付されており、必要に応じて重複説明は省略されている。
　上述したプログラムは、コンピュータに読み込まれた場合に、実施形態で説明された１又はそれ以上の機能をコンピュータに行わせるための命令群（又はソフトウェアコード）を含む。プログラムは、非一時的なコンピュータ可読媒体又は実体のある記憶媒体に格納されてもよい。限定ではなく例として、コンピュータ可読媒体又は実体のある記憶媒体は、random-access memory（RAM）、read-only memory（ROM）、フラッシュメモリ、solid-state drive（SSD）又はその他のメモリ技術、CD-ROM、digital versatile disc（DVD）、Blu-ray（登録商標）ディスク又はその他の光ディスクストレージ、磁気カセット、磁気テープ、磁気ディスクストレージ又はその他の磁気ストレージデバイスを含む。プログラムは、一時的なコンピュータ可読媒体又は通信媒体上で送信されてもよい。限定ではなく例として、一時的なコンピュータ可読媒体又は通信媒体は、電気的、光学的、音響的、またはその他の形式の伝搬信号を含む。
　実施の形態１
　図１に実施の形態１にかかるクーラント液管理装置１の概要を説明する概略図を示す。図１に示すように、実施の形態１にかかるクーラント液管理装置１は、制御ボックス１０とセンサボックス２０の２つのボックスが吊り下げワイヤＷｉ及び連結治具により連結される。また制御ボックス１０には表示部ＩＮＤが設けられる。また、クーラント液管理装置１は、滑車にかけられた吊り下げワイヤＷｉにより重りＷｔと連結される。そして、吊り下げワイヤＷｉを上下させることでクーラント液管理装置１は、クーラント槽内を上下する。
　センサボックス２０は、クーラント液に浸された状態でクーラント液の導電率等を測定するための各種センサが搭載される。なお、実施の形態１では、センサボックス２０には、導電率計が搭載される。なお、以下の説明では、導電率計を導電量計と称す。制御ボックス１０は、センサボックス２０に搭載されている各種センサを制御すると共に、外部と通信を行う通信インタフェース等が搭載される。また、制御ボックス１０は、センサボックス２０がクーラント液に浸された状態であても、クーラント液に浸からない位置に設けられる。
　続いて、実施の形態１にかかるクーラント液管理装置１の内部構成について説明する。そこで、図２に実施の形態１にかかるクーラント液管理装置１のブロック図を示す。図２に示すように、制御ボックス１０は、電源部１１、演算部１２、通信インタフェース１３、外部通信インタフェース１４、表示部ＩＮＤを有する。また、センサボックス２０は、電源部２１、相関ＢＲＩＸ算出部２２、制御信号生成部２３、導電量計２４、アナログデジタル変換回路２５、通信インタフェース２６を有する。
　電源部１１は、外部に設けられるＡＣアダプタから供給される電源を制御ボックス１０の各処理ブロックに分配する（不図示）。また、電源部１１は、センサボックス２０の電源部２１も電源を供給する。電源部２１は、センサボックス２０の各処理ブロックに供給された電源を分配する（不図示）。
　演算部１２は、例えば、マイクロコントローラのようなプログラムを実行可能な演算部である。演算部１２は、例えば、クーラント液管理装置１において実行されるクーラント液のＢＲＩＸ計測処理の計測シーケンスを制御する。また、演算部１２は、通信インタフェース１３、外部通信インタフェース１４を介した通信処理に関する制御、表示部ＩＮＤへの情報表示処理を行う。通信インタフェース１３は、センサボックス２０側の通信インタフェース２６を介して演算部１２と相関ＢＲＩＸ算出部２２が通信を行うための通信の具体的処理を行うインタフェース回路である。外部通信インタフェース１４は、上位システム（構内サーバー、クラウドサーバ等）と有線、或いは、無線により通信を行うためのインタフェース回路である。表示部ＩＮＤは、クーラント液管理装置１の動作状況、計測情報及びクーラント液の状態を表示するためのものである。
　相関ＢＲＩＸ算出部２２は、例えば、マイクロコントローラのようなプログラムを実行可能な演算部である。相関ＢＲＩＸ算出部２２は、例えば、導電量計２４から得られた導電率の値に基づきＢＲＩＸを算出するＢＲＩＸ計算処理、導電量計２４から得られた信号を具体的な数値に変換する導電率計算処理等を行う。このＢＲＩＸ計算処理の詳細は後述する。また、相関ＢＲＩＸ算出部２２は、制御信号生成部２３に導電量計２４を動作させるか否かの指示を与える。
　制御信号生成部２３は、相関ＢＲＩＸ算出部２２からの指示に基づき導電量計２４に動作の具体的な指示を与える。導電量計２４は、クーラント液に浸された状態でクーラント液の導電率を測定する。アナログデジタル変換回路２５は導電量計２４から出力されたアナログ値を有する出力信号をデジタル値を有する信号に変換する。通信インタフェース２６は、相関ＢＲＩＸ算出部２２が演算部１２と通信を行うための具体的な処理を行うインタフェース回路である。
　続いて、実施の形態１の導電量計２４について詳細に説明する。そこで、図３に実施の形態１にかかるクーラント液管理装置の導電量計のブロック図を示す。図３に示すように導電量計２４は、正弦波生成部３０、アンプ３１、一次側トロイダルトランス３２、二次側トロイダルトランス３３、ローノイズアンプ３４、信号振幅検出器３５を有する。
　正弦波生成部３０は、制御信号生成部２３からの指示を受けて、一次側トロイダルトランス３２を動作させるための正弦波信号を出力する。アンプ３１は、アンプ３１から与えられた正弦波信号を互いに位相が１８０度異なる差動信号に変換して一次側トロイダルトランス３２を駆動する。一次側トロイダルトランス３２は、アンプ３１から与えられた差動信号に基づき電界を発生させる。そして、導電量計２４では、一次側トロイダルトランス３２において発生した電界に基づきクーラント液に一次側トロイダルトランス３２及び二次側トロイダルトランス３３の内周を貫通して一巡する電流ループが発生する。
　そして、二次側トロイダルトランス３３は、一次側トロイダルトランス３２が発生させた電流ループに起因して発生した電界に応じた出力信号をローノイズアンプ３４に出力する。ローノイズアンプ３４は、二次側トロイダルトランス３３の出力信号を増幅して、信号振幅検出器３５に与える。信号振幅検出器３５は、ローノイズアンプ３４から与えられた信号の振幅を検出して、当該検出信号を測定結果としてアナログデジタル変換回路２５に出力する。
　実施の形態１にかかるクーラント液管理装置１では、導電量計２４によって取得された導電率を用いてＢＲＩＸ値を計算により算出する。そこで、実施の形態１におけるＢＲＩＸの計算方法について詳細に説明する。
　まず、ＢＲＩＸについて説明する。ＢＲＩＸは、主に食品産業において、ショ糖、果糖、転化糖、ブドウ糖など、いわゆる糖の含有量を測るために、糖度として用いられる物理量である。ＢＲＩＸ値は、２０℃のショ糖溶液の質量百分率に相当する値で定められている。ショ糖１ｇのみを溶質として含む水溶液１００ｇをＢＲＩＸ屈折計で測定したときその示度ＢＲＩＸ値が１％となる。
　一方、導電率は、上記屈折率計で測定した値とは全く異なるものであるため、実施の形態１にかかるクーラント液管理装置１では、相関ＢＲＩＸ算出部２２を用いて、導電率からＢＲＩＸに相当する相関ＢＲＩＸ値を算出する。そこで、相関ＢＲＩＸ算出部２２について以下で詳細に説明する。なお、以下の説明では、相関ＢＲＩＸ算出部２２の機能のうち相関ＢＲＩＸ値に関する部分を特に詳細に説明し、導電量計２４の制御に関する部分については説明を省略する。
　実施の形態１にかかる相関ＢＲＩＸ算出部２２では、更液（クーラント液の全交換）した直後の未使用状態のクーラント液の導電率とＢＲＩＸ値とを測定して導電率とＢＲＩＸの関係を示す回帰関数を作成して、この回帰関数から相関ＢＲＩＸ値Ｂｄの算出に用いる変換式（例えば、ＢＲＩＸ近似関数）に適用されるＢＲＩＸ近似係数を算出する。そして、算出したＢＲＩＸ近似係数を適用したＢＲＩＸ近似関数を用いて導電率から相関ＢＲＩＸ値Ｂｄを算出する。
　その後、相関ＢＲＩＸ算出部２２はクーラント液の相関ＢＲＩＸ値Ｂｄを出力する計測処理を開始する前準備として計測で利用するパラメータ情報を収集する。このパラメータ情報には、クーラント槽にクーラント液を入れた更液日時Ｄ０、測定開始前などの特定の基準日時に導電量計２４で計測したクーラント槽内のクーラント液の導電率である基準導電率Ｃｒ、同時に別途準備したＢＲＩＸ計で測定したクーラント槽内のクーラント液のＢＲＩＸである基準ＢＲＩＸ値Ｂｒ、基準導電率Ｃｒ及び基準ＢＲＩＸ値Ｂｒを測定した基準測定日時Ｄｒ１、クーラント液の希釈に用いる希釈水の導電率である希釈水導電率Ｃｗが含まれる。そして、相関ＢＲＩＸ算出部２２は、測定開始前に、これらパラメータ情報を用いて経時変化係数ｋを算出する。
　そして、測定開始後は、相関ＢＲＩＸ算出部２２は、導電量計２４でクーラント液の導電率測定値Ｃｄを取得し、取得した導電率測定値Ｃｄに対して経時変化係数ｋを適用した補正処理により補正導電率Ｃｄｎを生成し、この補正導電率をＢＲＩＸ近似関数に適用して算出される相関ＢＲＩＸ値Ｂｄをクーラント液のＢＲＩＸ値として出力する。
　上記処理を行う相関ＢＲＩＸ算出部２２について、さらに詳細に説明する。そこで、図４に実施の形態１にかかる相関ＢＲＩＸ算出部２２のブロック図を示す。図４に示すように、相関ＢＲＩＸ算出部２２は、データ格納部４１、パラメータ保存部４２、内蔵時計４３、経過時間カウンタ４４、経時変化係数算出部４５、導電率変換部４６、相関ＢＲＩＸ変換部４７を有する。
　相関ＢＲＩＸ算出部２２は、例えば、マイクロコンピュータのようなプログラムを実行可能な演算部である。そして、データ格納部４１及びパラメータ保存部４２は、演算部内のハードディスク、或いは、ＳＳＤ（Solid-State Drive）のような不揮発性の記憶装置、又は、ＤＲＡＭ等の揮発性メモリを利用可能である。一方、その他の構成ブロックについてはプログラム（例えば、ＢＲＩＸ計測プログラム）により実現可能である。
　データ格納部４１は、クーラント液の情報をライブラリ化して保持する記憶部である。データ格納部４１には、事前に調査したクーラント液のＢＲＩＸ近似係数が、クーラント液のメーカー名及び型番（或いは、品番）と紐付けられた情報が保持される。なお、ＢＲＩＸ近似係数は、事前の検証により、別途準備した導電量計とＢＲＩＸ計とを用いて作成する相関ＢＲＩＸ近似関数に含まれる係数である。実施の形態１にかかる相関ＢＲＩＸ算出部２２では、相関ＢＲＩＸ近似関数として（１）式で示される二次関数を用いる。
Figure JPOXMLDOC01-appb-M000001
 ここで、Ｂｎは測定したクーラント液の相関ＢＲＩＸ値であり、Ｃｎは測定したクーラント液の導電率測定値であり、ａｎ、ｂｎ及びｃｎはＢＲＩＸ近似係数である。ここで、図５に実施の形態１にかかる相関ＢＲＩＸ算出部で利用される変換式を用いた導電率とＢＲＩＸとの関係を説明するグラフを示す。ＢＲＩＸ近似係数は、測定時に計測される複数の計測点から得られるＢＲＩＸ近似関数の曲線に沿うように（１）式が沿うように決定する。なお、図５に示した曲線の始点は、希釈液の相関ＢＲＩＸ値のものであり、相関ＢＲＩＸ関数は、希釈水導電率Ｃｗに対して相関ＢＲＩＸ値Ｂｄとして０を出力する。なお、クーラント液の利用期間が長くなると、希釈水の蒸発等によりクーラント液の導電率を変化させる電解成分の濃度が高くなるためＢＲＩＸ値に変化がなくても相関ＢＲＩＸ関数の傾きが小さくなるような変化が生じる（図５の時間経過後のＢＲＩＸ近似関数）。そこで、相関ＢＲＩＸ算出部２２では経時変化係数ｋを適用して（１）式から得られる導電率と相関ＢＲＩＸ値Ｂｄの関係が一定になるように維持する。
　パラメータ保存部４２は、経過時間カウンタ４４及び導電率変換部４６が利用するパラメータを保存する。このパラメータには、クーラント情報、更液日時Ｄ０、基準測定日時Ｄｒ１、基準導電率Ｃｒ、希釈水導電率Ｃｗ、基準ＢＲＩＸ値Ｂｒ、経時変化定数Ｄ１、Ｄ２が保存される。これらパラメータは、導電量計２４、別途準備するＢＲＩＸ計、データ格納部４１に格納されているクーラント情報、内蔵時計４３から得られる情報を自動的に収集するものであっても、制御ボックス１０側の表示部ＩＮＤ等の入力装置（不図示）から与えられるものであってもよい。なお、更液日時Ｄ０及び経時変化定数Ｄ１、Ｄ２は、表示部ＩＮＤから入力される情報である。
　また、基準測定日時Ｄｒ１は、基準導電率Ｃｒ、基準ＢＲＩＸ値Ｂｒを測定した日時であり、実施の形態１にかかるクーラント液管理装置１では、基準測定日時を基準日とする。基準導電率Ｃｒは、クーラント槽に注液されたクーラント液について基準測定日時Ｄｒ１に導電量計２４で測定された導電率である。基準ＢＲＩＸ値Ｂｒは、クーラント槽に注液されたクーラント液について基準測定日時Ｄｒ１に外部のＢＲＩＸ計により測定されたＢＲＩＸ値である。また、希釈水導電率Ｃｗは、クーラント槽に注液されたクーラント液を希釈した際に注液した希釈水について導電量計２４を用いて計測した導電率である。経時変化定数Ｄ１、Ｄ２は、システム利用時に設定する値であり、例えば、経時変化定数Ｄ１は５０日程度の値であり、経時変化定数Ｄ２は３００日程度の値であり、導電率が２倍程度に変化する期間を指定することが好ましい。
　内蔵時計４３は、日時をカウントして現在時刻情報Ｄｎを出力する。この現在時刻情報Ｄｎの１部は基準測定日時Ｄｒ１としても利用可能である。経過時間カウンタ４４は、経過時間カウンタ４４が出力する基準経過時間Ｄｒ＝Ｄｒ１－Ｄ０で値がセットされ、その後内蔵時計４３からの時刻情報に基づき１秒ごとにカウントアップされる経過時間Ｄを出力する。
　経時変化係数算出部４５は、クーラント液の導電率の経時変化の測定結果から予め作成した経時変化関数を用いて基準日から現在までの前記導電率の変化量を補正する経時変化係数ｋを算出する。経時変化係数算出部４５は、基準導電率Ｃｒと、基準ＢＲＩＸ値Ｂｒと、希釈水導電率Ｃｗと、基準測定日時Ｄｒ１と、を経時変化関数に適用して経時変化係数ｋを算出する。
　ここで、相関ＢＲＩＸ算出部２２では、経時変化関数として導電率測定値の変化を非線形近似した関数を用いる。そこで、この経時変化関数について図６を参照して説明する。図６に実施の形態１にかかる相関ＢＲＩＸ算出部で利用される経時変化関数を説明するグラフを示す。図６に示すように、経時変化関数α（Ｄ）は、更液日時Ｄ０からの経過日数をＤとした関数である。経時変化関数α（Ｄ）は、経過日数Ｄが小さい期間の変化割合が大きく、経過日数Ｄが大きくなると変化の割合が小さくなる。
　この経時変化関数α（Ｄ）についてさらに詳しく説明する。まず、経過時間Ｄが経過した後の導電率測定値Ｃｄは、（２）式で表わされ、この（２）式からα（Ｄ）は（３）式にて表わすことができる。
Figure JPOXMLDOC01-appb-M000002
Figure JPOXMLDOC01-appb-M000003
 （２）式及び（３）式には導電率変換部４６が出力する補正導電率Ｃｄｎが含まれる。この補正導電率Ｃｄｎは、（１）式導出時に測定したクーラント液の導電率に相当する値である。導電率測定値Ｃｄは、経過時間Ｄにおいて測定したクーラント液の導電率である。そして、（３）式は、経時変化係数ｋ、経時変化定数Ｄ１、Ｄ２を用いて（４）式の近似式とすることができる。
Figure JPOXMLDOC01-appb-M000004
　そして、経時変化係数ｋは、（１）式～（３）式及び基準測定日時Ｄｒ１において測定された基準導電率Ｃｒと、基準ＢＲＩＸ値Ｂｒ、及び基準測定日時までの経過日時Ｄｒ（＝Ｄｒ１－Ｄ０）により（５）式のように表わされる。
Figure JPOXMLDOC01-appb-M000005
 （５）式のＣｒｎは、基準初期導電率であり、基準ＢＲＩＸ値Ｂｒを（１）式に適用し、かつ、（１）式の解の公式から算出できる。基準初期導電率Ｃｒｎは（６）式で算出できる。
Figure JPOXMLDOC01-appb-M000006
 このように、経時変化係数算出部４５では、パラメータ保存部４２から読み出した基準導電率Ｃｒ、希釈水導電率Ｃｗ、経時変化定数Ｄ１、Ｄ２及び基準ＢＲＩＸ値Ｂｒを（５）式及び（６）式に適用して経時変化係数ｋを算出する。
　導電率変換部４６は、クーラント液の導電率を計測する導電量計２４から導電率測定値Ｃｄを取得して、導電率測定値Ｃｄに経時変化係数ｋを適用して、導電率測定値Ｃｄから経時的な電解成分濃度変化の影響を取り除いた補正導電率Ｃｄｎを算出する。具体的には、導電率変換部４６は以下のような計算を行うことで補正導電率Ｃｄｎを算出する。
　補正導電率Ｃｄｎは、（２）式を変形し、かつ（５）式を用いて（７）式で表わすことが出来る。
Figure JPOXMLDOC01-appb-M000007
　相関ＢＲＩＸ変換部４７は、補正導電率Ｃｄｎを予め決定した変換式に適用して補正導電率Ｃｄｎに対応する相関ＢＲＩＸ値Ｂｄを出力する。ここで予め決定した変換式は事前に算出したＢＲＩＸ近似係数を（１）式に適用したものであり、（８）式で表わされる。
Figure JPOXMLDOC01-appb-M000008
　上記説明より、実施の形態１にかかる相関ＢＲＩＸ算出部２２では、クーラント液の経時変化から導き出され、かつ、時間とともに値が変化する経時変化係数ｋを用いて導電率測定値Ｃｄから電解成分の濃度変化の影響を除去する。これにより、相関ＢＲＩＸ算出部２２では、事前に準備した相関ＢＲＩＸ近似関数を変更することなく、クーラント液の実際のＢＲＩＸとのズレが少なく精度の高い相関ＢＲＩＸ値Ｂｄを出力することができる。
　また、相関ＢＲＩＸ算出部２２において算出される経時変化係数ｋは、時間の関数であり、計測時間をカウントすることで時間とともに自動的に変化させていくことができる。つまり、実施の形態１にかかる相関ＢＲＩＸ算出部２２を用いることで、パラメータの調整を行うことなくクーラント液のＢＲＩＸの測定を連続して行うことができる。
　実施の形態２
　実施の形態２では、実施の形態１で説明した相関ＢＲＩＸ算出部２２の別の形態となる相関ＢＲＩＸ算出部２２ａについて説明する。なお、実施の形態２の説明において、実施の形態１で説明した構成要素と同じ構成要素については、実施の形態１と同じ符号を付して説明を省略する。
　実施の形態１にかかる相関ＢＲＩＸ算出部２２では、ＢＲＩＸ近似関数の算出において、未使用のクーラント液を用いたが、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、ＢＲＩＸ近似関数を運用中のクーラント液を基準に作成する。また、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、基準日をクーラント液管理装置１の設置日等のクーラント液が既に利用状態にある日、かつ、クーラント液の導電率やＢＲＩＸを測定した日とする。
　図７に実施の形態２にかかる相関ＢＲＩＸ算出部２２ａのブロック図を示す。図７に示すように、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａは、パラメータ保存部５２、内蔵時計４３、経過時間カウンタ５４、経時変化係数算出部５５、導電率変換部５６、相関ＢＲＩＸ変換部５７を有する。
　パラメータ保存部５２は、経時変化係数算出部５５が計算で利用するパラメータを保存する。パラメータ保存部５２に保存されるパラメータには、第１基準測定日時Ｄｒ１、第１基準導電率Ｃｒ１、希釈後導電率Ｃｒ１’、第１基準ＢＲＩＸ値Ｂｒ１、希釈後ＢＲＩＸ値Ｂｒ１’、第２基準測定日時Ｄｒ２、第２基準導電率Ｃｒ２、第２基準ＢＲＩＸ値Ｂｒ２、希釈水導電率Ｃｗが含まれる。
　実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、基準日をクーラント液管理装置１の設置日等のクーラント液が既に利用状態にある日、かつ、クーラント液の導電率やＢＲＩＸを測定した日とする。そして、この測定日の日時を第１基準測定日時Ｄｒ１とする。測定日に測定した導電率測定値Ｃｄを第１基準導電率Ｃｒ１とする。希釈後導電率Ｃｒ１’は、第１基準導電率Ｃｒ１を測定したクーラント液をさらに予め設定された濃度（例えば、１／２の濃度）にまで希釈を進めた希釈クーラント液の導電率である。第１基準ＢＲＩＸ値Ｂｒ１は、基準日に測定したクーラント液のＢＲＩＸ値である。希釈後ＢＲＩＸ値Ｂｒ１’は、希釈後導電率Ｃｒ１’を測定したクーラント液のＢＲＩＸ値である。
　また、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、第１基準測定日時Ｄｒ１より後の日であって、相関ＢＲＩＸ値Ｂｄと実際のＢＲＩＸ値とのズレが大きくなったと判断した日に再度各種測定を行う。この再測定日が基準測定日時Ｄｒ２である。第２基準導電率Ｃｒ２は、基準測定日時Ｄｒ２において測定した導電率測定値Ｃｄである。第２基準ＢＲＩＸ値Ｂｒ２は、基準測定日時Ｄｒ２において測定したＢＲＩＸである。また、希釈水導電率Ｃｗは、クーラント液を希釈する際に利用した希釈液の導電率である。
　経過時間カウンタ５４は、第１基準測定日時Ｄｒ１からの経過した時間の長さを経過時間Ｄとして出力する。経時変化係数算出部５５は、第２基準導電率Ｃｒ２と、第２基準ＢＲＩＸ値Ｂｒ２と、希釈水導電率Ｃｗと、基準測定日時Ｄｒ１と、基準測定日時Ｄｒ２と、を経時変化関数α（Ｄ）に適用して経時変化係数ｋを算出する。ここで、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、経過時間ＤがＤｒ２までの経過時間よりも小さいときには経時変化係数を０、つまり、経時変化関数α（Ｄ）を０とする。また、経時変化係数算出部５５は、相関ＢＲＩＸ変換部５７で利用されるＢＲＩＸ近似係数ａｒ、ｂｒ、ｃｒの算出も行うものとする。
　ここで、経時変化係数算出部５５におけるＢＲＩＸ近似係数ａｒ、ｂｒ、ｃｒ及び経時変化係数ｋの算出方法について説明する。相関ＢＲＩＸ算出部２２ａでは、運用開始時に相関ＢＲＩＸ変換部５７が相関ＢＲＩＸ値Ｂｄを算出するためにＢＲＩＸ近似係数ａｒ、ｂｒ、ｃｒが必要になる。そこで、まずＢＲＩＸ近似係数ａｒ、ｂｒ、ｃｒの算出方法について説明する。
　相関ＢＲＩＸ変換部５７で用いられるＢＲＩＸ近似関数は、（９）式で表わされる。
Figure JPOXMLDOC01-appb-M000009
 ここで、（９）式のＢｒは、第１基準測定日時Ｄｒ１におけるクーラント液のＢＲＩＸであり、Ｃｒは第１基準測定日時Ｄｒ１におけるクーラント液の導電率である。そして、（９）式のＢＲＩＸ近似係数ａｒ、ｂｒ、ｃｒは、第１基準測定日時Ｄｒ１に測定された第１基準導電率Ｃｒ１、希釈後導電率Ｃｒ１’、第１基準ＢＲＩＸ値Ｂｒ１、希釈後ＢＲＩＸ値Ｂｒ１’、希釈水導電率Ｃｗ、希釈水のＢＲＩＸ値（＝０）を用いて（１０）式～（１２）式で求められる。
Figure JPOXMLDOC01-appb-M000010
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
　ここで、（９）式で示したＢＲＩＸ近似関数により得られる導電率とＢＲＩＸ値との関係について説明する。そこで、図８に実施の形態２にかかる相関ＢＲＩＸ算出部２２ａで利用される変換式を用いた導電率とＢＲＩＸとの関係を説明するグラフを示す。図８に示すように、実施の形態２におけるＢＲＩＸ近似関数は、希釈水導電率Ｃｗを始点として、第１基準導電率Ｃｒ１と第１基準ＢＲＩＸ値Ｂｒ１により決まる第１点と、希釈後導電率Ｃｒ１’と希釈後ＢＲＩＸ値Ｂｒ１’により決まる第２点とを結ぶ曲線に沿うように設定される。そして、時間が経過すると、導電率測定値Ｃｄと相関ＢＲＩＸ値Ｂｄとにより決定されるＢＲＩＸ近似関数は、基準測定日時Ｄｒ１において作成したＢＲＩＸ近似関数よりも傾きが緩やかになる。そこで、相関ＢＲＩＸ算出部２２ａでは、経時変化係数ｋを用いて導電率測定値Ｃｄを補正導電率Ｃｄｒに変換することで導電率測定値Ｃｄと相関ＢＲＩＸ値Ｂｄとの関係が基準測定日時Ｄｒ１におけるＢＲＩＸ近似関数の関係になるような補正処理を行う。
　続いて、実施の形態２における経時変化係数ｋの算出方法について説明する。まず、実施の形態２では、経過時間Ｄが第２基準測定日時Ｄｒ２までの経過時間よりも小さい期間はｋ＝０とする。そして、経過時間Ｄが第２基準測定日時Ｄｒ２までの経過時間以上となったときに、０ではない経時変化係数ｋを用いる。以下の説明では、０ではないｋの算出方法について説明する。
　まず、経過時間Ｄの時の導電率測定値Ｃｄは、（１３）式で表わされる。
Figure JPOXMLDOC01-appb-M000013
 また、（１３）式から経時変化係数α（Ｄ）は（１４）式で表わされる。
Figure JPOXMLDOC01-appb-M000014
 （１３）式及び（１４）式には導電率変換部５６が出力する補正導電率Ｃｄｒが含まれる。この補正導電率Ｃｄｒは、（９）式導出時に測定したクーラント液の導電率（例えば第１基準導電率Ｃｒ１）に相当する値である。導電率測定値Ｃｄは、経過時間Ｄにおいて測定したクーラント液の導電率である。そして、（１４）式は、経時変化係数ｋを用いて（１５）式の近似式とすることができる。
Figure JPOXMLDOC01-appb-M000015
　そして、経時変化係数ｋは、（９）式、（１４）式及び（１５）式及び第２基準測定日時Ｄｒ２において測定された第２基準導電率Ｃｒ２、及び第２基準測定日時までの経過時間Ｄｒ=Ｄｒ２-Ｄｒ１として（１６）式のように表わされる。
Figure JPOXMLDOC01-appb-M000016
 （１６）式のＣｒ２ｎは、第２基準日初期導電率であり、第２基準ＢＲＩＸ値Ｂｒ２を（９）式に適用し、かつ、（９）式の解の公式から算出できる。基準初期導電率Ｃｒ２ｎは（１７）式で算出できる。
Figure JPOXMLDOC01-appb-M000017
 このように、経時変化係数算出部５５では、パラメータ保存部５２から読み出したパラメータを（１６）式及び（１７）式に適用して経時変化係数ｋを算出する。
　ここで、実施の形態２にかかる経時変化関数についてグラフを用いて説明する。図９に実施の形態２にかかる相関ＢＲＩＸ算出部２２ａで利用される経時変化関数を説明するグラフを示す。図９に示すように、実施の形態２にかかる経時変化係数α（Ｄ）は、計測を開始する経過時間Ｄが第１基準測定日時Ｄｒ１から第２基準測定日時Ｄｒ２に達するまでの期間に有る場合は０になる。一方、経過時間Ｄが第２基準測定日時Ｄｒ２を越えると傾きが経時変化係数ｋとなる直線的に時間の経過とともに上昇する。
　導電率変換部５６は、クーラント液の導電率を計測する導電量計から導電率測定値Ｃｄを取得して、導電率測定値Ｃｄに経時変化係数ｋを適用して、導電率測定値Ｃｄから経時的な電解成分濃度変化の影響を取り除いた補正導電率Ｃｄｒを算出する。具体的には、導電率変換部５６は以下のような計算を行うことで補正導電率Ｃｄｒを算出する。
　補正導電率Ｃｄｒは、（１３）式を変形し、かつ（１６）式を用いて（１８）式で表わすことが出来る。
Figure JPOXMLDOC01-appb-M000018
 ここで、経時変化係数ｋは、経過時間がＤ＜Ｄｒ２のときは０となる。
　相関ＢＲＩＸ変換部５７は、補正導電率Ｃｄｒを予め決定した変換式に適用して補正導電率Ｃｄｒに対応する相関ＢＲＩＸ値Ｂｄを出力する。ここで予め決定した変換式は事前に算出したＢＲＩＸ近似係数を（９）式に適用したものであり、（１９）式で表わされる。
Figure JPOXMLDOC01-appb-M000019
　上記説明より、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａでは、クーラント液の経時変化から導き出され、かつ、時間とともに値が変化する経時変化係数ｋを用いて導電率測定値Ｃｄから電解成分の濃度変化の影響を除去する。これにより、相関ＢＲＩＸ算出部２２ａでは、設置日等に作成した相関ＢＲＩＸ近似関数を変更することなく、クーラント液の実際のＢＲＩＸ値とのズレが少なく精度の高い相関ＢＲＩＸ値Ｂｄを出力することができる。なお、実施の形態２にかかる経時変化関数α（Ｄ）は、直線近似した経時変化特性から経時変化係数ｋを算出しているため、実施の形態２にかかる相関ＢＲＩＸ算出部２２ａが算出する相関ＢＲＩＸ値Ｂｄは、実施の形態１にかかる相関ＢＲＩＸ算出部２２よりも精度が劣る。
　また、相関ＢＲＩＸ算出部２２ａでは、ＢＲＩＸ近似係数の算出、経時変化係数ｋの算出を利用中のクーラント液から取得できる導電率及びＢＲＩＸ値を用いて算出する。これにより、相関ＢＲＩＸ算出部２２ａを用いたクーラント液管理装置は、実施の形態１のように利用中のクーラント液の未使用時のＢＲＩＸ近似特性を事前に準備する必要がなく、実施の形態１にかかるクーラント液管理装置よりも設置にかかる手間を大幅に簡素化できる。
　実施の形態３
　実施の形態３では、実施の形態１で説明した相関ＢＲＩＸ算出部２２の別の形態となる相関ＢＲＩＸ算出部２２ｂについて説明する。なお、実施の形態３の説明において、実施の形態１で説明した構成要素と同じ構成要素については、実施の形態１と同じ符号を付して説明を省略する。
　図１０に実施の形態３にかかる相関ＢＲＩＸ算出部２２ｂのブロック図を示す。図１０に示すように、実施の形態３にかかる相関ＢＲＩＸ算出部２２ｂは、実施の形態１にかかる相関ＢＲＩＸ算出部２２に希釈水補正部６８を追加したものである。また、相関ＢＲＩＸ算出部２２ｂでは、希釈水補正部６８の追加に伴い、パラメータ保存部４２、経時変化係数算出部４５、相関ＢＲＩＸ変換部４７に代えてパラメータ保存部６２、経時変化係数算出部６５、相関ＢＲＩＸ変換部６７を有する。
　パラメータ保存部６２には、保存されるパラメータとして、クーラント情報、更液日時Ｄ０、基準測定日時Ｄｒ１、基準導電率Ｃｒ、希釈水導電率Ｃｗ、基準ＢＲＩＸ値Ｂｒ、経時変化定数Ｄ１、Ｄ２、指定導電率Ｃ１が含まれる。経時変化係数算出部６５は、経時変化係数ｋの算出に加えて基準導電率Ｃ０の算出を行う。希釈水補正部６８は、相関ＢＲＩＸ変換部の変換式（例えば、（８）式）を決定した際にクーラント液を希釈した希釈水の導電率を示す希釈水導電率を基準導電率Ｃ０と読み替え、実際にクーラント槽に注液したクーラント液を希釈した希釈水の導電率を希釈水導電率Ｃｗとして、その差を補正する。希釈水補正部６８は補正後の導電率として希釈水補正導電率Ｃｄｎｗを出力する。相関ＢＲＩＸ変換部６７は、（８）式で示したＢＲＩＸ近似関数の補正導電率Ｃｄｎとして希釈水補正導電率Ｃｄｎｗを用いて基準ＢＲＩＸ値Ｂｄを算出する。
　ここで、希釈水補正部６８で補正に利用する入力パラメータである基準導電率Ｃ０、指定導電率Ｃ１、希釈水導電率Ｃｗについて説明する。まず、基準導電率Ｃ０は、（８）式を作成する際に計測したクーラント液を希釈した希釈水の導電率である。また、指定導電率Ｃ１は、固定値として表示部ＩＮＤ等の入力部から入力する値である。この指定導電率Ｃ１は、ＢＲＩＸ近似関数の傾きを調整するパラメータであり、希釈水補正部６８による補正後のＢＲＩＸ近似関数を検証して決定する。
　基準導電率Ｃ０は、ＢＲＩＸ近似係数ａｎ、ｂｎ、ｃｎを用いて（８）式の解の公式から導き出されるものであり（２０）式に基づき算出される。この基準導電率Ｃ０は、（８）式で算出される相関ＢＲＩＸ値Ｂｄが０となる導電率である。
Figure JPOXMLDOC01-appb-M000020
　希釈水補正部６８では、基準導電率Ｃ０と希釈水導電率Ｃｗとの差をなくすような補正を行う。具体的には、希釈水補正部６８では、（２１）式に基づき補正処理を行い希釈水補正導電率Ｃｄｎｗを出力する。
Figure JPOXMLDOC01-appb-M000021
 一方、経時変化係数ｋの算出に関しても補正が必要である。（５）式を参照して、
Figure JPOXMLDOC01-appb-M000022
 ここで（２２）式のＣｒｎｗは、（８）式に於ける基準初期導電率Ｃｒｎに相当するものであるが、基準導電率Ｃ０と希釈水導電率Ｃｗとの差をなくすために、以下の補正を行う。
Figure JPOXMLDOC01-appb-M000023
 ここで（２３）式のＣｒｎは、（６）式にて求めることができる。
（２１）式（２２）式、及び（２３）式に基づく補正を行う事で、（８）式の算出時に利用した希釈水と、クーラント槽内のクーラント液を希釈した希釈水との間に導電率の差があっても、ＢＲＩＸ近似関数から得られる曲線は、そのまま利用することができる。
　上記説明より、実施の形態３にかかる相関ＢＲＩＸ算出部２２ｂを用いることで、希釈水補正部６８で利用される変換式作成時と、クーラント槽内のクーラント液の希釈水との間に導電率の差がある場合であっても、導電率から算出したＢＲＩＸ値の精度を高く維持することができる。
　実施の形態４
　実施の形態４では、測定値の温度依存性をキャンセルする測定値温度係数補正部について説明する。この測定値温度係数補正部は、電気伝導方式によって計測される溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値の温度係数を補正するものであり、導電量計２４から得られる導電率測定値に限らず様々な測定値に適用出来る。このような測定値を出力する測定器としては、例えば、導電量計、塩分濃度計、ＴＤＳ（総溶解固形物）計などがある。以下の説明では、実施の形態１～３で扱った導電率測定値を測定値温度係数補正部が補正対象とする測定値の例として説明を行う。
　実施の形態４では、クーラント液の温度に起因する導電率のズレを補正する。そこで、図１１に実施の形態４にかかるクーラント液管理装置４のブロック図を示す。図１１に示すように、実施の形態４にかかるクーラント液管理装置４は、実施の形態１にかかるクーラント液管理装置１に温度センサＴＳを追加している。この温度センサＴＳはクーラント槽内のクーラント液の液温を測定する。また、実施の形態４にかかるクーラント液管理装置４は、実施の形態１にかかるクーラント液管理装置１の相関ＢＲＩＸ算出部２２を相関ＢＲＩＸ算出部２２ｃに置き換えている。
　図１２に実施の形態４にかかる相関ＢＲＩＸ算出部２２ｃのブロック図を示す。図１２に示すように、相関ＢＲＩＸ算出部２２ｃは、実施の形態１にかかる相関ＢＲＩＸ算出部２２に測定値温度係数補正部７１を追加したものである。相関ＢＲＩＸ算出部２２ｃでは、測定値温度係数補正部７１が、導電率測定値Ｃｄをクーラント液の温度に応じて補正して、温度補正済みの導電率測定値Ｃｄを導電率変換部４６に出力する。つまり、測定値温度係数補正部７１は、測定値温度係数補正部７１を用いて導電率測定値Ｃｄに含まれる温度変動成分を除去する温度補正処理を行う。
　図１３に実施の形態４にかかる測定値温度係数補正部７１のブロック図を示す。この図１３を参照して測定値温度係数補正部７１について詳細を説明する。図１３に示すように、測定値温度係数補正部７１は、測定値メモリ８１、温度係数算出パラメータ保存部８２、内蔵時計８３、温度係数算出部８４、行列メモリ８５、マトリクスメモリ８６、導電率温度補正部８７を有する。なお、以下の説明では、温度補正前の導電率測定値の符号をＣ、温度補正後の導電率測定値の符号をＣｄとする。
　測定値メモリ８１は、クーラント液の温度を計測する温度センサから取得した温度測定値Ｔと導電率測定値Ｃとを関連付けて時系列に記憶する。また、測定値メモリ８１には、温度測定値Ｔと導電率測定値Ｃとともに導電率の変動量に関連するデータが保存される。そこで、測定値メモリ８１に保存されるデータについて説明する。
　図１４に実施の形態４にかかる測定値温度係数補正部の測定値メモリに蓄積されるデータのフォーマットを説明する図を示す。図１４に示すように、測定値メモリ８１には測定時刻と温度測定値Ｔと、温度補正前の導電率測定値Ｃとが関連付けて１つのエントリとして保存される。また、各エントリには、１つ前の時刻のエントリと現時刻のエントリとの導電率測定値の差が導電率差分として保存される。また、測定値メモリ８１には、温度係数算出部８４から指定される変動幅Ａにより指定される範囲内に導電率差分が入っている場合に１となり、変動率差分が範囲外である場合に０となる変動幅フラグＦが保存される。図１４では、変動率フラグは変動幅Ａが±１．０％のみを示したが、変動幅Ａは、他の値も設定される。また、測定値メモリ８１には、変動幅フラグＦが連続して１となる間カウントアップされたカウント値が記録される安定期間カウンタＷの値が保存される。このカウント値は、温度係数算出部８４で算出されるものとする。なお、図１４に示す例では、１０分間隔で各エントリが保存される例を示した。
　温度係数算出パラメータ保存部８２は、外部から作業者が入力する温度係数算出パラメータが保存される。この温度係数算出パラメータには、溶液基準温度Ｔｒｅｆ、初期温度係数α０、温度係数範囲αｍｉｎ／αｍａｘ、計算期間Ｄｍｉｎ／Ｄｍａｘが含まれる。
　内蔵時計８３は、測定値メモリ８１の計測時間、温度係数算出部８４における期間計算で利用される現在時刻情報Ｄｎを出力する。
　温度係数算出部８４は、測定値メモリに記憶された導電率測定値Ｃと温度測定値Ｔとに対して行う統計処理から導き出した温度係数αを算出する。この統計処理では、行列メモリ８５とマトリクスメモリ８６を用いる。そこで、この統計処理について詳細に説明する。
　温度係数算出部８４は、変動幅Ａと安定期間Ｔｃの安定期間カウント値Ｗの値の組み合わせを切替えながら測定値メモリ８１に蓄積されているデータを抽出する。ここで、変動幅Ａは、前エントリの導電率測定値Ｃと現エントリの導電率測定値Ｃとの変化の割合を示すものであり、例えば、０．１％、０．２％、０．３％、０．５％、０．７％、１％、２％、３％程度の大きさを設定することができる。また、安定期間カウント値Ｗは、例えば、２、４、８、１２、２４、４８、７２、９６（時間）分の値を設置することができる。
　例えば、変動幅Ａを０．２％、安定期間カウント値Ｗを１２時間とした場合、温度係数算出部８４は、測定値メモリ８１を参照して、変動幅Ａが０．２％以下となり、安定期間カウント値Ｗが１２時間分以上カウント値のカウントアップが連続する部分の導電率測定値Ｃと温度測定値Ｔを安定期間Ｔｃのデータとして抽出する。そして、温度係数算出部８４は、抽出したデータに対して統計処理を施し、統計計算の結果を行列メモリ８５に書き込む。なお、温度係数算出部８４は、変動幅Ａと安定期間カウント値Ｗの組み合わせ毎にデータテーブルを生成する。
　そこで、図１５に実施の形態４にかかる測定値温度係数補正部７１の行列メモリ８５に蓄積されるデータのフォーマットを説明する図を示す。図１５に示すように、行列メモリ８５には行列上にデータが格納される。行列メモリ８５では、行方向に安定期間Ｔｃ毎にデータを格納する。行列メモリ８５では、列方向に各安定期間Ｔｃの統計処理の結果が格納される。
　図１５に示す例では、ＣＬ１列目に安定期間Ｔｃの番号が格納される。ＣＬ２列目には、安定期間Ｔｃが開始される測定値メモリ８１のエントリ番号、ＣＬ３列目には安定期間Ｔｃが終了する測定値メモリ８１のエントリ番号が格納される。ＣＬ４列目には安定期間Ｔｃの長さｗが格納される。ＣＬ５列目には各安定期間Ｔｃの開始時刻が格納される。ＣＬ６行目には各安定期間Ｔｃの平均温度Ｘａｖが格納される。ＣＬ７列目には各安定期間Ｔｃの平均導電率Ｙａｖが格納される。
　ＣＬ８列目には各安定期間Ｔｃの温度偏差Ｓｘが格納される。ＣＬ９列目には各安定期間Ｔｃの導電率偏差Ｓｙが格納される。ＣＬ１０列目には各安定期間Ｔｃの共分散Ｓｘｙが格納される。ＣＬ１１列目には各安定期間Ｔｃの相関係数Ｒ２が格納される。
　ＣＬ１２列目には１次係数ａが格納される。１次係数ａはａ＝Ｓｘｙ／Ｓｘ２により算出される値である。ＣＬ１３列目には０次係数ｂが格納される。０次係数ｂはｂ＝Ｙａｖ－ａ＊Ｘａｖにより算出される値である。ＣＬ１４列目にはゼロ温度ｔ０が格納される。ゼロ温度ｔ０はｔ０＝ｂ／ａにより算出される値である。ＣＬ１５列目には温度係数候補値αｔｃが格納される。温度係数候補値αｔｃはαｔｃ＝１／（Ｘａｖ－ｔ０）により算出される値である。ＣＬ１６列目には温度係数制限Ｚが格納される。温度係数制限Ｚは、温度係数範囲αｍｉｎ／αｍａｘの範囲内に温度係数候補値αｔｃがあれば１となり、温度係数候補値αｔｃが温度係数範囲αｍｉｎ／αｍａｘの範囲外であれば０となるフラグである。温度係数制限Ｚが０であるエントリは、後述するマトリクスメモリ８６を用いた処理の対象外とされる。
　ＣＬ１７列目には相関係数平均Ｒ２ａｖが格納される。ＣＬ１８列目には共分散平均Ｓｘｙａｖ／Ａが格納される。ＣＬ１９列目には加重係数Ｋが格納される。加重係数ＫはＫ＝ｗ＊Ｓｘｙにより算出される値である。ＣＬ２０列目には温度係数加重平均αａｖが格納される。温度係数加重平均αａｖはαａｖ＝Σαｔｃ＊Ｋ／ΣＫにより算出される値である。なお、行列メモリ８５では、下に位置する行に格納される相関係数平均Ｒ２ａｖ、共分散平均Ｓｘｙａｖ／Ａ、温度係数加重平均αａｖは、それより上の行のエントリのうち温度係数制限Ｚが１となったエントリの値を含む平均値であり、最も下のエントリの相関係数平均Ｒ２ａｖ、共分散平均Ｓｘｙａｖ／Ａ、温度係数加重平均αａｖが統計処理対象とした変動幅Ａと安定期間カウント値Ｗについての統計結果となる。
　続いて、温度係数算出部８４は、行列メモリ８５に格納されているテーブルデータを統計処理の対象としている変動幅Ａと安定期間カウント値Ｗとについて生成した後に、マトリクスメモリ８６に統計データを集計したデータを書き込み、その後行列メモリ８５に格納したデータをクリアして次の変動幅Ａと安定期間カウント値Ｗの組み合わせに利用する。
　温度係数算出部８４は、行列メモリ８５に格納した統計データの集計処理においては、行列メモリ８５に格納されたデータの最も下の行のエントリの相関係数平均Ｒ２ａｖ、共分散平均Ｓｘｙａｖ／Ａ、温度係数加重平均αａｖと有効エントリ数（温度係数制限Ｚのフラグ値が１となっているエントリ数、有効Ｔｃ数）と全エントリ数（全Ｔｃ数）を出力する。
　ここで、図１６に実施の形態４にかかる測定値温度係数補正部７１のマトリクスメモリ８６に蓄積されるデータのフォーマットを説明する図を示す。図１６に示すようにマトリクスメモリ８６には、行方向に変動幅Ａ毎のエントリが並べられ、列方向に安定期間カウント値Ｗ毎のエントリが並べられる。そして、温度係数算出部８４は、行列メモリ８５を用いた統計処理の対象とした変動幅Ａと安定期間カウント値Ｗに該当するエントリに、集計した相関係数平均Ｒ２ａｖ、共分散平均Ｓｘｙａｖ／Ａ、温度係数加重平均αａｖと有効エントリ数（温度係数制限Ｚのフラグ値が１となっているエントリ数、有効Ｔｃ数）と全エントリ数（全Ｔｃ数）を格納する。
　そして、温度係数算出部８４は、全ての変動幅Ａの候補と安定期間カウント値Ｗの候補との組み合わせについて統計処理と集計処理を行い、マトリクスメモリ８６のエントリを埋める。その後、温度係数算出部８４は、予め設定された有効係数条件を満たすエントリを探索する。有効係数条件は、例えば、以下の５条件を満たすエントリである。
第１条件：有効Ｔｃ数が２以上
第２条件：相関平均が０．５以上
第３条件：共分散平均が０．５以上
第４条件：共分散平均が最大となるエントリ
第５条件：第１条件～第４条件を満たすエントリが複数有った場合、相関平均が最も大きいエントリを選出する。
温度係数算出部８４は、上記５条件をエントリ毎に検証し、全ての条件を満たすエントリに含まれる温度係数加重平均αａｖを温度係数αとして決定し、導電率温度補正部８７に出力する。一方、温度係数算出部８４は、上記５条件を満たすエントリがない場合、温度係数算出パラメータ保存部８２に格納されている初期温度係数α０を温度係数αとして出力する。また、温度係数算出部８４は、温度係数αとともに溶液基準温度Ｔｒｅｆを温度係数算出パラメータ保存部８２から読み出して導電率温度補正部８７に与える。
　導電率温度補正部８７は、温度係数αと溶液基準温度Ｔｒｅｆとを用いて導電率測定値Ｃを補正して導電率測定値Ｃｄを出力する。導電率温度補正部８７は、例えば、（２４）式に基づき温度補正を行う。
Figure JPOXMLDOC01-appb-M000024
　上記説明より、実施の形態４にかかる相関ＢＲＩＸ算出部２２ｃでは、測定値温度係数補正部７１において導電率測定値Ｃｄの温度変動成分をキャンセルする補正を行う。これにより、実施の形態４にかかる相関ＢＲＩＸ算出部２２ｃが出力される相関ＢＲＩＸ値Ｂｄには、温度による誤差が含まれない。つまり、実施の形態４にかかる相関ＢＲＩＸ算出部２２ｃを用いることで相関ＢＲＩＸ値Ｂｄの精度を高めることができる。
　なお、本発明は上記実施の形態に限られたものではなく、趣旨を逸脱しない範囲で適宜変更することが可能である。
　　　（付記１）
　クーラント液の導電率の経時変化の測定結果から予め作成した経時変化関数を用いて基準日から現在までの前記導電率の変化量を補正する経時変化係数を算出する経時変化係数算出部と、
　クーラント液の導電率を計測する導電率計から導電率測定値を取得して、前記導電率測定値に前記経時変化係数を適用して、前記導電率測定値から経時的な電解成分濃度変化の影響を取り除いた補正導電率を算出する導電率変換部と、
　前記補正導電率を予め決定した変換式に適用して前記補正導電率に対応する相関ＢＲＩＸ値を算出する相関ＢＲＩＸ変換部と、
　を有するＢＲＩＸ測定装置。
　　　（付記２）
　前記基準日と、前記基準日に測定された前記クーラント液の前記導電率である基準導電率と、前記基準日に測定された前記クーラント液のＢＲＩＸである基準ＢＲＩＸ値と、前記クーラント液に混ぜられている希釈液の導電率である希釈水の前記導電率である希釈水導電率と、を保持するパラメータ保存部をさらに有し、
　前記経時変化関数は、前記導電率測定値の変化を非線形近似した関数であり、
　前記経時変化係数算出部は、前記基準導電率と、前記基準ＢＲＩＸ値と、前記希釈水導電率と、前記基準日と、を前記経時変化関数に適用して前記経時変化係数を算出する付記１に記載のＢＲＩＸ測定装置。
　　　（付記３）
　前記基準日と、前記基準日に測定された前記クーラント液の導電率である第１基準導電率と、前記基準日に測定された前記クーラント液のＢＲＩＸである第１基準ＢＲＩＸ値と、前記第１基準導電率を測定した前記クーラント液をさらに予め設定された濃度にまで希釈を進めた希釈クーラント液の導電率である希釈後導電率と、前記希釈クーラント液のＢＲＩＸである希釈後ＢＲＩＸ値と、前記クーラント液に混ぜられている希釈液の導電率である希釈水導電率と、前記基準日より後の計測日と、前記計測日に測定された前記クーラント液の前記導電率である第２基準導電率と、前記計測日に測定された前記クーラント液のＢＲＩＸである第２基準ＢＲＩＸ値と、を保持するパラメータ保存部をさらに有し、
　前記経時変化関数は、前記導電率測定値の変化を線形近似し、所定の期間で傾きが切り替わる関数であり、
　前記経時変化係数算出部は、前記第２基準導電率と、前記第２基準ＢＲＩＸ値と、前記希釈水導電率と、前記基準日と、前記計測日を前記経時変化関数に適用して前記経時変化係数を算出する付記１に記載のＢＲＩＸ測定装置。
　　　（付記４）
　前記相関ＢＲＩＸ変換部は、前記第１基準導電率、前記第１基準ＢＲＩＸ値、前記希釈後導電率、前記希釈後ＢＲＩＸ値、を用いて前記変換式における係数を決定する付記３に記載のＢＲＩＸ測定装置。
　　　（付記５）
　前記相関ＢＲＩＸ変換部の前記変換式を決定した際に前記クーラント液を希釈した希釈水の導電率を示す希釈水導電率と、クーラント槽に注液したクーラント液を希釈した希釈水の導電率との差を補正する希釈水補正部をさらに有する付記１乃至４のいずれか１項に記載のＢＲＩＸ測定装置。
　　　（付記６）
　前記導電率測定値を前記クーラント液の温度に応じて補正して、温度補正済みの前記導電率測定値を前記導電率変換部に出力する測定値温度係数補正部をさらに有し、
　前記測定値温度係数補正部は、
　前記クーラント液の温度を計測する温度センサから取得した温度測定値と前記導電率測定値とを関連付けて時系列に記憶する測定値メモリと、
　前記測定値メモリに記憶された前記導電率測定値と前記温度測定値とに対して行う統計処理から導き出した温度係数を算出する温度係数算出部と、
　前記温度係数を用いて前記導電率測定値を補正して、前記導電率変換部に温度補正済みの前記導電率測定値を出力する導電率温度補正部と、を有する付記１乃至５のいずれか１項に記載のＢＲＩＸ測定装置。
　　　（付記７）
　前記クーラント液の温度の基準値を示す溶液基準温度と、前記温度係数の初期値である初期温度係数と、前記温度係数の上限と下限の範囲を示す温度係数範囲と、計算対象とする温度計測の期間を示す計算期間と、を温度係数算出パラメータとして記憶する温度係数算出パラメータ記憶部をさらに有し、
　前記温度係数算出部は、
　前記測定値メモリを参照し、時系列で連続する２つの前記導電率測定値の差が予め設定された変動幅以下になる安定期間の長さを異なる前記変動幅毎に算出し、
　前記安定期間毎に前記安定期間に含まれる前記導電率測定値及び前記温度測定値の平均値、偏差、共分散、相関係数を用いて温度係数加重平均を算出し、
　前記共分散の平均値と前記相関係数の平均値がともに既定値以上にある前記安定期間の長さと前記変動幅の大きさの組み合わせのうち前記共分散の平均値が最も大きな組み合わせにおいて算出される前記温度係数加重平均を前記温度係数として算出する付記６に記載のＢＲＩＸ測定装置。
　　　（付記８）
　前記温度係数算出部は、
　前記温度係数加重平均の算出に用いられる温度係数候補値が前記温度係数範囲から外れた当該温度係数候補値が含まれる前記安定期間のデータを無効とし、温度係数候補値が前記温度係数範囲内の当該温度係数候補値が含まれる前記安定期間のデータを有効とし、
　有効な前記安定期間のデータから算出された前記温度係数加重平均を前記温度係数の算出対象とする付記７に記載のＢＲＩＸ測定装置。
　　　（付記９）
　前記温度係数算出部は、
　前記共分散の平均値と前記相関係数の平均値がともに既定値以上にある前記安定期間の長さと前記変動幅の大きさの組み合わせがない場合、前記初期温度係数を前記温度係数として出力する付記７又は８に記載のＢＲＩＸ測定装置。
　　　（付記１０）
　演算部で実行されることで、クーラント液の導電率を計測する導電率計から取得した導電率測定値から前記クーラント液のＢＲＩＸ値を算出するＢＲＩＸ測定プログラムであって、
　予め前記クーラント液の導電率の経時変化の測定結果から作成した経時変化関数を用いて基準日から現在までの前記導電率の変化量を補正する経時変化係数を算出する経時変化係数算出処理と、
　前記導電率測定値に前記経時変化係数を適用して、前記導電率測定値から経時的な電解成分濃度変化の影響を補正により取り除いた補正導電率を算出する導電率変換処理と、
　前記補正導電率を予め決定した変換式に適用して前記補正導電率に対応する相関ＢＲＩＸ値を算出する相関ＢＲＩＸ変換処理と、
　を行うＢＲＩＸ測定プログラム。
　１、４　クーラント液管理装置
　１０　制御ボックス
　１１　電源部
　１２　演算部
　１３　通信インタフェース
　１４　外部通信インタフェース
　１５　制御信号生成部
　２０　センサボックス
　２１　電源部
　２２　相関ＢＲＩＸ算出部
　２３　制御信号生成部
　２４　導電量計
　２５　アナログデジタル変換回路
　２６　通信インタフェース
　３０　正弦波生成部
　３１　アンプ
　３２　一次側トロイダルトランス
　３３　二次側トロイダルトランス
　３４　ローノイズアンプ
　３５　信号振幅検出器
　４１　データ格納部
　４２、５２、６２　パラメータ保存部
　４３　内蔵時計
　４４、５４　経過時間カウンタ
　４５、５５、６５　経時変化係数算出部
　４６、５６　導電率変換部
　４７、５７、６７　相関ＢＲＩＸ変換部
　６８　希釈水補正部
　７１　測定値温度係数補正部
　８１　測定値メモリ
　８２　温度係数算出パラメータ保存部
　８３　内蔵時計
　８４　温度係数算出部
　８５　行列メモリ
　８６　マトリクスメモリ
　８７　導電率温度補正部
　ＴＳ　温度センサ
　ＩＮＤ　表示部

Claims (5)

  1. 　溶液の温度を計測する温度センサから取得した温度測定値と、電気伝導方式によって計測される前記溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値と、を関連付けて時系列に記憶する測定値メモリと、
    　前記測定値メモリに記憶された前記溶液測定値と前記温度測定値とについて統計処理から導き出した温度係数を算出する温度係数算出部と、
    　前記温度係数を用いて前記溶液測定値を補正した温度補正済みの前記溶液測定値を出力する測定値温度補正部と、を有する測定値温度係数補正装置。
  2. 　前記溶液の温度の基準値を示す溶液基準温度と、前記温度係数の初期値である初期温度係数と、前記温度係数の上限と下限の範囲を示す温度係数範囲と、計算対象とする温度計測の期間を示す計算期間と、を温度係数算出パラメータとして記憶する温度係数算出パラメータ記憶部をさらに有し、
    　前記温度係数算出部は、
    　前記測定値メモリを参照し、時系列で連続する２つの前記溶液測定値の差が予め設定された変動幅以下になる安定期間の長さを異なる前記変動幅毎に算出し、
    　前記安定期間毎に前記安定期間に含まれる前記溶液測定値及び前記温度測定値の平均値、偏差、共分散、相関係数を用いて温度係数加重平均を算出し、
    　前記共分散の平均値と前記相関係数の平均値がともに既定値以上にある前記安定期間の長さと前記変動幅の大きさの組み合わせのうち前記共分散の平均値が最も大きな組み合わせにおいて算出される前記温度係数加重平均を前記温度係数として算出する請求項１に記載の測定値温度係数補正装置。
  3. 　前記温度係数算出部は、
    　前記温度係数加重平均の算出に用いられる温度係数候補値が前記温度係数範囲から外れた当該温度係数候補値が含まれる前記安定期間のデータを無効とし、温度係数候補値が前記温度係数範囲内の当該温度係数候補値が含まれる前記安定期間のデータを有効とし、
    　有効な前記安定期間のデータから算出された前記温度係数加重平均を前記温度係数の算出対象とする請求項２に記載の測定値温度係数補正装置。
  4. 　前記温度係数算出部は、
    　前記共分散の平均値と前記相関係数の平均値がともに既定値以上にある前記安定期間の長さと前記変動幅の大きさの組み合わせがない場合、前記初期温度係数を前記温度係数として出力する請求項２に記載の測定値温度係数補正装置。
  5. 　演算部で実行されることで、温度係数を用いて測定値を補正する測定値温度係数補正プログラムであって、
    　溶液の温度を計測する温度センサから取得した温度測定値と、電気伝導方式によって計測される前記溶液の特性値であって温度によって値が変化する温度依存性を有する溶液測定値と、を関連付けて時系列に測定値メモリに記憶する測定値記憶処理と、
    　前記測定値メモリに記憶された前記溶液測定値と前記温度測定値とについて統計処理から導き出した温度係数を算出する温度係数算出処理と、
    　前記温度係数を用いて前記溶液測定値を補正した温度補正済みの前記溶液測定値を出力する測定値温度補正処理と、
    　を前記演算部に実行させる測定値温度係数補正プログラム。
PCT/JP2022/035234 2021-10-15 2022-09-21 測定値温度係数補正装置及び測定値温度係数補正プログラム WO2023063044A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2023555064A JPWO2023063044A1 (ja) 2021-10-15 2022-09-21

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2021169463 2021-10-15
JP2021-169463 2021-10-15

Publications (1)

Publication Number Publication Date
WO2023063044A1 true WO2023063044A1 (ja) 2023-04-20

Family

ID=85987669

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/035234 WO2023063044A1 (ja) 2021-10-15 2022-09-21 測定値温度係数補正装置及び測定値温度係数補正プログラム

Country Status (2)

Country Link
JP (1) JPWO2023063044A1 (ja)
WO (1) WO2023063044A1 (ja)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040085072A1 (en) * 2002-11-05 2004-05-06 The Furukawa Battery Co., Ltd. Method and system for monitoring state of lead acid battery
WO2019151295A1 (ja) * 2018-01-31 2019-08-08 センスプロ株式会社 オイル状態判定システム、オイル状態判定方法、及び、オイル状態判定プログラム
JP2021041487A (ja) * 2019-09-10 2021-03-18 ディー・クルー・テクノロジーズ株式会社 クーラント液管理装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040085072A1 (en) * 2002-11-05 2004-05-06 The Furukawa Battery Co., Ltd. Method and system for monitoring state of lead acid battery
WO2019151295A1 (ja) * 2018-01-31 2019-08-08 センスプロ株式会社 オイル状態判定システム、オイル状態判定方法、及び、オイル状態判定プログラム
JP2021041487A (ja) * 2019-09-10 2021-03-18 ディー・クルー・テクノロジーズ株式会社 クーラント液管理装置

Also Published As

Publication number Publication date
JPWO2023063044A1 (ja) 2023-04-20

Similar Documents

Publication Publication Date Title
JP6153528B2 (ja) 電池監視のためのシステムおよび方法
JP5827426B1 (ja) 予兆診断システム及び予兆診断方法
JP5827425B1 (ja) 予兆診断システム及び予兆診断方法
US20160169585A1 (en) Dynamic monitoring, diagnosis, and control of cooling tower systems
JP2018169161A (ja) 電池の劣化診断装置、劣化診断方法、及び劣化診断システム
BR112015001689A2 (pt) desenvolvimento de projeto e implementação de sistema e algoritmo de controle baseados em analisador
EP2884354A1 (en) Model-based predictive controller with steady-state model adaptation
CN105866504A (zh) 一种基于卡尔曼滤波的光纤电流互感器温度补偿方法
WO2016147722A1 (ja) 推定装置、推定方法およびプログラム
KR101737968B1 (ko) 발전플랜트 성능데이터 학습에 의한 발전기 출력 추정 방법
WO2023063044A1 (ja) 測定値温度係数補正装置及び測定値温度係数補正プログラム
JP5939439B2 (ja) 異常診断装置
CN117169638B (zh) 变压器绝缘纸老化程度检测方法、系统、设备及存储介质
KR101433724B1 (ko) 하위 구성요소의 상태가 반영된 플랜트의 건강지수 측정방법 및 그 방법을 수행하기 위한 프로그램이 저장된 컴퓨터 판독 가능한 저장매체
JP6715800B2 (ja) 燃料削減率出力システム、燃料削減率出力方法、および燃料削減率出力プログラム
KR101325638B1 (ko) 하위 구성요소의 상태가 반영된 플랜트의 건강지수 측정방법 및 그 방법을 수행하기 위한 프로그램이 저장된 컴퓨터 판독 가능한 저장매체
JP5772464B2 (ja) 寿命推定装置、寿命推定方法及びプログラム
CN110020000B (zh) 判定异常风速数据的方法和装置
WO2023184237A1 (zh) 计算电子系统的剩余使用寿命的方法、装置及计算机介质
He et al. Modeling and detecting approach for the change point of electronic product infant failure rate
CN112685449B (zh) 一种热力系统性能计算方法、装置及设备
JP2016045853A (ja) 異常診断装置及び異常診断方法
US20190243349A1 (en) Anomaly analysis method, program, and system
CN112825104B (zh) 基于温湿度因子的叠层片式电感可靠性预计修正方法及介质
WO2012048734A1 (en) Method for controlling an industrial process

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22880740

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2023555064

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE