CN109165480B - Method and system for predicting pit water inflow based on clay stress path - Google Patents
Method and system for predicting pit water inflow based on clay stress path Download PDFInfo
- Publication number
- CN109165480B CN109165480B CN201811297493.4A CN201811297493A CN109165480B CN 109165480 B CN109165480 B CN 109165480B CN 201811297493 A CN201811297493 A CN 201811297493A CN 109165480 B CN109165480 B CN 109165480B
- Authority
- CN
- China
- Prior art keywords
- parameters
- clay
- model
- dimensional
- water inflow
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 68
- 239000004927 clay Substances 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 23
- 239000002689 soil Substances 0.000 claims abstract description 38
- 230000006835 compression Effects 0.000 claims abstract description 37
- 238000007906 compression Methods 0.000 claims abstract description 37
- 238000007596 consolidation process Methods 0.000 claims abstract description 27
- 230000008878 coupling Effects 0.000 claims abstract description 6
- 238000010168 coupling process Methods 0.000 claims abstract description 6
- 238000005859 coupling reaction Methods 0.000 claims abstract description 6
- 238000012544 monitoring process Methods 0.000 claims description 10
- 238000004062 sedimentation Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 abstract description 6
- 230000002265 prevention Effects 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The embodiment of the invention discloses a method and a system for predicting water inflow of a pit based on a clay stress path, wherein the method comprises the following steps: establishing a clay consolidation model according to pre-acquired clay compression parameters, wherein the clay compression parameters are parameters obtained after the clay in the range of a preset project area is pre-treated; establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area; coupling the cohesive soil consolidation model and the three-dimensional hydrogeologic model, and predicting the water inflow of the pit within the range of the preset project area. The method not only considers the water inflow situation of the pit predicted by the traditional technology, but also breaks the water balance of the viscous soil layer due to external pressure or external reasons, and releases water flow. The water inflow of the pit predicted by the traditional technology and the water flow released by the clay layer are taken into consideration at the same time, so that the water inflow of the pit can be predicted more accurately based on the stress path of the clay, and a reliable basis is provided for mine water prevention and control.
Description
Technical Field
The invention relates to the technical field of geological environment prediction, in particular to a method and a system for predicting pit water inflow based on a cohesive soil stress path.
Background
Because the clay has good water-resisting property, the clay is treated according to the water-resisting layer in the process of hydrogeological exploration of mines and prediction of water inflow of pits in general. In the practical research, the property of the clay soil for consolidation and water release under certain conditions is ignored by both the prior textbook and the specification and the technical staff, so that a large error exists in the prediction of the water inflow of the pit, and even the prediction error is directly caused.
Disclosure of Invention
The embodiment of the invention aims to provide a method and a system for predicting pit water inflow based on a clay stress path, which are used for solving the problems that the existing prediction of the pit water inflow has errors and even the prediction is wrong.
In order to achieve the above object, an embodiment of the present invention provides a method for predicting pit water inflow based on a clay stress path, the method comprising:
establishing a clay consolidation model according to pre-acquired clay compression parameters, wherein the clay compression parameters are parameters obtained after the clay in the range of a preset project area is pre-treated;
establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area;
coupling the cohesive soil consolidation model and the three-dimensional hydrogeologic model, and predicting the water inflow of the pit within the range of the preset project area.
The embodiment of the invention has the following advantages: and establishing a cohesive soil consolidation model according to the pre-acquired cohesive soil compression parameters. The clay consolidation model can simulate the water release amount of clay under different time or different pressure conditions. And a three-dimensional hydrogeologic model can be established according to geological parameters in the range of the preset project area, and the water flow condition when the viscous soil is assumed to be used as the water-resisting layer and does not release water flow can be predicted through the model. The two models are coupled, so that the water release condition of the clay soil and the common geology under different time periods and different pressures can be seen, and the total water inflow of the pit can be accurately predicted. The method is convenient for subsequent staff to take effective measures to be applied to the aspect of water control, and personal safety is ensured.
To achieve the above object, an embodiment of the present invention provides a system for predicting pit water inflow based on a clay stress path, the system comprising:
the processing unit is used for establishing a cohesive soil consolidation model according to pre-acquired cohesive soil compression parameters, wherein the cohesive soil compression parameters are parameters acquired after the cohesive soil in the preset project area is pre-processed;
establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area;
and the prediction unit is used for coupling the cohesive soil consolidation model and the three-dimensional hydrogeologic model and predicting the water inflow of the pit within the range of the preset project area.
The embodiment of the invention has the following advantages: and establishing a cohesive soil consolidation model according to the pre-acquired cohesive soil compression parameters. The clay consolidation model can simulate the water release amount of clay under different time or different pressure conditions. And a three-dimensional hydrogeologic model can be established according to geological parameters in the range of the preset project area, and the water flow condition when the viscous soil is assumed to be used as the water-resisting layer and does not release water flow can be predicted through the model. The two models are coupled, so that the water release condition of the clay soil and the common geology under different time periods and different pressures can be seen, and the total water inflow of the pit can be accurately predicted. The method is convenient for subsequent staff to take effective measures to be applied to the aspect of water control, and personal safety is ensured.
Drawings
FIG. 1 is a schematic flow chart of a method for predicting pit water inflow based on a clay stress path according to embodiment 1 of the present invention;
fig. 2 is a schematic diagram of a system for predicting pit water inflow based on a clay stress path according to embodiment 2 of the present invention.
Detailed Description
Further advantages and effects of the present invention will become apparent to those skilled in the art from the disclosure of the present invention, which is described by the following specific examples.
It should be understood that the structures, proportions, sizes, etc. shown in the drawings are for illustration purposes only and should not be construed as limiting the invention to the extent that it can be practiced, since modifications, changes in the proportions, or otherwise, used in the practice of the invention, are not intended to be critical to the essential characteristics of the invention, but are intended to fall within the spirit and scope of the invention. Also, the terms such as "upper", "lower", "left", "right", "middle", etc. are used herein for descriptive purposes only and are not intended to limit the scope of the invention for which the invention may be practiced or for which the relative relationships may be altered or modified without materially altering the technical context.
Example 1
The embodiment 1 of the invention provides a method for predicting pit water inflow based on a clay stress path, and specifically as shown in fig. 1, fig. 1 is a flow chart of a method for predicting pit water inflow based on a clay stress path, which is provided by the embodiment of the invention, and comprises the following steps:
and 110, building a cohesive soil consolidation model according to the pre-acquired cohesive soil compression parameters.
The clay compression parameters are parameters obtained by preprocessing clay in a preset project area range.
Alternatively, the pre-acquired parameters may include at least a compression coefficient and a compression modulus. In a specific example, the clay compaction parameters may be obtained by a high pressure consolidation test.
The high-pressure consolidation test is a test method for measuring the drainage, exhaust and bubble compression properties of a soil body under the action of external force. In general, the soil body is subjected to three principal stresses. The purpose of the test was to determine the relationship between deformation and pressure, deformation and time, void ratio and pressure of the test sample under confined and axial drainage conditions to draw a compression curve. Meanwhile, the compression coefficient and the compression modulus can be obtained. And establishing a cohesive soil consolidation model according to the compression coefficient and the compression model. And drawing a compression curve, and observing deformation conditions of the clay in different time periods. The deformation of clay under different pressures, the void ratio between clay particles under different pressures, etc. And establishing a clay consolidation model according to the pre-acquired clay compression parameters, wherein the clay compression parameters are substituted into preset software to establish the clay consolidation model. I.e. the clay consolidation model has an initial morphology, and these parameters need to be substituted into the initial model, such as the Feflow software.
And 120, establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the pre-acquired preset project area.
Optionally, the staff can obtain the geological parameters according to the existing or newly added drilling construction, geophysical prospecting survey and other means. The geologic parameters include engineering geologic parameters and hydrogeologic parameters. Such as formation lithology, surface elevation, formation thickness, permeability coefficient, and the like.
And the three-dimensional hydrogeologic model can be built by using professional simulation software such as Visual Modflow or Feflow. That is, the geologic parameters are input into the initial three-dimensional hydrogeologic model of the corresponding simulation software, and a three-dimensional hydrogeologic model actually corresponding to the present embodiment is established.
And 130, coupling the cohesive soil consolidation model and the three-dimensional hydrogeological model, and predicting the pit water inflow in the range of the preset project area.
And establishing a cohesive soil consolidation model according to the pre-acquired cohesive soil compression parameters. The clay consolidation model can simulate the water release amount of clay under different time or different pressure conditions. And a three-dimensional hydrogeologic model can be established according to geological parameters in the range of the preset project area, and the water flow condition when the viscous soil is assumed to be used as the water-resisting layer and does not release water flow can be predicted through the model. The two models are coupled, so that the water release condition of the clay soil and the common geology under different time periods and different pressures can be seen, and the total water inflow of the pit can be accurately predicted.
Optionally, the present embodiment may further include: monitoring different stratum within the range of the preset project area, and respectively determining the sedimentation degree of the different stratum;
and when the sedimentation degree of any stratum in different strata exceeds a preset threshold, sending out alarm information.
The step of monitoring may be performed at a fixed time period before, during, or after steps 110-130, or may be real-time monitoring. The specific setting rule can be set according to actual conditions. The method is mainly used for detecting different stratum in the range of the preset project area and determining the sedimentation degree of the different stratum. The concrete way to be adopted is to install a settlement meter in different strata to monitor the settlement degree. Meanwhile, non-settlement monitoring points can be respectively set outside the project area range to play a role in reference. The project area range can be understood as subsidence area monitoring points, ground elevation of each monitoring point is measured regularly, and the subsidence degree of the ground is monitored.
Once it is determined that the degree of settlement of any one of the different formations exceeds a preset threshold, an alarm message needs to be sent. To prompt the staff to take safety precautions in time. Unnecessary loss is avoided, and even casualties occur.
The method for predicting the pit water inflow based on the cohesive soil stress path provided by the embodiment of the invention not only considers the situation of the pit water inflow predicted by the traditional technology, but also releases water flow after the cohesive soil layer breaks the water balance due to external pressure or external reasons. The water inflow of the pit predicted by the traditional technology and the water flow released by the clay layer are taken into consideration at the same time, so that the water inflow of the pit can be predicted more accurately based on the stress path of the clay, and a reliable basis is provided for mine water prevention and control.
Example 2
Corresponding to the above embodiment 1, embodiment 2 of the present invention provides a system for predicting pit water inflow based on a clay stress path, specifically as shown in fig. 2, fig. 2 is a schematic structural diagram of a system for predicting pit water inflow based on a clay stress path, where the system includes: a processing unit 201 and a prediction unit 202.
The processing unit 201 is configured to establish a clay consolidation model according to pre-acquired clay compression parameters, where the clay compression parameters are parameters acquired after the clay within a preset project area is pre-processed;
establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area;
and a prediction unit 202, configured to couple the clay consolidation model and the three-dimensional hydrogeologic model, and predict a pit water inflow within a range of a preset project area.
Optionally, the clay compression parameters include at least: compression coefficient and compression modulus.
Optionally, the geological parameters include at least: engineering geologic parameters and hydrogeologic parameters; the processing unit 201 is specifically configured to substitute the engineering geological parameter and the hydrogeological parameter into preset three-dimensional model building software to build a three-dimensional hydrogeological model.
Optionally, the system further comprises: the monitoring unit 203 is configured to monitor different strata within a preset project area, and determine different stratum sedimentation degrees respectively;
and the alarm unit 204 is used for sending alarm information when the sedimentation degree of any stratum in different strata exceeds a preset threshold value.
The functions performed by each component in the system for predicting the water inflow of a pit based on the clay stress path provided in the embodiment of the present invention are described in detail in the method for predicting the water inflow of a pit based on the clay stress path provided in embodiment 1, and will not be described in detail here.
The system for predicting the pit water inflow based on the cohesive soil stress path provided by the embodiment of the invention not only considers the situation of the pit water inflow predicted by the traditional technology, but also releases water flow after the cohesive soil layer breaks the water balance due to external pressure or external reasons. The water inflow of the pit predicted by the traditional technology and the water flow released by the clay layer are taken into consideration at the same time, so that the water inflow of the pit can be predicted more accurately based on the stress path of the clay, and a reliable basis is provided for mine water prevention and control.
While the invention has been described in detail in the foregoing general description and specific examples, it will be apparent to those skilled in the art that modifications and improvements can be made thereto. Accordingly, such modifications or improvements may be made without departing from the spirit of the invention and are intended to be within the scope of the invention as claimed.
Claims (4)
1. A method of predicting pit water inflow based on a clay stress path, the method comprising:
establishing a clay consolidation model according to pre-acquired clay compression parameters, wherein the clay compression parameters are parameters obtained by pre-treating clay in a preset project area range;
establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area;
coupling the cohesive soil consolidation model and the three-dimensional hydrogeologic model, and predicting pit water inflow in the range of the preset project area;
the clay compression parameters at least comprise: compression coefficient and compression modulus;
the geological parameters at least comprise: engineering geologic parameters and hydrogeologic parameters; the method for establishing the three-dimensional hydrogeologic model according to the geological parameters in the range of the pre-acquired preset project area specifically comprises the following steps:
substituting the engineering geological parameters and the hydrogeological parameters into preset three-dimensional model building software to build a three-dimensional hydrogeological model.
2. The method according to claim 1, wherein the method further comprises: monitoring different stratum within the preset project area range, and respectively determining the sedimentation degree of the different stratum;
and when the sedimentation degree of any stratum in the different strata exceeds a preset threshold, sending out alarm information.
3. A system for predicting pit water inflow based on a clay stress path, the system comprising:
the processing unit is used for establishing a cohesive soil consolidation model according to pre-acquired cohesive soil compression parameters, wherein the cohesive soil compression parameters are parameters acquired after the cohesive soil in the preset project area is pre-processed;
establishing a three-dimensional hydrogeologic model according to the geological parameters in the range of the preset project area;
the prediction unit is used for coupling the cohesive soil consolidation model and the three-dimensional hydrogeologic model and predicting the pit water inflow in the range of the preset project area;
the clay compression parameters at least comprise: compression coefficient and compression modulus;
the geological parameters at least comprise: engineering geologic parameters and hydrogeologic parameters; the processing unit is specifically used for substituting the engineering geological parameters and the hydrogeological parameters into preset three-dimensional model building software to build a three-dimensional hydrogeological model.
4. A system according to claim 3, wherein the system further comprises: the monitoring unit is used for monitoring different stratum in the preset project area range and respectively determining the sedimentation degree of the different stratum;
and the alarm unit is used for sending alarm information when determining that the sedimentation degree of any stratum in different strata exceeds a preset threshold value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811297493.4A CN109165480B (en) | 2018-11-01 | 2018-11-01 | Method and system for predicting pit water inflow based on clay stress path |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811297493.4A CN109165480B (en) | 2018-11-01 | 2018-11-01 | Method and system for predicting pit water inflow based on clay stress path |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109165480A CN109165480A (en) | 2019-01-08 |
| CN109165480B true CN109165480B (en) | 2023-08-11 |
Family
ID=64876650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811297493.4A Active CN109165480B (en) | 2018-11-01 | 2018-11-01 | Method and system for predicting pit water inflow based on clay stress path |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109165480B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010047859A1 (en) * | 2008-10-20 | 2010-04-29 | Exxonmobil Upstream Research Company | Method for modeling deformation in subsurface strata |
| CN107480349A (en) * | 2017-07-26 | 2017-12-15 | 中国铁路总公司 | Surface subsidence Forecasting Methodology along high ferro based on three-dimensional geological model and back analysis |
| CN108643246A (en) * | 2018-04-02 | 2018-10-12 | 上海交通大学 | The prediction technique that artesian water precipitation influences surrounding enviroment in foundation pit construction |
-
2018
- 2018-11-01 CN CN201811297493.4A patent/CN109165480B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010047859A1 (en) * | 2008-10-20 | 2010-04-29 | Exxonmobil Upstream Research Company | Method for modeling deformation in subsurface strata |
| CN107480349A (en) * | 2017-07-26 | 2017-12-15 | 中国铁路总公司 | Surface subsidence Forecasting Methodology along high ferro based on three-dimensional geological model and back analysis |
| CN108643246A (en) * | 2018-04-02 | 2018-10-12 | 上海交通大学 | The prediction technique that artesian water precipitation influences surrounding enviroment in foundation pit construction |
Non-Patent Citations (1)
| Title |
|---|
| 三维可视化建模与矿坑涌水量预测研究;吴烨等;《四川地质学报》;20100331(第01期);第61-64页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109165480A (en) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Long-term performance of metro tunnels induced by adjacent large deep excavation and protective measures in Nanjing silty clay | |
| Robertson | The James K. Mitchell Lecture: Interpretation of in-situ tests–some insights | |
| WO2017035371A1 (en) | Method for estimating stress magnitude | |
| KR101096817B1 (en) | Method for managing knowledge-based injecting construction using ags and gin | |
| WO2021136077A1 (en) | Analysis system and method for designing and correcting water curtain system of underground water-sealed oil storage | |
| CN107478357A (en) | Stress field of the surrounding rock crack field Integrated Monitoring System and quantitatively determine method | |
| Meng et al. | In situ investigation and numerical simulation of the failure depth of an inclined coal seam floor: a case study | |
| Xu et al. | Predicting the excavation damaged zone within brittle surrounding rock masses of deep underground caverns using a comprehensive approach integrating in situ measurements and numerical analysis | |
| CN108952657A (en) | A kind of horizontal well platform fracturing fracture length determining method and device | |
| CN116564065A (en) | Mine disaster prediction and early warning method combining on-site monitoring and numerical simulation | |
| CN104727815A (en) | Real-time well drilling formation correction method and device | |
| CN109165480B (en) | Method and system for predicting pit water inflow based on clay stress path | |
| Zhu et al. | Calculation model of overburden rock failure volume in mined-out area based on Brillouin optical time-domain analysis technology | |
| CN109709609B (en) | Cement consolidation grouting quality inspection and evaluation method based on rock mass quality index | |
| Farkas et al. | Effect of foliation and fluid viscosity on hydraulic fracturing tests in mica schists investigated using distinct element modeling and field data | |
| CN111680896A (en) | Coal mine underground reservoir safety distance determination method | |
| Wen et al. | The use of Hoek Brown failure criterion on determination of the geo-mechanical parameters of a grouting consolidation body | |
| Wang et al. | Robust optimization design of bolt-shotcrete support structure in tunnel | |
| CN112946778A (en) | Method for early warning karst collapse based on underground water turbidity monitoring | |
| EP3341562A1 (en) | Method for estimating stress magnitude | |
| CN114647913A (en) | Maximum level principal ground stress acquisition method based on Mohr-Coulomb criterion | |
| Yi et al. | Application of Formation Pressure While Drilling Prediction Technology in Offshore HTHP Wells | |
| Zhang et al. | Detection and evaluation of crack development near the fault zone under the influence of coal mining | |
| CN110439614B (en) | Coal mining bottom plate dredging water simulation device and use method | |
| Yu et al. | On how water injection may induce fault reactivation and slippage: A numerical method |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Zhixue Inventor after: Li Yuanhang Inventor after: Li Hongzhao Inventor after: Zhu Xiaoqin Inventor after: Guo Qiushan Inventor after: Liu Xiaolong Inventor after: Wang Guodong Inventor after: Yan Pengyang Inventor after: Hao Yunyang Inventor after: Zhang Tao Inventor after: Yuan Bing Inventor after: Sun Binchao Inventor after: Liu Zhou Inventor before: Li Zhixue Inventor before: Liu Zhou Inventor before: Ji Lipu Inventor before: Li Hongzhao Inventor before: Zhu Xiaoqin Inventor before: Guo Qiushan Inventor before: Liu Xiaolong Inventor before: Wang Guodong Inventor before: Yan Pengyang Inventor before: Hao Yunyang Inventor before: Xi Dapeng Inventor before: Zhang Tao Inventor before: Yuan Bing Inventor before: Sun Binchao |