CN107870034B - Underwater acoustic velocity measurement method based on phase difference - Google Patents
Underwater acoustic velocity measurement method based on phase difference Download PDFInfo
- Publication number
- CN107870034B CN107870034B CN201711000856.9A CN201711000856A CN107870034B CN 107870034 B CN107870034 B CN 107870034B CN 201711000856 A CN201711000856 A CN 201711000856A CN 107870034 B CN107870034 B CN 107870034B
- Authority
- CN
- China
- Prior art keywords
- signal
- discrete
- sequence
- cross
- wave signal
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses underwater sound based on phase differenceA method for fast measurement includes using dual AD sampling modules to synchronously and respectively pre-sample emission signals x1(t) and pre-sampled received signal x2(t) collecting to obtain two paths of discrete sampling sequences x1(n) and x2(n), then two discrete sampling sequences x1(n) and x2(n) performing cross-spectrum processing and cross-correlation calculation, and decomposing the measurement of the propagation time of the acoustic wave signal in the distance d into the sum of the phase difference in a single period and the whole period number, so that the measurement accuracy of the propagation time of the acoustic wave signal can be realized by improving the calculation accuracy of the phase difference in the single period; the method has the advantages that the propagation time of the sound wave signal is obtained by accurately calculating the phase difference and the whole period number of the sound wave signal in a single period, and on the basis of time-guaranteeing measurement, the method is simple in measurement process, high in measurement precision and low in cost.
Description
Technical Field
The invention relates to an underwater acoustic velocity measurement method, in particular to an underwater acoustic velocity measurement method based on phase difference.
Background
At present, the conventional underwater sound velocity measurement method mainly comprises an empirical formula method and a direct measurement method. The empirical formula method is used for calculating based on parameters obtained by measuring the temperature, salinity and depth in the marine environment to obtain the acoustic velocity of water sound. Although the empirical formula method has high accuracy in measuring the sound velocity, in the method, corresponding instruments are firstly adopted to obtain the temperature, the salinity and the depth in the marine environment, and then the parameters are obtained from the instruments to perform calculation, so that the measurement process is complex, the time for obtaining the sound velocity is long, the hysteresis ratio is obvious, and the underwater sound velocity of the current sea area cannot be obtained in real time.
The direct measurement method mainly utilizes the propagation characteristic of sound waves in seawater, sound wave signals are acquired through an analog-to-digital converter and then directly transmitted to a processor, the processor calculates the propagation time of the sound waves, and then the underwater sound velocity is obtained through calculation. The method has the advantages of simple measurement process, short time for obtaining the sound velocity, negligible hysteresis and capability of obtaining the underwater sound velocity of the current sea area in real time. However, in this method, the accuracy of the digital measurement of the propagation time directly restricts the accuracy of the sound velocity. Although increasing the acquisition frequency of the analog-to-digital converter can improve the measurement accuracy of the propagation time to some extent, the analog-to-digital converter with higher acquisition frequency has higher cost, and the acquisition frequency of the analog-to-digital converter in the state of the art is relatively limited, so that the measurement accuracy of improving the sound velocity in this way is relatively limited.
Disclosure of Invention
The invention aims to solve the technical problem of providing a phase difference-based underwater acoustic velocity measurement method which is simple in measurement process, high in measurement precision and low in cost on the basis of time-guaranteeing measurement.
The technical scheme adopted by the invention for solving the technical problems is as follows: a phase difference-based underwater acoustic velocity measurement method is characterized by comprising the following steps:
placing a transmitting transducer and a receiving transducer in the measured seawater at intervals, and recording the distance between the transmitting transducer and the receiving transducer as d, wherein d is more than or equal to 10cm and less than or equal to 20 cm;
second, using signal generator to generate frequency f0Sine wave signal of f01MHz, the sine wave signal is denoted as se (t), where t is the time domain time;
amplifying the power of the sine wave signal by 4-8 times through a power amplifier to obtain a driving signal, and recording the driving signal as s (t);
fourthly, attenuating the driving signal s (t) through an attenuation circuit, wherein the attenuation multiple is 8-16 to obtain a pre-sampling emission signal, and recording the pre-sampling emission signal as x1(t); and driving the transmitting transducer by the driving signal s (t), wherein the transmitting transducer emits a sound wave signal, the receiving transducer receives the sound wave signal and converts the sound wave signal into a receiving electrical signal, the receiving electrical signal is amplified by 4-8 times and then is subjected to filtering treatment to obtain a pre-sampling receiving signal, and the pre-sampling receiving signal is recorded as x2(t);
Adopting two-way AD sampling module to respectively pre-sample emission signal x1(t) and pre-sampled received signal x2(t) collecting, wherein the sampling frequency of the two-way AD sampling module is fs,20MHz≤fsLess than or equal to 80MHz, and transmitting signal x by pre-sampling1(t) taking a sample to obtain a discrete sample sequence denoted as x1(n) for the pre-sampled received signal x2(t) taking a sample to obtain a discrete sample sequence denoted as x2(N), N is a discrete serial number, N is 0,1,2,3, …, N-1, N is an integer of 2000 or more and 2500 or less;
sixthly, discrete sampling sequence x1(n) and a discrete sampling sequence x2(n) cross-spectrum processing, discrete sampling sequence x1The frequency spectrum of (n) is denoted X1(k) Discrete sampling sequence x2The frequency spectrum of (n) is denoted X2(k) Discrete sample sequence x1(n) and a discrete sampling sequence x2The cross spectrum of (n) is marked as Y (k), and k is a discrete sampling sequence x1(n) and a discrete sampling sequence x2(N) a sequence number of a frequency domain, where k is 0,1,2,3, …, N ' -1, where N ' is a number of points of a discrete sampling sequence participating in FFT operation, and a value of N ' is 128 or 256; the specific process is as follows:
sixthly, 1, respectively adopting FFT calculation to obtain discrete sampling sequence x1(n) frequency spectrum X1(k) And a discrete sample sequence x2(n) frequency spectrum X2(k);
Sixthly-2 adopts the formulaCalculating to obtain x1(n) and x2(n) a cross-spectrum Y (k), whereinIs X1(k) The symbol x is a conjugate operator, and the symbol x is an inner product operator;
seventhly, calculating cross spectrum Y (k) at frequency point f0Phase angle ofThe specific process is as follows:
seventhly-1 sets an intermediate parameter, which is marked as k0By the formula k0=[f0/fs(N'-1)]Calculating cross spectrum Y (k) middle frequency point f0Corresponding k0The term "[ 2 ]]"is the operator of rounding, symbol/is the operator of dividing;
seventhly-2 adopting arctangent function to solve phase angle phi0I.e. byIm(Y(k0) Is Y (k)0) Imaginary part of Re (Y (k)0) Is Y (k)0) The real part of (a);
calculating the whole period number M of sound wave signal propagation, and the specific process is as follows:
[ 1 ] discrete sampling sequence x1(n) and a discrete sampling sequence x2The cross-correlation function of (n) is denoted as R (r) and adopts the formulaCalculating to obtain corresponding cross-correlation function values, wherein L is an integer which is more than or equal to 100 and less than or equal to 200, r is 0,1,2,3, …, and N-L;
(viii) the cross-correlation function with the largest cross-correlation function value calculated in the step (viii) -2 is denoted as r (q), and q is an integer which is not less than 0 and not more than N-L;
-3 according to the formula M ═ qf0/fs]Calculating to obtain M, symbol "[ alpha ]]"is the operator of rounding, symbol/is the operator of dividing;
ninthly, recording the time of the acoustic signal passing through the distance d as T, and adopting a formulaCalculating to obtain T;
and (c) calculating the sound speed c by using the formula c ═ d/T.
Compared with the prior art, the method has the advantages that the pre-sampling transmitting signals x are synchronously and respectively subjected to pre-sampling by adopting the double-path AD sampling module1(t) and Pre-samplingReceived signal x2(t) collecting to obtain two paths of discrete sampling sequences x1(n) and x2(n), then two discrete sampling sequences x1(n) and x2(n) performing cross-spectrum processing and cross-correlation calculation, and decomposing the measurement of the propagation time of the sound wave signal in the distance d into the sum of the phase difference and the whole period number in a single period, thereby improving the calculation accuracy of the phase difference in the single period and measuring the propagation time of the sound wave signal.
Drawings
FIG. 1 is an error curve for measuring absolute errors of travel times at different SNR according to the method of the present invention;
FIG. 2 is an error curve of the absolute error of sound velocity measured by the method of the present invention under different SNR.
Detailed Description
The invention is described in further detail below with reference to the accompanying examples.
Example (b): a phase difference-based underwater acoustic velocity measurement method comprises the following steps:
placing a transmitting transducer and a receiving transducer in the measured seawater at intervals, and recording the distance between the transmitting transducer and the receiving transducer as d, wherein d is more than or equal to 10cm and less than or equal to 20 cm;
second, using signal generator to generate frequency f0Sine wave signal of f01MHz, the sine wave signal is denoted as se (t), where t is the time domain time;
amplifying the power of the sine wave signal by 4-8 times through a power amplifier to obtain a driving signal, and recording the driving signal as s (t);
fourthly, attenuating the driving signal s (t) through an attenuation circuit, wherein the attenuation multiple is 8-16 to obtain a pre-sampling emission signal, and recording the pre-sampling emission signal as x1(t); and driving the transmitting transducer by using the driving signal s (t), wherein the transmitting transducer transmitsGenerating a sound wave signal, receiving the sound wave signal by the receiving transducer, converting the sound wave signal into a receiving electric signal, amplifying the receiving electric signal by 4-8 times, filtering the wave to obtain a pre-sampling receiving signal, and recording the pre-sampling receiving signal as x2(t);
Adopting two-way AD sampling module to respectively pre-sample emission signal x1(t) and pre-sampled received signal x2(t) collecting, wherein the sampling frequency of the two-way AD sampling module is fs,20MHz≤fsLess than or equal to 80MHz, and transmitting signal x by pre-sampling1(t) taking a sample to obtain a discrete sample sequence denoted as x1(n) for the pre-sampled received signal x2(t) taking a sample to obtain a discrete sample sequence denoted as x2(N), N is a discrete serial number, N is 0,1,2,3, …, N-1, N is an integer of 2000 or more and 2500 or less;
sixthly, discrete sampling sequence x1(n) and a discrete sampling sequence x2(n) cross-spectrum processing, discrete sampling sequence x1The frequency spectrum of (n) is denoted X1(k) Discrete sampling sequence x2The frequency spectrum of (n) is denoted X2(k) Discrete sample sequence x1(n) and a discrete sampling sequence x2The cross spectrum of (n) is marked as Y (k), and k is a discrete sampling sequence x1(n) and a discrete sampling sequence x2(N) a sequence number of a frequency domain, where k is 0,1,2,3, …, N ' -1, where N ' is a number of points of a discrete sampling sequence participating in FFT operation, and a value of N ' is 128 or 256; the specific process is as follows:
sixthly, 1, respectively adopting FFT calculation to obtain discrete sampling sequence x1(n) frequency spectrum X1(k) And a discrete sample sequence x2(n) frequency spectrum X2(k);
Sixthly-2 adopts the formulaCalculating to obtain x1(n) and x2(n) a cross-spectrum Y (k), whereinIs X1(k) Of a conjugated complex sequence of (A), a symbolSymbol x is the conjugate operator, symbol x is the inner product operator;
seventhly, calculating cross spectrum Y (k) at frequency point f0Phase angle ofThe specific process is as follows:
seventhly-1 sets an intermediate parameter, which is marked as k0By the formula k0=[f0/fs(N'-1)]Calculating cross spectrum Y (k) middle frequency point f0Corresponding k0The term "[ 2 ]]"is the operator of rounding, symbol/is the operator of dividing;
seventhly-2 adopting arctangent function to solve phase angle phi0I.e. byIm(Y(k0) Is Y (k)0) Imaginary part of Re (Y (k)0) Is Y (k)0) The real part of (a);
calculating the whole period number M of sound wave signal propagation, and the specific process is as follows:
[ 1 ] discrete sampling sequence x1(n) and a discrete sampling sequence x2The cross-correlation function of (n) is denoted as R (r) and adopts the formulaCalculating to obtain corresponding cross-correlation function values, wherein L is an integer which is more than or equal to 100 and less than or equal to 200, r is 0,1,2,3, …, and N-L;
(viii) the cross-correlation function with the largest cross-correlation function value calculated in the step (viii) -2 is denoted as r (q), and q is an integer which is not less than 0 and not more than N-L;
-3 according to the formula M ═ qf0/fs]Calculating to obtain M, symbol "[ alpha ]]"is the operator of rounding, symbol/is the operator of dividing;
ninthly, recording the time of the acoustic signal passing through the distance d as T, and adopting a formulaCalculating to obtain T;
and (c) calculating the sound speed c by using the formula c ═ d/T.
The method of the invention is simulated in a Matlab platform, the sound velocity is assumed to be 1500 m/s, the method is adopted for measurement when the SNR is respectively 20dB, 30dB and 40dB, the absolute error of the measured propagation time is shown in figure 1, and the absolute error of the measured sound velocity is shown in figure 2.
Analyzing fig. 1 and 2 may conclude that: along with the improvement of the SNR, the measurement error is obviously reduced, and the precision is improved. When the SNR is 40dB, the absolute error of the measurement of the propagation time is less than 0.674 multiplied by 10-9Seconds, i.e., 0.674 nanoseconds; and the corresponding sound velocity measurement absolute error is less than 0.012 m/s. In actual measurement, the propagation distance of the acoustic signal is only 10 cm-20 cm, and the SNR (signal to noise ratio) of the environment is>40dB is easy to meet, and high-precision sound velocity measurement can be obtained by measuring by the method.
Claims (1)
1. A phase difference-based underwater acoustic velocity measurement method is characterized by comprising the following steps:
placing a transmitting transducer and a receiving transducer in the measured seawater at intervals, and recording the distance between the transmitting transducer and the receiving transducer as d, wherein d is more than or equal to 10cm and less than or equal to 20 cm;
second, using signal generator to generate frequency f0Sine wave signal of f01MHz, the sine wave signal is denoted as se (t), where t is the time domain time;
amplifying the power of the sine wave signal by 4-8 times through a power amplifier to obtain a driving signal, and recording the driving signal as s (t);
fourthly, attenuating the driving signal s (t) through an attenuation circuit, wherein the attenuation multiple is 8-16 to obtain a pre-sampling emission signal, and recording the pre-sampling emission signal as x1(t); and driving the transmitting transducer by the driving signal s (t), wherein the transmitting transducer transmits a sound wave signal, the receiving transducer receives the sound wave signal and converts the sound wave signal into a receiving electrical signal, and the receiving electrical signal is amplified 4Filtering after 8 times to obtain pre-sampled received signal, and recording the pre-sampled received signal as x2(t);
Adopting two-way AD sampling module to respectively pre-sample emission signal x1(t) and pre-sampled received signal x2(t) collecting, wherein the sampling frequency of the two-way AD sampling module is fs,20MHz≤fsLess than or equal to 80MHz, and transmitting signal x by pre-sampling1(t) taking a sample to obtain a discrete sample sequence denoted as x1(n) for the pre-sampled received signal x2(t) taking a sample to obtain a discrete sample sequence denoted as x2(N), N is a discrete serial number, N is 0,1,2,3, …, N-1, N is an integer of 2000 or more and 2500 or less;
sixthly, discrete sampling sequence x1(n) and a discrete sampling sequence x2(n) cross-spectrum processing, discrete sampling sequence x1The frequency spectrum of (n) is denoted X1(k) Discrete sampling sequence x2The frequency spectrum of (n) is denoted X2(k) Discrete sample sequence x1(n) and a discrete sampling sequence x2The cross spectrum of (n) is marked as Y (k), and k is a discrete sampling sequence x1(n) and a discrete sampling sequence x2(N) a sequence number of a frequency domain, where k is 0,1,2,3, …, N ' -1, where N ' is a number of points of a discrete sampling sequence participating in FFT operation, and a value of N ' is 128 or 256; the specific process is as follows:
sixthly, 1, respectively adopting FFT calculation to obtain discrete sampling sequence x1(n) frequency spectrum X1(k) And a discrete sample sequence x2(n) frequency spectrum X2(k);
Sixthly-2 adopts the formulaCalculating to obtain x1(n) and x2(n) a cross-spectrum Y (k), whereinIs X1(k) The symbol x is a conjugate operator, and the symbol x is an inner product operator;
seventhly, calculating cross spectrum Y (k) at frequency point f0Phase angle ofThe specific process is as follows:
seventhly-1 sets an intermediate parameter, which is marked as k0By the formula k0=[f0/fs(N'-1)]Calculating cross spectrum Y (k) middle frequency point f0Corresponding k0The term "[ 2 ]]"is the operator of rounding, symbol/is the operator of dividing;
seventhly-2 adopting arctangent function to solve phase angle phi0I.e. byIm(Y(k0) Is Y (k)0) Imaginary part of Re (Y (k)0) Is Y (k)0) The real part of (a);
calculating the whole period number M of sound wave signal propagation, and the specific process is as follows:
[ 1 ] discrete sampling sequence x1(n) and a discrete sampling sequence x2The cross-correlation function of (n) is denoted as R (r) and adopts the formulaCalculating to obtain corresponding cross-correlation function values, wherein L is an integer which is more than or equal to 100 and less than or equal to 200, r is 0,1,2,3, …, and N-L;
(viii) the cross-correlation function with the largest cross-correlation function value calculated in the step (viii) -2 is denoted as r (q), and q is an integer which is not less than 0 and not more than N-L;
-3 according to the formula M ═ qf0/fs]Calculating to obtain M, symbol "[ alpha ]]"is the operator of rounding, symbol/is the operator of dividing;
ninthly, recording the time of the acoustic signal passing through the distance d as T, and adopting a formulaCalculating to obtain T;
and (c) calculating the sound speed c by using the formula c ═ d/T.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711000856.9A CN107870034B (en) | 2017-10-24 | 2017-10-24 | Underwater acoustic velocity measurement method based on phase difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711000856.9A CN107870034B (en) | 2017-10-24 | 2017-10-24 | Underwater acoustic velocity measurement method based on phase difference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107870034A CN107870034A (en) | 2018-04-03 |
CN107870034B true CN107870034B (en) | 2019-12-24 |
Family
ID=61753218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711000856.9A Active CN107870034B (en) | 2017-10-24 | 2017-10-24 | Underwater acoustic velocity measurement method based on phase difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107870034B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111352099B (en) * | 2018-12-20 | 2022-05-10 | 宁波大学科学技术学院 | Time delay estimation method based on cross-correlation signal phase decomposition |
CN112067301B (en) * | 2020-09-10 | 2022-05-06 | 宁波大学科学技术学院 | S-shaped water turbine combination mode comprehensive performance test experimental device |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6296878A (en) * | 1985-10-24 | 1987-05-06 | Tech Res & Dev Inst Of Japan Def Agency | Sonic wave propagation simulating system |
CN104634878A (en) * | 2014-12-16 | 2015-05-20 | 北京林业大学 | Wood damage monitoring method based on acoustic emission technique |
CN105572650A (en) * | 2015-12-15 | 2016-05-11 | 宁波大学 | Broadband multiple correlation flow velocity measurement method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8305257B2 (en) * | 2009-09-02 | 2012-11-06 | Trizna Dennis B | Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents |
-
2017
- 2017-10-24 CN CN201711000856.9A patent/CN107870034B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6296878A (en) * | 1985-10-24 | 1987-05-06 | Tech Res & Dev Inst Of Japan Def Agency | Sonic wave propagation simulating system |
CN104634878A (en) * | 2014-12-16 | 2015-05-20 | 北京林业大学 | Wood damage monitoring method based on acoustic emission technique |
CN105572650A (en) * | 2015-12-15 | 2016-05-11 | 宁波大学 | Broadband multiple correlation flow velocity measurement method |
Non-Patent Citations (1)
Title |
---|
基于AD8302的相位差法测声速;花添雨 等;《物理实验》;20170831;第37卷(第8期);10-13 * |
Also Published As
Publication number | Publication date |
---|---|
CN107870034A (en) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7535797B2 (en) | High-resolution ultrasound displacement measurement apparatus and method | |
CN101813673B (en) | Acoustic signal processing device and method for detecting concentration of trace binary-component gas | |
Ricci et al. | FPGA-based system for in-line measurement of velocity profiles of fluids in industrial pipe flow | |
CN106461437B (en) | Extended range ADC flow meter | |
CN112254802B (en) | Voice reconstruction method based on millimeter wave radar phase ranging | |
CN107870034B (en) | Underwater acoustic velocity measurement method based on phase difference | |
CN102353957B (en) | Method for processing multi-beam bathymetric data based on variable-bandwidth filter | |
CN105997147A (en) | Ultrasonic pulse Doppler imaging method and device | |
CN114167423A (en) | Radar sea wave parameter measuring method based on depth regression network | |
CN103845080A (en) | Ultrasonic umbilical cord blood measuring system and method based on linear frequency modulation coding | |
CN103654861A (en) | Ultrasonic umbilical cord blood flow measuring system and ultrasonic umbilical cord blood flow measuring method both based on pseudorandom codes | |
Sun et al. | Measurement of sound absorption by underwater acoustic material using pulse-separation method | |
CN105232088A (en) | Ultrasonic imaging processing method and system based on RF data | |
RU2460093C1 (en) | Method of measuring distance using sonar | |
Kalmykov et al. | A FMCW—Interferometry approach for ultrasonic flow meters | |
CN101856242A (en) | Doppler imaging method and Doppler imaging device for pulse waves | |
CN108020818A (en) | A kind of sine pulse signal distance measuring method under noise based on sliding DFT | |
CN106992820B (en) | A kind of the underwater sound signal design and detection method of oceanographic instrumentation safety beacon | |
RU2593622C1 (en) | Method of measuring radial velocity of object at its noise emission | |
CN108761465A (en) | A kind of ternary array passive ranging method based on distance dimension Wave beam forming | |
JP3881078B2 (en) | Frequency estimation method, frequency estimation device, Doppler sonar and tidal meter | |
Ashraf et al. | High accuracy time of flight measurement using digital signal processing techniques for subsea applications | |
JP4964344B2 (en) | Doppler measuring instrument and tide meter | |
RU2467350C2 (en) | Method and device of signal detection with alternating doppler effect present | |
RU126921U1 (en) | DIGITAL ULTRASONIC DEVICE FOR MEASURING THE FREQUENCY OF THE HEART CONTRACTIONS OF THE FRUIT |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221010 Address after: No.19 Yunkai Road, Beibei District, Chongqing 400700 Patentee after: Bihu Jiaxiao photoelectric technology (Chongqing) Co.,Ltd. Address before: 315212 No. 505 Yuxiu Road, Zhuangshi Street, Zhenhai District, Ningbo City, Zhejiang Province Patentee before: COLLEGE OF SCIENCE & TECHNOLOGY NINGBO University |
|
TR01 | Transfer of patent right |