CN108981835A - Utilize the method for the indirect metering gas usage amount of gaseous-pressure - Google Patents
Utilize the method for the indirect metering gas usage amount of gaseous-pressure Download PDFInfo
- Publication number
- CN108981835A CN108981835A CN201810659291.3A CN201810659291A CN108981835A CN 108981835 A CN108981835 A CN 108981835A CN 201810659291 A CN201810659291 A CN 201810659291A CN 108981835 A CN108981835 A CN 108981835A
- Authority
- CN
- China
- Prior art keywords
- peak
- point
- value
- pressure
- initial period
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000007789 gas Substances 0.000 claims abstract description 37
- 239000000567 combustion gas Substances 0.000 claims abstract description 16
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 230000011218 segmentation Effects 0.000 claims abstract description 6
- 239000012634 fragment Substances 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 239000002737 fuel gas Substances 0.000 claims 1
- 239000000284 extract Substances 0.000 abstract description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000003325 tomography Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000011897 real-time detection Methods 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
- G01F1/86—Indirect mass flowmeters, e.g. measuring volume flow and density, temperature or pressure
- G01F1/88—Indirect mass flowmeters, e.g. measuring volume flow and density, temperature or pressure with differential-pressure measurement to determine the volume flow
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The present invention provides a kind of method using the indirect metering gas usage amount of gaseous-pressure, which solve existing electronic gas flowmeters it is at high cost, structure is complicated the technical issues of.The following steps are included: step 1, acquires pressure value；Step 2, judge firepower and data sectional: the change point of firepower being judged to be segmented according to the difference of mean value between the difference and adjacent segment between pressure data consecutive points, the firepower size of every segmentation is identical；Step 3, extract peak-to-valley value: extracting cycle peak-to-valley value, it is denoted as a peak-to-valley value, second extraction is carried out to a peak-to-valley value, be denoted as secondary peak-to-valley value: step 4, calculating cycle length: the first secondary peak point taken out in current fragment data is combustion gas period starting point, finds a next adjacent peak point after which, the spacing between the two peak points is calculated, as the initial period；Step 5, flow rate calculation.The present invention is widely used in gas metering technical field.
Description
Technical field
The present invention relates to a kind of methods of metering gas usage amount, are measured indirectly more particularly to a kind of using gaseous-pressure
The method of combustion gas usage amount.
Background technique
With the development of economy, more and more home-use to have gone up natural gas from rural area to city.Most of family uses
Diaphragm gas meter be the instrument that natural gas volume metering is carried out using flexible thin-walled chamber, although this metering method calculates metering
Accurately, but combustion gas usage amount can only go the number on observation gas meter, flow meter by people, and cannot be real-time transmitted to gas company.
Although most physical structure is complicated, expensive in addition, existing market has electronic gas flowmeter.
Summary of the invention
The present invention for the technical issues of existing electronic gas flowmeter is at high cost, structure is complicated, provide it is a kind of it is at low cost,
The small method using the indirect metering gas usage amount of gaseous-pressure of operand.
For this purpose, the technical scheme is that, comprising the following steps:
Step 1, pressure value is acquired；
Step 2, judge firepower and data sectional:
According to the difference of mean value between the difference and adjacent segment between pressure data consecutive points judge the change point of firepower to be segmented,
The firepower size of every segmentation is identical；
Step 3, peak-to-valley value is extracted:
Extracting cycle peak-to-valley value is denoted as a peak-to-valley value, carries out second extraction to a peak-to-valley value, is denoted as secondary peak valley
Value:
Step 4, calculating cycle length:
The first secondary peak point taken out in current fragment data is combustion gas period starting point, is found after which next
An a adjacent peak point, calculates the spacing between the two peak points, as the initial period；
Step 5, flow rate calculation:
(1) flow measurement is carried out since the starting point in initial period, if metering failure twice in succession, is needed again really
Fixed cycle continues flow measurement, until calculating to the last one peak point；
(2) the combustion gas service life of accumulation calculating, so that calculating combustion gas uses volume.
Preferably, it carries out filtering out the fluctuation of gas pipeline pipe network before step 3: every segmentation is smoothed in terms of respectively
Calculate and filter out the pipe network fluctuation of gas pipeline.
Preferably, threshold value house point processing is carried out after a peak extraction to it.
Preferably, judge whether the initial period is reasonable before step 5:
If the secondary peak point as starting point continuously can snap into peak point plus the initial period three times, determine just
Period selection begin rationally；Otherwise, it is determined that the initial period is chosen unreasonable, next secondary peak point and next after it is found
An a adjacent peak point, chooses the initial period, until finding the reasonable initial period again.
Preferably, in step 5 flow measurement method are as follows:
Peak point is snapped into whenever starting point adds an initial period, starting point is replaced with into the peak point, as new week
The starting point of phase, and gas quantity counts and adds one.
The invention has the advantages that being extracted using the periodically variable relationship of gas outlet pressure by signal processing
The period of change of gas pressure out, the final metering realized to combustion gas volume is used.Method is simple, and operand is small, calculated result
Accurately, robustness is good, being capable of real-time detection domestic gas service condition.
Detailed description of the invention
Fig. 1 is the method for the present invention flow chart；
Fig. 2 is that the present invention judges firepower and data sectional flow chart；
Fig. 3 is the flow chart that the present invention filters out the fluctuation of gas pipeline pipe network；
Fig. 4 is the flow chart of extracting cycle peak-to-valley value of the present invention；
Fig. 5 is calculating cycle length flow chart of the present invention.
Specific embodiment
The present invention is described further below with reference to embodiment.
The device of this method is equipped with single-chip microcontroller and temperature and pressure transmitter, the diaphragm gas meter of temperature and pressure transmitter acquisition
The gaseous-pressure data of outlet are stored into the RAM of single-chip microcontroller, and a combustion gas usage amount is calculated when storage arrives a certain amount of.
As shown in Figure 1, calculate combustion gas usage amount method the following steps are included:
Step 1, pressure value is acquired:
Diaphragm gas meter gas outlet gas pressure force value is acquired by the temperature and pressure transmitter in metering device valve body；
Step 2, judge firepower and data sectional:
After accumulating certain pressure value data volume, according to the difference of mean value between the difference and adjacent segment between pressure data consecutive points
Value judges the change point of firepower to be segmented, and the firepower size of every segmentation at this time is identical；
Step 3, the fluctuation of gas pipeline pipe network is filtered out:
Every segmentation is smoothed respectively to calculate and filter out the fluctuation of the pipe network of gas pipeline；
Step 4, extracting cycle peak-to-valley value:
Peak-to-valley value is extracted by the method that consecutive points compare, is denoted as a peak-to-valley value, threshold is carried out to the peak-to-valley value after extraction
Value house point processing, retains valid data to reduce redundancy；Peak is extracted again by the method that consecutive points compare to a peak-to-valley value
Valley is denoted as secondary peak-to-valley value；
Step 5, calculating cycle length:
(1) the first secondary peak point taken out in current fragment data is combustion gas period starting point, is found after which
A next adjacent peak point, calculates the spacing between the two peak points, as the initial period；
(2) if the secondary peak point as starting point continuously can snap into peak point plus the initial period three times, sentence
Determine the initial period to choose rationally；Otherwise, it is determined that the initial period choose it is unreasonable, find next secondary peak point and its after
A next adjacent peak point, chooses the initial period, until finding the reasonable initial period again.
Step 6, flow rate calculation:
(1) after initial period selection is correct, flow measurement is carried out since the starting point (secondary peak point) in initial period,
If metering failure twice in succession, it is likely encountered the fiery point of small change, needs to return to step 5 and redefines the period, continue flow
Metering, until calculating to the last one peak point；
(2) the combustion gas service life of accumulation calculating, so that calculating combustion gas uses volume.
As shown in Fig. 2, step 2 judges comprising the concrete steps that for firepower and data sectional:
Firstly, input raw pressure signal；Then, (increase by 5 every time from the 51st point to the 51st point traversal reciprocal
Point), if the difference of the pressure value of i-th point of pressure value and the i-th+6 point is greater than 50, and before i-th point 50 points average value
It is greater than 30 with the average value difference of 50 points after the i-th+6 point, then has change fire, i-th point is tomography point, and the position of i-th point of record exists
In array cp；Otherwise, as without becoming fire, first point is denoted as tomography point, and there are in array cp；Each tomography point it is next
Point is to being set as one section between next tomography point, after waiting all data sectionals complete, end judges firepower and data sectional part.
As shown in figure 3, step 3 filters out comprising the concrete steps that for gas pipeline pipe network fluctuation:
(1) the average value bid of the pressure value of preceding 20 points is taken；
(2) first points of pressure value and average value bid makes the difference；If absolute value of the difference is more than or equal to 7, average value
Bid is updated to the value of the point and the mean value of average value bid, and then the pressure value of the point is updated to average value bid；If difference
Absolute value is less than 7, then average value bid does not update, the pressure value of the point is updated to average value bid；It loops through all
Sampled point；
(3) fixed length of window is 50 progress curve smoothings:
A. preceding 25 points: the value of nth point use the mean value of the value of its adjacent each n-1 point in front and back and the value of itself as its
New value；
B. the point (removing the point of preceding 25 points and rear 25 points) among: with the value of its adjacent each 25 point in front and back and its
The mean value of the value of the value new as its itself；
C. after 25 points: the value of nth point use the mean value of the value of its adjacent each 25-n point in front and back and the value of itself as
Its new value；
(4) step (3) are repeated twice, obtains the smooth pipe network curve of cyclical fluctuations；
(5) original waveform subtracts pipe network perturbation waveform, the waveform after as filtering out pipe network fluctuation.
As shown in figure 4, step 4 extracting cycle peak-to-valley value comprises the concrete steps that:
(1) array after filtering out pipe network fluctuation is set as Temp [], and all points obtain Temp1 [] after taking opposite number；
(2) a kind of step-length N is set, is begun stepping through from first point of Temp [] array, if a value is more adjacent than it
Value is all denoted as greatly peak value, and spacing of the next peak point of peak value apart from a upper peak point is at least N；
(3) a kind of step-length N is set, is begun stepping through from first point of Temp1 [] array, if a value is more adjacent than it
Value is all denoted as greatly peak value, and spacing of the next peak point of peak value apart from a upper peak point is at least N, then takes opposite number
The valley point of as former array；
(4) position of all peak points is stored in array locs [], numerical value is stored in array value []；By institute
There is the position of valley point to be placed in array locsd [], numerical value is stored in array valued [].
After peak-to-valley value extracts, carries out threshold value and gives up point:
(1) peak point extracted for previous step is begun stepping through from the 1st point in locs [], checks corresponding position
Value [] in value, if value is greater than+6Pa or value is less than -6Pa, the numerical value for retaining the point of original position is constant；It is no
Then, the numerical value of original position point is set to 0；
(2) location information of the non-zero point of all values in value [] and numerical information are covered into original array locs []
With array value []；
(3) identical as step (1) (2) to the processing of valley point.
As shown in figure 5, step 5 comprises the concrete steps that:
(1) peak value is extracted to array value [] again, is denoted as two minor peaks, records the position of first secondary peak point
It sets and size, as initial period starting point ND；It is final that week is begun look for from the array value [] of peak point of storage
Stop NE, the next peak point after taking ND is as cycle arrest point NE；
(2) the distance between ND and NE nd_length are calculated and illustrates the cycle arrest got if nd_length < 20
Point NE is miscellaneous point, then a peak point after using the point covers current point as NE；
(3) three are found near these three points of ND+nd_length, ND+2*nd_length and ND+3*nd_length
It whether there is peak point in point, if continuously all existed three times, it is correct to illustrate that initial period nd_length chooses, is flowed
Amount calculates；Otherwise, it is incorrect to illustrate that initial period nd_length chooses, next secondary peak point is found after current ND and is made
For period starting point ND, correct initial period nd_length is chosen again, carries out flow rate calculation.
Flow rate calculation: the initial period choose it is correct after, whenever ND+nd_length snaps into peak point, by ND replacement peak thus
It is worth point, as the starting point in new period, and gas quantity counts and adds one.
The single-chip microcontroller of the present embodiment uses 32 ARM microcontrol processor STM32F103ZET6 of ST company production；Gas
Temperature and pressure transmitter uses the DPS310 sensor of Infineon company production.
The present embodiment utilizes the periodically variable relationship of gas outlet pressure, extracts gas pressure by signal processing
Period of change, the final metering realized to combustion gas volume is used.Method is simple, and operand is small, and calculated result is accurate, robustness
It is good, it being capable of real-time detection domestic gas service condition.Hardware configuration only needs microcontrol processor, pressure sensor, structure letter
It is single, it is at low cost.
The only above person, only specific embodiments of the present invention, when cannot be limited the scope of implementation of the present invention with this, therefore
The displacement of its equivalent assemblies, or according to equivalent changes and modifications made by the invention patent protection scope, should still belong to power of the present invention
The scope that sharp claim is covered.
Claims (5)
1. a kind of method using the indirect metering gas usage amount of gaseous-pressure, characterized in that the following steps are included:
Step 1, pressure value is acquired；
Step 2, judge firepower and data sectional:
According to the difference of mean value between the difference and adjacent segment between pressure data consecutive points judge the change point of firepower to be segmented, every point
The firepower size of section is identical；
Step 3, peak-to-valley value is extracted:
Extracting cycle peak-to-valley value is denoted as a peak-to-valley value, carries out second extraction to a peak-to-valley value, is denoted as secondary peak-to-valley value:
Step 4, calculating cycle length:
The first secondary peak point taken out in current fragment data is combustion gas period starting point, finds next phase after which
An adjacent peak point, calculates the spacing between the two peak points, as the initial period；
Step 5, flow rate calculation:
(1) flow measurement is carried out since the starting point in initial period, if metering failure twice in succession, needs to redefine week
Phase continues flow measurement, until calculating to the last one peak point；
(2) the combustion gas service life of accumulation calculating, so that calculating combustion gas uses volume.
2. the method according to claim 1 using the indirect metering gas usage amount of gaseous-pressure, which is characterized in that in institute
It carries out filtering out the fluctuation of gas pipeline pipe network before stating step 3: every segmentation is smoothed respectively to calculate and filter out fuel gas conduit
The pipe network in road fluctuates.
3. the method according to claim 1 using the indirect metering gas usage amount of gaseous-pressure, which is characterized in that described
Threshold value house point processing is carried out after peak extraction to it.
4. the method according to claim 1 using the indirect metering gas usage amount of gaseous-pressure, which is characterized in that in institute
Judge whether the initial period is reasonable before stating step 5:
If the secondary peak point as starting point continuously can snap into peak point plus the initial period three times, initial week is determined
Phase is chosen reasonable；Otherwise, it is determined that the initial period chooses unreasonable, the next secondary peak point of searching and next phase after it
An adjacent peak point, chooses the initial period, until finding the reasonable initial period again.
5. the method according to claim 1 using the indirect metering gas usage amount of gaseous-pressure, which is characterized in that described
The method of flow measurement in step 5 are as follows:
Peak point is snapped into whenever starting point adds an initial period, starting point is replaced with into the peak point, as the new period
Starting point, and gas quantity counts and adds one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810659291.3A CN108981835B (en) | 2018-06-25 | 2018-06-25 | Method for indirectly measuring gas usage amount by gas pressure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810659291.3A CN108981835B (en) | 2018-06-25 | 2018-06-25 | Method for indirectly measuring gas usage amount by gas pressure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108981835A true CN108981835A (en) | 2018-12-11 |
| CN108981835B CN108981835B (en) | 2020-05-12 |
Family
ID=64538585
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810659291.3A Active CN108981835B (en) | 2018-06-25 | 2018-06-25 | Method for indirectly measuring gas usage amount by gas pressure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108981835B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113984134A (en) * | 2021-10-22 | 2022-01-28 | 上海观燃智能科技有限公司 | Method and system for measuring natural gas consumption by adopting sink constant pressure |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104374444A (en) * | 2014-11-11 | 2015-02-25 | 大庆市日上仪器制造有限公司 | Indirect metering method and device for samples injected without passing mass flow meter |
| US20160341620A1 (en) * | 2015-05-21 | 2016-11-24 | Airgas, Inc. | Method and system for monitoring pressure in a gas containment unit |
| CN106322120A (en) * | 2016-08-12 | 2017-01-11 | 姚俊 | Gas metering device capable of preventing gas leakage |
| CN207486871U (en) * | 2017-11-29 | 2018-06-12 | 宁波方太厨具有限公司 | Firepower displaying device and application have the gas-cooker of the firepower displaying device |
-
2018
- 2018-06-25 CN CN201810659291.3A patent/CN108981835B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104374444A (en) * | 2014-11-11 | 2015-02-25 | 大庆市日上仪器制造有限公司 | Indirect metering method and device for samples injected without passing mass flow meter |
| US20160341620A1 (en) * | 2015-05-21 | 2016-11-24 | Airgas, Inc. | Method and system for monitoring pressure in a gas containment unit |
| CN106322120A (en) * | 2016-08-12 | 2017-01-11 | 姚俊 | Gas metering device capable of preventing gas leakage |
| CN207486871U (en) * | 2017-11-29 | 2018-06-12 | 宁波方太厨具有限公司 | Firepower displaying device and application have the gas-cooker of the firepower displaying device |
Non-Patent Citations (1)
| Title |
|---|
| 商荣辉: "膜式燃气表内部流动数值模拟研究", 《中国计量学院学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113984134A (en) * | 2021-10-22 | 2022-01-28 | 上海观燃智能科技有限公司 | Method and system for measuring natural gas consumption by adopting sink constant pressure |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108981835B (en) | 2020-05-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cardell-Oliver et al. | Smart meter analytics to pinpoint opportunities for reducing household water use | |
| CN109858728B (en) | Load prediction method based on industry-divided electricity utilization characteristic analysis | |
| WO2012082802A9 (en) | Methods and system for nonintrusive load monitoring | |
| CN104930353B (en) | A kind of pressure regulation of public supply mains subtracts leakage assessment system | |
| CN104035458A (en) | Intelligent temperature control method, device and system | |
| CN109739905A (en) | A kind of pipe gallery fire alarm method and system | |
| CN108981835A (en) | Utilize the method for the indirect metering gas usage amount of gaseous-pressure | |
| CN106844290B (en) | A kind of time series data processing method based on curve matching | |
| CN111859292A (en) | Water supply leakage monitoring method for night water use active cell | |
| CN108153266B (en) | Calibration method of pump station information intelligent monitoring terminal based on power detection | |
| CN1651871A (en) | Equal water quantity correcting and regulating method of telemetering water meter accuracy | |
| CN203037739U (en) | Accurate acquisition device of power grid frequency in audio forensics | |
| CN111042779B (en) | Method for calculating apparent water absorption index based on ground indicator diagram | |
| CN201908645U (en) | Natural gas-liquid two-phase on-line metering device with orifice plate differential pressure | |
| CN114857498B (en) | Method for confirming resistance characteristics of water supply network | |
| CN109990875A (en) | A kind of gas station's the measurement and calibration system | |
| CN110688362A (en) | Data sectional type storage method based on time stamp | |
| CN115310852A (en) | Improved river ecological runoff evaluation method | |
| CN113407797A (en) | Data mining method for gas stealing behavior by utilizing fractal calculation | |
| CN109323702A (en) | A kind of navigation route planning method and device | |
| CN114963535A (en) | Water heater scale detection method, device, server, storage medium and product | |
| CN112132442A (en) | Method for evaluating load identification effect under intermittent start-stop electric heating working condition | |
| CN116104479A (en) | Oil pumping well production metering method based on edge calculation | |
| CN114033357B (en) | Method and system for measuring dynamic change of load of oil pumping well comprehensive tester | |
| CN108956160A (en) | Vehicular amusement apparatus caution sound measurement method of parameters, system, device and computer readable storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: Room b221, innovation and entrepreneurship base, No. 213-3, Torch Road, torch high tech Industrial Development Zone, Weihai City, Shandong Province Applicant after: Shandong Zhuocheng Intelligent Technology Co.,Ltd. Address before: 264209, Weihai torch high tech Industrial Development Zone, Shandong Province, Torch Road, -2, -213 innovation and entrepreneurship base, room 525, block A Applicant before: WEIHAI ZHUOCHENG GAS SAFETY EQUIPMENT Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210720 Address after: 264200 603, block B, innovation and entrepreneurship base, torch high tech Industrial Development Zone, Weihai City, Shandong Province Patentee after: Weihai Shengtai Intelligent Instrument Co.,Ltd. Address before: Room b221, innovation and entrepreneurship base, 213-3 Torch Road, torch hi tech Industrial Development Zone, Weihai City, Shandong Province, 264200 Patentee before: Shandong Zhuocheng Intelligent Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240423 Address after: 264209 room b221, innovation and entrepreneurship base, No. 213-3, Torch Road, torch high tech Industrial Development Zone, Weihai City, Shandong Province Patentee after: Shandong Zhuocheng Intelligent Technology Co.,Ltd. Country or region after: China Address before: 264200 603, block B, innovation and entrepreneurship base, torch high tech Industrial Development Zone, Weihai City, Shandong Province Patentee before: Weihai Shengtai Intelligent Instrument Co.,Ltd. Country or region before: China |
|
| TR01 | Transfer of patent right |