CN106092159B - 具解调结果可信度自诊断功能的光纤光栅传感方法 - Google Patents
具解调结果可信度自诊断功能的光纤光栅传感方法 Download PDFInfo
- Publication number
- CN106092159B CN106092159B CN201610363698.2A CN201610363698A CN106092159B CN 106092159 B CN106092159 B CN 106092159B CN 201610363698 A CN201610363698 A CN 201610363698A CN 106092159 B CN106092159 B CN 106092159B
- Authority
- CN
- China
- Prior art keywords
- demodulation result
- fiber grating
- spectrum
- confidence level
- diagnostic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000013307 optical fiber Substances 0.000 title description 9
- 238000001228 spectrum Methods 0.000 claims abstract description 42
- 239000000835 fiber Substances 0.000 claims abstract description 40
- 230000003595 spectral effect Effects 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000000985 reflectance spectrum Methods 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 3
- 238000004092 self-diagnosis Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 230000002542 deteriorative effect Effects 0.000 description 2
- 238000002405 diagnostic procedure Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Transform (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种具解调结果可信度自诊断功能的光纤光栅传感系统，该系统通过在常规的光纤光栅传感系统中嵌入一个特殊的光谱诊断软件模块，以实现对传感解调结果可信度的在线自诊断。光谱诊断软件模块包括以下内容：1.对传感系统获得的光纤传感器光谱进行预处理；2.计算光谱特征指标及其相对变化；3.根据光谱特征指标及其相对变化对光谱信号恶化造成的解调误差进行诊断分析；4.利用分析所得的解调误差对解调结果可信度进行标记并存储。利用此诊断软件模块，光纤光栅传感系统可在线自动地排除错误的解调结果，确保解调结果的可信度。
Description
技术领域
本发明属于光栅传感器领域，具体涉及一种具解调结果可信度自诊断功能的光纤光栅传感系统及方法。
背景技术
光纤光栅具有抗电磁干扰、电绝缘性好、体积小、重量轻、可串联复用等特点，利用其制作的传感器在民用设施及军用设备上都有广泛的应用。光纤光栅在通过设计不同结构对其进行封装之后，可以实现温度、应力/应变、压力、湿度、折射率、流速、流量、粘度等物理量的检测。但是无论进行怎样封装，其最基本的传感原理都是利用光纤光栅光谱中心波长的偏移实现对外界物理量变化的传感。光纤光栅传感系统一般主要包含传感器网络，光谱采集子模块和光谱解调子模块。实现光谱采集功能的硬件组合形式主要有扫描激光器与光电探测器的组合，以及宽带光源与光谱仪的组合，但是无论通过什么方式对光谱进行采集，最终实现光谱信号解调的原理都是检测光纤光栅传感器光谱中心波长的偏移。而实现光纤光栅传感器光谱的偏移解调的方法有很多，例如比较传统的峰值检测法，质心法，拟合法等，以及比较新颖的相关法，相位法，补偿法等。
然而，在实际应用中，尤其在进行长期的传感测量时，光纤光栅传感光谱往往会出现逐步恶化的情况。光谱的恶化既有可能含有信噪比的降低，也有可能含有光谱形状的变形。光谱的恶化非常容易造成解调精度的下降，甚至解调的错误。在目前的传感系统中，不论采用哪种解调方法，总是能够输出解调信息，但是，解调出的信息是真实反映了被测物理量的变化，还是由于光纤光栅传感器的性能退化造成的，这在实际应用中是非常难以判断的。因此，为了尽量避免实际应用中出现解调错误的情况，对解调信息的可信度进行实时在线的判断是非常必要的，然而目前的解调系统均缺少这样的判断环节。
发明内容
鉴于此，本发明的目的是提供一种具解调结果可信度自诊断功能的光纤光栅传感系统及方法。
本发明的目的之一是通过以下技术方案实现的，一种具解调结果可信度自诊断功能的光纤光栅传感系统，包括光纤光栅光谱自诊断解调子系统，所述光纤光栅光谱自诊断解调子系统包括常规的光纤光栅传感解调软件模块，所述光纤光栅光谱自诊断解调子系统还包括可对传感光谱恶化造成的解调误差进行分析的光谱诊断软件模块，所述光谱诊断软件模块与常规的光纤解调软件模块相结合。
本发明的目的之二是通过以下技术方案实现的，一种具解调结果可信度自诊断功能的方法，包括以下步骤：
S1.对传感系统采集的光纤传感器光谱进行预处理；
S2.计算光谱特征指标及其相对变化；
S3.根据光谱特征及其相对变化分析光谱恶化造成的解调误差；
S4.利用分析诊断的解调误差对解调结果可信度进行标记并存储。
进一步，光谱特征指标及其相对变化包括A1、A2、A3、A4中的任意一个或几个的组合；
A1＝FWHM/FWHM0，A2＝SMR，A3＝ILA/ILA0，A4＝IOA，其中，
FWHM表示光谱带宽，FWHM＝ΔλL+ΔλR，其中ΔλL及ΔλR分别表示光纤光栅反射光谱的左半宽和右半宽度，FWHM0为光纤光栅初始光谱的带宽。
SMR表示旁瓣主瓣比，计算公式为SMR＝Rsl/Rmax，Rmax及Rsl分别表示光纤光栅反射光谱主瓣和旁瓣的最大值。
ILA表示局部非对称性指标，计算公式为ILA0表示光纤光栅初始光谱的局部非对称性指标；
IOA表示整体非对称性指标，计算公式为IOA＝min[ρa(τ)]，其中非对称程度函数ρa(τ)的计算式为
上式中R(λ)为光纤光栅的反射光谱。
进一步，所述光谱特征指标及其相对变化为A1、A2、A3、A4，步骤S3和S4中诊断分析方法为：判断光谱特征是否满足条件Case1，若满足则诊断标记为解调结果可信；若不满足，则判断光谱特征是否满足条件Case2，若满足则诊断标记为解调结果不可信；若不满足，则比较Error的值与5％的大小，若Error<5％，则诊断标记为解调结果可信，若Error≥5％，则诊断标记为解调结果不可信；
Case 1:A1<1.618and A2<0.200and A3<1.618and A4<0.020；
Case 2:A1>2.000and A2>0.618and A3>2.000and A4>0.060；
其中 w1,w2,w3,w4表示权重。
进一步，w1+w2+w3+w4＝10％。
由于采用了上述技术方案，本发明具有如下的优点：
本发明直接利用传感器的光谱特征及其相对变化对解调结果的可信度状态进行诊断的方法非常适于现场实时在线的监测，并且此方法非常简单实用，通常只需要在现有传感系统中的解调软件上嵌入光谱诊断软件模块即可实现具有可信度标记或解调误差估计的解调结果，利用诊断信息可以及时地排除错误的解调结果，确保解调结果的可信度。
附图说明
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步的详细描述，其中：
图1为常规光纤光栅传感系统基本结构；
图2为在线自诊断型光纤光栅传感系统结构；
图3为光纤光栅光谱的特征刻画参数示意图；
图4为在线光谱自诊断软件模块的一般流程。
具体实施方式
以下将结合附图，对本发明的优选实施例进行详细的描述；应当理解，优选实施例仅为了说明本发明，而不是为了限制本发明的保护范围。
常规的光纤光栅传感系统基本结构如图1所示，而利用本发明实现对光纤光栅传感解调结果可信度在线实时判断的传感系统基本结构如图2所示。主要的特征在于在常规的光纤光栅解调子系统中增加了一个光谱诊断软件模块，并能输出解调结果可信度的诊断标记。
由于在实际的光纤光栅传感应用中，后端能够获的信息只有传感器的光谱信息，因此就只能通过对光谱特征的分析来实现对解调结果可信度实时在线的诊断。
在线光谱自诊断的一般流程如图4所示，在进行光谱特征指标及其相对变化的计算之前，通常需要对光谱进行一定的预处理，预处理主要包括为特定的光纤光栅传感光谱确定一定的光谱计算范围及进行一定的平滑处理。倘若光纤光栅传感器进行串联复，光谱计算范围内应当仅含一支传感器的光谱；光谱平滑的方式很多，其主要目的是消除随机噪声对计算光纤光栅光谱特征的影响。
一种具解调结果可信度自诊断功能的方法，包括以下步骤：S1.对传感系统采集的光纤传感器光谱进行预处理；S2.计算光谱特征指标及其相对变化；S3.根据光谱特征及其相对变化分析光谱恶化造成的解调误差；S4.利用分析诊断的解调误差对解调结果可信度进行标记并存储。
光谱特征指标及其相对变化包括A1、A2、A3、A4中的任意一个或几个的组合；
A1＝FWHM/FWHM0，A2＝SMR，A3＝ILA/ILA0，A4＝IOA，
光谱带宽(FWHM)，反映主瓣的宽度，计算公式为FWHM＝ΔλL+ΔλR，其中ΔλL及ΔλR表示光纤光栅反射光谱的左半宽和右半宽度；
旁瓣主瓣比(Sidelobe-to-Mainlobe Ratio，SMR)，反映光谱主瓣与旁瓣的关系，计算公式为SMR＝Rsl/Rmax，Rmax及Rsl分别表示光纤光栅反射光谱主瓣和旁瓣的最大值，与常用的边模抑制比SMSR(Side mode suppression ratio)的关系为SMSR＝10log(1/SMR)；
局部非对称性指标(Index of Local Asymmetry，ILA)，反映光谱主瓣的非对称性，计算公式为ILA0表示光纤光栅初始光谱的局部非对称性指标；整体非对称性指标(Index of Overall Asymmetry，IOA)，反映用于计算的光谱范围[λmin,λmax]内整体非对称性，计算公式为IOA＝min[ρa(τ)]，其中非对称程度函数ρa(τ)的计算式为
上式中R(λ)为光纤光栅的反射光谱。
假设光纤光栅传感器原始光谱的带宽及局部非对称性指标分别为FWHM0和ILA0，原始光谱一般可选择为光纤光栅传感系统安装校正之后的光谱。然后利用A1＝FWHM/FWHM0，A2＝SMR，A3＝ILA/ILA0，A4＝IOA可以对解调结果的可信度进行诊断，具体的诊断方法如表1所示：
表1光纤光栅解调结果可信度的诊断方法
关于表1中权重(w1,w2,w3,w4)的设置：
权重满足关系w1+w2+w3+w4＝10％；系统中采用的解调算法不同，权重的分配可具有差异，例如对于峰值检测法，对其影响较大的是带宽及局部非对称性的变化，因此w1，w3的权重应当较大，可取为(w1＝4.5％，w2＝0.5％，w3＝4.5％,w4＝0.5％)；对于质心法，旁瓣及主瓣对其解调结果的影响都较大，因此w1，w2，w3的权重应当较大，可取为(w1＝3％，w2＝3％， w3＝3％，w4＝1％)；对于高斯拟合法，谱形整体的非对称性其解调结果的影响较大，因此w4 的权重应当较大，可取为(w1＝2％，w2＝2％，w3＝2％，w3＝4％)。
本发明直接利用传感器的光谱特征及其相对变化对解调结果可信度进行诊断的方法非常适于现场实时在线的监测，并且此方法非常简单实用，通常只需要在现有传感系统中的解调软件上嵌入光谱诊断软件模块即可实现具有可信度标记或误差估计的解调结果，利用诊断信息可以及时地排除错误的解调结果，确保解调结果的可信度。
在本发明中，并不一定限制于综合采用利用四个光谱指标及其相对变化，也可以单独采用某个指标或其相对变化，以及某些指标的组合进行判断(优选实施例中判断依据采用了黄金比例0.618和1.618)。
以上所述仅为本发明的优选实施例，并不用于限制本发明，显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种具解调结果可信度自诊断功能的方法，其特征在于：包括以下步骤：
S1.对传感系统采集的光纤传感器光谱进行预处理；
S2.计算光谱特征指标及其相对变化；
S3.根据光谱特征及其相对变化分析光谱恶化造成的解调误差；
S4.利用分析诊断的解调误差对解调结果可信度进行标记并存储；
光谱特征指标及其相对变化包括A1、A2、A3、A4中的任意一个或几个的组合；
A1＝FWHM/FWHM0，A2＝SMR，A3＝ILA/ILA0，A4＝IOA，
其中，FWHM表示光谱带宽，FWHM＝ΔλL+ΔλR，其中ΔλL及ΔλR分别表示光纤光栅反射光谱的左半宽和右半宽度，FWHM0为光纤光栅初始光谱的带宽；SMR表示旁瓣主瓣比，计算公式为SMR＝Rsl/Rmax，Rmax及Rsl分别表示光纤光栅反射光谱主瓣和旁瓣的最大值；
ILA表示局部非对称性指标，计算公式为：
ILA0表示光纤光栅初始光谱的局部非对称性指标；
IOA表示整体非对称性指标，计算公式为IOA＝min[ρa(τ)]，其中非对称程度函数ρa(τ)的计算式为：
上式中R(λ)为光纤光栅的反射光谱。
2.根据权利要求1所述的具解调结果可信度自诊断功能的方法，其特征在于：所述光谱特征指标及其相对变化为A1、A2、A3、A4，步骤S3和S4中诊断分析方法为：
判断光谱特征及其相对变化是否满足条件Case1，若满足则诊断结果为解调结果可信；若不满足，则判断光谱特征及其相对变化是否满足条件Case2，
若满足则诊断结果为解调结果不可信；若不满足，则比较Error的值与5％的大小，若Error<5％，则诊断标记为解调结果可信，若Error≥5％，则诊断标记为解调结果不可信；
Case1:A1<1.618 and A2<0.200 and A3<1.618 and A4<0.020；
Case2:A1>2.000 and A2>0.618 and A3>2.000 and A4>0.060；
其中，w1,w2,w3,w4表示权重。
3.根据权利要求2所述的具解调结果可信度自诊断功能的方法，其特征在于：w1+w2+w3+w4＝10％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610363698.2A CN106092159B (zh) | 2016-05-27 | 2016-05-27 | 具解调结果可信度自诊断功能的光纤光栅传感方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610363698.2A CN106092159B (zh) | 2016-05-27 | 2016-05-27 | 具解调结果可信度自诊断功能的光纤光栅传感方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106092159A CN106092159A (zh) | 2016-11-09 |
| CN106092159B true CN106092159B (zh) | 2018-12-28 |
Family
ID=57230047
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610363698.2A Active CN106092159B (zh) | 2016-05-27 | 2016-05-27 | 具解调结果可信度自诊断功能的光纤光栅传感方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106092159B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107796421B (zh) * | 2017-10-13 | 2018-08-21 | 北京卫星环境工程研究所 | 快速确定fbg反射光谱异常模式的方法 |
| JP6709241B2 (ja) * | 2018-02-26 | 2020-06-10 | 株式会社Ｓｕｂａｒｕ | 診断装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6573990B1 (en) * | 2000-02-21 | 2003-06-03 | Tektronix, Inc. | Optical system providing concurrent detection of a calibration signal and a test signal in an optical spectrum analyzer |
| JP2007178258A (ja) * | 2005-12-28 | 2007-07-12 | Yokogawa Electric Corp | 光スペクトラムアナライザおよびこれを用いた光ファイバセンサシステム |
| CN102419348B (zh) * | 2011-08-19 | 2013-05-08 | 北京航空航天大学 | 一种基于光纤布拉格光栅的声发射信号功率型无损检测方法 |
| US9417103B2 (en) * | 2011-09-20 | 2016-08-16 | Schlumberger Technology Corporation | Multiple spectrum channel, multiple sensor fiber optic monitoring system |
| CN104316090B (zh) * | 2014-11-14 | 2016-08-31 | 山东大学 | 一种温度自补偿高分辨率高频光纤光栅解调系统的温度自补偿方法 |
| CN105444922B (zh) * | 2015-11-13 | 2016-09-14 | 济南大学 | 光纤光栅温度传感器波长偏移校正方法及测温装置 |
-
2016
- 2016-05-27 CN CN201610363698.2A patent/CN106092159B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106092159A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10260964B2 (en) | Optical fiber sensor system for detecting temperature changes in an aircraft | |
| EP2726836B1 (en) | Otdr trace analysis in pon systems | |
| CN108871607A (zh) | 一种面向分布式光纤拉曼传感器的高精度温度解调方法 | |
| RU2601669C2 (ru) | Система мониторинга измерительной схемы турбореактивного двигателя | |
| EP2989705B1 (en) | Method and apparatus for defect pre-warning of power device | |
| CN104111082A (zh) | 一种高精度fbg传感信号寻峰方法 | |
| CN106092159B (zh) | 具解调结果可信度自诊断功能的光纤光栅传感方法 | |
| CN108696368B (zh) | 一种网元健康状态的检测方法及设备 | |
| CN104221475A (zh) | 用于照明控制的光传感器的故障检测、定位及性能监控 | |
| US6786635B2 (en) | Turbine blade (bucket) health monitoring and prognosis using neural network based diagnostic techniques in conjunction with pyrometer signals | |
| CN110324081A (zh) | 分布式多节点协同的光纤布拉格光栅传感器故障定位方法 | |
| EP3126858B1 (en) | System and method for health monitoring of electrical systems | |
| EP2506103B1 (fr) | Procédé et dispositif de détection automatique d'une donnée air erronée sur un aéronef | |
| CN107077134A (zh) | 用于对监控飞机发动机的系统进行验证的工具 | |
| CN116955091B (zh) | 基于机器学习的数据中心故障检测系统 | |
| US8285514B2 (en) | Sensor fault detection systems and methods thereof | |
| EP3259571B1 (en) | Method and apparatus for detecting shared risk link groups | |
| KR102575917B1 (ko) | 클라우드 기반 가상 센서를 이용한 IoT 센서 이상 진단 방법 및 시스템 | |
| Abdelli et al. | A BiLSTM-CNN based multitask learning approach for fiber fault diagnosis | |
| CN108254000A (zh) | 一种基于改进模拟退火法的fbg传感网光谱解调方法 | |
| CN109307482B (zh) | 温度不敏感的光纤形状重构系统 | |
| JP7415654B2 (ja) | 電池監視システム | |
| CN108917636B (zh) | 一种基于子集窗自适应算法的分布式光纤应变解调方法 | |
| KR102032653B1 (ko) | 결함 평가를 이용한 기계 상태 모니터링 장치 및 방법 | |
| KR20160039078A (ko) | 발전소 배관 운전신호 분석 방법 |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210720 Address after: 710077 718, block a, Haixing city square, Keji Road, high tech Zone, Xi'an City, Shaanxi Province Patentee after: Xi'an zhicaiquan Technology Transfer Center Co.,Ltd. Address before: 400044 No. 174 Shapingba street, Shapingba District, Chongqing Patentee before: Chongqing University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211103 Address after: 257000 Room 308, building 3, Dongying Software Park, No. 228, Nanyi Road, development zone, Dongying City, Shandong Province Patentee after: Dongkai Shuke (Shandong) Industrial Park Co.,Ltd. Address before: 710077 718, block a, Haixing city square, Keji Road, high tech Zone, Xi'an City, Shaanxi Province Patentee before: Xi'an zhicaiquan Technology Transfer Center Co.,Ltd. |