CN113008413B - Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network - Google Patents
Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network Download PDFInfo
- Publication number
- CN113008413B CN113008413B CN202110290463.6A CN202110290463A CN113008413B CN 113008413 B CN113008413 B CN 113008413B CN 202110290463 A CN202110290463 A CN 202110290463A CN 113008413 B CN113008413 B CN 113008413B
- Authority
- CN
- China
- Prior art keywords
- temperature
- heating
- sensor
- point
- ijk
- 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
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/10—Thermometers specially adapted for specific purposes for measuring temperature within piled or stacked materials
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Aiming at the horizontal warehouse adopted by most grains stored in China, the invention provides a method for accurately solving the target point position of heating grains by utilizing the existing temperature measuring network by combining the existing real-time dynamic grain temperature monitoring network among the horizontal warehouse grains and utilizing the temperature measuring point data and the change thereof acquired by the temperature measuring network. The method specifically comprises the following steps: s1, setting a temperature measurement network in a flat warehouse; s2, determining a neighborhood of a heating point in the horizontal warehouse; s3, combining a temperature measuring network, and calculating according to the neighborhood of the heating point to obtain a target point position of the heating grain. According to the method, the position of the starting point can be accurately calculated, so that ventilation or bin turnover treatment can be carried out on a heating area, rapid cooling is realized, and the grain storage safety is ensured.
Description
[ field of technology ]
The invention belongs to the technical field of grain storage, and particularly relates to a method for accurately reversely solving a heating grain target point position by using an existing temperature measuring network.
[ background Art ]
Most of the grains stored in China adopt a horizontal warehouse or a tall horizontal warehouse, and in the grain storage process, measures such as grain storage moisture control, grain insect killing, grain storage storehouse ventilation and cooling and the like are needed to achieve a systematic grain storage environment and maintain the grain storage environment as much as possible. However, even under the systematic grain environment, the grain seeds serve as organisms, and the moisture rising and insect pest growing caused by grain biological activities (grain respiration and the like) occur, and macroscopic appearance is that some randomly dispersed heating points always exist in the stored grains, and the heating points can cause the temperature rise of local grain piles, so that the quality degradation of the stored grains is caused, and the stored grains are lost. In order to control grain storage losses such as grain browning, mildew and the like caused by overhigh local temperature rise, ventilation or bin turnover treatment is generally required for the heating grain stack. If the heating position of the grain pile can be accurately positioned, ventilation or bin turnover can be carried out point to point, and rapid cooling can be realized on the premise of energy conservation, so that a method for rapidly positioning the heating position of the grain pile is urgently needed.
[ invention ]
In order to solve the problems, the invention provides a method for accurately solving the target point position of heating grains by using the existing temperature measuring network by combining the existing real-time dynamic grain temperature monitoring network between the warehouse grains in the flat house and utilizing the temperature measuring point data and the change thereof acquired by the temperature measuring network aiming at the flat house warehouse adopted by most grains stored in China.
The invention is realized by the following technical scheme, and provides a method for accurately reversely solving a heating grain target point position by using the existing temperature measuring network, which comprises the following steps:
s1, setting a temperature measurement network in a flat warehouse;
s2, determining a neighborhood of a heating point in the horizontal warehouse;
s3, combining a temperature measuring network, and calculating according to the neighborhood of the heating point to obtain a target point position of the heating grain.
In particular, the step S1 is specifically implemented according to the following method:
the sensors in the storehouse of the square storehouse are uniformly distributed according to rows, columns and layers, and the distributed sensors are numbered and are designated as a ijk The physical position of each sensor is marked by a Cartesian coordinate system, and the sensors are temperature measuring sensors.
In particular, the step S2 is specifically implemented according to the following method:
and (3) carrying out inspection on a temperature measuring network in the flat warehouse in real time, if one sensor reports a temperature abnormal value, taking each temperature measuring point with the associated row, column and layer numbers of +1 and-1 and 27 temperature points in total of the sensor reporting the temperature abnormal value as calculation basis when the heating point is in the neighborhood of the sensor reporting the temperature abnormal value.
In particular, in the step S2, if the plurality of sensors report abnormal temperature values, the sensors reporting abnormal temperature values and the sensors with row, column and layer numbers +1 and-1 associated with the sensors report abnormal temperature values are taken as calculation basis for a total of 27 temperature points.
In particular, the temperature in the flat warehouse is diffused around by taking the target heating point O as the center of a circle and is attenuated linearly, and the sensor a is used for ijk The temperature is highest, so sensor a ijk Nearest to the target heating point O, when the point O is equal to a ijk When the distance is relatively short, the distance d between the isotherms and the sensor is known from the geometric relationship of circumference segmentation x 、d y 、d z Aliquoting, wherein the step S3 is specifically implemented according to the following method:
s31 sensor a for reporting abnormal temperature and highest temperature ijk Is (x) 1 ,y 1 ,z 1 ) The display temperature value is K 1 The reported coordinate point positions of the second and third high temperature sensors are (x) 2 ,y 2 ,z 2 )、(x 3 ,y 3 ,z 3 ) The display temperature value is K 2 、K 3 Defining target heating point as O (a, b, c) and temperature as K 0 The coordinate points of the 3 sensors are respectively positioned at the centers of O (a, b, c) and the radii are respectively R 1 、R 2 、R 3 And Wen Qiumian to give the following formula:
in the formula (1), ρ x 、ρ y 、ρ z The attenuation rates of the temperature along the linearity of each direction are respectively shown in the unit of ℃/m;
s32 temperature is diffused around the target heating point O as the center of circle and linearly decays, because of the sensor a ijk The temperature is highest, so sensor a ijk Nearest to the target heating point O, to sensor a ijk Adjacent sensors a in the same column and layer as the center i-1jk 、a ijk 、a i+1jk Temperature change is affected by ρ x Due to the outward diffusion of temperature around the heating target point, sensor a i-1jk And a i+1jk The middle temperature is far away from the target heating point and is connected with the sensor a ijk At the same time, the two are positioned on the same side of the conduction region of the heating point, and the temperature between the two is linearly attenuated along the x direction, so that the sensor a is selected i-1jk And a i+1jk The lower temperature K shown in x To calculate the basis, the attenuation rate ρ along the x-axis is obtained according to the following equation x ：
Similarly, with sensor a ijk Selecting adjacent sensors a in the same row and layer as the center ij-1k 、a ij+1k The medium temperature shows a smaller temperature K y For calculation basis, the attenuation rate ρ along the y-axis can be obtained y ：
Similarly, with sensor a ijk Selecting adjacent sensors a in the same row and column as the center ijk-1 、a ijk+1 The medium temperature shows a smaller temperature K z For calculation basis, the attenuation rate ρ along the z-axis can be obtained z ：
S33 reporting the highest abnormal temperature sensor a ijk The point position coordinates of the same plane as the target heating point O, which are right under or right above the point position, are O' (x) 1 、y 1 C) defining the temperature at this point as K' 0 The relative positional relationship and the temperature conduction relationship can be known:
s34 target heating spot O (a, b, c) conducting temperature to O' (x) 1 、y 1 C) satisfying the following formula:
the coordinate positions of the target heating points can be calculated and obtained by combining the formulas (1), (2), (3), (4), (5) and (6) to be O (a, b and c), and the target point positions of the heating grains can be obtained.
The invention provides a method for accurately solving the target point position of heating grains by utilizing the existing temperature measuring network, which can rapidly calculate the target point position of the heating grains, thereby carrying out ventilation or bin turnover treatment on the heating grains according to the calculation result, rapidly and accurately realizing the cooling of the heating point, avoiding the problems of storage loss and the like caused by the deterioration of the quality of stored grains due to the fact that the heating grains cannot be timely treated, and simultaneously avoiding unnecessary operation caused by inaccurate calculation result and increasing the labor cost.
[ description of the drawings ]
FIG. 1 is a diagram of a storage warehouse and its storage state in a method of accurately solving a target point of a heated grain by using an existing temperature measurement network, taking a tall bungalow as an example;
FIG. 2 is a diagram showing a network arrangement for temperature measurement between the storage and the storehouse in a method for accurately reversely solving a target point position of heating grains by using an existing temperature measurement network according to the present invention;
FIG. 3 is a diagram showing the correlation of heating points obtained by the temperature measurement network in the method of precisely solving the heating grain target points by using the existing temperature measurement network;
FIG. 4 is a diagram Wen Qiumian showing an abnormal temperature sensor, a second high temperature sensor, a third high temperature sensor and a target heating point in a method for accurately solving the target heating point by using the existing temperature measuring network;
FIG. 5 is a graph showing the gradient of temperature field along x, y and z directions in a method of precisely solving the target point position of heating grain by using the existing temperature measurement network;
FIG. 6 is a graph showing the gradient of temperature field along x and y directions in a method of precisely solving the target point position of heating grain by using the existing temperature measuring network;
FIG. 7 is a graph showing the gradient of temperature field along y and z directions in a method of precisely solving the target point position of heating grain by using the existing temperature measuring network;
FIG. 8 is a schematic diagram of calculation of a target heating spot in a method for accurately solving a heating grain target spot by using an existing temperature measurement network according to the present invention.
[ detailed description ] of the invention
The present invention will be described in further detail with reference to the following examples in order to make the objects, technical solutions and advantages of the present invention more apparent.
Referring to fig. 1-2, fig. 1 is a grain storehouse and grain storage state diagram thereof, in fig. 2, each layer of temperature measuring network consists of 6 rows, 8 columns and 4 layers of sensor temperature measuring points, the positions and the positions are 1.5m from the storehouse wall, and the distance between the temperature measuring sensors between the rows and the columns is d x ＝d y For 4m, the surface layer sensor is 0.5m away from the grain surface, the lower layer sensor is 0.3m away from the bottom of the bin, the upper layer, the middle layer and the middle layer are equally divided, and for the convenience of calculation, the layer spacing d is set z For normal temperature grain storage bin, normal grain storage temperature is required to be within 28 ℃ (within 15 ℃ of low temperature bin) when the grain stacking height is 6.2 m. Generally, if the grain storage temperature dynamically measured by a certain temperature measuring sensor exceeds 28 ℃, but does not exceed 30 ℃, the grain storage temperature can be treated by a ventilation cooling mode, and if the grain storage temperature exceeds 30 ℃, the grain storage temperature needs to be treated by a bin turning operation. The invention provides a method for accurately reversely solving a heating grain target point position by utilizing an existing temperature measuring network, which comprises the following specific implementation steps:
s1, numbering each temperature measuring sensor on the three-dimensional temperature measuring network, and setting the number as a ijk I, j and k are numbers of 6 rows, 8 columns and 4 layers respectively, a Cartesian coordinate system is established for the selected grain storehouse, each a ijk A Cartesian coordinate system (x, y, z) is used for representing the definite physical position;
s2, through temperature inspection, wherein the coordinate of the highest temperature point reported in real time is (13.5,9.5,3.9) and the number is a 432 The sensor shows a temperature K 1 =32 ℃, as shown in FIG. 3, take the point and ANDA total of 27 temperature measuring points in the related field are taken as observation areas;
s3, assuming that the attenuation rates of the internal temperature of the grain pile along the x, y and z axis directions are ρ respectively x 、ρ y 、ρ z As shown in fig. 3, the coordinates of the second high temperature point among 27 high temperature points in the observation area are (13.5,9.5,5.7), number a 431 The sensor shows a temperature K 2 The coordinates of the third highest temperature point are (13.5,5.5,3.9), number a =30.9℃ 422 The sensor shows a temperature K 3 The specific calculation steps of the heating target point coordinates are as follows:
s31, the target heating point position is O (a, b, c), the temperature is unknown, and the assumption is K 0 As shown in fig. 4, sensor a 432 、a 431 、a 422 The coordinate points of (a) and (b) are respectively calculated by adopting a formula (1) on the spherical surfaces with the same temperature by taking O (a, b and c) as the spherical surfaces with the same temperature, so as to obtain the following formula:
s32 as shown in FIG. 8, with sensor a 432 Adjacent sensors a in the same column and layer as the center 332 、a 432 、a 532 Temperature change is affected by ρ x Sensor a) 332 Is 26.7 ℃, sensor a) 532 The temperature of (a) is 29.1 ℃, and the sensor a is used for the outward diffusion of the temperature around the heating target point 332 And sensor a 432 On the same side of the conduction region of the heating point, the attenuation rate ρ along the x-axis is obtained by calculation according to the formula (2) x ：
Similarly, with sensor a 432 As the center, adjacent sensors a of the same row and layer 422 、a 442 The temperatures were 30.5℃and 25.8℃respectively. Selecting 25.8deg.C with smaller temperature as calculation basis, and calculating according to formula (3) to obtain attenuation rate along y-axisρ y ：
Similarly, with sensor a 432 Centered, adjacent sensors a in the same row and column 431 、a 433 The temperature is 30.9 ℃ and 28.8 ℃ respectively, 28.8 ℃ with smaller temperature is selected as the calculation basis, and the attenuation rate rho along the z axis can be obtained by calculating according to the formula (4) z ：
S33 as shown in FIG. 8, the sensor shows the highest temperature a 432 The point position coordinates of the same plane as the target heating point O and just above the point position are O '(13.5, 9.5, c), and the temperature of the point position is assumed to be K' 0 According to the relative position relation and the temperature conduction relation, the following formula is obtained by combining the formula (5):
s34 the target heating spot O (a, b, c) conducts the temperature to O' (13.5, 9.5, c), the following formula is obtained in combination with formula (6):
the above formulas are combined, and the calculation can be obtained: k (K) 0 ＝34.8℃、K' 0 ＝33.1℃、ρ x ＝1.33℃/m、ρ y ＝1.55℃/m、ρ z Because the target heating point coordinates were calculated to be O (14.4,8,4.5), the average temperature measured in the grain depot at the location of the pair of coordinates was 34.8 ℃ within a radius of 0.25m, the temperatures within the area of 0.25m at the point were all approximately 34.8 ℃ and the temperatures outside the area of 0.25m at the point were all less than 34 ℃, indicating thatThe coordinates of the heating points obtained by the method are positioned in the actual heat source field, and can reflect the actual situation.
In summary, the method provided by the invention can rapidly calculate the target point position of the heating grain, so that ventilation or bin turnover treatment can be performed on the heating grain according to the calculation result, the temperature of the heating point can be rapidly and accurately reduced, and the problems of storage loss and the like caused by deterioration of the quality of the stored grain due to the fact that the heating grain cannot be timely treated are avoided.
Claims (1)
1. The method for accurately reversely solving the heating grain target point position by utilizing the existing temperature measurement network is characterized by comprising the following steps of:
s1, setting a temperature measuring network in a flat warehouse, wherein the step S1 is implemented according to the following method:
the sensors in the storehouse of the square storehouse are uniformly distributed according to rows, columns and layers, and the distributed sensors are numbered and are designated as a ijk The physical position of each sensor is marked by a Cartesian coordinate system, and the sensors are temperature measuring sensors;
s2, determining a neighborhood of heating points in a flat warehouse, wherein the step S2 is implemented specifically according to the following method:
the method comprises the steps of carrying out inspection on a temperature measuring network in a flat warehouse in real time, if one sensor reports a temperature abnormal value, heating points are in the neighborhood of the sensor reporting the temperature abnormal value, 27 temperature points are taken as calculation basis in total by the associated row, column and layer numbers of the temperature measuring points +1 and-1 and the sensor reporting the temperature abnormal value, and step S2 takes 27 temperature points as calculation basis by the sensor reporting the highest temperature abnormal value and the sensor associated row, column and layer numbers of +1 and-1 if a plurality of sensors report the temperature abnormal value;
s3, combining a temperature measuring network, calculating according to the neighborhood of the heating point to obtain a target point position of the heating grain, wherein the temperature in the flat warehouse is diffused and linearly attenuated around the target heating point O as the center of a circle, and the sensor a is used for ijk The temperature is highest, so sensor a ijk Nearest to the target heating point O, when the point O is equal to a ijk When the distance is relatively short, the distance d between the isotherms and the sensor is known from the geometric relationship of circumference segmentation x 、d y 、d z Aliquoting, wherein the step S3 is specifically implemented according to the following method:
s31 sensor a for reporting abnormal temperature and highest temperature ijk Is (x) 1 ,y 1 ,z 1 ) The display temperature value is K 1 The reported coordinate point positions of the second and third high temperature sensors are (x) 2 ,y 2 ,z 2 )、(x 3 ,y 3 ,z 3 ) The display temperature value is K 2 、K 3 Defining target heating point as O (a, b, c) and temperature as K 0 The coordinate points of the 3 sensors are respectively positioned at the centers of O (a, b, c) and the radii are respectively R 1 、R 2 、R 3 And Wen Qiumian to give the following formula:
in the formula (1), ρ x 、ρ y 、ρ z The attenuation rates of the temperature along the linearity of each direction are respectively shown in the unit of ℃/m;
s32 with sensor a ijk Adjacent sensors a in the same column and layer as the center i-1jk 、a ijk 、a i+1jk Temperature change is affected by ρ x Due to the outward diffusion of temperature around the heating target point, sensor a i-1jk And a i+1jk The middle temperature is far away from the target heating point and is connected with the sensor a ijk At the same time, the two are positioned on the same side of the conduction region of the heating point, and the temperature between the two is linearly attenuated along the x direction, so that the sensor a is selected i-1jk And a i+1jk The lower temperature K shown in x To calculate the basis, the attenuation rate ρ along the x-axis is obtained according to the following equation x ：
Similarly, with sensor a ijk Selecting adjacent sensors a in the same row and layer as the center ij-1k 、a ij+1k The medium temperature shows a smaller temperature K y For calculation basis, the attenuation rate ρ along the y-axis can be obtained y ：
Similarly, with sensor a ijk Selecting adjacent sensors a in the same row and column as the center ijk-1 、a ijk+1 The medium temperature shows a smaller temperature K z For calculation basis, the attenuation rate ρ along the z-axis can be obtained z ：
S33 reporting the highest abnormal temperature sensor a ijk The point position coordinates of the same plane as the target heating point O, which are right under or right above the point position, are O' (x) 1 、y 1 C) defining the temperature at this point as K' 0 The relative positional relationship and the temperature conduction relationship can be known:
s34 target heating spot O (a, b, c) conducting temperature to O' (x) 1 、y 1 C) satisfying the following formula:
the coordinate positions of the target heating points can be calculated and obtained by combining the formulas (1), (2), (3), (4), (5) and (6) to be O (a, b and c), and then the target point positions of the heating grains can be accurately obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110290463.6A CN113008413B (en) | 2021-03-18 | 2021-03-18 | Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110290463.6A CN113008413B (en) | 2021-03-18 | 2021-03-18 | Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113008413A CN113008413A (en) | 2021-06-22 |
| CN113008413B true CN113008413B (en) | 2023-10-10 |
Family
ID=76409635
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110290463.6A Active CN113008413B (en) | 2021-03-18 | 2021-03-18 | Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113008413B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101848246A (en) * | 2010-07-05 | 2010-09-29 | 河南工业大学 | Grain cabin capacity information wireless monitoring system, monitoring method and networking method |
| CN102752578A (en) * | 2012-06-05 | 2012-10-24 | 深圳市粮食集团有限公司 | Unusual monitoring system and method for grain storage |
| CN104132735A (en) * | 2014-06-27 | 2014-11-05 | 河海大学常州校区 | Indoor-temperature 3D visualized monitoring system |
| CN204255890U (en) * | 2014-12-12 | 2015-04-08 | 河南工业大学 | In a kind of novel grain heap, temperature transmits test unit |
| CN108805501A (en) * | 2018-06-07 | 2018-11-13 | 吉林大学 | A kind of granary storage state monitoring and managing method based on correlation statistics |
| CN110810027A (en) * | 2019-12-02 | 2020-02-21 | 吉林大学 | Grain temperature-based circulation temperature control method in granary |
-
2021
- 2021-03-18 CN CN202110290463.6A patent/CN113008413B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101848246A (en) * | 2010-07-05 | 2010-09-29 | 河南工业大学 | Grain cabin capacity information wireless monitoring system, monitoring method and networking method |
| CN102752578A (en) * | 2012-06-05 | 2012-10-24 | 深圳市粮食集团有限公司 | Unusual monitoring system and method for grain storage |
| CN104132735A (en) * | 2014-06-27 | 2014-11-05 | 河海大学常州校区 | Indoor-temperature 3D visualized monitoring system |
| CN204255890U (en) * | 2014-12-12 | 2015-04-08 | 河南工业大学 | In a kind of novel grain heap, temperature transmits test unit |
| CN108805501A (en) * | 2018-06-07 | 2018-11-13 | 吉林大学 | A kind of granary storage state monitoring and managing method based on correlation statistics |
| CN110810027A (en) * | 2019-12-02 | 2020-02-21 | 吉林大学 | Grain temperature-based circulation temperature control method in granary |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113008413A (en) | 2021-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lee et al. | Optimal sensor placement for monitoring and controlling greenhouse internal environments | |
| CN108061877B (en) | A kind of passive track-corelation direction cross positioning method based on angle information | |
| CN108805501B (en) | Granary grain storage state supervision method based on correlation statistics | |
| CN101587182A (en) | Locating method for RFID indoor locating system | |
| US20070032979A1 (en) | Method and apparatus for three-dimensional measurements | |
| CN113008413B (en) | Method for accurately reversely solving target point position of heating grain by utilizing existing temperature measuring network | |
| CN103512662B (en) | Temperature detection range determining device and method | |
| CN109000810A (en) | A kind of multidimensional grain temperature display systems and method | |
| CN110810027B (en) | Grain temperature-based circulation temperature control method in granary | |
| CN113408206B (en) | Indoor natural illuminance modeling method | |
| CN113063524B (en) | Method for constructing temperature field of high-temperature abnormal region of horizontal warehouse by using limited temperature measuring points | |
| CN108802860B (en) | Data correction method and data correction device | |
| CN107014329A (en) | 3D Calibration Fields for portable articulated coordinate machine error calibration | |
| CN107340614A (en) | Detection or correction position system and its method | |
| CN111579086A (en) | Remote infrared temperature measurement precision correction method based on distance compensation | |
| CN106937298B (en) | A kind of improved wireless sensor network 3-D positioning method | |
| CN104684083A (en) | AP (access point) selecting method based on clustering idea | |
| CN104640206B (en) | A kind of multipoint positioning method of centering type star cloth station system | |
| CN116125370A (en) | Multi-platform direction finding rapid fusion positioning method | |
| CN108307399A (en) | Mangrove growing environment wireless monitor system | |
| CN110631721B (en) | Granary heat insulation judgment method based on grain condition big data | |
| CN106204525B (en) | The detection unified approach of tree crown Light interception-rate | |
| CN208109039U (en) | A kind of flatness checking device of liquid crystal display panel | |
| CN111412888A (en) | Building size measuring system and measuring method | |
| Shi et al. | Secure Bearing-Based Target Localization for Multi-Agent Networks Against Malicious Agents |
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 |