CN108399148A - A kind of data processing control method, computer program based on intercept method - Google Patents
A kind of data processing control method, computer program based on intercept method Download PDFInfo
- Publication number
- CN108399148A CN108399148A CN201810087806.7A CN201810087806A CN108399148A CN 108399148 A CN108399148 A CN 108399148A CN 201810087806 A CN201810087806 A CN 201810087806A CN 108399148 A CN108399148 A CN 108399148A
- Authority
- CN
- China
- Prior art keywords
- ratio
- intercept
- data processing
- mea
- processing control
- Prior art date
- Legal status (The legal status 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 status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention belongs to isotopic geochronology studying technological domains, disclose a kind of data processing control method, computer program based on intercept method, degrade signal conic fitting, take intercept；The standard specimen initial ratio R of data calibration carries out linear fit and obtains the initial ratio r of virtual standard specimen interpolation R ' correcting samples；T is obtained according to decay equation calculation；For each measuring point, the data of corresponding isotopic ratio are taken, using least square method, conic fitting is carried out with the time, carries out error analysis.The operation interface of the present invention is beautiful, easy to operate, has higher accuracy and accuracy.The present invention can preferably serve the isotopic geochronology research work in field of earth sciences.
Description
Technical field
The invention belongs to data processing software more particularly to a kind of data processing control method based on intercept method, calculate
Machine program.
Background technology
Based on U-Th-Pb decay systems, using the equipment (LA- that laser ablation system and plasma mass spectrograph are mating
ICP-MS) in-situ micro area U-Th-Pb isotope chronology test analysis can be carried out to rich U mineral.However, due in laser pair
There are element fractionation effects (to refer to caused by one or more factors, the test of different nucleic during sample is degraded
The ratio between counting and concentration have differences), it is therefore desirable to off-line correction is carried out to test data using special technology software.Currently,
Common data processing software has the Glitter softwares of Australian Macquarie universities exploitation and China University of Geosciences (military
The Chinese) exploitation ICPMSDataCal softwares.Since element fractionation effect was not only related with sample matrices, but with laser ablation parameter and
It is related to degrade the time.Therefore, if being intended to effectively correct element fractionation effect, in addition to needing using the matrix with sample to be tested
Other than the standard sample to match is calibrated, the entire instantaneous signal for degrading the period is also detected in real time, using suitable " plan
Slightly " it is corrected in conjunction with external standard method principle.Currently used " strategy " is " method of average ", i.e., is degraded in the period with interception and instantaneously believed
Number one section average (mean value of isotope signals ratio or the ratio of isotope signals mean value), substitute into calculation formula
(1), the calibration for carrying out isotope ratio calculates.The Glitter softwares and Chinese Geological of Australian Macquarie universities exploitation are big
" strategy " for learning the ICPMSDataCal softwares use of (Wuhan) professor Liu Yongsheng exploitation is " method of average ".However, for member
Element is fractionated serious sample, and the calibration result of " method of average " will be largely dependent upon the selection of instantaneous signal section,
To can also vary with each individual, this will directly affect accuracy, accuracy and the objectivity of data.Paul J.Sylvester and
Majid Ghaderi (1997) once attempted to calibrate strategy using " intercept method " solve scheelite and man-made glass standard sample it
Between the unmatched problem of matrix, principle is will to degrade instantaneous signal in the period to carry out linear fit, with the value of " zero moment "
(i.e. intercept) replaces the counting in formula (1) or counting ratio.The uncertainty of " intercept method " calibration strategy is mainly and linear fit
Goodness it is related, not by calibration executor's subjectivity idea interfered.It is, in principle, that " intercept method " calibration strategy is in processing element
When the more serious sample of fractionating degree, effect can be than " method of average " advantageously, and result is more objective.
Ratioij=(Ratioijs*Cjs*Ci)/(Cis*Cj) (1)
Ratioij:The ratio of sample to be tested Radionuclide i, j after correction
Ratioijs:The ratio of standard sample Radionuclide i, j
Ci、Cj:The actual measurement of sample to be tested Radionuclide i, j count
Cis、Cjs:The actual measurement of standard sample Radionuclide i, j count
In conclusion problem of the existing technology is：It is existing to element fractionation effect carry out off-line correction based on " flat
The data processing software of equal method " is largely dependent upon the selection of instantaneous signal section, and to can also vary with each individual, this will be straight
Connecing influences accuracy, accuracy and the objectivity of data.It is fixed for the serious mineral of element fractionation effect (such as baddeleyite) U-Pb
Year test, satisfactory calibration result can not be then reached using existing software.Technical method that the present invention uses and establishment
Software does not depend on the selection of mass spectrograph instantaneous signal section, and " intercept method " dexterously avoids element test signal and divide at any time
The deviation of test result precision caused by evaporating is determined Nian Yineng for the serious mineral of the elements such as baddeleyite fractionation and is obtained accurately
Correction.
Invention content
The data processing control method that in view of the problems of the existing technology, the present invention provides a kind of based on intercept method,
Computer program.
It is described based on intercept method the invention is realized in this way a kind of data processing control method based on intercept method
Data processing control method includes：Signal conic fitting is degraded, intercept is taken；The standard specimen initial ratio R of data calibration is carried out
Linear fit obtains the initial ratio r of virtual standard specimen interpolation R ' correcting samples；T is obtained according to decay equation calculation；For each
Measuring point takes the data of corresponding isotopic ratio, using least square method, carries out conic fitting with the time, carries out error point
Analysis.
2, the data processing control method based on intercept method as described in claim 1, which is characterized in that signal will be degraded
Conic fitting is carried out, intercept is taken；The standard specimen initial ratio R of be useful for data calibration is carried out linear fit to obtain virtually
Standard specimen interpolation R ' is used for the initial ratio r of correcting sample, and t is obtained according to decay equation calculation；Intercept is initial time element ratio
Value.
Further, it for each measuring point, takes in 40 second data acquisition time of corresponding isotopic ratio between 5 seconds to 35 seconds
Data are carried out conic fitting with the time, are obtained following formula using least square method：
R=at2+bt+c；
Intercept c is obtained, is the initial ratio of the initial ratio Rmea and the corresponding isotope ratio of sample of the corresponding isotope ratio of standard specimen
Value rmea；Using least square method, Rmea and time are subjected to linear fit：
Rmea=at+b；
Obtain the correction factor R ' on sample corresponding time point；
Obtain rcal：
rcal=rmea*Rsta/R'；
rcal：The corresponding isotopic ratio of sample after correction；
rmea：Sample surveys corresponding isotopic ratio；
Rsta：The corresponding isotopic ratio of standard sample；
R'：Correction factor or virtual standard ratio, standard ratio interpolation；
206Pb/238U and207Pb/235U age t are obtained by following formula：
T=ln (rcal+1)/λ；
λ：235U、238The decay coefficient of U；
207Pb/206Pb ages t is obtained by following formula using iterative method：
λ1：238The decay coefficient of U；
λ2：235The decay coefficient of U.
Further, error analysis：
For206Pb/238U and207Pb/235U age t have：
Due to rcal=rmea*Rsta/ R', then have lnrcal=lnrmea+lnRsta- lnR', then：
Wherein, RstaWith σ RstaFor given value, rcal、rmeaIt is obtained by calculation with R'；Sample spot measured value rmenBy secondary
What curve matching obtained：
R=at2+bt+c；
Being converted to linear formula is：
R=T+c；
T=at2+bt；
Then have：
Wherein：
Separately have：
：
Correction factor R ' carries out linear fit with the time by actual measurement standard specimen ratio R mea and obtains, then：
It obtains206Pb/238U and207Pb/235The error at U ages；
207Pb/206The computational methods of Pb age errors are：
r:207Pb/206Pb
t:207Pb/206Pb age
u:235U/238U
λ1:238The decay coefficient of U
λ2:235The decay coefficient of U.
Another object of the present invention is to provide a kind of meters of the data processing control method based on intercept method described in realize
Calculation machine program.
Another object of the present invention is to provide a kind of information data processing terminals equipped with the computer program.
Another object of the present invention is to provide a kind of computer readable storage mediums, including instruction, when it is in computer
When upper operation so that computer executes the data processing control method based on intercept method.
The present invention is provided by the principle of " intercept method " element Fractionation Correction and the design method and partial code of software
One kind being based on " intercept method " element Fractionation Correction principle, and the data processing worked out using Visual Basic programming softwares is soft
Part.Beautiful interface, it is easy to operate, compare existing software, efficiency can be improved ten times, this by substantially reduce researcher it is valuable when
Between cost.Meanwhile serious mineral are fractionated for elements such as processing baddeleyites and determine annual data with essence more higher than existing software
Exactness and accuracy, due to baddeleyite be it is optimal in mafic, ultrabasic rock determine year mineral, the present invention is in addition in tradition
Rich U, Th mineral, which determine year field, to be had other than broad prospect of application, determines the year larger geologic body of difficulty in mafic, ultrabasic rock etc.
It will also be played an irreplaceable role in Chronologic Study.
Description of the drawings
Fig. 1 is the data processing control method flow chart provided in an embodiment of the present invention based on intercept method.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
The application principle of the present invention is explained in detail below in conjunction with the accompanying drawings.
As shown in Figure 1, the data processing control method provided in an embodiment of the present invention based on intercept method includes the following steps：
S101：Signal conic fitting is degraded, intercept is taken；The standard specimen initial ratio R of be useful for data calibration is carried out
Linear fit obtains the initial ratio r of virtual standard specimen interpolation R ' correcting samples, and t is obtained according to decay equation calculation；
S102：For each measuring point, the data of corresponding isotopic ratio are taken, using least square method, two are carried out with the time
Secondary curve matching；
S103：Carry out error analysis.
With reference to concrete analysis, the invention will be further described.
The preparation method of data processing software provided in an embodiment of the present invention based on intercept method is based on " intercept method " element
The principle of Fractionation Correction, calculating process and error propagation strategy are as follows：
1) correction principle：Based on degrading, start time is substantially not present element fractionation and element is fractionated approximation and is at any time
The fact that conic section, will degrade signal and carry out conic fitting, takes intercept, i.e. initial time element ratios.To own again
Standard specimen initial ratio R for data calibration carries out linear fit and obtains virtual standard specimen interpolation R ', is used for the initial of correcting sample
Ratio r, and then t is obtained according to decay equation calculation.
2) calculating process：
For each measuring point, the data between 5 seconds to 35 seconds in 40 second data acquisition time of corresponding isotopic ratio are taken, are adopted
With least square method, conic fitting is carried out with the time, obtains following formula：
R=at2+bt+c；
Obtain intercept c, the i.e. initial ratio of the initial ratio Rmea of the corresponding isotope ratio of standard specimen and the corresponding isotope ratio of sample
Value rmea.Using least square method, Rmea and time are subjected to linear fit:
Rmea=at+b；
Obtain correction factor (or virtual standard ratio, standard ratio interpolation) R ' on sample corresponding time point.According to following formula
Obtain rcal：
rcal=rmea*Rsta/R'
rcal：The corresponding isotopic ratio of sample after correction
rmea：Sample surveys corresponding isotopic ratio
Rsta：The corresponding isotopic ratio of standard sample
R'：Correction factor (or virtual standard ratio, standard ratio interpolation)
206Pb/238U and207Pb/235U age t are obtained by following formula：
T=ln (rcal+1)/λ
λ：235U、238The decay coefficient of U
207Pb/206Pb ages t can be obtained by following formula using iterative method：
λ1：238The decay coefficient of U
λ2：235The decay coefficient of U
3) error analysis：
For generic function y=f (x), haveSimilarly, for206Pb/238U and207Pb/235U ages t is then
Have：
Due to rcaL=rmea*Rsta/ R', then have lnrcal=lnrmea+lnRsta- lnR', then：
Wherein, RstaWith σ RstaFor given value, rcal、rmeaIt can be obtained by calculation with R', computational methods are shown in that " (2) calculate
Process ", then need to only obtain σ rmenWith σ R'.
Sample spot measured value rmenIt is obtained by conic fitting：
R=at2+bt+c
Converting it into linear formula is：
R=T+c
T=at2+bt
Then have：
Wherein：
Separately have：
(3) formula of substitution obtains：
In addition, being obtained since correction factor R ' carries out linear fit by actual measurement standard specimen ratio R mea with the time, then：
(5), (6) formula are substituted into (3) formula, then (3) formula substitution (2) formula can be obtained206Pb/238U and207Pb/235The U ages
Error.
207Pb/206The computational methods of Pb age errors are：
r:207Pb/206Pb
t:207Pb/206Pb age
u:235U/238U
λ1:238The decay coefficient of U
λ2:235The decay coefficient of U
The application principle of the present invention is further described with reference to specific embodiment.
Embodiment 1：
In order to realize the efficient calculating and processing of Elemental analysis data, while interface is beautiful as possible, and operation is simple as possible
Just, software have the function of following menu bar, window and：
1) File menu, containing functions such as importing, export data, preservations
The data of importing are EXCEL, and suffix is the data file of .csv formats (Excel file).No less than 300 can be imported
A such file.
Derived file is EXCEL file format.
Preservation refers to the operation and setting for preserving user to standard specimen, data-signal.
2) actions menu, containing functions such as standard specimen setting, signal behaviors.
Outer target standard value may be present in .txt texts, and can be read by software.External standard in one group of test data is answered
It can be automatically recognized, internal standard can set input field and be manually entered.Standard specimen is corrected calculating after setting.
There is the line chart of data-signal in .csv files in signal behavior window, and selected containing sample selection, element (ratio),
The functions such as the selection of signal section, filter peak.
It is realized according to data-oriented regression criterion limit value at filter peak.
Signal behavior window can be opened multiple.
3) result is at figure
Age result is respectively displayed in positive and negative harmonious figure, and age range can be arranged；Rare earth element content results are aobvious
Show in partition of rare earth element figure there are more set a reference values available.
Portion of software code is extracted：
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
All any modification, equivalent and improvement etc., should all be included in the protection scope of the present invention made by within refreshing and principle.
Claims (7)
1. a kind of data processing control method based on intercept method, which is characterized in that the data processing control based on intercept method
Method processed includes：Signal conic fitting is degraded, intercept is taken；The standard specimen initial ratio R of data calibration carries out linear fit and obtains
To the initial ratio r of virtual standard specimen interpolation R ' correcting samples；T is obtained according to decay equation calculation；For each measuring point, take corresponding
The data of isotopic ratio are carried out conic fitting with the time, are carried out error analysis using least square method.
2. the data processing control method based on intercept method as described in claim 1, which is characterized in that signal progress will be degraded
Conic fitting takes intercept；The standard specimen initial ratio R of be useful for data calibration is subjected to linear fit and obtains virtual standard specimen
Interpolation R ' is used for the initial ratio r of correcting sample, and t is obtained according to decay equation calculation；Intercept is initial time element ratios.
3. the data processing control method based on intercept method as described in claim 1, which is characterized in that for each measuring point,
Take the data between 5 seconds to 35 seconds in 40 second data acquisition time of corresponding isotopic ratio, using least square method, with the time into
Row conic fitting, obtains following formula：
R=at2+bt+c；
Intercept c is obtained, is the initial ratio of the initial ratio Rmea and the corresponding isotope ratio of sample of the corresponding isotope ratio of standard specimen
rmea；Using least square method, Rmea and time are subjected to linear fit：
Rmea=at+b；
Obtain the correction factor R ' on sample corresponding time point；
Obtain rcal：
rcal=rmea*Rsta/R'；
rcal：The corresponding isotopic ratio of sample after correction；
rmea：Sample surveys corresponding isotopic ratio；
Rsta：The corresponding isotopic ratio of standard sample；
R'：Correction factor or virtual standard ratio, standard ratio interpolation；
206Pb/238U and207Pb/235U age t are obtained by following formula：
T=ln (rcal+1)/λ；
λ：235U、238The decay coefficient of U；
207Pb/206Pb ages t is obtained by following formula using iterative method：
λ1：238The decay coefficient of U；
λ2：235The decay coefficient of U.
4. the data processing control method based on intercept method as described in claim 1, which is characterized in that error analysis：
For206Pb/238U and207Pb/235U age t have：
Due to rcal=rmea*Rsta/ R' then has ln rcal=ln rmea+ln Rsta- ln R', then：
Wherein, RstaWith σ RstaFor given value, rcal、rmeaIt is obtained by calculation with R'；
Sample spot measured value rmenIt is obtained by conic fitting：
R=at2+bt+c；
Being converted to linear formula is：
R=T+c；
T=at2+bt；
Then have：
Wherein：
Separately have：
：
Correction factor R ' carries out linear fit with the time by actual measurement standard specimen ratio R mea and obtains, then：
It obtains206Pb/238U and207Pb/235The error at U ages；
207Pb/206The computational methods of Pb age errors are：
r:207Pb/206Pb
t:207Pb/206Pb age
u:235U/238U
λ1:238The decay coefficient of U
λ2:235The decay coefficient of U.
5. a kind of computer journey for realizing the data processing control method based on intercept method described in Claims 1 to 4 any one
Sequence.
6. a kind of information data processing terminal equipped with computer program described in claim 5.
7. a kind of computer readable storage medium, including instruction, when run on a computer so that computer is executed as weighed
Profit requires the data processing control method based on intercept method described in 1~4 any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810087806.7A CN108399148A (en) | 2018-01-30 | 2018-01-30 | A kind of data processing control method, computer program based on intercept method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810087806.7A CN108399148A (en) | 2018-01-30 | 2018-01-30 | A kind of data processing control method, computer program based on intercept method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108399148A true CN108399148A (en) | 2018-08-14 |
Family
ID=63095764
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810087806.7A Pending CN108399148A (en) | 2018-01-30 | 2018-01-30 | A kind of data processing control method, computer program based on intercept method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108399148A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105510427A (en) * | 2014-09-25 | 2016-04-20 | 中国科学院地球化学研究所 | Numerical treatment method for multi-element isotope double spike technology |
| CN106483189A (en) * | 2016-09-26 | 2017-03-08 | 中国地质调查局西安地质调查中心 | A kind of baddeleyite U Pb based on LA ICP (Q) MS determines year method |
-
2018
- 2018-01-30 CN CN201810087806.7A patent/CN108399148A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105510427A (en) * | 2014-09-25 | 2016-04-20 | 中国科学院地球化学研究所 | Numerical treatment method for multi-element isotope double spike technology |
| CN106483189A (en) * | 2016-09-26 | 2017-03-08 | 中国地质调查局西安地质调查中心 | A kind of baddeleyite U Pb based on LA ICP (Q) MS determines year method |
Non-Patent Citations (2)
| Title |
|---|
| CARL SPANDLERA 等: "《A new titanite standard for in situ analysis of Sm–Nd isotopes and U–Pb geochronology》", 《CHEMICAL GEOLOGY》 * |
| 李艳广: "《斜锆石LA-ICP-MS U-Pb定年方法及应用》", 《地质学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bernasconi et al. | InterCarb: A community effort to improve interlaboratory standardization of the carbonate clumped isotope thermometer using carbonate standards | |
| Sundell et al. | Rapid U‐Pb geochronology by laser ablation multi‐collector ICP‐MS | |
| Cheng et al. | Improvements in 230Th dating, 230Th and 234U half-life values, and U–Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry | |
| Halekas | Seasonal variability of the hydrogen exosphere of Mars | |
| Vervoort et al. | Isotopic composition of Yb and the determination of Lu concentrations and Lu/Hf ratios by isotope dilution using MC‐ICPMS | |
| McLean et al. | An algorithm for U‐Pb isotope dilution data reduction and uncertainty propagation | |
| Rudge et al. | The double spike toolbox | |
| Gehrels et al. | Detrital zircon geochronology by laser-ablation multicollector ICPMS at the Arizona LaserChron Center | |
| Uemura et al. | Ranges of moisture-source temperature estimated from Antarctic ice cores stable isotope records over glacial–interglacial cycles | |
| Bracciali et al. | UPb LA-(MC)-ICP-MS dating of rutile: New reference materials and applications to sedimentary provenance | |
| US10768159B2 (en) | Method for tracing wine origin based on multielements and stable isotopes | |
| John | Optimizing sample and spike concentrations for isotopic analysis by double-spike ICPMS | |
| CN102141539A (en) | Method for uranium lead dating of baddeleyite by using secondary ion mass spectroscopy | |
| Coble et al. | Calibration of Nu-Instruments Noblesse multicollector mass spectrometers for argon isotopic measurements using a newly developed reference gas | |
| Harmon et al. | Limits of U-series analysis: Phase 1 results of the Uranium-Series Intercomparison Project | |
| CN111141812A (en) | U-Pb dating data correcting method, storage medium and correcting equipment for sphene | |
| Jenkins et al. | A comprehensive global oceanic dataset of helium isotope and tritium measurements | |
| Stauffer et al. | A post‐2013 dropoff in total ozone at a third of global ozonesonde stations: Electrochemical concentration cell instrument artifacts? | |
| Fowler et al. | In‐situ measurements of ion density in the Martian ionosphere: Underlying structure and variability observed by the MAVEN‐STATIC instrument | |
| CN114168906A (en) | Mapping geographic information data acquisition system based on cloud computing | |
| Burger et al. | KiDS-1000 cosmology: Constraints from density split statistics | |
| CN106442599A (en) | Method and device for determining oil content of rocks | |
| CN105717190A (en) | Value determining method for arsenic content in tobacco quality control sample | |
| Pourmand et al. | Introducing a comprehensive data reduction and uncertainty propagation algorithm for U‐Th geochronometry with extraction chromatography and isotope dilution MC‐ICP‐MS | |
| Barry et al. | High precision nitrogen isotope measurements in oceanic basalts using a static triple collection noble gas mass spectrometer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180814 |