CN106290231B - Nitrate esters explosive wastewater powder charge safe storage life predictor method - Google Patents
Nitrate esters explosive wastewater powder charge safe storage life predictor method Download PDFInfo
- Publication number
- CN106290231B CN106290231B CN201610614112.5A CN201610614112A CN106290231B CN 106290231 B CN106290231 B CN 106290231B CN 201610614112 A CN201610614112 A CN 201610614112A CN 106290231 B CN106290231 B CN 106290231B
- Authority
- CN
- China
- Prior art keywords
- sample
- temperature
- infrared
- nitrate esters
- storage life
- 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.)
- Active
Links
- 239000000843 powder Substances 0.000 title claims abstract description 105
- 238000003860 storage Methods 0.000 title claims abstract description 96
- 239000002360 explosive Substances 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 76
- 150000002823 nitrates Chemical class 0.000 title claims abstract description 72
- 239000002351 wastewater Substances 0.000 title claims abstract description 57
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 81
- 230000032683 aging Effects 0.000 claims abstract description 80
- 230000001133 acceleration Effects 0.000 claims abstract description 70
- 230000006641 stabilisation Effects 0.000 claims abstract description 46
- 238000011105 stabilization Methods 0.000 claims abstract description 46
- 230000015556 catabolic process Effects 0.000 claims abstract description 31
- 238000006731 degradation reaction Methods 0.000 claims abstract description 31
- 238000005086 pumping Methods 0.000 claims abstract description 22
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000003333 near-infrared imaging Methods 0.000 claims abstract description 9
- 238000010521 absorption reaction Methods 0.000 claims abstract description 3
- 239000003381 stabilizer Substances 0.000 claims description 73
- 239000003380 propellant Substances 0.000 claims description 52
- 239000000126 substance Substances 0.000 claims description 42
- 239000010410 layer Substances 0.000 claims description 35
- 238000002329 infrared spectrum Methods 0.000 claims description 30
- 238000004448 titration Methods 0.000 claims description 24
- 238000010606 normalization Methods 0.000 claims description 15
- 238000009795 derivation Methods 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 10
- 239000011229 interlayer Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000011068 loading method Methods 0.000 claims description 7
- 239000003814 drug Substances 0.000 claims description 6
- 238000013213 extrapolation Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000004128 high performance liquid chromatography Methods 0.000 claims description 5
- LNTHITQWFMADLM-UHFFFAOYSA-N gallic acid Chemical compound OC(=O)C1=CC(O)=C(O)C(O)=C1 LNTHITQWFMADLM-UHFFFAOYSA-N 0.000 claims description 4
- 238000003556 assay Methods 0.000 claims description 3
- 238000004587 chromatography analysis Methods 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 229940079593 drug Drugs 0.000 claims description 2
- 150000002148 esters Chemical class 0.000 claims description 2
- 229940074391 gallic acid Drugs 0.000 claims description 2
- 235000004515 gallic acid Nutrition 0.000 claims description 2
- 239000007791 liquid phase Substances 0.000 claims 1
- 230000009466 transformation Effects 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 104
- PZIMIYVOZBTARW-UHFFFAOYSA-N centralite Chemical compound C=1C=CC=CC=1N(CC)C(=O)N(CC)C1=CC=CC=C1 PZIMIYVOZBTARW-UHFFFAOYSA-N 0.000 description 29
- 229910002651 NO3 Inorganic materials 0.000 description 18
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 description 18
- 101100133184 Caenorhabditis elegans nep-1 gene Proteins 0.000 description 16
- 150000001875 compounds Chemical class 0.000 description 16
- AAOVKJBEBIDNHE-UHFFFAOYSA-N diazepam Chemical compound N=1CC(=O)N(C)C2=CC=C(Cl)C=C2C=1C1=CC=CC=C1 AAOVKJBEBIDNHE-UHFFFAOYSA-N 0.000 description 16
- 229960003529 diazepam Drugs 0.000 description 16
- 102100038782 Carbohydrate sulfotransferase 1 Human genes 0.000 description 14
- 101710158459 Carbohydrate sulfotransferase 1 Proteins 0.000 description 14
- 108010060966 silk gland factor-1 Proteins 0.000 description 14
- 239000002131 composite material Substances 0.000 description 12
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 10
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 9
- 230000003595 spectral effect Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 8
- 239000007789 gas Substances 0.000 description 8
- 238000012937 correction Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 239000012086 standard solution Substances 0.000 description 7
- 238000002512 chemotherapy Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- AKHNMLFCWUSKQB-UHFFFAOYSA-L sodium thiosulfate Chemical compound [Na+].[Na+].[O-]S([O-])(=O)=S AKHNMLFCWUSKQB-UHFFFAOYSA-L 0.000 description 5
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000010561 standard procedure Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000005844 autocatalytic reaction Methods 0.000 description 3
- 230000001066 destructive effect Effects 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000010238 partial least squares regression Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 235000019345 sodium thiosulphate Nutrition 0.000 description 3
- XWNSFEAWWGGSKJ-UHFFFAOYSA-N 4-acetyl-4-methylheptanedinitrile Chemical compound N#CCCC(C)(C(=O)C)CCC#N XWNSFEAWWGGSKJ-UHFFFAOYSA-N 0.000 description 2
- CPELXLSAUQHCOX-UHFFFAOYSA-M Bromide Chemical compound [Br-] CPELXLSAUQHCOX-UHFFFAOYSA-M 0.000 description 2
- VEXZGXHMUGYJMC-UHFFFAOYSA-N Hydrochloric acid Chemical compound Cl VEXZGXHMUGYJMC-UHFFFAOYSA-N 0.000 description 2
- SNIOPGDIGTZGOP-UHFFFAOYSA-N Nitroglycerin Chemical compound [O-][N+](=O)OCC(O[N+]([O-])=O)CO[N+]([O-])=O SNIOPGDIGTZGOP-UHFFFAOYSA-N 0.000 description 2
- 239000004153 Potassium bromate Substances 0.000 description 2
- QAOWNCQODCNURD-UHFFFAOYSA-N Sulfuric acid Chemical compound OS(O)(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-N 0.000 description 2
- 238000009835 boiling Methods 0.000 description 2
- 230000031709 bromination Effects 0.000 description 2
- 238000005893 bromination reaction Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 235000019441 ethanol Nutrition 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 229960003711 glyceryl trinitrate Drugs 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 238000004811 liquid chromatography Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 229920001220 nitrocellulos Polymers 0.000 description 2
- 229940094037 potassium bromate Drugs 0.000 description 2
- 235000019396 potassium bromate Nutrition 0.000 description 2
- POCJOGNVFHPZNS-ZJUUUORDSA-N (6S,7R)-2-azaspiro[5.5]undecan-7-ol Chemical class O[C@@H]1CCCC[C@]11CNCCC1 POCJOGNVFHPZNS-ZJUUUORDSA-N 0.000 description 1
- WKBOTKDWSSQWDR-UHFFFAOYSA-N Bromine atom Chemical compound [Br] WKBOTKDWSSQWDR-UHFFFAOYSA-N 0.000 description 1
- 241001411320 Eriogonum inflatum Species 0.000 description 1
- 238000004566 IR spectroscopy Methods 0.000 description 1
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 1
- 229920002472 Starch Polymers 0.000 description 1
- 239000012496 blank sample Substances 0.000 description 1
- GDTBXPJZTBHREO-UHFFFAOYSA-N bromine Substances BrBr GDTBXPJZTBHREO-UHFFFAOYSA-N 0.000 description 1
- 229910052794 bromium Inorganic materials 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical class C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 239000012259 ether extract Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- XMBWDFGMSWQBCA-UHFFFAOYSA-N hydrogen iodide Chemical compound I XMBWDFGMSWQBCA-UHFFFAOYSA-N 0.000 description 1
- 238000004643 material aging Methods 0.000 description 1
- DOTMOQHOJINYBL-UHFFFAOYSA-N molecular nitrogen;molecular oxygen Chemical compound N#N.O=O DOTMOQHOJINYBL-UHFFFAOYSA-N 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000036632 reaction speed Effects 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000012488 sample solution Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000011734 sodium Substances 0.000 description 1
- 229910052708 sodium Inorganic materials 0.000 description 1
- PANBYUAFMMOFOV-UHFFFAOYSA-N sodium;sulfuric acid Chemical compound [Na].OS(O)(=O)=O PANBYUAFMMOFOV-UHFFFAOYSA-N 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 239000008107 starch Substances 0.000 description 1
- 125000000446 sulfanediyl group Chemical group *S* 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
- G01N31/16—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using titration
Abstract
The present invention relates to a kind of nitrate esters explosive wastewater powder charge safe storage life predictor methods, prepare stabilization agent effective content near-infrared quantitative model modeling sample using thermal acceleration aging process, establish stabilization agent effective content near-infrared quantitative model；Stabilization agent effective content transformation temperature coefficient is obtained by thermal acceleration degradation, establishes safe storage life predicting equation；Critical quantity is consumed as critical pumping rate using stabilization agent effective content, the single temperature thermal acceleration ageing process stabilization agent effective content of Near-Infrared Absorption Method tracing detection powder charge, the crash time that effective content reaches critical pumping rate consumption is obtained, safe storage life at a temperature of nitrate esters explosive wastewater powder charge storage environment is estimated by safe storage life predicting equation.This method, which only needs near-infrared method to track a powder charge sample difference ageing time stabilization agent effective content at a temperature, can obtain safe storage life, and sample size reduces 96%, safety, low cost.
Description
Technical field
The invention belongs to explosive wastewater Life Assessment Technology fields, and it is pre- to be related to a kind of explosive wastewater storage process safe storage life
Estimate method.Especially it is a kind of using near-infrared diffuse spectrometry detection the single temperature thermal acceleration aging sample of nitrate esters explosive wastewater
Product effective stabilizer content obtains the crash time that effective stabilizer content reaches critical pumping rate consumption, is store safely according to powder charge
Life prediction equation is deposited, the method for estimating nitrate esters explosive wastewater safe storage life at a temperature of storage environment.
Background technique
Explosive wastewater safe storage life is under storage requirement, and the pervious storage of autocatalysis accelerated decomposition does not occur for explosive wastewater
Deposit the time.
Nitrocotton, nitroglycerine are the primary energy content in nitrate esters explosive wastewater, nitrocotton and nitroglycerine molecule
All contain-C-ONO in structure2Structure, ammonal contain-N-NO2Structure, nitrate esters fire is fried in heated or storage process
Medicine each component occurs slowly to decompose, O-NO2With-N-NO2Key fracture releases energy, while generating the nitrogen oxygen with autocatalysis
Compound, nitrogen oxides are further catalyzed energy ingredient decomposition, continue slow heat release and generate heat accumulation, reduce chemical stability,
Lead to burning, explosion, influences storing stability.Therefore, it needs to be added stabilization agent in nitrate esters explosive wastewater as stabilization agent, uses
The nitrogen oxide that nitrate esters explosive wastewater decomposes to give off is absorbed, so that nitrogen oxides be inhibited to decompose nitrate esters explosive wastewaters
Autocatalysis, while generating still there is stable effect stabilization agent derivative, delay energy ingredient to decompose, raising nitrate
Class explosive wastewater safe storage life.Stabilization agent and its derivative total content are known as effective stabilizer content or stabilization agent effectively contains
It measures, the time undergone under storage requirement when effective stabilizer content drops to critical pumping rate is that nitrate esters explosive wastewater is store safely
Deposit the service life.
Using thermal acceleration aging process tracking effective stabilizer content, (GJB770B-2005 method 506.1 estimated storage in the past
Service life thermal acceleration aging process) nitrate esters explosive wastewater safe storage life is estimated, commonly using chemical titration, (bromine is reacted with stabilization agent
Generate bromide chemical method) or chromatography determination stabilization agent effective content be damage detection method, using aether backflow 24
The hour even longer time extracts and crushes sample stabilization agent effective content, and low boiling point solvent ether is inflammable and explosive, passes through machinery
Safety is poor when crushing nitrate esters explosive wastewater aging sample preparation extracted by ether sample, cumbersome, time-consuming and chemical drop
The bromide for determining method generation brings pollution.
It is regular by nitrate esters explosive wastewater material aging sample effective stabilizer changes of contents under the different temperatures of acquisition,
When estimating safe storage life with Bei Seluote (Berthelot) equation, due to thermal acceleration degradation at 50 DEG C~95 DEG C extremely
4 temperature spots are used less, and each temperature spot sampling number is at least 6 times, and the aged samples amount needed is big, duration of test runs 6
It is a month or more, extremely time-consuming, and long-time high temperature thermal acceleration aging and pass through mechanical cutting nitrate esters explosive wastewater aging sample
Risk is high when preparing extraction sample.
Hot product can be generated for the nitrate esters explosive wastewater powder charge thermal acceleration degradation thermal histories with certain size
It is tired, on the one hand lead to thermal field uneven distribution, cause stabilization agent content distribution uneven, another aspect heap, nitrate esters fire is fried
Easily because of heat accumulation burning and explosion accident occur for medicine powder charge multi-temperature point degradation, accordingly, it is difficult to pass through the multiple temperature of powder charge
Degree point degradation simultaneously obtains stable agent content using chemical titration or chromatography, carry out newly formula long-life powder charge design and
It estimates safe storage life in storage process to estimate, it is therefore desirable to use new method.
Near-infrared spectrum technique be frequency multiplication, sum of fundamental frequencies information based on X-H (X-C, N, O) in organic molecule to substance into
The spectral technique of row qualitative and quantitative analysis has the advantages that convenient, lossless, green.The nitrate esters fire established according to this method
Explosive loading safe storage life predicting equation, it is effective using stabilization agent in near-infrared method detection powder charge thermal acceleration ageing process
Content method, can safe, time saving, inexpensive, environmentally friendly acquisition powder charge safe storage life.
Summary of the invention
For defect or deficiency existing for above-mentioned existing nitrate esters explosive wastewater safe storage life technology, mesh of the invention
Be, it is pre- to provide a kind of nitrate esters explosive wastewater powder charge safe storage life based near infrared detection effective stabilizer content
Estimate method, this method can safe, time saving, inexpensive, environmentally friendly acquisition nitrate esters explosive wastewater powder charge safe storage life.
In order to realize above-mentioned task, the invention adopts the following technical scheme:
A kind of nitrate esters explosive wastewater powder charge safe storage life predictor method, which is characterized in that this method is added using heat
Fast aging process prepares nitrate esters explosive wastewater effective stabilizer content near-infrared quantitative model modeling sample, using chemical titration
Or high performance liquid chromatography obtain effective stabilizer content chemical score, using near-infrared diffuse spectrometer acquisition spectrum, pass through
Chemical Measurement establishes effective stabilizer content near-infrared quantitative model；Critical quantity is consumed as critical using effective stabilizer content
Criterion carries out powder charge thermal acceleration degradation, obtains temperature coefficient by effective stabilizer changes of contents rule, establishes powder charge peace
Full storage life predicting equation；Sample is effectively stable in the single temperature thermal acceleration ageing process of the lossless tracing detection of near-infrared method
Agent content obtains the crash time that effective stabilizer content reaches critical pumping rate consumption, passes through powder charge safe life predicting equation
The safe storage life of nitrate esters explosive wastewater powder charge at a temperature of acquisition storage environment.Specifically follow these steps to carry out:
(1) thermal acceleration aging process prepares near-infrared method detection effective stabilizer content modeling sample
Nitrate esters explosive wastewater material be fitted into capillary exhaust pipe ground decrement bottle in (semi-enclosed condition) under certain temperature
Thermal acceleration aging, interval time sampling, prepares the nitrate esters explosive wastewater sample of different effective stabilizer contents, is divided into internal school
Positive collection sample and external certificate collection sample.
(2) effective stabilizer near-infrared quantitative model
1. obtaining nitrate esters explosive wastewater sample effective stabilizer content using chemical titration or high performance liquid chromatography
Chemical score determines content range.
2. using near infrared spectrometer to internal calibration set sample and external certificate collection sample.Acquire near infrared spectrum.It adopts
Sample loading mode is diffusing reflection, and sampling wave band is 4000cm-1~12800cm-1, optimize spectral scan parameter, determine optimum resolution,
Scanning times and sample duplicate measurements number.Collected near infrared spectrum is transmitted to computer.
3. carrying out pretreatment and Linear Quasi to internal calibrations collection sample near infrared spectrum is divided into using chemometrics method
It closes.Key band is selected, in normalization, baseline is smooth, first derivation, second order derivation, multiplicative scatter correction Pretreated spectra side
Selected in method alone or in combination mode acquisition spectrum is pre-processed, then returned using Partial Least Squares, that is, PLS method
Return fitting, establish calibration model, external certificate is carried out to model built with external certificate collection sample, the internal inspection according to model
The appraisement system combined with external certificate, successive optimization model.
According to the effective stabilizer content near-infrared quantitative model of foundation, using near infrared spectrometer to selected sample into
The worst error of row effective stabilizer assay, acquisition near-infrared predicted value, chemical score and near-infrared predicted value should be less than marking
Quasi- method repeatability error.
(3) safe storage life predicting equation
Nitrate esters and nitramine class energetic material decompose the change for causing nitrate esters explosive wastewater effective stabilizer changes of contents
Reaction is learned, the Arrhenius Equation is obeyed, indicates the accurate empirical formula of k-T relationship are as follows:
The temperature coefficient r for defining reaction speed is the every variation multiplying power for rising 10 DEG C of rate constants of temperature:
R=kT+10/kT (1)
In formula: r indicates that the temperature difference is 10 DEG C of reaction rate temperature coefficient；
kTRate constant at T DEG C of expression；
kT+10Indicate rate constant when (T+10) DEG C.
Due in T1And T2At two temperature, same reaction, aA+bB+ ... → lL+ ..., all from same initial dense are carried out respectively
Degree proceeds to identical conversion ratio, and required time is respectively t1And t2, k1With k2Respectively T1And T2When rate constant, then have
t2/t1=k1/k2 (2)
If: Tn-Tn-1=10 DEG C
Powder charge safe storage life predicting equation can be obtained by formula (1) (2):
In formula: r indicates that the temperature difference is 10 DEG C of reaction rate temperature coefficient；
t0Indicate storage environment temperature safety storage life, d；
tnIndicate high temperature thermal acceleration ageing time, d；
TnIndicate high temperature thermal acceleration degradation temperature, DEG C；
T0Indicate storage environment temperature, DEG C.
(4) temperature coefficient r
Logarithm is taken to formula (3) both sides, is obtained after arrangement:
Tn=A+Blgtn (4)
In formula: A, B indicates coefficient, A=T0+10×lgt0/ lgr, B=-10/lgr
Temperature coefficient r=10-10/B, sample different temperatures T is tracked using near-infrared methodniLower thermal acceleration degradation is not
With the effective stabilizer content of ageing time, effective stabilizer content-time changing curve is obtained, is disappeared with effective stabilizer content
Consumption 50% is used as critical pumping rate, obtains different temperatures TniCorresponding crash time tni, substituted into (4) formula progress unitary and linearly returned
Return, coefficient B is obtained using least square method, and then temperature coefficient r is calculated.
(5) thermal acceleration degradation and life prediction
By nitrate esters explosive wastewater powder charge, T at Yu Danyi temperaturenLower thermal acceleration degradation, near infrared spectrum scanning is not
Central area and grain surface with ageing time lower section powder column superposed surfaces consume most fast region with stabilization agent effective content and make
Temperature T is obtained using the position Guan Chong stabilization agent effective content consumption 50% as critical pumping rate for the position Guan ChongnLower stabilization agent disappears
The crash time t of consumptionn, according to the powder charge safe storage life at a temperature of formula (3) extrapolation storage environment.
Nitrate esters explosive wastewater safe storage life near-infrared rapid Estimation method of the invention, bring Advantageous effect
Fruit is embodied in the following aspects:
1. during stabilization agent plays stable effect in nitrate esters explosive wastewater, not only stabilization agent itself has stable effect
Fruit, stabilization agent derivative still have stable effect, are the important component of effective stabilizer content, therefore, nitrate esters
Stable agent content in explosive wastewater ageing process is effective stabilizer content.Since the derivative species of stabilization agent are more, structure is complicated,
It is difficult to obtain, the method for obtaining near-infrared effective stabilizer content modeling sample by degradation overcomes a variety of stabilization agents
Derivative preparation and proportion problem, have avoided the high nitrate esters explosive wastewater manufacturing process of complicated risk, and quick, convenient,
Low cost, environmental protection.
2. the near-infrared effective stabilizer content modeling sample of thermal acceleration aging process preparation is obtained by control sample time
Different effective stabilizer contents, reach and are evenly distributed in variation range, solve technique manufacture bring low content sample point
The non-uniform defect of cloth.
3. chemical titration and liquid chromatography are that stable agent content damages detection method, material can only be stabilized
The detection of agent effective content, can not detect powder charge stabilization agent effective content.Near infrared light has very strong penetration capacity,
In test sample, the scanning of near infrared spectrum is directly done to sample, this non-destructive testing technology can without destroying powder charge
Directly measure each position stabilization agent effective content of powder charge.In addition, progress stabilization agent effectively contains after near-infrared stabilizes agent content modeling
When amount detection, testing result can be respectively provided in 2~5min, has avoided effective stabilizer in chemical titration and liquid chromatography
Sample needs the inflammable process of refluxing extraction for 24 hours in low boiling point ether solvent when separation, solves that time-consuming, low efficiency and dirt
The problems such as dye, easy to operate, low consumption, environmental protection, safety, quickly.
4. when the previous temperature coefficient for obtaining typical nitrate esters explosive wastewater, carrying out destructive inspection using chemical titration
It surveys, different temperatures difference ageing time typical case nitrate esters explosive wastewater effective stabilizer content is tracked, to obtain at least four temperature
The corresponding relationship of lower 6 ageing times and effective stabilizer content, and then be fitted and obtain temperature coefficient, it is seen then that using chemistry drop
Determine method and obtain temperature coefficient, at least needs 25 or more sample measurement effective stabilizer contents.Near-infrared effective stabilizer content
After modeling, when tracking 4 temperature thermal acceleration aging sample stabilization agent changes of contents, only 4 samples is needed to carry out non-damaged data,
Chemical titration is avoided and has needed repeatedly to take out the destructive detection of aging sample progress, sample size reduces 84%, while avoiding
Repeatedly (each temperature spot at least 6 times) samples the risky operation process that is mechanically pulverized, more quickly, convenient, realizes essence
Cost is significantly reduced while safe.
5. the near-infrared tracking single temperature speed degradation difference ageing time charge center region of nitrate explosive wastewater has
Stable agent content is imitated, using the consumption 50% of effective stabilizer content as critical pumping rate, using predicting equationPowder charge safe storage life is obtained, the nitrate esters explosive wastewater powder charge safe storage longevity is established
Order predictor method.It is this using single temperature thermal acceleration aging and lossless to simple sample progress effective stabilizer content near-infrared
The method of detection, compared with GJB770B-2005 method 506.1 estimates storage life thermal acceleration aging process, without carrying out multiple temperature
Spend thermal acceleration aging and repeatedly sampling carry out the destructive detection of stabilization agent effective content titration, it is only necessary to near-infrared method with
A powder charge sample difference ageing time stabilization agent effective content at one temperature of track can obtain safe storage life, sample
Amount reduces 96%, and the safe storage life that safer, inexpensive, timesaving realizes nitrate explosive wastewater powder charge is estimated.
Detailed description of the invention
Fig. 1 capillary exhaust pipe ground is reduced bottle .1- capillary exhaust pipe, and 2- ground is reduced bottle cap, and 3- ground is reduced bottle body
Fig. 2 capillary exhaust pipe is reduced case .4- capillary exhaust pipe, and 5- is reduced case lid, and 6- is reduced cabinet
Fig. 3 nitrate esters propellant and nitric-ester-compound explosive herbal cake sandwich style interlayer structure sample .7- upper layer tablet 8-
Sandwich of layers 9- lower layer tablet
Fig. 4 nitrate esters propellant and the upper layer nitric-ester-compound explosive powder charge two-part Grain structure sample .10- powder column, 11-
Central area, 12- lower layer powder column superposed surfaces, 13- lower layer powder column
Fig. 5 triple-base propellant SGF-1 internal calibrations collection sample near infrared spectrum
Fig. 6 triple-base propellant SGF-1 different temperatures effective stabilizer content changes over time relation curve
Fig. 7 composite modified double-base propellant GST-1 internal calibrations collection sample atlas of near infrared spectra
Fig. 8 composite modified double-base propellant GST-1 different temperatures effective stabilizer content changes over time relation curve
Fig. 9 nitric-ester-compound explosive NEP-1 internal calibrations collection sample atlas of near infrared spectra
Figure 10 nitric-ester-compound explosive NEP-1 different temperatures effective stabilizer content changes over time relation curve
Specific embodiment
Nitrate esters explosive wastewater powder charge safe storage life predictor method of the invention, which is characterized in that this method utilizes
Thermal acceleration aging process prepares effective stabilizer content near-infrared modeling sample, is obtained using chemical titration or high performance liquid chromatography
Stabilization agent effective content chemical score is taken, is diffused spectrometer acquisition spectrum, is established by Chemical Measurement effective using near-infrared
Stable agent content near-infrared quantitative model；Using effective stabilizer content consumption 50% as safe storage life critical pumping rate,
Nitrate esters explosive wastewater thermal acceleration degradation is carried out, temperature coefficient is obtained by effective stabilizer changes of contents rule, is established
Powder charge safe storage life predicting equation；The single temperature thermal acceleration ageing process nitrate esters explosive wastewater of Near-Infrared Absorption Method tracing detection
Powder charge effective stabilizer content obtains the crash time that effective stabilizer content reaches critical pumping rate consumption, passes through powder charge safety
Storage life predicting equation obtains safe storage life at a temperature of nitrate esters explosive wastewater powder charge storage environment.
(1) thermal acceleration aging process prepares near-infrared method detection effective stabilizer content modeling sample
Esters of gallic acid propellant powder is fitted into capillary gas vent ground decrement bottle to be added at (semi-enclosed condition, Fig. 1) Yu Danyi temperature
Fast degradation, interval time sampling, prepares the propellant powder sample of different effective stabilizer contents, is divided into internal calibrations collection sample
With external certificate collection sample.
Nitrate esters propellant herbal cake sandwich style interlayer structure sample, by upper layer tablet 7, sandwich of layers 8, lower layer's tablet 9
It is sequentially overlapped composition herbal cake (Fig. 3), herbal cake is fitted into (semi-enclosed condition, Fig. 1) in the ground decrement bottle with capillary exhaust pipe, in
Thermal acceleration aging at single temperature, interval time take out sandwich of layers 8, obtain different effective stabilizer content propellant samples, point
For internal calibrations collection sample and external certificate collection sample.
Nitric-ester-compound explosive herbal cake sandwich style interlayer structure sample, by upper layer tablet 7, sandwich of layers 8, lower layer's tablet 9 according to
Secondary superposition forms herbal cake (Fig. 3), and herbal cake is fitted into (semi-enclosed condition, Fig. 1) in the ground decrement bottle with capillary exhaust pipe, Yu Dan
Thermal acceleration aging at a temperature of one, interval time take out sandwich of layers 8, obtain the explosive sample of different effective stabilizer contents, be divided into
Internal calibrations collection sample and external certificate collection sample (2) effective stabilizer content near-infrared quantitative model
1. obtaining nitrate esters explosive wastewater sample effective stabilizer content using chemical titration or high performance liquid chromatography
Chemical score determines content range.
2. acquiring near infrared spectrum to internal calibration set sample and verifying collection sample using near infrared spectrometer.Sample mode
For diffusing reflection, sampling wave band is 4000cm-1~12800cm-1, optimize spectral scan parameter, determine optimum resolution, scanning time
Several and sample duplicate measurements number.Collected near infrared spectrum is transmitted to computer.
3. carrying out pretreatment and linear fit to internal calibration set sample near infrared spectrum using chemometrics method.Choosing
Key band is selected, in normalization, baseline is smooth, first derivation, second order derivation, selects in multiplicative scatter correction preprocessing procedures
It selects mode alone or in combination to pre-process acquisition spectrum, regression fit is then carried out using Partial Least Squares, that is, PLS method,
Calibration model is established, external certificate is carried out to model built with external certificate collection sample, internal inspection and outside according to model
Verify the appraisement system combined, successive optimization model.
According to the effective stabilizer content near-infrared quantitative model of foundation, using near infrared spectrometer to selected sample into
The assay of row effective stabilizer obtains near-infrared predicted value, and chemical score and the worst error of near-infrared predicted value should be less than
The repeatability error of standard method.
(3) safe storage life
Nitrate esters explosive wastewater powder charge safe storage life predicting equation
In formula: r indicates that the temperature difference is 10 DEG C of reaction rate temperature coefficient；
t0Indicate storage environment temperature safety storage life, d；
tnIndicate high temperature thermal acceleration ageing time, d；
TnIndicate high temperature thermal acceleration degradation temperature, DEG C；
T0Indicate storage environment temperature, DEG C.
(4) temperature coefficient r
1. debugging safety-type constant temperature oil bath baking oven, keep it constant in 4 different temperatures, temperature interval is 10 DEG C.
2. nitrate esters explosive wastewater is placed on capillary gas vent ground decrement bottle (semi-enclosed condition, Fig. 1), it is respectively placed in 4
Thermal acceleration aging is carried out in a safety-type constant temperature oil bath baking oven of different temperatures, it is fixed according to the effective stabilizer content near-infrared of foundation
Measure model, using near infrared spectrometer to effective stabilizer content carry out tracking and measuring, using effective stabilizer consumption 50% as
Safe storage life critical pumping rate obtains 4 groups of temperature and high temperature thermal acceleration aging crash time corresponding relationship, using formula (4),
Fitting obtains temperature coefficient:
R=10-10/B
(5) thermal acceleration degradation and life prediction
Nitrate esters propellant powder is stacked in capillary exhaust pipe decrement case (Fig. 2), nitrate esters hair is prepared
Penetrate medicine powder charge sample, T at Yu Danyi temperaturenLower thermal acceleration degradation, Selection Center region carry out near infrared spectrum scanning, obtain
The effective stabilizer content for obtaining different ageing times obtains temperature T using effective content consumption 50% as critical pumping rateniLower peace
Determine the time t of agent consumptionn, according to the safe storage life at a temperature of formula (4) extrapolation nitrate esters gun propellant charge storage environment
t0。
Nitrate esters propellant charge is radially divided into two sections, vertically with lower section powder column 13 by upper section powder column 10
Lower superposition is fabricated to two-part powder column (see Fig. 4), T at Yu Danyi temperaturenIt is old to obtain powder charge two-part for lower thermal acceleration degradation
Change sample, the central area 11 of 13 superposed surfaces 12 of near infrared spectrum scanning difference ageing time lower section powder column and grain surface, with
Stabilization agent effective content consumes most fast region and is used as the position Guan Chong, using the position Guan Chong stabilization agent effective content consume 50% as
Critical pumping rate obtains temperature TnThe crash time t of lower stabilization agent consumptionn, according to nitrate at a temperature of formula (3) extrapolation storage environment
The safe storage life t of class propellant loading0。
Nitric-ester-compound explosive powder column is radially divided into two sections, by upper section powder column 10 and lower section powder column 13 vertically under
Superposition is fabricated to two-part powder column (see Fig. 4), T at Yu Danyi temperaturenLower thermal acceleration degradation obtains the aging of powder charge two-part
Sample, the central area 11 of 13 superposed surfaces 12 of near infrared spectrum scanning difference ageing time lower section powder column and grain surface, with peace
Determine agent effective content and consume most fast region as the position Guan Chong, consumes 50% as facing using the position Guan Chong stabilization agent effective content
Boundary's criterion obtains temperature TnThe crash time t of lower stabilization agent consumptionn, extrapolate at a temperature of storage environment according to formula (3) by nitrate
The safe storage life t of class explosive loading0。
Embodiment 1: nitrate esters gun propellant charge safe storage life
It is illustrated using typical triple-base propellant SGF-1 as representative.
(1) thermal acceleration aging process prepares No. II centralite effective content near-infrared modeling sample
Typical triple-base propellant SGF-1 is packed into (semi-hermetic in capillary gas vent ground decrement bottle (Φ=4cm, H=8cm)
Condition, Fig. 1) the thermal acceleration aging at 95 DEG C, it samples within every 4 hours, prepares the triple-base propellant sample of different stable agent contents, receive
Collect 60 SGF-1 aged samples, wherein 29 samples (number 1-29) are used as internal calibrations collection, carries out linear fit and inside
Cross validation, 31 samples integrate that (number is used for external certificate as 30-60) as external certificate.
(2) No. II centralite effective content near-infrared quantitative models
1. the device used, which includes that match is silent, flies the scientific and technological near infrared spectrometer Antaris II of generation that, computer, chemical titration
The device of No. II centralite effective content is measured, and using silent winged scientific and technological chemo metric software TQAnalyst of generation that of match etc..
2. No. II centralite effective content chemical score is obtained according to National Military Standard in triple-base propellant SGF-1
(GJB770B-2005 method 210.1) " centralite bromination method " measures No. II centralite effective content of all samples, range
It is 1.38%~0.20%.
3. close red to triple-base propellant internal calibrations collection sample and the acquisition of external certificate collection sample using near infrared spectrometer
External spectrum.Sample mode is diffusing reflection, and sampling wave band is 4000cm-1-12800cm-1, resolution ratio 8cm-1, scanning 64 times, often
A sample duplicate measurements 5 times.Collected internal calibrations collection near infrared spectrum (Fig. 5) is transmitted to computer by USB data line.
Using your science and technology chemo metric software TQAnalyst of the silent winged generation of match to internal calibration set sample near infrared spectrum into
Row pretreatment and linear fit.Key band is selected, using normalization method, baseline is smooth, first derivation and multiplicative scatter correction
The methods of carry out preprocessing procedures, establish calibration model, external certificate carried out to model built with external certificate collection sample,
The appraisement system that internal inspection and external certificate according to model combine, successive optimization model.Table 1 gives several different pre-
The parameter of processing method, the preferable spectral model of different dimensions.
The parameter of the preferable spectral model of table 1
| Preprocess method | Dimension | R | RMSEP |
| Min-max normalization | 5 | 0.9994 | 0.104 |
| Multiplicative scatter correction | 5 | 0.9991 | 0.095 |
| Normalization+baseline is smooth+first derivation | 5 | 0.9996 | 0.081 |
| First derivative+MSC | 6 | 0.9995 | 0.215 |
| First derivative+SNC | 6 | 0.9987 | 0.097 |
As can be seen from Table 1,4000cm-1~7500cm-1Use normalization, baseline smooth in range and first derivation group
It is optimal to close pretreated model.The preprocess method is finally used, the Partial Least-Squares Regression Model established is when dimension selects 5
Most preferably, evaluation index R=0.9996, the RMSEP=0.081 of model, the regression model is as in triple-base propellant SGF-1 II
Determine agent effective content near-infrared quantitative model.
It is right near infrared spectrometer Antaris II according to the triple-base propellant SGF-1 near-infrared Quantitative Prediction Model of foundation
10 samples for accurately measuring No. II centralite effective content with chemical titration are predicted that detailed results are shown in Table 2, chemistry
The worst error of value and near-infrared predicted value should be less than the repeatability error of standard method.
The table centralite effective content prediction result of sample to be tested II of 2 triple-base propellant SGF-1 10
(3) temperature coefficient r
1. debugging safety-type oil bath baking oven, keep its temperature constant at 95 ± 1 DEG C, 85 ± 1 DEG C, 75 ± 1 DEG C and 65 ± 1 DEG C
2. triple-base propellant SGF-1 is packed into (semi-hermetic in capillary gas vent ground decrement bottle (Φ=4cm, H=8cm)
Condition, Fig. 1), it is respectively placed in 95 ± 1 DEG C, 85 ± 1 DEG C, 75 ± 1 DEG C and 65 ± 1 DEG C safety-type oil bath baking ovens and carries out thermal acceleration
Aging effectively contains the sample of different ageing times using No. II centralite of near infrared spectrometer and triple-base propellant SGF-1
Measure near-infrared quantitative model, No. II centralite effective content variation (being shown in Table 3) of tracking and measuring.Obtain ageing time at 4 temperature
With the corresponding relationship of No. II centralite content.
Ageing time and No. II centralite effective content at table 3 triple-base propellant SGF-1,4 temperature
Fig. 4 gives using No. II centralite effective content consumption 50% as the critical pumping rate of ultimate stage, obtains 4 groups of temperature
With thermal acceleration aging crash time corresponding relationship.
4 different temperatures of table and thermal acceleration aging crash time corresponding relationship
| Temperature/DEG C | 95 | 85 | 75 | 65 |
| Time/d | 5.2395 | 19.1281 | 64.6583 | 248.4874 |
It is fitted, obtains using equation (5):
T=107.9954-17.9331lgt (R2=0.9996)
Temperature coefficient:
R=3.61
(4) thermal acceleration degradation and life prediction
Triple-base propellant SGF-1 powder charge safe storage life predicting equation
Triple-base propellant SGF-1 is stacked in capillary exhaust pipe decrement case (25cm × 25cm × 25cm)
(Fig. 2) prepares triple-base propellant accumulation sample, the thermal acceleration aging at 71 DEG C, in near-infrared method inspection center region II number
Determine agent effective content and is consumed to the acquisition of critical pumping rate 50% high temperature thermal acceleration aging crash time tn=105d, outside according to formula (5)
Push away the safe storage life 55.6a of triple-base propellant SGF-1 powder charge at 30 DEG C of storage environment temperature.
Embodiment 2: nitrate esters propellant loading safe storage life
It is illustrated using exemplary complex modified double base propellant GST-1 as representative.
(1) thermal acceleration aging process prepares compound diazepam agent effective content near-infrared modeling sample
Exemplary complex modified double base propellant GST-1 herbal cake sandwich style interlayer structure sample (Φ=30mm, upper layer tablet
7H1=5mm, sandwich of layers 8H2=10mm, lower layer tablet 9H3=5mm, Fig. 3) it is packed into the ground decrement bottle with capillary gas vent
In (semi-enclosed condition, Fig. 1) thermal acceleration aging at 90 DEG C, take within every 2 hours sandwich of layers 8, obtain answering for different stable agent contents
Composite-modified double-base propellants sample collects 60 GST-1 aged samples, wherein 29 samples (number 1-29) are used as internal school
Positive collection, carries out linear fit and cross-validation, and 31 samples integrate as external certificate (number as 30-60), for external
Verifying.
(2) compound diazepam agent effective content near-infrared quantitative model
1. the device used, which includes that match is silent, flies the scientific and technological near infrared spectrometer Antaris II of generation that, computer, chemical titration
The device of No. II centralite effective content is measured, and using silent winged scientific and technological chemo metric software TQAnalyst of generation that of match etc..
2. chemical titration obtains compound diazepam agent effective content chemistry in exemplary complex modified double base propellant GST-1
Value.
It is 2mm~3mm powder, extracted by ether, extraction process that compound diazepam agent effective content near-infrared modeling sample, which crushes,
A drop diphenylamines sulfuric acid solution is added dropwise to react with residue, terminates and extracts when no blue occurs.
The stuffed conical flask for filling ether extract is placed in not higher than 50 DEG C water-baths of temperature and boils off ether, 50ml is added
Ethyl alcohol；The potassium bromate standard solution 25ml for being accurately added that concentration is 0.2mol/L into stuffed conical flask with brown buret, 8 DEG C
10ml hydrochloric acid solution will accurately be added in stuffed conical flask after keeping the temperature 10~15min in~15 DEG C of water-baths, uniformly rocks 30s, then
The liquor kalii iodide of 10ml mass fraction 15% is added, stoppers bottle stopper, is rapidly 0.1mol/L thio with concentration after gently shaking up
2~3ml starch solution that concentration is 5g/L is added when solution is in yellow in sulfuric acid sodium standard solution titration solution rapid titration,
Continue to be titrated to blue disappearance.Blank test twice is carried out under similarity condition, that is, takes 50ml ethyl alcohol, and 25ml potassium bromate mark is added
Quasi- solution is carried out by above procedure.The volume differences of the consumption sodium thiosulfate standard solution titration solution of blank test twice do not surpass
Cross 0.2ml.
Compound diazepam agent effective content is calculated by formula (6):
In formula:
W-compound diazepam agent effective content (mass percent), %；
V0Sodium thiosulfate standard solution volume, ml consumed by-blank sample；
Sodium thiosulfate standard solution volume consumed by V-titration sample solution, ml；
C-Concentration of Sodium Thiohsulfate Standard Solution, mol/L；
M-and 1.0ml sodium thiosulfate standard solution [c (Na2S2O3)=1.000mol/L] comparable average compound diazepam
Agent molal weight, g/mmol；
M-sample mass, g.
The compound diazepam agent effective content for measuring all samples according to the method described above, in the range of 1.48%~0.30%.
3. using near infrared spectrometer to composite modified double-base propellant powder charge internal calibrations collection sample and external certificate collection
Sample acquires near infrared spectrum.Sample mode is diffusing reflection, and sampling wave band is 4000cm-1-12800cm-1, resolution ratio 8cm-1,
Scanning 64 times, each sample duplicate measurements 5 times.Collected near infrared spectrum (Fig. 7) is transmitted to calculating by USB data line
Machine.
Using your science and technology chemo metric software TQAnalyst of the silent winged generation of match to internal calibration set sample near infrared spectrum into
Row pretreatment and linear fit.Key band is selected, using normalization method, baseline is smooth, first derivation and multiplicative scatter correction
The methods of carry out preprocessing procedures, linear regression is carried out using Partial Least Squares (PLS), calibration model is established, with outer
Portion's verifying collection sample carries out external certificate to model built, the evaluation body that the internal inspection and external certificate according to model combine
System, successive optimization model.Table 5 gives the parameter of preferable spectral model when several different pretreatments methods, different dimensions.
The parameter of the preferable spectral model of table 5
| Preprocess method | Dimension | R | RMSEP |
| Min-max normalization | 5 | 0.9988 | 0.123 |
| Multiplicative scatter correction | 5 | 0.9991 | 0.085 |
| Normalization+baseline is smooth+first derivation | 5 | 0.9996 | 0.085 |
| First derivative+MSC | 7 | 0.9995 | 0.124 |
| First derivative+SNC | 7 | 0.9994 | 0.256 |
As can be seen from Table 5,4000cm-1~7000cm-1The interior use normalization of range, baseline is smooth and first derivation is pre-
The model of processing is optimal.The preprocess method is finally used, it is most that dimension, which selects the Partial Least-Squares Regression Model established when 5,
Good, evaluation index R=0.9996, the RMSEP=0.085 of model, the regression model is as composite modified double-base propellant GST-1
Compound diazepam agent effective content near-infrared quantitative model.
According to the composite modified double-base propellant powder charge GST-1 near-infrared Quantitative Prediction Model of foundation, near infrared spectrum is used
Instrument Antaris II predicts 10 samples for accurately measuring compound diazepam agent effective content with chemical titration, in detail
It the results are shown in Table 6, the worst error of chemical score and near-infrared predicted value is not more than -0.04%, chemical score and near-infrared predicted value
Worst error is less than standard method repeatability error.
10 sample to be tested compound diazepam agent effective contents of 6 composite modified double-base propellant powder charge GST-1 of table predict knot
Fruit
| Sample number into spectrum | Chemical score/% | Near-infrared predicted value/% | Error/% |
| 1 | 1.48 | 1.50 | 0.02 |
| 2 | 1.37 | 1.37 | 0.00 |
| 3 | 1.26 | 1.25 | -0.01 |
| 4 | 1.14 | 1.12 | -0.02 |
| 5 | 1.05 | 1.03 | -0.02 |
| 6 | 0.89 | 0.90 | 0.01 |
| 7 | 0.75 | 0.74 | -0.01 |
| 8 | 0.68 | 0.70 | 0.02 |
| 9 | 0.55 | 0.53 | -0.02 |
| 10 | 0.41 | 0.42 | 0.01 |
(3) temperature coefficient r
1. debugging safety-type low temperature oil bath baking oven, keep its temperature constant at 95 ± 1 DEG C, 85 ± 1 DEG C, 75 ± 1 DEG C and 65 ± 1
℃
2. composite modified double-base propellant powder charge GST-1 herbal cake sandwich style interlayer structure sample (Φ=30mm, upper layer medicine
Piece 7H1=5mm, sandwich of layers 8H2=10mm, lower layer tablet 9H3=5mm, Fig. 3), it is fitted into the ground decrement bottle of capillary gas vent
(semi-enclosed condition, Fig. 1) is respectively placed in 95 ± 1 DEG C, 85 ± 1 DEG C, 75 ± 1 DEG C and 65 ± 1 DEG C safety-type low temperature oil bath baking ovens
Thermal acceleration aging is carried out, for the sample of different ageing times, using composite modified double-base propellant powder charge GST-1 compound diazepam
The variation (being shown in Table 7) of 8 compound diazepam agent effective content of agent effective content near-infrared quantitative model tracking and measuring sandwich of layers obtains 4
The corresponding relationship (Fig. 8) of a temperature lower time and compound diazepam agent effective content.
Ageing time and compound diazepam agent effective content at 74 temperature of composite modified double-base propellant powder charge GST-1 of table
Fig. 8 give using compound diazepam agent consumption 50% as safe storage life critical pumping rate, obtain 4 temperature and
Thermal acceleration aging crash time corresponding relationship (table 8).
8 different temperatures of table and thermal acceleration aging crash time corresponding relationship
| Temperature/DEG C | 95 | 85 | 75 | 65 |
| Safe storage life/d | 5.4276 | 13.8333 | 58.8182 | 172.0000 |
It is fitted, obtains using equation (4):
T=108.47-19.3660lgt (R2=0.9937)
Temperature coefficient:
R=3.28
(4) thermal acceleration degradation and life prediction
Composite modified double-base propellant powder charge GST-1 safe storage life predicting equation
Composite modified double-base propellant powder charge GST-1 powder column is radially divided into two sections, by upper section powder column 10 and lower section
Lower superposition is fabricated to two-part powder column (see Fig. 4) to powder column 13 vertically, and the thermal acceleration degradation at 71 DEG C obtains powder charge two
Segmentation aging sample, the central area 11 of 13 superposed surfaces 12 of near infrared spectrum scanning difference ageing time lower section powder column and powder column table
Face consumes most fast central area 11 as the position Guan Chong using stabilization agent effective content, with the position Guan Chong stabilization agent effective content,
Using effective content consumption 50% as critical pumping rate, the crash time t of 71 DEG C of stabilization agents of temperature consumption is obtainedn=83d, according to formula
(7) extrapolate storage environment at a temperature of composite modified double-base propellant powder charge safe storage life 29.7a.
Embodiment 3: nitric-ester-compound explosive powder charge safe storage life
It is illustrated using typical nitrate explosive NEP-1 as representative.
(1) thermal acceleration aging process prepares No. II centralite effective content modeling sample
Typical nitrate explosive NEP-1 herbal cake sandwich style interlayer structure sample (Φ=30mm, upper layer tablet 7H1=
5mm, sandwich of layers 8H2=10mm, lower layer tablet 9H3=5mm, Fig. 3) it is fitted into (half in the ground decrement bottle with capillary gas vent
Confined condition, Fig. 1) thermal acceleration aging is carried out at 85 DEG C, it takes sandwich of layers 8 within every 2.5 hours, obtains No. II centralite effective content
Different calibration set samples and totally 60, external certificate collection sample, wherein 29 samples (number 1-29) are modeling internal calibrations
Collect sample, participates in linear fit and cross-validation, 31 samples (number 30-60) are used as external certificate collection, for outer
Portion's verifying.
(2) No. II centralite effective content near-infrared modelings
1. the device used, which includes that match is silent, flies the desk-top near infrared spectrometer of generation that science and technology II type of ANTARIS, computer, chemistry
The device of No. II centralite effective content of titration measuring, and using the silent winged scientific and technological chemo metric software TQ of generation that of match
Analyst etc..
2. the chemical score of No. II centralite effective content obtains military according to country in typical nitrate explosive NEP-1 sample
Standard (GJB770B-2005 method 210.1) " centralite bromination method ", measures No. II centralite effective content of all samples,
Range is 0.75%~0.20%.
3. being acquired in nitrate explosive NEP-1 using the silent winged generation that desk-top near infrared spectrometer of science and technology II type of ANTARIS of match
Portion's calibration set sample and external certificate collection sample acquisition atlas of near infrared spectra (see Fig. 9).Sample mode is diffusing reflection, samples wave band
For 4000cm-1~12800cm-1, resolution ratio 8cm-1, scanning times 64, each sample duplicate measurements 5 times.It is collected close
Infrared spectroscopy is transmitted to computer by USB data line.
4. by the silent winged generation that science and technology chemo metric software TQ Analyst of match to internal calibration set sample near infrared light
Spectrum carries out pretreatment and linear fit.Key band is selected, using normalization method, baseline is smooth and the methods of first derivation is to light
Spectrum is pre-processed, and Partial Least Squares (PLS) carries out linear regression, is established No. II centralite effective content near-infrared and is quantified mould
Type, the system combined using cross-validation collection external certificate carry out successive optimization to model.Table 9 give it is several not
The parameter of preferable spectral model when same preprocess method, different dimensions.
The parameter of 9 typical case's preferable spectral model of nitrate explosive NEP-1 of table
| Preprocess method | Dimension | R | RMSEP |
| Min-max normalization | 5 | 0.9965 | 0.215 |
| Multiplicative scatter correction | 5 | 0.9992 | 0.125 |
| Normalization+baseline is smooth+first derivation | 5 | 0.9996 | 0.048 |
| First derivative+MSC | 8 | 0.9995 | 0.218 |
| First derivative+SNC | 8 | 0.9994 | 0.226 |
As can be seen from Table 9,4000cm-1~7000cm-1Use normalization, baseline smooth in range and first derivation group
It is optimal to close pretreated model.The preprocess method is finally used, the Partial Least-Squares Regression Model established is when dimension selects 5
Most preferably, the evaluation index R=0.9996, RMSEP=0.048 of the regression model.The regression model is as nitrate explosive NEP-1
In No. II centralite effective content near-infrared quantitative model.
According to above-mentioned model, it is different that No. II centralite effective content is accurately measured to chemical titration near infrared spectrometer
10 samples to be tested predicted that detailed results are shown in Table 10, and chemical score and the worst error of near-infrared predicted value are not more than
0.03%, the worst error of chemical score and near-infrared predicted value should be less than the repeatability error of standard method.
No. II centralite effective content chemical score of 10 samples to be tested of the typical nitrate explosive NEP-1 of table 10 and prediction
Value
| Sample number into spectrum | Chemical score/% | Near-infrared predicted value/% | Deviation/% |
| 1 | 0.75 | 0.77 | 0.02 |
| 2 | 0.71 | 0.68 | -0.03 |
| 3 | 0.66 | 0.67 | 0.01 |
| 4 | 0.63 | 0.62 | -0.01 |
| 5 | 0.59 | 0.57 | -0.02 |
| 6 | 0.55 | 0.58 | 0.03 |
| 7 | 0.50 | 0.49 | -0.01 |
| 8 | 0.45 | 0.46 | 0.01 |
| 9 | 0.41 | 0.40 | -0.01 |
| 10 | 0.37 | 0.38 | 0.01 |
(3) temperature coefficient r
1. debugging safety-type low temperature oil bath baking oven, keep its temperature constant at 85 ± 1 DEG C, 75 ± 1 DEG C, 65 ± 1 DEG C and 55 ± 1
℃
2. nitrate explosive NEP-1 herbal cake sandwich style interlayer structure sample (Φ=30mm, upper layer tablet 7H1=5mm,
Sandwich of layers 8H2=10mm, lower layer tablet 9H3=5mm, Fig. 3), it is fitted into (semi-hermetic item in the ground decrement bottle of capillary gas vent
Part, Fig. 2), it is respectively placed in 85 ± 1 DEG C, 75 ± 1 DEG C, 65 ± 1 DEG C and 55 ± 1 DEG C safety-type low temperature oil bath baking ovens and carries out hot add
Fast aging, it is quantitative using No. II centralite effective content near-infrared of nitrate explosive NEP-1 for the sample of different ageing times
Model following measures the variation (being shown in Table 11) of effective stabilizer content, obtains 4 temperature lower times and effectively contains with No. II centralite
The corresponding relationship (Figure 10) of amount.
Ageing time and No. II centralite effective content at 11 4 temperature of nitrate explosive NEP-1 of table
Figure 10 gives using No. II centralite effective content consumption 50% as the critical pumping rate of safe storage life, obtains
4 temperature and thermal acceleration aging crash time corresponding relationship (table 12).
12 nitrate explosive NEP-1 different temperatures of table and thermal acceleration aging crash time corresponding relationship
| Temperature/DEG C | 85 | 75 | 65 | 55 |
| Safe storage life/d | 5.4110 | 14.0299 | 55.8201 | 163.3997 |
It is fitted, obtains using equation (4):
T=98.8203-19.7129lgt (R2=0.9950)
Temperature coefficient:
R=3.22
(4) thermal acceleration degradation and life prediction
Nitrate explosive NEP-1 safe storage life predicting equation
Nitrate explosive NEP-1 powder column is radially divided into two sections, by upper section powder column 10 and lower section powder column 13 along axial direction
Superposition is fabricated to two-part powder column (see Fig. 4) up and down, the thermal acceleration degradation at 71 DEG C, obtains powder charge two-part aging examination
Sample, the central area 11 of 13 superposed surfaces 12 of near infrared spectrum scanning difference ageing time lower section powder column and grain surface, with stable
It is the position Guan Chong that agent effective content, which consumes most fast central area to make 11, is made with the position Guan Chong stabilization agent effective content consumption 50%
For critical pumping rate, the crash time t of 71 DEG C of stabilization agents of temperature consumption is obtainedn=24d, according to formula (8) extrapolation storage environment temperature
The safe storage life 8.0a of nitrate explosive NEP-1 powder charge at 30 DEG C.
Claims (1)
1. a kind of nitrate esters explosive wastewater powder charge safe storage life predictor method, which is characterized in that this method utilizes thermal acceleration
Aging process prepares nitrate esters explosive wastewater effective stabilizer content near-infrared modeling sample, using chemical titration or efficient liquid phase
Chromatography obtain stabilization agent effective content chemical score, using near-infrared diffuse spectrometer acquisition spectrum, pass through Chemical Measurement
Establish effective stabilizer content near-infrared quantitative model；Consumed for 50% facing as safe storage life using effective stabilizer content
Boundary's criterion carries out nitrate esters explosive wastewater accelerated aging test, obtains temperature coefficient by effective stabilizer changes of contents rule,
Establish powder charge safe storage life predicting equation；Carry out single temperature powder charge thermal acceleration degradation, Near-Infrared Absorption Method tracing detection
Ageing process nitrate esters explosive wastewater powder charge effective stabilizer content obtains effective stabilizer content and reaches critical pumping rate consumption
Crash time passes through safety at a temperature of powder charge safe storage life predicting equation acquisition nitrate esters explosive wastewater powder charge storage environment
Storage life specifically follows these steps to carry out:
(1) it is fried to prepare nitrate esters explosive wastewater effective stabilizer content near-infrared modeling sample nitrate esters fire for thermal acceleration aging process
Medicine is fitted into the ground decrement bottle of the semi-enclosed condition with capillary exhaust pipe, thermal acceleration aging at Yu Danyi temperature, when interval
Between sample, prepare the nitrate esters explosive wastewater sample of different effective stabilizer contents, be divided into internal calibrations collection sample and outside is tested
Card collection sample；
Esters of gallic acid propellant powder is fitted into capillary gas vent ground decrement bottle with powder sample；
Nitrate esters propellant is fitted into capillary exhaust pipe ground decrement bottle with herbal cake sample, and herbal cake is by upper layer tablet (7), sandwich
Layer (8), lower layer's tablet (9) are sequentially overlapped composition sandwich style interlayer structure；
Nitric-ester-compound explosive is fitted into capillary exhaust pipe ground decrement bottle with herbal cake sample, and herbal cake is by upper layer tablet (7), sandwich of layers
(8), lower layer's tablet (9) is sequentially overlapped composition sandwich style interlayer structure；
(2) effective stabilizer content near-infrared quantitative model
1. obtaining nitrate esters explosive wastewater sample effective stabilizer content chemistry using chemical titration or high performance liquid chromatography
Value, determines content range；
2. acquiring near infrared spectrum, sample mode to internal calibration set sample and external certificate collection sample using near infrared spectrometer
For diffusing reflection, sampling wave band is 4000cm-1~12800cm-1, resolution ratio 8cm-1, scan 64 times, each sample duplicate measurements 5
Secondary, collected near infrared spectrum is transmitted to computer；
3. carrying out pretreatment and linear fit, 4000cm to internal calibration set sample near infrared spectrum using chemometrics method-1~7500cm-1It uses normalization, baseline smooth in range and first derivation combined method progress Pretreated spectra, then uses
Partial Least Squares, that is, PLS method carries out regression fit, establishes calibration model, is carried out with external certificate collection sample to model built outer
Portion's verifying, the appraisement system that internal inspection and external certificate according to model combine, successive optimization model obtain it and effectively pacify
Determine agent content near-infrared quantitative model；
According to the effective stabilizer content near-infrared quantitative model of foundation, had using near infrared spectrometer to selected sample
The assay of stabilization agent is imitated, near-infrared predicted value is obtained, chemical score and the worst error of near-infrared predicted value should be less than standard
The repeatability error of method；
(3) safe storage life predicting equation
Nitrate esters explosive wastewater powder charge safe storage life predicting equation
In formula: r indicates that the temperature difference is 10 DEG C of reaction rate temperature coefficient；
t0Indicate storage environment temperature safety storage life, d；
tnIndicate high temperature thermal acceleration ageing time, d；
TnIndicate high temperature thermal acceleration degradation temperature, DEG C；
T0Indicate storage environment temperature, DEG C；
(4) temperature coefficient r
Logarithm is taken to formula (1) both sides, is obtained after arrangement:
Tn=A+Blgtn (2)
In formula: A, B indicates coefficient, A=T0+10×lgt0/ lgr, B=-10/lgr
Temperature coefficient r=10-10/B, sample different temperatures T is tracked using near-infrared methodniLower thermal acceleration degradation difference aging
The effective stabilizer content of time obtains effective stabilizer content-time changing curve, with effective stabilizer content consumption 50%
As critical pumping rate, different temperatures T is obtainedniCorresponding crash time tni, substituted into (2) formula and carry out one-variable linear regression, used
Least square method obtains coefficient B, and then temperature coefficient r is calculated；
(5) thermal acceleration degradation and life prediction
Nitrate esters propellant powder is stacked in capillary exhaust pipe decrement case, nitrate esters gun propellant charge heap is prepared
Accumulate sample, T at Yu Danyi temperaturenLower thermal acceleration degradation, Selection Center region carry out near infrared spectrum scanning, are had according to it
Stable agent content near-infrared quantitative model is imitated, the effective stabilizer content of different ageing times is obtained, is consumed with effective content
50% is used as critical pumping rate, obtains temperature TniThe time t of lower stabilization agent consumptionn, according to formula (1) extrapolation nitrate esters propellant powder dress
Safe storage life t at a temperature of medicine storage environment0；
Nitrate esters propellant charge is radially divided into two sections, vertically by upper section powder column (10) and lower section powder column (13)
Lower superposition is fabricated to two-part powder column, T at Yu Danyi temperaturenLower thermal acceleration degradation obtains powder charge two-part aging sample,
The central area (11) of near infrared spectrum scanning difference ageing time lower section powder column (13) superposed surfaces (12) and grain surface, to have
It imitates stable agent content and consumes most fast region and be used as the position Guan Chong, consume 50% as facing using the position Guan Chong effective stabilizer content
Boundary's criterion obtains temperature TnThe crash time t of lower stabilization agent consumptionn, according to nitrate esters at a temperature of formula (1) extrapolation storage environment
The safe storage life t of propellant loading0；
Nitric-ester-compound explosive powder column is radially divided into two sections, by upper section powder column (10) and lower section powder column (13) vertically under
Superposition is fabricated to two-part powder column, T at Yu Danyi temperaturenLower thermal acceleration degradation obtains powder charge two-part aging sample, closely
The central area (11) of IR spectrum scanning difference ageing time lower section powder column (13) superposed surfaces (12) and grain surface, with effective
Stable agent content consumes most fast region and is used as the position Guan Chong, consumes 50% as critical using the position Guan Chong effective stabilizer content
Criterion obtains temperature TnThe crash time t of lower stabilization agent consumptionn, extrapolate at a temperature of storage environment according to formula (1) by nitrate esters
The safe storage life t of explosive loading0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610614112.5A CN106290231B (en) | 2016-07-29 | 2016-07-29 | Nitrate esters explosive wastewater powder charge safe storage life predictor method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610614112.5A CN106290231B (en) | 2016-07-29 | 2016-07-29 | Nitrate esters explosive wastewater powder charge safe storage life predictor method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106290231A CN106290231A (en) | 2017-01-04 |
| CN106290231B true CN106290231B (en) | 2019-04-19 |
Family
ID=57663134
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610614112.5A Active CN106290231B (en) | 2016-07-29 | 2016-07-29 | Nitrate esters explosive wastewater powder charge safe storage life predictor method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106290231B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110108744B (en) * | 2019-05-08 | 2021-10-08 | 西安近代化学研究所 | Explosive classification method based on thermal accelerated aging test |
| CN112595656B (en) * | 2020-12-09 | 2022-05-27 | 中国兵器工业第五九研究所 | Testing device and evaluation method for adaptability of explosive device long-storage environment for bomb |
| CN113591282B (en) * | 2021-07-15 | 2023-08-11 | 西安近代化学研究所 | Laboratory assessment method for explosive charge life |
| CN114550846B (en) * | 2022-04-22 | 2022-07-26 | 北京理工大学 | Propellant life estimation method based on single-temperature high-temperature accelerated aging test |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792964B (en) * | 2015-05-04 | 2016-12-07 | 西安近代化学研究所 | A kind of based on the predictor method of NEPE propellant storage life under deadweight induction pressure condition |
| CN105928901B (en) * | 2016-07-11 | 2019-06-07 | 上海创和亿电子科技发展有限公司 | A kind of near-infrared quantitative model construction method that qualitative, quantitative combines |
-
2016
- 2016-07-29 CN CN201610614112.5A patent/CN106290231B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106290231A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106290236B (en) | Nitrate esters explosive wastewater powder charge ultimate stage rapid Estimation method | |
| CN106290231B (en) | Nitrate esters explosive wastewater powder charge safe storage life predictor method | |
| CN103018195B (en) | Method for determination of PCTFE content in PBX explosive by near infrared spectrum | |
| Wang et al. | Soil pH value, organic matter and macronutrients contents prediction using optical diffuse reflectance spectroscopy | |
| CN106290237B (en) | Composite modified double-base propellant powder charge ultimate stage rapid Estimation method | |
| CN102288572A (en) | Method for quickly detecting content of index ingredient of traditional Chinese medicinal material by utilizing near infrared spectrum technique | |
| CN104316489A (en) | Method of detecting adulteration of ganoderma lucidum extract product by near infrared spectroscopy | |
| CN102445428B (en) | Analytical method of tetravalent uranium | |
| CN106525761A (en) | Nitrite detection method based on terahertz spectroscopy scanning | |
| CN106290232B (en) | Accumulation storage single-base gun propellant safe storage life predictor method | |
| CN106290222B (en) | Composite modified double-base propellant powder charge safe storage life predictor method | |
| CN102954946B (en) | By the method for infrared spectrum measurement sulfur content in crude oil | |
| CN103115889A (en) | Method for predicating sulphur content of crude oil by infrared transmittance spectroscopy | |
| CN103063599B (en) | The method of prediction oil density is composed by transmitted infrared light | |
| Xie et al. | An approach to detecting diphenylamine content and assessing chemical stability of single-base propellants by near-infrared reflectance spectroscopy | |
| CN103134770A (en) | Method for eliminating influence on infrared spectrum detection of soil total nitrogen content from moisture | |
| Liu et al. | A feasibility quantitative analysis of free fatty acids in polished rice by fourier transform near‐infrared spectroscopy and chemometrics | |
| CN106290233B (en) | Double base propellant powder charge ultimate stage rapid Estimation method | |
| CN106093292B (en) | A kind of NEPE propellant loadings safe storage life predictor method | |
| CN106018333B (en) | Accumulate double-base propellant ultimate stage rapid Estimation method | |
| CN106290172B (en) | Accumulate triple-base propellant ultimate stage rapid Estimation method | |
| CN103134764B (en) | The method of prediction true boiling point curve of crude oil is composed by transmitted infrared light | |
| CN106226266B (en) | A kind of double base propellant powder charge safe storage life predictor method | |
| CN106290223B (en) | Accumulation storage triple-base propellant safe storage life predictor method | |
| CN106290235B (en) | Accumulate single-base gun propellant ultimate stage rapid Estimation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |