CN118013895A - 一种水库深孔前淤积厚度的计算方法 - Google Patents
一种水库深孔前淤积厚度的计算方法 Download PDFInfo
- Publication number
- CN118013895A CN118013895A CN202410425704.7A CN202410425704A CN118013895A CN 118013895 A CN118013895 A CN 118013895A CN 202410425704 A CN202410425704 A CN 202410425704A CN 118013895 A CN118013895 A CN 118013895A
- Authority
- CN
- China
- Prior art keywords
- section
- formula
- reservoir
- deep hole
- sand
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 35
- 238000004062 sedimentation Methods 0.000 title claims abstract description 19
- 239000004576 sand Substances 0.000 claims abstract description 95
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 47
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 17
- 239000013049 sediment Substances 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 22
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 239000002245 particle Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 3
- 230000000903 blocking effect Effects 0.000 description 3
- 230000008021 deposition Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 239000010802 sludge Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水库深孔前淤积厚度的计算方法，属于水利工程领域中的水库调度研究技术领域，设定水库上游进口处为第一断面，下游出口处为第二断面，根据水流挟沙力的影响因素建立进口第一断面和出口第二断面的水流挟沙力关系式；选用合适的水流挟沙力计算公式，建立出口处断面含沙量计算公式；结合河床变形方程与含沙量计算公式提出一种水库深孔前淤积厚度计算公式。本发明建立的含沙量计算模型可以较好地反映挟沙力的沿程变化规律，综合河床变形公式得到了水库深孔前沿淤积厚度的计算公式，弥补了尚无有效预测坝前深孔泥沙淤积厚度的不足。
Description
技术领域
本发明涉及一种水库深孔前淤积厚度的计算方法，属于水利工程领域中的水库调度研究技术领域。
背景技术
在现代水利工程中，由于各地区水资源分布与需求存在较大差异，所以常常因地制宜选用不同形式的水工建筑物来调蓄水流以满足需要。水库因其具有防洪、蓄水、发电、养殖等功能在水利工程中得到了广泛的应用，水库作为主要的防洪措施可以拦蓄洪水，削减洪峰流量，但多沙水流带来的大量泥沙受拦水坝和库尾回水的影响常常在水库尾部淤积，水库库容减小，拦蓄洪水的能力下降，水库使用寿命缩减。在实际工程中，为了缓解水库淤积的情况常常采用蓄清排浊、人工清淤等手段进行干涉以保持库容，但在此之前必须对库区泥沙淤积情况进行分析，预测库区泥沙分布与含沙量差异，确定泥沙淤积强度，为水库水沙调度工作提供理论支持。
目前，国内外学者针对水库坝前深孔淤积的研究多基于模型试验或现场测量，较少有合适的方法能快速计算淤积深度，给水库运用的效果带来不利影响。
因此，需要有一种水库深孔前淤积厚度的计算方法，基于河床变形方程预测深孔前的泥沙淤积厚度，有效弥补了关于深孔前泥沙淤积厚度尚无有效计算方法的不足。
发明内容
本发明要解决的技术问题是：为了克服现有技术的不足，提供以挟沙力的影响因素为出发点，结合泥沙连续方程推导出可以较好地反映挟沙力沿程变化的含沙量计算公式，将含沙量计算公式与河床变形方程结合的一种水库深孔前淤积厚度的计算方法。
本发明解决上述问题所采用的技术方案为：一种水库深孔前淤积厚度的计算方法，包括如下步骤：
步骤1，设定水库上游进口处为第一断面，第一断面即为进口断面，下游深孔位置为第二断面，第二断面即为出口断面，测量进口断面至深孔位置间的距离为L；
步骤2，统计水库河床坡度J、进口断面单宽流量q、进口含沙量；
步骤3，测量入库悬移质泥沙的中值粒径D 50，计算该粒径的沉速；
步骤4，测量进口断面的平均水深；
步骤5，根据一维泥沙连续方程推导水库下游深孔位置断面的含沙量公式；
具体推导方式如下：
公式一、一维泥沙连续方程式：
公式一中，x为从上游到下游方向上的距离；为任意位置的含沙量；/>为任意位置的水流挟沙力；/>为恢复饱和系数；/>为任意位置的断面单宽流量，水库上游和下游流量恒定，即得出公式二：
则公式一转换为公式三：
对公式三沿x方向积分，积分范围为进口断面至下游深孔位置断面处，距离为L，得出公式四：
公式四进一步化简得到公式五：
公式五中，为进口断面含沙量，/>为进口断面的水流挟沙力，/>为深孔位置断面的含沙量、/>为深孔位置断面的水流挟沙力；
水库河床坡度为J，则进口断面至下游距离为x处断面的水深公式为公式六：
公式六中，为进口断面的水深；
由于上游和下游流量恒定，则进口断面至下游距离为x处断面位置挟沙力与进口断面挟沙力的关系公式为公式七：
公式七转换为公式八，得出：
将公式五中按照泰勒展开变为：
约去高阶项，简化为公式十：
将公式八和公式十代入公式五，然后进行积分，得到公式十一：
公式十一即为水库下游深孔位置断面含沙量的计算公式；
步骤6，计算含沙量的一阶衰减量：
其中，为水库进口含沙量的一阶衰减量，/>为进口断面的平均水深，单位为m；为恢复饱和系数；/>为悬移质泥沙沉速，单位为m/s；L为进口断面至深孔位置间的距离，单位为m；q为进口断面单宽流量，单位为m2/s；J为水库河床坡度；
步骤7，计算含沙量的二阶衰减量：
其中，为水库进口含沙量的二阶衰减量，其他参数含义同前；
步骤8，计算含沙量的三阶衰减量：
其中，为水库进口含沙量的三阶衰减量，其他参数含义同前；
步骤9，计算含沙量的局部增量：
其中，为含沙量的局部增量，其他参数含义同前；
步骤10，将步骤6、步骤7、步骤8和步骤9的结果代入公式十一，可以得到第二断面（坝前深孔位置断面）的平均含沙量：
步骤11，计算进口断面和深孔断面的水流挟沙力；
步骤12，结合河床变形方程与步骤10的含沙量计算公式，得到水库深孔前淤积厚度计算公式；
式中，为水库深孔前淤积厚度，单位m；
进一步的，所述步骤3泥沙沉速的计算方法为：
；
其中，为悬移质泥沙沉速，单位为m/s； D 50为悬沙泥沙的粒径，单位为m；/>为运动粘性系数，取0.000001，单位为m2/s；
进一步的，所述恢复饱和系数：冲刷计算取1.0，淤积计算取0.25；
进一步的，所述步骤10中进口断面和深孔断面的水流挟沙力计算方法为：
其中，为出口断面的水流挟沙力，单位为kg/m3；V断面平均流速，单位m/s；g为重力加速度，单位为m/s2；/>为悬移质泥沙沉速，单位为m/s；h为断面平均水深，单位为m。
与现有技术相比，本发明的优点在于：
本发明一种水库深孔前淤积厚度的计算方法，考虑到传统方法计算含沙量时只能反映较短距离内的含沙量变化情况，河道较长时需要分多段计算，连续性较差，所以本发明充分考虑了河道长度、比降等因素，以泥沙连续方程为基础，结合了河床变形公式，本发明可以较好地反映长距离水库深孔前淤积厚度情况，连续性得到增强；同时本发明根据不同的水库特征可以灵活修改相关参数，实用性较强，可以适用不同情况的水库深孔前淤积厚度预测。
附图说明
图1为本发明一种水库深孔前淤积厚度的计算方法的流程图。
具体实施方式
一种水库深孔前淤积厚度的计算方法，包括如下步骤：
步骤1，设定水库上游进口处为第一断面，第一断面即为进口断面，下游深孔位置为第二断面，第二断面即为出口断面，测量进口断面至深孔位置间的距离为L；
步骤2，统计水库河床坡度J、进口断面单宽流量q、进口含沙量；
步骤3，测量入库悬移质泥沙的中值粒径D 50，计算该粒径的沉速；
；
步骤4，测量进口断面的平均水深；
步骤5，根据一维泥沙连续方程推导水库下游深孔位置断面的含沙量公式；
具体推导方式如下：
公式一、一维泥沙连续方程式：
公式一中，x为从上游到下游方向上的距离；为任意位置的含沙量；/>为任意位置的水流挟沙力；/>为恢复饱和系数；/>为任意位置的断面单宽流量，水库上游和下游流量恒定，即得出公式二：
则公式一转换为公式三：
对公式三沿x方向积分，积分范围为进口断面至下游深孔位置断面处，距离为L，得出公式四：
公式四进一步化简得到公式五：；
公式五中，为进口断面含沙量，/>为进口断面的水流挟沙力，/>为深孔位置断面的含沙量、/>为深孔位置断面的水流挟沙力；
水库河床坡度为J，则进口断面至下游距离为x处断面的水深公式为公式六：
公式六中，为进口断面的水深；
由于上游和下游流量恒定，则进口断面至下游距离为x处断面位置挟沙力与进口断面挟沙力的关系公式为公式七：
公式七转换为公式八，得出：
将公式五中按照泰勒展开可变为：
约去高阶项，可简化为公式十：
将公式八和公式十代入公式五，然后进行积分，得到公式十一：
公式十一即为水库下游深孔位置断面含沙量的计算公式；
步骤6，计算含沙量的一阶衰减量：
其中，为水库进口含沙量的一阶衰减量，/>为进口断面的平均水深，单位为m；为恢复饱和系数；/>为悬移质泥沙沉速，单位为m/s；L为进口断面至深孔位置间的距离，单位为m；q为进口断面单宽流量，单位为m2/s；J为水库河床坡度；
步骤7，计算含沙量的二阶衰减量：
其中，为水库进口含沙量的二阶衰减量，其他参数含义同前；
步骤8，计算含沙量的三阶衰减量：
其中，为水库进口含沙量的三阶衰减量，其他参数含义同前；
步骤9，计算含沙量的局部增量：
其中，为含沙量的局部增量，其他参数含义同前；
步骤10，采用下式计算第二断面（坝前深孔位置断面）的平均含沙量：
步骤11，采用下式计算进口、出口断面水流挟沙力：
步骤12，结合河床变形方程与步骤9的含沙量计算公式，计算水库深孔前淤积厚度：
实施例1，如图1所示，本发明水库深孔前淤积厚度的计算方法的流程如下：
步骤1，设定水库上游进口处为第一断面，第一断面即为进口断面，下游出口处为第二断面，第二断面即为出口断面，测量进口断面至深孔位置间的距离为L=110000m；
步骤2，统计水库河床坡度J为0.00108、进口河道单宽流量q为100m2/s、进口含沙量为10kg/m3；
步骤3，测量得到入库悬移质泥沙的中值粒径D 50=0.00003m，计算得到泥沙沉速为：0.00058m/s；
步骤4，测量进口断面的平均水深为25m；
步骤5，根据一维泥沙连续方程推导水库下游深孔位置断面的含沙量公式；
具体推导方式如下：
公式一、一维泥沙连续方程式：
公式一中，x为从上游到下游方向上的距离；为任意位置的含沙量；/>为任意位置的水流挟沙力；/>为恢复饱和系数；/>为任意位置的断面单宽流量，水库上游和下游流量恒定，即得出公式二：
则公式一转换为公式三：
对公式三沿x方向积分，积分范围为进口断面至下游深孔位置断面处，距离为L，得出公式四：
公式四进一步化简得到公式五：
公式五中，为进口断面含沙量，/>为进口断面的水流挟沙力，/>为深孔位置断面的含沙量、/>为深孔位置断面的水流挟沙力；
水库河床坡度为J，则进口断面至下游距离为x处断面的水深公式为公式六：
公式六中，为进口断面的水深；
由于上游和下游流量恒定，则进口断面至下游距离为x处断面位置挟沙力与进口断面挟沙力的关系公式为公式七：
公式七转换为公式八，得出：
将公式五中按照泰勒展开变为：
约去高阶项，简化为公式十：
将公式八和公式十代入公式五，然后进行积分，得到公式十一：
公式十一即为水库下游深孔位置断面含沙量的计算公式；
步骤6，计算含沙量的一阶衰减量：
恢复饱和系数，为淤积条件，取值为0.25；计算得/>=0.001071；
步骤7，计算含沙量的二阶衰减量：
恢复饱和系数取值为0.25；计算得/>=0.0000682；
步骤8，计算含沙量的三阶衰减量：
恢复饱和系数取值为0.25；计算得/>=0.0000057；
步骤9，计算含沙量的局部增量：
恢复饱和系数取值为0.25；计算得/>=1.011376；
步骤10，采用下式计算第二断面的平均含沙量：
经过计算得到=8.70kg⁄m3；
步骤11，计算进口、出口断面水流挟沙力分别为18.01kg⁄m3和0.016 kg⁄m3；
步骤12，采用水库深孔前淤积厚度计算公式计算得水库深孔前3个月淤积强度=3.6m。
综上所述，考虑到传统方法计算含沙量时只能反映较短距离内的含沙量变化情况，河道较长时需要分多段计算，连续性较差，所以本发明充分考虑了河道长度、比降等因素，以泥沙连续方程为基础，结合了河床变形公式，本发明可以较好地反映长距离水库深孔前淤积厚度情况，连续性得到增强；同时本发明根据不同的水库特征可以灵活修改相关参数，实用性较强，可以适用不同情况的水库深孔前淤积厚度预测。
除上述实施例外，本发明还包括有其他实施方式，凡采用等同变换或者等效替换方式形成的技术方案，均应落入本发明权利要求的保护范围之内。
Claims (5)
1.一种水库深孔前淤积厚度的计算方法，其特征在于：包括如下步骤：
步骤1，设定水库上游进口处为第一断面，第一断面即为进口断面，下游深孔位置为第二断面，第二断面即为出口断面，测量进口断面至深孔位置间的距离为L；
步骤2，统计水库河床坡度J、进口断面单宽流量q、进口含沙量；
步骤3，测量入库悬移质泥沙的中值粒径D 50，计算该粒径的沉速；
步骤4，测量进口断面的平均水深；
步骤5，根据一维泥沙连续方程推导水库下游深孔位置断面的含沙量公式；
步骤6，计算含沙量的一阶衰减量：
；
其中， 为水库进口含沙量的一阶衰减量， />为进口断面的平均水深，单位为m；/>为恢复饱和系数；/>为悬移质泥沙沉速，单位为m/s；L为进口断面至深孔位置间的距离，单位为m；q为进口断面单宽流量，单位为m2/s；J为水库河床坡度；
步骤7，计算含沙量的二阶衰减量：
；
其中，为水库进口含沙量的二阶衰减量，其他参数含义同前；
步骤8，计算含沙量的三阶衰减量：
；
其中，为水库进口含沙量的三阶衰减量，其他参数含义同前；
步骤9，计算含沙量的局部增量：
；
其中，为含沙量的局部增量，其他参数含义同前；
步骤10，采用下式计算第二断面的平均含沙量：
；
步骤11，计算进口断面和深孔断面的水流挟沙力；
步骤12，结合河床变形方程与步骤10的含沙量计算公式，得到水库深孔前淤积厚度计算公式；
；
式中，为水库深孔前淤积厚度，单位m。
2.根据权利要求1所述的一种水库深孔前淤积厚度的计算方法，其特征在于：所述步骤3泥沙沉速的计算方法为：
；
其中，为悬移质泥沙沉速，单位为m/s； D 50为悬沙泥沙的粒径，单位为m；/>为运动粘性系数，取0.000001，单位为m2/s。
3.根据权利要求1所述的一种水库深孔前淤积厚度的计算方法，其特征在于：所述恢复饱和系数，冲刷计算取1.0，淤积计算取0.25。
4.根据权利要求1所述的一种水库深孔前淤积厚度的计算方法，其特征在于：所述步骤9中进口断面和深孔断面的水流挟沙力计算方法为：
；其中，/>为出口断面的水流挟沙力，单位为kg/m3；V断面平均流速，单位m/s；g为重力加速度，单位为m/s2；/>为悬移质泥沙沉速，单位为m/s；h为断面平均水深，单位为m。
5.根据权利要求1所述的一种水库深孔前淤积厚度的计算方法，其特征在于：所述步骤5具体推导方式如下：
公式一、一维泥沙连续方程式：
；
公式一中，x为从上游到下游方向上的距离；为任意位置的含沙量；/>为任意位置的水流挟沙力；/>为恢复饱和系数；/>为任意位置的断面单宽流量，水库上游和下游流量恒定，即得出公式二：
；
则公式一转换为公式三：
；
对公式三沿x方向积分，积分范围为进口断面至下游深孔位置断面处，距离为L，得出公式四：
；
公式四进一步化简得到公式五：
；
公式五中，为进口断面含沙量，/>为进口断面的水流挟沙力，/>为深孔位置断面的含沙量、/>为深孔位置断面的水流挟沙力；
水库河床坡度为J，则进口断面至下游距离为x处断面的水深公式为公式六：
；
公式六中，为进口断面的水深；
由于上游和下游流量恒定，则进口断面至下游距离为x处断面位置挟沙力与进口断面挟沙力的关系公式为公式七：
；
公式七转换为公式八，得出：
；
将公式五中按照泰勒展开变为：
；
约去高阶项，简化为公式十：
；
将公式八和公式十代入公式五，然后进行积分，得到公式十一：
；
公式十一即为水库下游深孔位置断面含沙量的计算公式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410425704.7A CN118013895B (zh) | 2024-04-10 | 2024-04-10 | 一种水库深孔前淤积厚度的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410425704.7A CN118013895B (zh) | 2024-04-10 | 2024-04-10 | 一种水库深孔前淤积厚度的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118013895A true CN118013895A (zh) | 2024-05-10 |
| CN118013895B CN118013895B (zh) | 2024-06-11 |
Family
ID=90960408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410425704.7A Active CN118013895B (zh) | 2024-04-10 | 2024-04-10 | 一种水库深孔前淤积厚度的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118013895B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111241467A (zh) * | 2020-01-07 | 2020-06-05 | 河海大学 | 一种多沙河流水库排沙比的计算方法 |
| CN113591411A (zh) * | 2021-08-09 | 2021-11-02 | 长江水利委员会水文局长江三峡水文水资源勘测局 | 一种水库泥沙淤积预测预报方法 |
| CN114091359A (zh) * | 2022-01-21 | 2022-02-25 | 中国长江三峡集团有限公司 | 一种水库水沙预测模型训练、水库水沙预测方法及装置 |
-
2024
- 2024-04-10 CN CN202410425704.7A patent/CN118013895B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111241467A (zh) * | 2020-01-07 | 2020-06-05 | 河海大学 | 一种多沙河流水库排沙比的计算方法 |
| CN113591411A (zh) * | 2021-08-09 | 2021-11-02 | 长江水利委员会水文局长江三峡水文水资源勘测局 | 一种水库泥沙淤积预测预报方法 |
| CN114091359A (zh) * | 2022-01-21 | 2022-02-25 | 中国长江三峡集团有限公司 | 一种水库水沙预测模型训练、水库水沙预测方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| ZHIJING LI;HONGLU QIAN;ZHIXIAN CAO;HUAIHAN LIU;GARETH PENDER;PENGHUI HU;: "Enhanced bed load sediment transport by unsteady flows in a degrading channel", INTERNATIONAL JOURNAL OF SEDIMENT RESEARCH, no. 03, 15 September 2018 (2018-09-15), pages 111 - 123 * |
| 李文杰等: "基于原型沙的三峡库区细颗粒泥沙冲淤模型试验研究", 《泥沙研究》, no. 3, 15 June 2016 (2016-06-15), pages 17 - 23 * |
| 梁艳洁;谢慰;赵正伟;罗秋实;付健;: "东庄水库运用方式对渭河下游减淤作用研究", 人民黄河, no. 10, 20 October 2016 (2016-10-20), pages 131 - 136 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118013895B (zh) | 2024-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108629055B (zh) | 一种基于饱和输沙原理的沙质内河航道回淤量预报方法 | |
| CN113887087B (zh) | 潮汐河段底沙输移引起的航道淤积量计算方法及系统 | |
| Yang et al. | Developing a comprehensive evaluation method for Interconnected River System Network assessment: A case study in Tangxun Lake group | |
| CN114548511A (zh) | 一种低水头水电站群下游水库日平均入库流量的预测方法 | |
| CN118013895B (zh) | 一种水库深孔前淤积厚度的计算方法 | |
| CN110472367B (zh) | 一种多沙河流干支流水沙全交互模拟方法及系统 | |
| CN106320255B (zh) | 一种漫滩水流滩槽过流量的计算方法 | |
| CN111428915A (zh) | 一种基于大数据的江河入海口深水航道回淤量预测方法 | |
| CN111241467A (zh) | 一种多沙河流水库排沙比的计算方法 | |
| CN106383997A (zh) | 一种反推三峡水库区间入流过程的计算方法 | |
| CN206693373U (zh) | 一种截流限流系统 | |
| CN116663716A (zh) | 一种梯级水电站水量平衡计算的水位点预测方法及系统 | |
| CN101418564B (zh) | 一种稳定输水渠道 | |
| CN116956763A (zh) | 一种基于调控下游退水闸的输水渡槽降水位增流方法 | |
| CN115545460A (zh) | 一种水库泥沙淤积风险二维评估方法 | |
| Jiao et al. | Morphodynamic characteristics and medium-term simulation of the north—south passage under the impact of the Yangtze Estuary deepwater navigation channel project | |
| CN114595596A (zh) | 一种水库淤积物密实沉降速率的模拟方法 | |
| CN109898461A (zh) | 一种计算排沙漏斗排沙效率的方法 | |
| CN110188310B (zh) | 梯级水库多流态复合的排沙预测方法及装置 | |
| CN111539153B (zh) | 一种基于预构泥沙信息库的水沙联合优化调度方法 | |
| Saleh et al. | Analysis of the Euphrates River’s movement within Al-Hindiya, Karbala, relative to steady flow conditions using the HEC-RAS model | |
| Tilai et al. | Analysis of sediment deposition downstream tidal sluice of estuary | |
| CN103870685A (zh) | 一种典型淤地坝坝控流域水土保持措施的评价方法 | |
| CN214530485U (zh) | 一种适用于多沙河流水电站的泥沙管理结构 | |
| CN209603104U (zh) | 用于下游水位变幅较大的厂房尾水布置结构 |
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 |