CN107179062B - 一种实验水底地形观测方法 - Google Patents
一种实验水底地形观测方法 Download PDFInfo
- Publication number
- CN107179062B CN107179062B CN201710254619.9A CN201710254619A CN107179062B CN 107179062 B CN107179062 B CN 107179062B CN 201710254619 A CN201710254619 A CN 201710254619A CN 107179062 B CN107179062 B CN 107179062B
- Authority
- CN
- China
- Prior art keywords
- topographic profile
- sink
- vector
- water
- topographic
- 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
- 238000002474 experimental method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 20
- 239000013598 vector Substances 0.000 claims abstract description 73
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 18
- 238000000605 extraction Methods 0.000 claims description 6
- 239000011521 glass Substances 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 239000000725 suspension Substances 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/06—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring contours or curvatures
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本专利涉及一种实验水底地形观测方法,在实验过程中停止水流,能使水体逐渐澄清,便于测量并记录水槽中部地形剖面;当求出两组基向量之间的变换矩阵后,便可以在任意时间点上,根据从水槽侧壁观测到的地形剖面推算出水槽中部地形剖面,实现对水槽中部的地形剖面的间接观测,因为水槽侧壁处地形剖面便于观察记录,而水槽中部的地形剖面因在水流中有悬浮的泥沙而难以被测出,因此本专利通过水槽侧壁处地形剖面间接观测水槽中部地形剖面,实验效率高。
Description
技术领域
本专利涉及水流-泥沙实验领域,具体涉及一种实验水底地形观测方法。
背景技术
为了研究水流通过局部侵蚀和淤积对水底地形造成的变化,常在实验室水槽中进行实验。因为水槽中部的泥沙受其他无关变量影响较小,例如水槽侧壁对水流有回旋或反弹等作用,而在水槽中部这些作用的影响可最小化,因此,实验研究常基于水槽中部的地形测量数据来进行。水槽中部地形可通过超声波测深仪测量,但实验过程中粒径较小的泥沙会随着水流从水底起动而悬浮在水体中,使超声波测深难以精确地实时测量水底地形剖面的变化过程,虽然实验室水槽侧壁常装有高清玻璃,即使水体浑浊亦可通过侧壁的高清玻璃观测记录水底地形在侧壁的剖面,但是水槽侧壁对水流产生的侧壁效应会使水流呈现一定的三维效应,进而使水底地形也具有一定的三维形态,因此,通过水槽侧壁记录的侧壁处地形剖面和用超声波测深仪测量到的水槽中部地形剖面之间存在一定差别,难以通过水槽侧壁处地形剖面获得水槽中部地形剖面。
发明内容
为了克服现有技术中的缺陷,本专利提供一种实验水底地形观测方法,能够较好地通过便于观察记录的水槽侧壁处地形剖面获得难以测量得到的水槽中部地形剖面,极大地提高了实验效率。
针对上述技术问题,本专利是这样加以解决的:一种实验水底地形观测方法,其步骤如下:
在实验过程中一直测量并记录水槽侧壁地形剖面;
S1:实验开始,使水在水槽中进行流动;
S2:在某一个时间点停止水流,测量并记录此时水槽中部地形剖面;
S3:对同时刻的水槽中部地形剖面及水槽侧壁地形剖面进行同阶数的多项式拟合,分别得到水槽中部地形曲线及水槽侧壁地形曲线,并分别从得到的两个多项式中提取系数,组成中部地形剖面向量及侧壁地形剖面向量;
S4:重新打开水流,每隔一段的时间停止水流,测量并记录水槽中部地形剖面,再重复步骤S3,直至得到的中部地形剖面向量、侧壁地形剖面向量及多项式系数的数量均相同;
S5:若得到的多个中部地形剖面向量之间线性无关,多个侧壁地形剖面向量之间也线性无关,则进入下一步,若否,则返回步骤S2;
S6:将多个中部地形剖面向量及多个侧壁地形剖面向量转换为两组基向量,求出这两组基向量之间的变换矩阵;
S7:在任何时间点记录到的水槽侧壁地形剖面通过该变换矩阵可得到同时刻的水槽中部地形剖面。
在实验中,实验设备为水槽,水槽中装有泥沙,实验时让水在泥沙上流过,研究水流对水槽中水底地形的影响,本专利在实验过程中停止水流,能使水体逐渐澄清,便于测量并记录水槽中部地形剖面;当求出两组基向量之间的变换矩阵后,便可以在任意时间点上,根据从水槽侧壁观测到的地形剖面推算出水槽中部地形剖面,实现对水槽中部的地形剖面的间接观测,因为水槽侧壁处地形剖面便于观察测量并记录,而水槽中部的地形剖面因在水流中有悬浮的泥沙而难以被测出,因此本专利通过水槽侧壁处地形剖面间接观测水槽中部地形剖面,实验效率高。
进一步地,所述步骤S3中根据水槽侧壁地形剖面进行拟合后得到的多项式为y=a0+a1x+a2x2+…+am-1xm-1+amxm,根据水槽中部地形剖面进行拟合后得到的多项式为y=b0+b1x+b2x2+…+bm-1xm-1+bmxm,提取系数后得到的侧壁地形剖面向量为 中部地形剖面向量为
进一步地,设所述步骤S6中多个侧壁地形剖面向量得出的基向量为多个中部地形剖面向量得出的基向量为 那么存在即可求出这两组基向量之间的变换矩阵A。
进一步地,所述步骤S7的具体步骤为:
S71:对任一时刻观测到水槽侧壁地形剖面进行与之前步骤S2相同阶数的多项式拟合,提取该多项式的系数组成待处理向量;
S72:变换矩阵的逆矩阵乘待处理向量,得到目的向量;
S73:将目的向量提取为多项式的系数,得到水槽中部地形曲线的拟合多项式方程,最终可得到同时刻的水槽中部地形剖面。
因为两组基向量之间存在变换矩阵,所以此处利用线性代数中的坐标变换公式,将待处理向量转换为目的向量。
进一步地,通过超声波测深仪测量并记录水槽中部地形剖面。
进一步地,水槽侧壁为标有坐标轴及刻度的透明玻璃,通过相机拍摄水槽侧壁地形的照片,再从该照片中提取出水槽侧壁地形剖面。
进一步地,所述多项式拟合采用最小二乘法实现,并且均为相同阶数的拟合。
相比于现有技术,本专利的有益效果为:通过拟合获得多项式,提取系数形成向量,获取多组向量组成基向量,直到最后求出变换矩阵,求解出水槽侧壁处地形剖面与水槽中部地形剖面之间存在的理论关系,便可方便地从便于观察的水槽侧壁处地形剖面获取难以测量得到的水槽中部地形剖面,提高实验效率。
附图说明
图1是本专利的流程图。
图2是本专利对地形剖面的实测数据和多项式拟合曲线的示意图。
图3是本专利通过变换矩阵将水槽侧壁地形剖面转换为同时刻的水槽中部地形剖面的流程图。
具体实施方式
如图1所示的一种实验水底地形观测方法,其步骤如下:
水槽侧壁为标有坐标轴及刻度的透明玻璃,在实验过程中一直通过相机拍摄水槽侧壁地形的照片,再从该照片中提取出水槽侧壁地形剖面;
S1:实验开始,使水在水槽中进行流动;
S2:在某一个时间点停止水流,测量并记录此时水槽中部地形剖面;
S3:如图2所示,对同时刻的水槽中部地形剖面及水槽侧壁地形剖面进行同阶数的多项式拟合,分别得到水槽中部地形曲线及水槽侧壁地形曲线,并分别从得到的两个多项式中提取系数,组成中部地形剖面向量及侧壁地形剖面向量;
S4:重新打开水流,每隔一段的时间停止水流,测量并记录水槽中部地形剖面,再重复步骤S3,直至得到的中部地形剖面向量、侧壁地形剖面向量及多项式系数的数量均相同;
S5:若得到的多个中部地形剖面向量之间线性无关,多个侧壁地形剖面向量之间也线性无关,则进入下一步,若否,则返回步骤S2;
S6:将多个中部地形剖面向量及多个侧壁地形剖面向量转换为两组基向量,求出这两组基向量之间的变换矩阵;
S7:在任何时间点测量并记录到的水槽侧壁地形剖面通过该变换矩阵可得到同时刻的水槽中部地形剖面。
在实验中,实验设备为水槽,水槽中装有泥沙,实验时让水在泥沙上流过,研究水流对水槽中水底地形的影响,本专利在实验过程中停止水流,能使水体逐渐澄清,便于测量并记录水槽中部地形剖面;当求出两组基向量之间的变换矩阵后,便可以在任意时间点上,根据从水槽侧壁观测到的地形剖面推算出水槽中部地形剖面,实现对水槽中部的地形剖面的间接观测,因为水槽侧壁处地形剖面便于观察测量并记录,而水槽中部的地形剖面因在水流中有悬浮的泥沙而难以被测出,因此本专利通过水槽侧壁处地形剖面间接观测水槽中部地形剖面,实验效率高。
所述步骤S3中根据水槽侧壁地形剖面进行拟合后得到的多项式为y=a0+a1x+a2x2+…+am-1xm-1+amxm,根据水槽中部地形剖面进行拟合后得到的多项式为y=b0+b1x+b2x2+…+bm-1xm-1+bmxm,提取系数后得到的侧壁地形剖面向量为中部地形剖面向量为
所述步骤S6中得到n个侧壁地形剖面向量及中部地形剖面向量,其中n=m+1,那n个侧壁地形剖面向量得出的基向量为n个中部地形剖面向量得出的基向量为那么存在即可求出这两组基向量之间的变换矩阵A,也即
如图3所示,所述步骤S7的具体步骤为:
S71:对任一时刻观测到水槽侧壁地形剖面进行与之前步骤S2相同阶数的多项式拟合,得到y=c0+c1x+c2x2+…+cm-1xm-1+cmxm,提取该多项式的系数组成待处理向量
S72:记要得到目的向量那么向量可用下式求出:
其中即为之前所求变换矩阵A。
S73:将目的向量提取为多项式的系数,得到水槽中部地形剖面的拟合多项式方程y=d0+d1x+d2x2+…+dm-1xm-1+dmxm,最终可得到同时刻的水槽中部地形剖面。
在具体实施过程中,通过超声波测深仪测量并记录水槽中部地形剖面,其中的多项式拟合采用最小二乘法实现,并且均为相同阶数的拟合。
Claims (7)
1.一种实验水底地形观测方法,其特征在于,其步骤如下:
在实验过程中一直测量并记录水槽侧壁地形剖面;
S1:实验开始,使水在水槽中进行流动;
S2:在某一个时间点停止水流,测量并记录此时水槽中部地形剖面;
S3:对同时刻的水槽中部地形剖面及水槽侧壁地形剖面进行同阶数的多项式拟合,分别得到水槽中部地形曲线及水槽侧壁地形曲线,并分别从得到的两个多项式中提取系数,组成中部地形剖面向量及侧壁地形剖面向量;
S4:重新打开水流,每隔一段的时间停止水流,测量并记录水槽中部地形剖面,再重复步骤S3,直至得到的中部地形剖面向量、侧壁地形剖面向量及多项式系数的数量均相同;
S5:若得到的多个中部地形剖面向量之间线性无关,多个侧壁地形剖面向量之间也线性无关,则进入下一步,若否,则返回步骤S2;
S6:将多个中部地形剖面向量及多个侧壁地形剖面向量转换为两组基向量,求出这两组基向量之间的变换矩阵;
S7:在任何时间点记录到的水槽侧壁地形剖面通过该变换矩阵可得到同时刻的水槽中部地形剖面。
2.根据权利要求1所述的一种实验水底地形观测方法,其特征在于,所述步骤S3中根据水槽侧壁地形剖面进行拟合后得到的多项式为a0+a1x+a2x2+…+am-1xm-1+amxm,根据水槽中部地形剖面进行拟合后得到的多项式为b0+b1x+b2x2+…+bm-1xm-1+bmxm,提取系数后得到的侧壁地形剖面向量为中部地形剖面向量为所述m为非负整数。
3.根据权利要求1所述的一种实验水底地形观测方法,其特征在于,设所述步骤S6中多个侧壁地形剖面向量得出的基向量为 多个中部地形剖面向量得出的基向量为那么存在即可求出这两组基向量之间的变换矩阵A,所述n为非负整数。
4.根据权利要求1所述的一种实验水底地形观测方法,其特征在于,所述步骤S7的具体步骤为:
S71:对任一时刻观测到水槽侧壁地形剖面进行与之前步骤S2相同阶数的多项式拟合,提取该多项式的系数组成待处理向量;
S72:变换矩阵的逆矩阵乘待处理向量,得到目的向量;
S73:将目的向量提取为多项式的系数,得到水槽中部地形曲线的拟合多项式方程,最终可得到同时刻的水槽中部地形剖面。
5.根据权利要求1至4任一项所述的一种实验水底地形观测方法,其特征在于,通过超声波测深仪测量并记录水槽中部地形剖面。
6.根据权利要求1至4任一项所述的一种实验水底地形观测方法,其特征在于,水槽侧壁为标有坐标轴及刻度的透明玻璃,通过相机拍摄水槽侧壁地形的照片,再从该照片中提取出水槽侧壁地形剖面。
7.根据权利要求1至4任一项所述的一种实验水底地形观测方法,其特征在于,所述多项式拟合采用最小二乘法实现,并且均为相同阶数的拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254619.9A CN107179062B (zh) | 2017-04-18 | 2017-04-18 | 一种实验水底地形观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254619.9A CN107179062B (zh) | 2017-04-18 | 2017-04-18 | 一种实验水底地形观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107179062A CN107179062A (zh) | 2017-09-19 |
CN107179062B true CN107179062B (zh) | 2019-06-21 |
Family
ID=59832001
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710254619.9A Active CN107179062B (zh) | 2017-04-18 | 2017-04-18 | 一种实验水底地形观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107179062B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109668511B (zh) * | 2019-01-30 | 2021-07-30 | 中国海洋大学 | 一种适用于水槽的人工鱼礁堆积形态测量装置及方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000233A (zh) * | 2007-01-17 | 2007-07-18 | 曲兆松 | 水下三维地形测量*** |
KR100936469B1 (ko) * | 2009-07-10 | 2010-01-20 | (주)새한지오텍 | 해저지형 및 지리정보에 대한 데이터 수집시스템 |
CN101709965A (zh) * | 2009-11-06 | 2010-05-19 | 天津大学 | 水槽三维地形自动测量装置 |
CN103591939A (zh) * | 2013-10-29 | 2014-02-19 | 中国科学院力学研究所 | 基于主动立体视觉技术的模拟海床地形测量方法及测量装置 |
CN203893831U (zh) * | 2014-06-17 | 2014-10-22 | 武汉大学 | 一种水下模型地形的测量与三维重建装置 |
CN104776977A (zh) * | 2015-04-28 | 2015-07-15 | 中国海洋大学 | 一种海岸工程泥沙物理模型试验底床动态综合观测方法 |
CN105301223A (zh) * | 2015-11-12 | 2016-02-03 | 重庆大学 | 一种便捷式实时监测水土流失模拟试验方法 |
CN105865421A (zh) * | 2016-04-07 | 2016-08-17 | 长沙理工大学 | 基于照相机图像处理技术的水槽三维地形测量装置 |
-
2017
- 2017-04-18 CN CN201710254619.9A patent/CN107179062B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000233A (zh) * | 2007-01-17 | 2007-07-18 | 曲兆松 | 水下三维地形测量*** |
KR100936469B1 (ko) * | 2009-07-10 | 2010-01-20 | (주)새한지오텍 | 해저지형 및 지리정보에 대한 데이터 수집시스템 |
CN101709965A (zh) * | 2009-11-06 | 2010-05-19 | 天津大学 | 水槽三维地形自动测量装置 |
CN103591939A (zh) * | 2013-10-29 | 2014-02-19 | 中国科学院力学研究所 | 基于主动立体视觉技术的模拟海床地形测量方法及测量装置 |
CN203893831U (zh) * | 2014-06-17 | 2014-10-22 | 武汉大学 | 一种水下模型地形的测量与三维重建装置 |
CN104776977A (zh) * | 2015-04-28 | 2015-07-15 | 中国海洋大学 | 一种海岸工程泥沙物理模型试验底床动态综合观测方法 |
CN105301223A (zh) * | 2015-11-12 | 2016-02-03 | 重庆大学 | 一种便捷式实时监测水土流失模拟试验方法 |
CN105865421A (zh) * | 2016-04-07 | 2016-08-17 | 长沙理工大学 | 基于照相机图像处理技术的水槽三维地形测量装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107179062A (zh) | 2017-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fressengeas et al. | Dislocation transport and intermittency in the plasticity of crystalline solids | |
CN107478267B (zh) | 基于光场相机的三维流场的温度场和速度场同步测试方法 | |
US8594939B2 (en) | Method for calibrating seismic imaging velocities | |
Rosenau et al. | Grounding line migration and high‐resolution calving dynamics of Jakobshavn Isbræ, West Greenland | |
CN111366936B (zh) | 多波束测深数据处理方法及装置 | |
CN107179062B (zh) | 一种实验水底地形观测方法 | |
CN103996171A (zh) | 基于时空图像的流体运动矢量估计方法 | |
CN109002418B (zh) | 基于体素生长和地面激光点云的树木胸径自动计算方法 | |
CN109633749B (zh) | 基于散射积分法的非线性菲涅尔体地震走时层析成像方法 | |
CN109085643A (zh) | 早至波的分步联合反演方法 | |
CN103425992B (zh) | 一种基于合成孔径声纳图像的海底底质分类方法及*** | |
CN106018198B (zh) | 一种气泡粒径的反演计算方法 | |
Mozumi et al. | 3D stereo imaging of abnormal waves in a wave basin | |
CN102538691B (zh) | 一种用图像识别技术测量高温钢管直径的方法 | |
CN105043390A (zh) | 基于泛克里金法的重力场插值方法 | |
CN102798380B (zh) | 线阵图像中目标运动参数的测量方法 | |
CN103591904A (zh) | 一种两步三维傅里叶变换测量物体内部三维变形场的方法 | |
CN106896105B (zh) | 一种检测单晶叶片小角度晶界取向差的方法 | |
CN100487357C (zh) | 地表微形态测量装置及其使用方法 | |
CN104704396A (zh) | 处理重力梯度数据的方法 | |
CN105931250B (zh) | 一种非球形气泡的等效粒径转换方法 | |
CN113533331B (zh) | 基于石盐包裹体进行深层古盐湖沉积温度确定的方法 | |
Moskalik et al. | Statistical analysis of Brepollen bathymetry as a key to determine average depths on a glacier foreland | |
CN105651212A (zh) | 一种土壤表面粗糙度检测装置及方法 | |
CN112379434B (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 | ||
GR01 | Patent grant |