CN107918722A - 一种平原水库井渗漏计算方法 - Google Patents

一种平原水库井渗漏计算方法 Download PDF

Info

Publication number
CN107918722A
CN107918722A CN201711328924.4A CN201711328924A CN107918722A CN 107918722 A CN107918722 A CN 107918722A CN 201711328924 A CN201711328924 A CN 201711328924A CN 107918722 A CN107918722 A CN 107918722A
Authority
CN
China
Prior art keywords
well
reservoir
leakage
wall
plain
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
Application number
CN201711328924.4A
Other languages
English (en)
Other versions
CN107918722B (zh
Inventor
李玉莹
吴先敏
高峰
单既连
刘英豪
李利红
李文义
吴建伟
牛超
杨士瑞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Water Conservancy Survey And Design Institute Co ltd
Original Assignee
Shandong Survey and Design Institute of Water Conservancy Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong Survey and Design Institute of Water Conservancy Co Ltd filed Critical Shandong Survey and Design Institute of Water Conservancy Co Ltd
Priority to CN201711328924.4A priority Critical patent/CN107918722B/zh
Publication of CN107918722A publication Critical patent/CN107918722A/zh
Application granted granted Critical
Publication of CN107918722B publication Critical patent/CN107918722B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Abstract

本发明公开了一种平原水库井渗漏计算方法。本发明针对库区内井对平原水库渗漏影响，提出了利用离散方法计算平原水库任意位置井对水库渗漏量影响计算公式，可应用于库区内任意位置井对水库渗漏影响的计算。离散单元法利用井在库区内距离坝轴线位置不同，根据计算精度要求离散成n个单元，首先计算每个单元的渗漏量，通过叠加计算整个井渗漏量，评估井对平原水库渗漏影响。依据本发明提供的计算公式，可以方便计算单个井对水库渗漏量影响，计算过程简明，结果准确。

Description

一种平原水库井渗漏计算方法
技术领域
本发明涉及土木工程技术领域，尤其涉及一种平原水库井渗漏计算方法，可应用于库区内任意位置井对水库渗漏影响的计算。
背景技术
平原水库是相对于山区水库而言的，其一般位于大江、大河下游冲击平原地区，这类地区地质岩性多数为黏土、粉质黏土、粉土、砂土、沙壤土、粉细砂等组成。平原水库所处位置如在井灌区时，原库底中有较多的农田灌溉和生活用水机井。部分水井是增加水库渗漏和加剧坝后农田浸没的因素，如处理不当会对今后水库运行带来影响。井的渗流计算采用轴对称来计算，一般计算方法是以井为中心，采用轴对称影响半径的方法对井渗流量进行计算分析。但对于平原水库库区内井来说，井的渗漏影响是空间三维问题，只能采用空间三维有限元程序进行计算，对分析人员和程序要求很高，且建模复杂，计算量巨大。
发明内容
本发明就是针对上述存在的缺陷而提供一种平原水库井渗漏计算方法。本发明针对库区内井对平原水库渗漏影响，提出了利用离散方法计算平原水库任意位置井对水库渗漏量影响计算公式，可应用于库区内任意位置井对水库渗漏影响的计算。离散单元法利用井在库区内距离坝轴线位置不同，根据计算精度要求离散成n个单元，首先计算每个单元的渗漏量，通过叠加计算整个井渗漏量，评估井对平原水库渗漏影响。
本发明的一种平原水库井渗漏计算方法技术方案为，用离散单元法计算平原水库任意位置井对水库渗漏量影响。
井在库区内任意位置，按照井和水库大坝轴线相对位置和计算精度要求，离散成n个单元，渗漏量分别为Q1……Qn，对应的轴对称半径分别为R1……Rn，井渗漏量Q井1利用传统轴对称计算，井对水库渗漏影响Q计算公式如下：
式中Q-井渗漏量(m3/d)；
Q井1—井轴对称计算渗漏量(m3/d)；
N—离散单元数。
本发明的有益效果为：本发明针对库区内井对平原水库渗漏影响，提出了利用离散方法计算平原水库任意位置井对水库渗漏量影响计算公式，可应用于库区内任意位置井对水库渗漏影响的计算。离散单元法利用井在库区内距离坝轴线位置不同，根据计算精度要求离散成n个单元，首先计算每个单元的渗漏量，通过叠加计算整个井渗漏量，评估井对平原水库渗漏影响。依据本发明提供的计算公式，可以方便计算单个井对水库渗漏量影响，计算过程简明，结果准确。
附图说明
图1所示为井在平原水库中的离散单元简图；
图2所示为单个离散单元计算简图。
图3所示为计算模型简图；
图中，1-坝轴线，2-坝体，3-坝基，4-截渗沟，5-坝基一，6-坝基二，7-坝基三，8-井。
具体实施方式：
为了更好地理解本发明，下面用具体实例来详细说明本发明的技术方案，但是本发明并不局限于此。
实施例1
某一小型平原水库坝轴线长513m，坝体为壤土均质坝，坝高10.0m，坝顶高程45.0m，库底高程35.0m，库内最高水位43.0m，相应库容8.39万m3，井位于库区内(见说明书附图图2)。坝体渗透系数为1×10-5cm/s。坝基1层为渗透系数1×10-4cm/s，坝基2层为强透水层,渗透系数5×10-3cm/s，坝基3层渗透系数为渗透系数2×10-5cm/s。正常运行时上游水位43.0m，下游截渗沟水位34.0m，井深17.0m，穿透坝基1层到达坝基2层强透水，计算模型见说明书附图图3。
水库根据计算精度要求离散为15度一个单元，共离散为24个半径，离散半径见下表1，根据计算时分别取有井和无井两种模型计算，经计算，无井时水库每天渗漏量36.82m3/d,年渗漏量1.34万m3/y，占水库最高水位库容16.03％；有井时水库每天渗漏量39.13m3/d,年渗漏量1.43万m3/y，占水库最高水位库容17.03％。因井增加渗漏量为844.61m3/y，占水库库容的1.0％。
表1 井渗漏离散半径和计算结果表

Claims (3)

1.一种平原水库井渗漏计算方法，其特征在于，利用离散单元法计算平原水库任意位置井对水库渗漏量影响。
2.根据权利要求1所述的一种平原水库井渗漏计算方法，其特征在于，井在库区内任意位置，按照井和水库大坝轴线相对位置和计算精度要求，离散成n个单元，渗漏量分别为Q1……Qn，对应的轴对称半径分别为R1……Rn，井渗漏量Q井1利用传统轴对称计算，井对水库渗漏影响Q。
3.根据权利要求2所述的一种平原水库井渗漏计算方法，其特征在于，井对水库渗漏影响Q计算公式为：
<mrow> <mi>Q</mi> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>Q</mi> <mi>i</mi> </msub> </mrow>
式中Q-井渗漏量(m3/d)；
Q井1—井轴对称计算渗漏量(m3/d)；
N—离散单元数。
CN201711328924.4A 2017-12-13 2017-12-13 一种平原水库井渗漏计算方法 Active CN107918722B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711328924.4A CN107918722B (zh) 2017-12-13 2017-12-13 一种平原水库井渗漏计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711328924.4A CN107918722B (zh) 2017-12-13 2017-12-13 一种平原水库井渗漏计算方法

Publications (2)

Publication Number Publication Date
CN107918722A true CN107918722A (zh) 2018-04-17
CN107918722B CN107918722B (zh) 2021-07-09

Family

ID=61893247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711328924.4A Active CN107918722B (zh) 2017-12-13 2017-12-13 一种平原水库井渗漏计算方法

Country Status (1)

Country Link
CN (1) CN107918722B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105676308A (zh) * 2016-01-14 2016-06-15 中国地质大学（武汉） 一种单井地下水渗流流速流向测量方法及测量仪
US20160208599A1 (en) * 2015-01-21 2016-07-21 Saudi Arabian Oil Company Measuring inter-reservoir cross flow rate through unintended leaks in zonal isolation cement sheaths in offset wells
CN106595506A (zh) * 2016-11-30 2017-04-26 山东金米尔仪器科技有限公司 一种水库坝体变形监测方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160208599A1 (en) * 2015-01-21 2016-07-21 Saudi Arabian Oil Company Measuring inter-reservoir cross flow rate through unintended leaks in zonal isolation cement sheaths in offset wells
CN105676308A (zh) * 2016-01-14 2016-06-15 中国地质大学（武汉） 一种单井地下水渗流流速流向测量方法及测量仪
CN106595506A (zh) * 2016-11-30 2017-04-26 山东金米尔仪器科技有限公司 一种水库坝体变形监测方法及系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
GEBREMEDHIN BERHANE等: ""Water leakage investigation of micro-dam reservoirs in Mesozoic sedimentary sequences in Northern Ethiopia"", 《JOURNAL OF AFRICAN EARTH SCIENCES》 *
康晓军: ""十三陵水库流域水文模拟及渗漏分析"", 《中国优秀硕士学位论文全文数据库·工程科技Ⅱ辑》 *
李伟: ""北塘水库对地下水的渗漏补给特性研究"", 《中国优秀硕士学位论文全文数据库·工程科技Ⅱ辑》 *
王俊鹏等: ""水库枢纽地下水环境影响评价与保护探讨——以凤凰湖调蓄水库工程为例"", 《地下水》 *
黄生文等: ""基于离散单元法的基坑涌水量数值模拟"", 《长沙理工大学学报(自然科学版)》 *

Also Published As

Publication number Publication date
CN107918722B (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN100487404C (zh) 一种测定河道草皮护坡抗冲性能的试验装置
CN108446489A (zh) 海岸线变迁对湿地地下水影响的测量方法及其处理装置
CN110046400A (zh) 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法
CN103650695B (zh) 一种根据土层厚度进行土地平整的方法
CN102322078B (zh) 一种承压水基坑开挖与降压耦合效应的突涌破坏计算处理方法
CN107229811B (zh) 一种饱和软土自重固结过程的预测方法
Bezuijen The influence of the leakage length on the initiation of backward erosion piping
CN107918722A (zh) 一种平原水库井渗漏计算方法
CN106759072A (zh) 兼顾床面稳定和快速排水的河工模型排水系统及设计方法
CN106202911A (zh) 一种河道生态流量地下水补给需水计算方法
CN101074899B (zh) 一种测定河道草皮护坡抗冲性能的方法
CN102930148B (zh) 一种基于随机起动的管涌渗透系数确定方法
CN106193139B (zh) 一种临河基坑施工对桩基影响的测试装置
Sun et al. Application of HEC-RAS for flood forecasting in perched river–A case study of hilly region, China
CN110472309B (zh) 一种二维渗流计算模型中考虑三维绕渗效应的方法、装置
CN206616513U (zh) 兼顾床面稳定和快速排水的河工模型排水系统
Bezuijen et al. Analytical groundwater flow calculations for understanding the flow and erosion in a coarse sand barrier
CN107038291B (zh) 高填方涵洞涵顶垂直土压力计算方法
CN105445177B (zh) 岸坡侵蚀量和侵蚀速率的监测方法
CN105040674B (zh) 一种高置换率散体材料桩复合地基固结度计算方法
Peterson Jr Intercepting drainage wells in artesian aquifer
CN113375751B (zh) 一种深厚覆盖层河床大坝渗漏量检测方法
Nagai et al. Pressures by breaking waves on composite-type breakwaters
Rantasewee et al. Groundwater flow modelling for the development of managed aquifer recharge scheme in irrigation projects, Thailand
Wang et al. Study on seepage control measures of plain reservoirs based on soil salinization after dam

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
TA01 Transfer of patent application right

Effective date of registration: 20210621

Address after: Room 501, G1 / F, 288 Ningxia road, Shinan District, Qingdao City, Shandong Province 266072

Applicant after: Qingdao Water Conservancy Survey and Design Institute Co.,Ltd.

Address before: No.121 Lishan Road, Lixia District, Jinan City, Shandong Province

Applicant before: SHANDONG SURVEY AND DESIGN INSTITUTE OF WATER CONSERVANCY

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant