CN106053701B - A kind of liquor classification method - Google Patents

A kind of liquor classification method Download PDF

Info

Publication number
CN106053701B
CN106053701B CN201610589133.6A CN201610589133A CN106053701B CN 106053701 B CN106053701 B CN 106053701B CN 201610589133 A CN201610589133 A CN 201610589133A CN 106053701 B CN106053701 B CN 106053701B
Authority
CN
China
Prior art keywords
sample
original producton
producton location
true
wine
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
Application number
CN201610589133.6A
Other languages
Chinese (zh)
Other versions
CN106053701A (en
Inventor
赵金松
苏占元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Chuanjiu Group Technology Development Co ltd
Sichuan University of Science and Engineering
Original Assignee
Sichuan University of Science and Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan University of Science and Engineering filed Critical Sichuan University of Science and Engineering
Priority to CN201610589133.6A priority Critical patent/CN106053701B/en
Publication of CN106053701A publication Critical patent/CN106053701A/en
Application granted granted Critical
Publication of CN106053701B publication Critical patent/CN106053701B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating 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/02Column chromatography
    • G01N30/88Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating 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/02Column chromatography
    • G01N30/04Preparation or injection of sample to be analysed
    • G01N30/06Preparation

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The present invention discloses a kind of liquor classification method, belongs to technical field of white spirit identification.The discrimination method is comprised the following steps:Select several true wine in certain original producton location, some wine samples are taken respectively in color-comparison tube, add a certain amount of inner mark solution, enter the gas chromatograph equipped with flame ionization ditector (FID) after well mixed to be detected, retention and standard specimen according to each group swarming on chromatogram contrast carry out it is qualitative, it is quantitative using areas of peak normalization method;The quantitative result that sample to be identified is obtained carries out Euclidean distance cluster analysis in the lump with the true wine database data in original producton location, judges whether sample is the true wine in original producton location by Dendrogram.

Description

A kind of liquor classification method
Technical field
The present invention relates to a kind of original producton location liquor classification method set up based on gas chromatographic technique, belong to liquor classification skill Art field.
Background technology
Geographical sign, is also called mark of origin,《AGREEMENT ON TRADE RELATED ASPECTS OF INTELLECTUAL PROPERTY RIGHTS》(abbreviation TRIPS agreements) The 1st section of 22 articles is defined as:" its mark goes out certain merchandise resources in certain member region, or the area in the region Or somewhere, the extra fine quality of the commodity, prestige or other features are mainly relevant with the geographic origin ".Alcohol product is important Agricultural product there is strong region feature, bigger protection is defined to drinks geographical sign in TRIPS agreements.
China white wine has great impetus as Chinese tradition alcoholic beverage to the development of national economy.In height It is also increasingly serious to white wine abuse under the ordering about of volume profit, the particularly geography symbol product with high added value.Personation Geography symbol product not only compromises the interests of consumer, and have impact on corporate reputation.Thus, China white wine authenticity Detection and adulterated detection have turned into the outstanding problem of Liquor-making industry concern.It is then desired to develop a kind of quick, reliable detection Method ensures the authenticity of white wine information.
At present, the Main Means to white wine detection and analysis are the detection of subjective appreciation, conventional chemical analysis method and sensor. But rely solely on the subjective experience of taster and judge the quality of white wine, subjectivity is very strong, be hardly formed one and objectively comment Price card is accurate.The domestic and international analysis method on original producton location food product is a lot, including ICP-MS, stable isotope mass spectrum, GC- MS, LC-MS, GC, HPLC, infrared spectrum and other analysis methods etc., but it is main to be applied particularly to the method for alcohol product Concentrate on GC-MS, IR, GC and HPLC etc..But no matter which kind of method used, have some limitations.Such as finger-print phase Close Y-factor method Y evaluation to evaluate differentiating that the uniformity of sample is based on fuzzy " similarity is high " or " correlation is high ", work as detected sample When the finger-print of product and the similarity of original producton location sample are between 0.8~0.9, it is difficult to make rational judgement.Simultaneously as The particularity that white wine quality is formed, same batch products by different storage time or storage temperature, in it material phase Mutually conversion, so as to cause the change of product original producton location characteristic (such as lipid, acids and aldehyde material) or forfeiture to make much to be based on Organic acid, alcohol, the analysis method (such as IR, GC-MS, GC and HPLC) of aldehyde composition are difficult to have breakthrough in temporal adaptation.And Based on element composition and content certain methods (such as ICP-MS, stable isotope mass spectrum) if by equipment price it is high, operate skill Art requirement is high, and fluctuation of each element in different batches sample is big, and (main mixing including raw material sources ground, industrial water receives season Section change etc.) etc. factor influence, there is not yet the successful identification case to carrying out counterfeit merchandise production under the conditions of original producton location, its space The intrinsic defect of method is there is also in adaptation.
The content of the invention
According to the deficiencies in the prior art, the present invention provides a kind of original producton location liquor classification set up based on gas chromatographic technique Method.
Technical scheme is as follows:
A kind of liquor classification method, comprises the following steps:
(1) the true wine sample in certain original producton location no less than 30 batches is selected, 3.0~8.0mL of pipette samples is in tool plug colorimetric respectively Guan Zhong, adds the inner mark solution of the n-butyl acetate that 0.05~0.3mL volume fractions are 2%, enters after being well mixed and is furnished with hydrogen The gas chromatograph of flame ionization detector (FID) detected, whole chromatographic peaks that sample analysis is obtained, containing ethanol and Internal standard peak, using areas of peak normalization method, the relative correction factor of each composition is not considered, tries to achieve each Component peak area percentages, And and standard control, the position of composition needed for being determined using retention time;
(2) each peak between ethanol and normal propyl alcohol is named as A areas, each peak name between ethyl lactate and ethyl hexanoate Be B areas, each peak between n-butanol and isoamyl alcohol is named as C areas, by required composition acetaldehyde, methyl alcohol, normal propyl alcohol, ethyl acetate, Isobutanol, acetal, n-butanol, isoamyl alcohol, ethyl butyrate, the degree of 11 compositions of ethyl lactate and ethyl hexanoate and Totally 14 peak area percents of parameters set up that to turn into this kind of original producton location true to the region of A, B, C tri- peak area sum degree Wine database, for subsequent data analysis;
(3) sample to be identified is detected according to above-mentioned processing method and identical chromatographic condition, is obtained 14 ginsengs Number, cluster analysis is carried out using SPSS softwares Euclidean distance, and such as sample to be identified is individually clustered in dendrogram bottom, And being more than or equal to 5 with the genetic distance of the true drinks group in original producton location, then the sample is the true wine in non-original producton location;As sample to be identified is gathered Class is in the middle of the true drinks in original producton location group's dendrogram or individually to gather be a class, but is less than 5 with the genetic distance of the true drinks group in original producton location, Then the sample is the true wine in original producton location.
Wherein, chromatographic condition described above is:Capillary column:AT-930 Liquor Analysis dedicated columns, column length 18m, internal diameter 0.53mm;Carrier gas is high-purity N2, flow velocity is 2.0mL/min;Sample size:1.0 μ L, split ratio is 20: 1, and injector temperature is 240 DEG C, dottle pin purging 6mL/min;Using temperature programming pattern, 50 DEG C of initial temperature keeps 6mL/min, then with 8 DEG C/min journeys After sequence is warming up to 140 DEG C, then 4min is kept to after 180 DEG C with 30 DEG C/min temperature programmings;Detector temperature is 250 DEG C, H2Stream Speed is 30mL/min, and air velocity is 300mL/min, and tail blows 25mL/min.
Described sample to be identified is no more than 3 months apart from the date of manufacture.
The effect of internal standard compound is mainly in the present invention:Internal standard compound and tested component are in same matrix, are eliminated every The secondary error due to caused by sample size difference, and then make subsequently to calculate required composition and internal standard compound using areas of peak normalization method Peak area percent it is more accurate.
Compared with prior art, advantages of the present invention is:The present invention is that the original producton location based on gas chromatographic technique foundation is white Wine discrimination method, by quantitative result and the true wine database data in original producton location that sample to be identified is obtained carry out in the lump Euclid away from From cluster analysis, judge whether sample is the true wine in original producton location by dendrogram;The discrimination method simple and fast, low cost is applicable Certain original producton location white wine in identification apart from the date of manufacture no more than 3 months, and accuracy has reached 90% or so.
Brief description of the drawings
Fig. 1 is the chromatogram of the true wine sample in certain original producton location;
Fig. 2 is the vertical icicle of the average connection (between group) obtained using SPSS software Euclidean distances;
Fig. 3 is the dendrogram of the average connection (between group) obtained using SPSS software Euclidean distances.
Specific embodiment
With reference to specific embodiment, the invention will be further described.
Embodiment 1
The 45 true wine samples in certain original producton location of batch 52%vol that acquisition in January, 2012 in June, 2013 produces, its production batch Number it is shown in Table 1.
The true wine lot number in certain original producton location of the 52%vol of table 1.
45 batch sample 5.0mL is drawn respectively in color-comparison tube, adds the acetic acid that 0.1mL volume fractions are 2% N-butyl inner mark solution, is well mixed the upper gas chromatograph for determination equipped with flame ionization ditector (FID).Wherein it is color Spectral condition is:Capillary column:AT-930 Liquor Analysis dedicated columns, column length 18m, internal diameter 0.53mm;Carrier gas is high-purity N2, flow velocity is 2.0mL/min;Sample size:1.0 μ L, split ratio is 20: 1, and injector temperature is 240 DEG C, dottle pin purging 6mL/min;Using journey Sequence heating mode, 50 DEG C of initial temperature keeps 6mL/min, then with 8 DEG C/min temperature programmings to after 140 DEG C, then with 30 DEG C/ Min temperature programmings keep 4min to after 180 DEG C;Detector temperature is 250 DEG C, H2Flow velocity is 30mL/min, and air velocity is 300mL/min, tail blows 25mL/min.
Whole chromatographic peaks (containing ethanol and n-butyl acetate internal standard) that sample analysis is obtained use areas of peak normalization method (not considering the relative correction factor of each composition) tries to achieve each Component peak area percentages;And and standard control, during using retaining Between determine needed for composition position.Each peak between ethanol and normal propyl alcohol is named as A areas;Between ethyl lactate and ethyl hexanoate Each peak be named as B areas;Each peak between n-butanol and isoamyl alcohol is named as C areas.By required composition acetaldehyde, methyl alcohol, normal propyl alcohol, Ethyl acetate, isobutanol, acetal, n-butanol, isoamyl alcohol, ethyl butyrate, ethyl lactate and ethyl hexanoate 11 the hundred of composition Point than content and the region of A, B, C tri- peak area sum degree, totally 14 peak area percents of parameters set up that turn into should The true wine database in original producton location is planted, for subsequent data analysis.
The result obtained by above-mentioned 45 batch sample is built into the true wine database (table 2) in this kind of original producton location of 52%vol.
The true wine database (unit in this kind of original producton location of 2. 52%vol of table:%)
Testing sample is drawn into 5.0mL in color-comparison tube, the n-butyl acetate that 0.1mL volume fractions are 2% is added Inner mark solution, is well mixed the upper gas chromatograph for determination equipped with flame ionization ditector (FID).Wherein it is chromatographic condition For:Capillary column:AT-930 Liquor Analysis dedicated columns, column length 18m, internal diameter 0.53mm;Carrier gas is high-purity N2, flow velocity is 2.0mL/ min;Sample size:1.0 μ L, split ratio is 20: 1, and injector temperature is 240 DEG C, dottle pin purging 6mL/min;Using temperature programming Pattern, 50 DEG C of initial temperature keeps 6mL/min, then with 8 DEG C/min temperature programmings to after 140 DEG C, then with 30 DEG C/min programs After being warming up to 180 DEG C, 4min is kept;Detector temperature is 250 DEG C, H2Flow velocity is 30mL/min, and air velocity is 300mL/ Min, tail blows 25mL/min.Whole chromatographic peaks (containing ethanol and n-butyl acetate internal standard) that sample analysis is obtained use peak area Normalization method (not considering the relative correction factor of each composition) tries to achieve each Component peak area percentages;And and standard control, profit The position of composition needed for being determined with retention time.
Each peak between ethanol and normal propyl alcohol is named as A areas;Each peak between ethyl lactate and ethyl hexanoate is named as B Area;Each peak between n-butanol and isoamyl alcohol is named as C areas, obtains required composition acetaldehyde, methyl alcohol, normal propyl alcohol, ethyl acetate, different Butanol, acetal, n-butanol, isoamyl alcohol, ethyl butyrate, the degree and A of 11 compositions of ethyl lactate and ethyl hexanoate, The region of B, C tri- peak area sum degree totally 14 peak area percents of parameter.
The wine sample analysis result (unit to be identified of table 3.:%)
Sample Z1 analysis results (being shown in Table 3) to be identified are added in table 2, wherein, the 2nd~46 row is respectively the 1st~45 Number genuine piece sample, the 2nd~15 row are respectively 14 peak area percent data of parameter, 14 parameters of sample Z1 to be identified The rows of data input Excel the 47th, and preserve.
Cluster analysis is carried out using SPSS softwares Euclidean distance, concrete operation step is as follows:
SPSS softwares are opened, is clicked on " file ", " opening text data ", in " file type " selection " Excel ", in catalogue In find the Excel file of above-mentioned preservation, open corresponding worksheet;Click on:" analysis " --- " classification " --- " system is gathered Class ";" V1 " is loaded into " mark case ", other all parameters are chosen, and are loaded into " variable ", and click on " drafting ", chooses " tree Shape figure ", is clicked on " method ", and " Eudidean distances " is selected on " interval " column, is ultimately generated cluster data and (is shown in Table 4, table 5 respectively And Fig. 2), it is drop-down to find dendrogram (see Fig. 3).
Table 4 is that the case treatment of the cluster obtained using SPSS software Euclidean distances is collected
Table 5 is the cluster table of the average connection (between group) obtained using SPSS software Euclidean distances
Cluster table
Observation dendrogram (see Fig. 3), sample Z1 to be identified is clustered in the middle of the true drinks group dendrogram in original producton location, can be sentenced The fixed sample is the true wine in this kind of original producton location of 52%vol.

Claims (2)

1. a kind of liquor classification method, it is characterised in that comprise the following steps:
(1) the true wine sample in certain original producton location no less than 30 batches is selected, 3.0~8.0mL of pipette samples is in color-comparison tube for difference, Add 0.05~0.3mL volume fractions be 2% n-butyl acetate inner mark solution, be well mixed after enter equipped with hydrogen flame from The gas chromatograph of sonization detector detected, whole chromatographic peaks that sample analysis is obtained are used into areas of peak normalization method, The relative correction factor of each composition is not considered, tries to achieve each Component peak area percentages, and and standard control, using retention time It is determined that the position of required composition;
(2) region between ethanol and normal propyl alcohol is named as A areas, the region between ethyl lactate and ethyl hexanoate is named as B Area, the region between n-butanol and isoamyl alcohol is named as C areas, by required composition acetaldehyde, methyl alcohol, normal propyl alcohol, ethyl acetate, isobutyl Alcohol, acetal, n-butanol, isoamyl alcohol, ethyl butyrate, the degree and A of 11 compositions of ethyl lactate and ethyl hexanoate, B, Set up as this kind of true wine number in original producton location by totally 14 peak area percents of parameter for C regional peak area sums degree According to storehouse, for subsequent data analysis;
(3) sample to be identified is detected according to above-mentioned processing method and identical chromatographic condition, is obtained 14 parameters, adopted Cluster analysis is carried out with SPSS softwares Euclidean distance, such as sample to be identified is individually clustered in dendrogram bottom, and with The genetic distance of the true drinks group in original producton location is more than or equal to 5, then the sample is the true wine in non-original producton location;As sample to be identified is clustered It is a class to gather in the middle of the true drinks group dendrogram in original producton location or individually, but is less than 5 with the genetic distance of the true drinks group in original producton location, then should Sample is the true wine in original producton location;
Above-mentioned chromatographic condition is:Capillary column:AT-930 Liquor Analysis dedicated columns, column length 18m, internal diameter 0.53mm;Carrier gas is height Pure N2, flow velocity is 2.0mL/min;Sample size:1.0 μ L, split ratio is 20: 1, and injector temperature is 240 DEG C, dottle pin purging 6mL/ min;Using temperature programming pattern, 50 DEG C of initial temperature keeps 6min, then with 8 DEG C/min temperature programmings to after 140 DEG C, then With 30 DEG C/min temperature programmings to after 180 DEG C, 4min is kept;Detector temperature is 250 DEG C, H2Flow velocity is 30mL/min, air Flow velocity is 300mL/min, and tail blows 25mL/min.
2. a kind of liquor classification method according to claim 1, it is characterised in that:Described sample to be identified distance production Date is no more than 3 months.
CN201610589133.6A 2016-07-25 2016-07-25 A kind of liquor classification method Active CN106053701B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610589133.6A CN106053701B (en) 2016-07-25 2016-07-25 A kind of liquor classification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610589133.6A CN106053701B (en) 2016-07-25 2016-07-25 A kind of liquor classification method

Publications (2)

Publication Number Publication Date
CN106053701A CN106053701A (en) 2016-10-26
CN106053701B true CN106053701B (en) 2017-06-20

Family

ID=57416918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610589133.6A Active CN106053701B (en) 2016-07-25 2016-07-25 A kind of liquor classification method

Country Status (1)

Country Link
CN (1) CN106053701B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108399433B (en) * 2018-03-15 2021-11-12 泸州老窖集团有限责任公司 Classification method based on fingerprint spectrum characteristics of white spirit
CN108445103A (en) * 2018-03-20 2018-08-24 江南大学 Method that is a kind of while detecting acetaldehyde and acetal
CN108765101A (en) * 2018-06-29 2018-11-06 四川理工学院 The flavour of food products self-sensing method of feedback is evaluated by terminal user based on internet
CN110780017A (en) * 2019-12-10 2020-02-11 天津达仁堂京万红药业有限公司 Fingerprint spectrum detection method for fen-flavor liquor
CN114846327A (en) * 2019-12-20 2022-08-02 日本电信电话株式会社 Determination system, verification device, and determination method
CN113917013B (en) * 2021-09-30 2024-02-02 山东省葡萄研究院 Method for measuring methanol and 5 higher alcohols in fruit distilled liquor by internal standard method
CN114062556B (en) * 2021-11-19 2023-10-31 泸州老窖集团有限责任公司 Carbon isotope composition detection method for white spirit flavor substances

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102023137B (en) * 2009-09-18 2014-05-14 贵州国台酒业有限公司 Method for identifying white spirits
EP2717675A4 (en) * 2011-06-10 2015-04-29 Univ Florida Hybrid tomatoes and methods of making hybrid tomatoes
CN103217487B (en) * 2012-12-14 2014-10-22 中华人民共和国中山出入境检验检疫局 Analysis method for tracing import Merlot grape wine production place by using dual finger print
CN105319308B (en) * 2014-07-07 2017-05-17 北京北分天普仪器技术有限公司 Gas chromatography/mass spectrometry analysis apparatus of various compositions of white spirit, and analysis method thereof
CN104297363B (en) * 2014-09-05 2016-05-18 浙江工业大学 A kind of flash distillation gas chromatography-mass spectrography is differentiated the method for fructus schisandrae and kadsura longepedunculata in conjunction with finger-print

Also Published As

Publication number Publication date
CN106053701A (en) 2016-10-26

Similar Documents

Publication Publication Date Title
CN106053701B (en) A kind of liquor classification method
Kuballa et al. Application of NMR for authentication of honey, beer and spices
Wu et al. Origin verification of French red wines using isotope and elemental analyses coupled with chemometrics
Jöbstl et al. Identification of the geographical origin of pumpkin seed oil by the use of rare earth elements and discriminant analysis
CN107478736B (en) Method for judging freshness of hops based on volatile component composition
Sun et al. Real wine or not? Protecting wine with traceability and authenticity for consumers: chemical and technical basis, technique applications, challenge, and perspectives
CN109781918B (en) Gas phase ion mobility spectrometry identification method for yellow rice wine produced by different enterprises
Liu et al. Combination of the 87Sr/86Sr ratio and light stable isotopic values (δ13C, δ15N and δD) for identifying the geographical origin of winter wheat in China
CN103645144B (en) Method for quantitatively analyzing components in white spirit
CN102645464B (en) Method for quickly detecting amino acids in yellow rice wine by using electronic tongue system
CN102967669B (en) Rapid measuring method for carbon acetate stable isotope in table vinegar
CN101487826B (en) Method for identifying aroma type and/or grade of white spirit
CN111044638A (en) Method for classifying and identifying different varieties of flue-cured tobacco leaves
CN111060642A (en) Method for classifying and identifying tobacco leaves of same variety and different producing areas
Liu et al. A review of recent compound-specific isotope analysis studies applied to food authentication
CN106770875A (en) The method that grape wine Production area recognition is carried out based on metabolism group research meanses
CN102507800A (en) Rapid aroma fingerprint identification method for geographical indication protection product of vinegar
CN111693625A (en) GC-IMS-based method for identifying adulterants in cubilose and application
Li et al. Exploring the regional typicality of Australian Shiraz wines using untargeted metabolomics
CN111337605B (en) Method for evaluating authenticity of lotus bee pollen
CN108195985A (en) It is a kind of that pseudo- method is mixed based on volatile ingredient combination PCA and CSA Model checking hops
CN101487824B (en) Method for recognizing beer flavor
CN113533587B (en) Method for identifying variety of chilli powder based on gas phase ion mobility spectrometry
CN115267038A (en) Method for identifying rape royal jelly and application
Yang et al. Qualitative analysis of age and brand of unblended brandy by electronic nose

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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 643000 Huidong learning street 180, Zigong, Sichuan

Patentee after: Sichuan University of Science & Engineering

Address before: 643000 Sichuan Province, Zigong City Hing Road Xueyuan Street No. 180

Patentee before: SICHUAN University OF SCIENCE & ENGINEERING

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210722

Address after: 610000 floor 1, No. 73, Dayuan south 2nd Street, high tech Zone, Chengdu, Sichuan

Patentee after: Sichuan Chuanjiu Group Technology Development Co.,Ltd.

Address before: 643000 Huidong learning street 180, Zigong, Sichuan

Patentee before: Sichuan University of Science & Engineering