The present invention relates to a data processing device, a data processing method and a recording medium that stores data processing program that are intended to process data acquired by gas chromatographic analysis, and for example, are used for data processing for predicting retention time at the time of temperature rising analysis on the basis of known data on a certain compound, which is acquired by isothermal analysis.
When a compound is qualified, for example, retention time of the target compound is measured by the gas chromatographic analysis, the measured retention time is converted into adjustment retention time or a retention index, and a database is searched for compounds having similar adjustment retention times or retention indexes to narrow down potential compounds.
However, the adjustment retention time and the retention index vary depending on temperature conditions at the time when the retention time is measured. More specifically, the adjustment retention time and the retention index vary depending on whether the retention time is measured by the isothermal analysis or the temperature rising analysis, and what temperature rising rate is set in the temperature rising analysis.
For this reason, in the case where an actually measured value and a known value in the database about a certain compound cannot be compared with each other by the gas chromatographic analysis due to varying temperature conditions, the adjustment retention time and the retention index in the actual temperature condition are predicted on the basis of the known adjustment retention time and retention index registered in the database (refer to Patent Literature 1).
There will be described a method of predicting the retention time to be measured by the temperature rising analysis at a predetermined temperature rising rate on the basis of the retention time measured by the isothermal analysis, which is registered in the database.
It is known that the relationship as represented by Expression 1 holds for the adjustment retention time measured on a certain compound by the isothermal analysis and the temperature of a column at the time of the isothermal analysis.
Where, tR′ represents the adjustment retention time and T represents the column temperature at the time of the isothermal analysis. That is, there is a linear relationship between a logarithm of the adjustment retention time at the time of the isothermal analysis and a reciprocal of the column temperature at the time of the isothermal analysis. Coefficients a and b can be determined by acquiring two adjustment retention times of the compound at different temperatures at the time of the isothermal analysis, which are registered in the database, and substituting the values into Expression 1 to solve the simultaneous equation. The retention time tR at any temperature T at the time of the isothermal analysis can be also found according to the Expression 1.
Following Expressions 2 to 4 are known as general expressions for predicting the retention time in the temperature rising analysis on the basis of the retention time in the isothermal analysis.
[Expression 2]
L=∫
0
ir
v(T)dt (2)
[Expression 3]
v=L/t
u (3)
[Expression 4]
T=rt+T, (4)
Where, L represents column length, v(T) represents compound moving rate in the column at temperature T, r represents setting temperature rising rate, t represents time, Ti represents initial temperature, tR represents retention time in the isothermal analysis, and tr represents predicted retention time in the temperature rising analysis.
Here, although the column length L is an unknown value, L is present in both the left-hand side and the right-hand side in the state where Expression 3 is substituted into Expression 2 and thus, can be deleted, resulting in that Expression 5 is solved.
That is, the time tr at the time when a definite integral of the right-hand side in Expression 5 becomes 1 is the retention time in the temperature rising analysis. The definite integral can be calculated by a sectional quadrature method of dividing integral interval into many sectors every short time period Δt as represented by Expression 6 by using a function between the retention time tR and the temperature T according to Expression 1 and a function between the temperature T and the time t according to Expression 4.
More specifically, as shown in
However, the predicted retention time in the temperature rising analysis, which is predicted based on known data in the isothermal analysis by using the above-mentioned method may largely deviated from the actually measured retention time actually measured when an oven is set such that the column temperature achieves the same temperature rising rate as the temperature rising rate used for prediction to perform the gas chromatographic analysis.
The Inventor earnestly examined this phenomenon, and found that the phenomenon occurred because the setting temperature rising rate set to the oven is deviated from the actual temperature rising rate that is the actual temperature rising rate in the column as represented in each graph in
That is, since the setting temperature rising rate set to the oven is only a target value, and the temperature in the column cannot be directly measured, the actual temperature rising rate in the column is likely to be late as shown in
Although deviation of the setting temperature rising rate from the actual temperature rising rate is assumed, since the temperature in the column cannot be actually measured, it is difficult to quantitatively assess the deviation by direct measurement. In addition, since the oven used for the gas chromatographic analysis varies, a deviation amount of the setting temperature rising rate from the actual temperature rising rate can differ according to a gas chromatograph that measures the retention time. Thus, it is difficult to formulate or generalize the deviation amount of the setting temperature rising rate from the actual temperature rising rate, and a method of quantitatively assess a temperature rising characteristic of each chromatograph has not been known yet.
Also in the isothermal analysis in addition to the above-mentioned temperature rising analysis, an error between the actually measured retention time and the predicted retention time calculated by Expression 1 and so on can occur. This is probably due to an error between a setting sustained temperature set to the oven and an actual sustained temperature in the column.
That is, due to the error between setting temperature conditions set to the oven such as the setting temperature rising rate and the setting sustained temperature, and actual temperature conditions such as the actual temperature rising rate and the actual sustained temperature in the column, the calculated predicted retention time may be largely deviated from the actually measured retention time.
In consideration of these matters, to further improve the accuracy of qualitative gas chromatographic analysis including the temperature rising analysis as well as the isothermal analysis, it is need to quantitatively assess the difference between the setting temperature condition set to the oven and the actual temperature condition in the column, which has not conventionally considered, or to use a similar condition to the actual temperature condition when the predicted retention time is calculated, thereby calculating the predicted retention time more accurately.
Patent Literature 1: International Publication No. WO2012/118209
The present invention is devised in consideration of the above-mentioned problem, and its object is to provide a data processing device for gas chromatograph, a data processing method, and recording medium that stores a data processing program that can estimate the deviation amount of the setting temperature condition from the actual temperature condition in the column, which has not conventionally considered, and improve the accuracy of predicting the predicted retention time and retention index from known data by the isothermal analysis on the basis of the estimated deviation amount.
That is, a data processing device for gas chromatograph according to the present invention includes an actually measured retention time storage part configured to store an actually measured retention time actually measured on a compound by isothermal analysis or temperature rising analysis in the state where a predetermined setting temperature condition is set to an oven, a predicted retention time storage part configured to store a predicted retention time calculated based on known data by the isothermal analysis on the compound and an assumed temperature condition assumed on a temperature condition in a column, and an actual temperature condition estimation part configured to estimate an actual temperature condition in the column on the basis of the actually measured retention time and the predicted retention time.
Here, the temperature condition conceptually contains sustained temperature in the isothermal analysis, the temperature rising rate in the temperature rising analysis, and an initial setting temperature. The retention time conceptually contains an adjustment retention time acquired by subtracting a void time from the retention time.
As such, since the actual temperature condition estimation part is configured to estimate the actual temperature condition in the column on the basis of the actually measured retention time and the predicted retention time, the actual temperature rising condition can be estimated without measuring the actual temperature condition in the column, which is difficult to be measure. Further, the predicted retention time can be calculated based on the estimated actual temperature rising condition more accurately, improving the estimation accuracy of the quality of an unknown sample.
To accurately estimate the temperature rising rate, which is especially difficult to be estimated by the gas chromatographic analysis, the actually measured retention time storage part may be configured to store the actually measured retention time actually measured on the compound by the temperature rising analysis in the state where a predetermined setting temperature rising rate that is the setting temperature condition is set to the oven, the predicted retention time storage part may be configured to store the predicted retention time calculated based on known data by the isothermal analysis on the compound, and an assumed temperature rising rate that is assumed as the assumed temperature condition and includes a rate correction value, and the actual temperature condition estimation part may be configured to estimate an actual temperature rising rate as the actual temperature condition in the column on the basis of the actually measured retention time and the predicted retention time.
Here, in “the assumed temperature rising rate acquired by adding the rate correction value to the setting temperature rising rate”, the rate correction value may be 0, and the assumed temperature rising rate may be the same as the setting temperature rising rate.
As such, since the actual temperature condition estimation part is configured to estimate the actual temperature rising rate that is the actual temperature rising rate in the column on the basis of the actually measured retention time and the predicted retention time, the actual temperature rising rate can be found without directly measuring the temperature in the column.
Next, a reason why the actual temperature rising rate in the column can be estimated according to the present invention will be qualitatively described.
It is likely that the actually measured retention time does not indicate a value measured when the temperature is raised at the setting temperature rising rate, but indicate a value measured when the temperature is raised at the actual temperature rising rate that is the actual temperature rising rate in the column. Thus, in the case where the value of the actually measured retention time matches or is close to the value of the predicted retention time, it is likely that the assumed temperature rising rate used to calculate the predicted retention time substantially matches the actual temperature rising rate.
Accordingly, for example, when an error between the actually measured retention time and the predicted retention time is assessed, and the error is smaller than a predetermined value, the assumed temperature rising rate used to calculate the predicted retention time can be estimated to be the actual temperature rising rate. Alternatively, for example, the rate correction values may be assigned in a predetermined range to calculate a plurality of predicted retention times, and the assumed temperature rising rate used to calculate the predicted retention time closest to the actually measured retention time can be estimated to be the actual temperature rising rate.
As described above, the actual temperature condition estimation part can estimate the actual temperature rising rate on the basis of the actually measured retention time and the predicted retention time.
To further improve the prediction accuracy of the calculated retention time in the temperature rising analysis by using the estimated actual temperature rising rate, a post-correction predicted retention time calculation part configured to calculate a post-correction predicted retention time on the basis of the actual temperature rising rate estimated by the actual temperature condition estimation part and the known data by the isothermal analysis may be further provided.
Even in the case where it is difficult to predict the value of the assumed temperature rising rate and therefore, the rate correction value that achieves a value close to the actual temperature rising rate cannot be set and only a random rate correction value can be set, to accurately estimate the actual temperature rising rate on the basis of the actually measured retention time and the predicted retention time, the predicted retention time storage part may be configured to store a plurality of predicted retention times calculated based on a plurality of rate correction values on the compound, and the actual temperature condition estimation part may be configured to include a relational expression identification part configured to identify a rate correction value-error relational expression that is a relational expression between the rate correction value and a retention time error on the basis of a plurality of combinations of the rate correction value and the retention time error that is an error between the actually measured retention time and the predicted retention time, a proper rate correction value calculation part configured to calculate a proper rate correction value that is the rate correction value at the time when the retention time error becomes a predetermined value on the basis of the rate correction value-error relational expression identified by the relational expression identification part, and an actual temperature rising rate calculation part configured to calculate the actual temperature rising rate on the basis of the proper rate correction value calculated by the proper rate correction value calculation part and the setting temperature rising rate.
As such, for example, the rate correction value-error relational expression indicating the relationship between the rate correction value and the retention time error that is the error between the actually measured retention time and the predicted retention time may be created from two predicted retention times calculated at the assumed temperature rising rates using two randomly-selected different rate correction values, and the proper rate correction value that is the rate correction value at which the error between the predicted retention time and the actually measured retention time is less than a predetermined error can be calculated by the rate correction value-error relational expression. Further, since the actually measured retention time substantially matches the predicted retention time, the temperature rising rate calculated by applying the proper rate correction value to the setting temperature rising rate becomes substantially same as the actual temperature rising rate.
Both in the case where the actual temperature rising rate is offset from the setting temperature rising rate in the time-axis direction, and where the error between the setting temperature rising rate and the actual temperature rising rate varies with time elapsed from start of temperature rising, to accurately estimate the actual temperature rising rate, the actually measured retention time storage part may be configured to store the actually measured retention time of each of the plurality of compounds, the predicted retention time storage part may be configured to store the plurality of predicted retention times calculated based on the assumed temperature rising rates of different rate correction values on each of the plurality of compounds, and the actual temperature condition estimation part may be configured to divide a column temperature rising section in which at least a column temperature is raised into a plurality of divided sections every predicted retention time calculated using the actually measured retention time or the setting temperature rising rate of each compound, and estimate the actual temperature rising rate in each divided section.
To estimate the actual temperature rising rate on the basis of the actually measured retention time and the predicted retention time on the compound, which reflect the temperature rising rate in each divided section the best to further improve the estimation accuracy, the relational expression identification part may be configured to identify the rate correction value-error relational expression of each compound on which the actually measured retention time is actually measured, the proper rate correction value calculation part may be configured to calculate the proper rate correction value on the basis of the rate correction value-error relational expression of each compound, and the actual temperature rising rate calculation part may be configured to calculate an intermediate corrected temperature rising rate that is a temperature rising rate acquired by adding a minimum error correction value calculated by the rate correct ion value-error relational expression of the compound corresponding to the actually measured retention time or the predicted retention time indicating an end of each divided section, to the setting temperature rising rate for each divided section, and calculate the actual temperature rising rate on the basis of the intermediate corrected temperature rising rate.
To evenly perform fitting of the actual temperature rising rate estimated for each divided section in almost all divided sections, not only in some divided sections to decrease the retention time error that is the error between the actually measured retention time and the predicted retention time, the actual temperature rising rate calculation part may be configured to calculate a final corrected temperature rising rate that is the temperature rising rate acquired by adding a rate recorrection value to the intermediate corrected temperature rising rate calculated for each divided section, and output the final corrected temperature rising rate at which an error between the predicted retention time of each compound, which is calculated using the final corrected temperature rising rate, and the actually measured retention time becomes minimum, as the actual temperature rising rate.
To accurately estimate the temperature condition in the column on the basis of an actual measurement result of the retention time of one or a plurality of compounds without directly measuring the temperature in the column, a data processing method for gas chromatograph including a retention time actually measuring step of actually measuring an actually measured retention time on a compound by isothermal analysis or temperature rising analysis in the state where a predetermined setting temperature condition is set to an oven, a retention time predicting step of calculating a predicted retention time on the basis of known data by the isothermal analysis on the compound and an assumed temperature condition assumed on a temperature condition in a column, and an actual temperature condition estimating step of estimating an actual temperature condition in the column on the basis of the actually measured retention time and the predicted retention time may be used.
To add a feature for estimating the actual temperature condition to the existing data processing device for gas chromatograph to improve the prediction accuracy of the predicted retention time and enable more accurate qualification of the compound, a data processing program for gas chromatograph including an actually measured retention time storage part configured to store an actually measured retention time actually measured on a compound by isothermal analysis or temperature rising analysis in the state where a predetermined setting temperature condition is set to an oven, a predicted retention time storage part configured to store a predicted retention time calculated based on known data by the isothermal analysis on the compound and an assumed temperature condition assumed on a temperature condition in a column, and an actual temperature condition estimation part configured to estimate an actual temperature condition in the column on the basis of the actually measured retention time and the predicted retention time, or a recording medium that records the data processing program for gas chromatograph may be used.
As described above, the data processing device for gas chromatograph, the data processing method, and the data processing program according to the present invention can indirectly estimate the actual temperature condition such as the actual temperature rising rate on the basis of the actually measured retention time and the predicted retention time predicted from known data without directly measuring the temperature in the column. Moreover, since the retention time error between the actually measured retention time and the predicted retention time reflects a difference between the assumed temperature condition including the setting temperature condition and the correction condition and the actual temperature condition, the actual temperature condition can be estimated based on the actually measured retention time and the predicted retention time.
An embodiment of the present invention will be described with reference to figures.
As shown in
As shown in
The user can set various temperature conditions including temperature rising analysis conditions and isothermal analysis conditions to the oven 14, at least initial oven temperature, setting temperature rising rate, temperature rising start time, and temperature rising end time can be set as the temperature conditions in the temperature rising analysis, and setting sustained temperature can be set as the temperature condition in the isothermal analysis. A temperature sensor (not shown) is provided in the oven 14 to perform temperature feedback control so as to minimize a difference between a target temperature at each time determined depending on the set initial oven temperature and setting temperature rising rate, and a measurement temperature measured by the temperature sensor. That is, such feedback control generates a slight error, such as a delay, between the setting temperature rising rate that is a target value set by the user and the actual temperature rising rate that is the actual temperature rising rate in the column 13.
The data processing device for gas chromatograph 100 is used to compare an adjustment retention time and a retention index on a compound, which are acquired by the gas chromatographic analysis, with the adjustment retention time and the retention index on the compound, which are recorded in the database or calculated based on known data to qualify the compound.
The data processing device for gas chromatograph 100 is configured such that, in the case where the gas chromatograph 101 performs the temperature rising analysis, by reflecting a temperature change occurring in the column 13 substantially as it is, a predicted retention time that is the retention time in the temperature rising analysis can be calculated from the known retention time in the isothermal analysis. In other words, the data processing device for gas chromatograph 100 is configured such that, to set the error between the predicted retention time and the actually measured retention time actually measured by the gas chromatograph 101 to a predetermined value or smaller, as needed, the temperature rising rate used to calculate the predicted retention time is corrected to a temperature rising rate other than the setting temperature rising rate set to the oven 14, to calculate the predicted retention time.
More specifically, the data processing device for gas chromatograph 100 is constituted of a so-called computer provided with a CPU, a memory, A/D, D/A converters, an input/output device, and a monitor, and is configured to, in such a manner that the components work together and execute the data processing program for gas chromatograph stored in the memory, fulfill functions of the components as shown in a functional block diagram in
That is, as shown in
Here, the retention index in the isothermal analysis is, as represented by Expression 7, a value calculated from a logarithm of the adjustment retention time of an unknown sample or compound and a logarithm of the adjustment retention time of a standard substance such as alkane. The retention index in the temperature rising analysis is, as represented by Expression 7′, a value calculated from the adjustment retention time of an unknown sample or a compound and the adjustment retention time of the standard substance such as alkane.
Where, I represents retention index, t′(X) represents adjustment retention time of unknown sample or compound, t′(PZ) represents adjustment retention time of alkane with carbon number Z, t′(PZ+1) represents adjustment retention time of alkane with carbon number Z+1, and Z represents carbon number of alkane.
When the temperature condition such as the temperature of the oven 14 in the gas chromatographic analysis varies, even the same compound has different retention indexes. However, the retention index can be theoretically prevented from being affected by a variation of the gas chromatograph 101 in the length of the column 13, the linear speed of carrier gas, and so on. Thus, the qualification part C1 can calculate the predicted retention index that reflects the same temperature condition as the temperature condition under which the retention index is actually measured, and compare the predicted retention index with the actually measured retention index.
That is, the qualification part C1 is configured such that, to compare the actually measured retention index with the predicted retention index, such that a value can be calculated from each of the actually measured retention time and the predicted retention time. The actually measured retention index is calculated from the retention time actually measured in the isothermal analysis or the temperature rising analysis on an unknown sample by the gas chromatograph 101, and the retention time of the standard substance. The predicted retention index is calculated the retention time of each compound converted from known data by the known data conversion part C2 so as to meet the temperature condition such as the isothermal analysis condition or the temperature rising analysis of the actually measured retention index, and the retention time of the standard substance. Then, the qualification part C1 compares the calculated actually measured retention index with the calculated predicted retention index, and when a compound having the predicted retention index matching the actually measured retention index is present, the qualification part C1 determines that the unknown sample corresponds to the matched compound, and outputs a name of the compound.
As shown in
Each of the parts will be described below.
The actually measured retention time storage part 2 serves to store the actually measured retention time that is the retention time actually measured when the gas chromatograph 101 performs the temperature rising analysis by using the setting temperature rising rate as the setting temperature condition on each of a plurality of compounds. In the present embodiment, for example, as shown in
As shown in
Where, tR represents retention time in isothermal analysis, to represents retention time (void time) of ingredients that are not held in the column 13, such as methane and noble gas, a and h represent coefficients unique to each compound, and T represents temperature. The function in Expression 8 will be hereinafter referred to as a RT enumeration function.
The isothermal-temperature rising conversion part 4 calculates the retention time in the temperature rising analysis from data known by the isothermal analysis. More specifically, the predicted retention time predicted as the retention time under any temperature rising analysis condition is calculated using the setting temperature rising rate, an initial temperature, and a rate correction value, which are set to the oven 14 by the user, by Expressions 8, 9, and 10. In other words, the isothermal-temperature rising conversion part 4 is configured to calculate the predicted retention time under the assumed temperature condition, which includes the setting temperature rising rate as a setting temperature rising condition and the rate correction value as a correction condition, and is assumed on the temperature condition in the column 13, by the Expressions 8, 9, and 10.
Where, tr represents predicted retention time in temperature rising analysis, tR(T) represents retention time in isothermal analysis at temperature T, which is found by Expression 8, r represents setting temperature rising rate, t represents time, Ti represents initial temperature, and h represents rate correction value.
Expression 10 is an expression indicating the temperature in the column 13 at the time of the temperature rising analysis. In the present embodiment, since it is assumed that the setting temperature rising rate set to the oven 14 is deviated from the actual temperature rising rate that is the actual temperature rising rate of the column 13, a value r+h of the assumed temperature rising rate acquired by adding the rate correction value h to the setting temperature rising rate r is used as a candidate for the actual temperature rising rate. As apparent from Expression 10, when the rate correction value is 0, the setting temperature rising rate set to the oven 14 is used as it is. Details of the method of determining the predicted retention time tr in the temperature rising analysis by Expression 9 are the same as those described in Background Art and thus, description thereof is omitted.
In the case where the RT enumeration function represented in Expression 8 cannot be directly acquired from the isothermal analysis known data storage part 3, the isothermal-temperature rising conversion part 4 calculates two adjustment retention times from two registered retention indexes of different temperatures in the isothermal analysis and then, identifies the RT enumeration function.
More specifically, the isothermal-temperature rising conversion part 4 finds the adjustment retention time by Expression 7 for calculating the retention index, and substitutes two known retention indexes of different temperatures measured on the same compound into Expression 11 and Expression 12 to convert the indexes into two adjustment retention times of different temperatures.
[Expression 11]
log t′T1=log t′T1(Pz)+(IT1−100×Z)(log t′T1(Pz+1)−log t′T1(Pz))/100 (11)
[Expression 12]
log t′T2=log t′T2(Pz)+(IT2−100×Z)(log t′T2(Pz+2)−log t′T2(Pz))/100 (12)
Where, IT1, IT2 represent known retention indexes at temperatures T1, T2 of compound, t′T1, t′T2 represent adjustment retention times at temperatures T1, T2 of compound, t′T1(PZ), t′T2(PZ) represent adjustment retention times at temperatures T1, T2 of alkane having carbon number of Z, t′T1(PZ+1), t′T2(PZ+1) represent adjustment retention times at known temperatures T1, T2 of alkane having carbon number of Z+1, and Z represents carbon number. In the case of using a Gauss symbol ∥ to represent the carbon number Z, Z is an integer selected by Z=[IT1] in Expression 11 and Z=[IT2] in Expression 12.
The isothermal-temperature rising conversion part 4 identifies the RT enumeration function from the two retention indexes as represented by Expression 13, by using the linear relationship between the logarithm of the adjustment retention time and the reciprocal of the temperature as represented by Expression 1 described in Background Art.
Where, T represents isothermal analysis temperature, and t′R represents adjustment retention time at isothermal analysis temperature T of compound.
In this manner, the isothermal-temperature rising conversion part 4 may identify the RT enumeration function for calculating the adjustment retention time in the isothermal analysis at any temperature T from the known retention index, and use the function in place of Expression 8.
The predicted retention time storage part stores the predicted retention time calculated based on known data by the isothermal analysis on the compound, and the assumed temperature condition including the setting temperature condition and the correction condition. More specifically, in the present embodiment, the predicted retention time storage part 5 stores the predicted retention time in the temperature rising analysis, which is converted by the isothermal-temperature rising conversion part 4, on the same compound as the compound subjected to the temperature rising analysis by the gas chromatograph 101. The predicted retention time storage part 5 stores the predicted retention time converted at the assumed temperature rising rate of different rate correction values that is the correction condition on each compound, which is set as the assumed temperature condition, the predicted retention time pairing with the rate correction value. In the present embodiment, the rate correction value is set to be −1 or more and 1 or less such that the assumed temperature rising rate is evenly distributed in a predetermined range centering on the setting temperature rising rate. The condition for the rate correction value is not limited to the above-mentioned example, and may be appropriately set according to a difference among the gas chromatographs 101.
As shown in
The actual temperature condition estimation part 6 shown in
For the above-mentioned estimation of the temperature rising rate, the actual temperature condition estimation part 6 is configured of, as shown in
The relational expression identification part 61 serves to identify a rate correction value-error relational expression that is a relational expression between the rate correction value and a retention time error on the basis of a plurality of combinations of the rate correction value and the retention time error that is an error between the actually measured retention time and the predicted retention time. In the present embodiment, the relational expression identification part 61 is configured of identify the rate correction value-error relational expression on each compound, the actually measured retention time of which is acquired. Here, using the setting temperature rising rate set to 5° C./min as a specific example,
As apparent from
When the rate correction value-error relational expression is identified on each compound, the proper rate correction value calculation part 62 shown in
Next, as shown in
More specifically, the actual temperature rising rate calculation part 63 is configured to first calculate an intermediate corrected temperature rising rate that is the temperature rising rate acquired by adding the minimum error correction value calculated by the rate correction value-error relational expression on the compound corresponding to the predicted retention time indicating the end of each divided section, to the setting temperature rising rate. That is, in the specific example shown in
In the present embodiment, as shown in
Next, the actual temperature rising rate calculation part 63 calculates a final corrected temperature rising rate that is the temperature rising rate acquired by adding a rate recorrection value to the intermediate corrected temperature rising rate as shown in
Finally, the post-correction predicted retention time calculation part 7 is configured to calculate the post-correction predicted retention time on the basis of the actual temperature rising rate estimated by the actual temperature condition estimation part 6 and known data by the isothermal analysis, which is stored in the isothermal analysis known data storage part 3, and output the predicted retention time on which the actual temperature rising rate reflects to represent a more correct value, to the qualification part C1.
As described above, the data processing device for gas chromatograph 100 in the present embodiment can estimate the actual temperature rising rate indicating the actual temperature change in the column 13 on the basis of the actually measured retention time and the predicted retention time.
As shown in
Then, the post-correction predicted retention time calculated using the actual temperature rising rate estimated by the actual temperature condition estimation part 6 is used by the qualification part C1 to calculate the predicted retention index, thereby further decreasing the error from the actual retention index.
Therefore, by comparing the actually measured retention index actually measured by the gas chromatograph 101 with the predicted retention index calculated from various known data, it is possible to determine a compound corresponding to an unknown sample with high accuracy.
Other embodiments will be described.
In the above-mentioned embodiment, the actual temperature condition estimation part is configured to output the intermediate corrected temperature rising rate recorrected according to the rate recorrection value, as the actual temperature rising rate and however, may be configured to output the intermediate corrected temperature rising rate as shown in
The assumed temperature rising rate is set by adding the rate correction value to the setting temperature rising rate and however, the assumed temperature rising rate may be set by multiplying the setting temperature rising rate by the rate correction value. Even when each parameter is determined with a random number without using the setting temperature rising condition and the correction condition, the assumed temperature condition assumed on the temperature condition in the column can be divided into the setting temperature rising condition and the correction condition. In other words, regardless of how to determine, the assumed temperature condition can be represented as the temperature condition including the setting temperature condition and the correction condition.
The actual temperature condition estimation part is configured to first identify the rate correction value-error relational expression and then, find the proper rate correction value by the identified expression and however, the rate correction values may be randomly assigned without identifying such rate correction value-error relational expression, and the rate correction value at which the calculated predicted retention time matches the actually measured retention time, or the retention time error that is the error between the calculated predicted retention time and the actually measured retention time is less than a predetermined value may be searched in a brute-force manner. In other words, the actual temperature condition estimation part may prepare a plurality of assumed temperature rising rates as candidates for the actual temperature rising rate on the basis of the setting temperature rising rate, compare the predicted retention time calculated at the assumed temperature rising rate with the actually measured retention time, and estimate the most proper assumed temperature rising rate as the actual temperature rising rate. In summary, by using the actually measured retention time and the predicted retention time, the actual temperature rising rate closer to the actual state than the setting temperature rising rate can be found.
The data processing device for gas chromatograph may be further provided with a receiving part for receiving an input from the user. More specifically, when the assumed temperature rising rate including the setting temperature rising rate and the rate correction value is set, the rate correction value may be received such that the user can set any desired rate correction value. The receiving part may be configured to accept user's selection on which of the actual temperature rising rates calculated using a plurality of correction conditions is used to calculate the predicted retention time.
In addition, in the embodiment, the temperature condition estimation part is configured to estimate the actual temperature rising rate in the column in the isothermal analysis and however, the temperature condition estimation part may be configured to estimate the actual temperature in the column with respect to the setting initial temperature set to keep the oven at a constant temperature, as another temperature condition, on the basis of the actually measured retention time and the predicted retention time.
The temperature condition estimation part may be configured to estimate an error between the setting temperature condition to the oven in the temperature rising analysis as well as the isothermal analysis and the actual temperature condition that is the actual temperature condition in the column, or estimate the actual temperature condition itself.
More specifically, in estimating the temperature condition at the time of the isothermal analysis, the setting temperature condition is the setting sustained temperature to the oven, the correction condition is correction temperature, and the assumed temperature condition is an assumed sustained temperature calculated from the setting sustained temperature and the plurality of correction temperatures. Then, the actual temperature condition estimation part may be configured to calculate the predicted retention time at a plurality of assumed sustained temperatures of different correction temperatures, and estimate the assumed sustained temperature that is deviated from the actually measured retention time by a predetermined value or less, as the actual sustained temperature that is the actual temperature condition.
As described above, according to the present invention, the error and relationship between the various setting temperature conditions (setting sustained temperature, setting initial temperature, setting temperature rising rate, and setting temperature rising time) set to the oven in the isothermal analysis or the temperature rising analysis, and the actual temperature conditions in the column (actual sustained temperature, actual initial temperature, actual temperature rising rate, and actual temperature rising time) can be estimated, thereby estimating the actual temperature condition.
In the embodiment, the actual temperature rising rate is estimated to accurately calculate the predicted retention index to be compared with the actually measured retention index and however, the actual temperature rising rate may be estimated to compare the adjustment retention times with each other, or for other purposes. The method of estimating the actual temperature rising rate according to the present invention may be used to recognize the characteristic of the oven and the characteristic of change in the temperature in the column more accurately.
In the case where the existing data processing device for gas chromatograph is provided with a feature for estimating the actual temperature rising rate, for example, the data processing program for gas chromatograph, which is recorded in a recording medium, may be installed into the data processing device
The embodiments may be modified in a various manner and combined with each other so as not to deviate from the subject matter of the present invention.
The present invention can provide the data processing device for gas chromatograph that can indirectly estimate the actual temperature condition such as the actual temperature rising rate without directly measuring the temperature in the column, and actually measuring various values with the gas chromatograph.
Number | Date | Country | Kind |
---|---|---|---|
2013-085302 | Apr 2013 | JP | national |