CN112749429A - 一种基于多元插值的海上风电固定式基础结构设计方法 - Google Patents
一种基于多元插值的海上风电固定式基础结构设计方法 Download PDFInfo
- Publication number
- CN112749429A CN112749429A CN202011386827.2A CN202011386827A CN112749429A CN 112749429 A CN112749429 A CN 112749429A CN 202011386827 A CN202011386827 A CN 202011386827A CN 112749429 A CN112749429 A CN 112749429A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- design
- wind power
- parameter
- environmental
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种基于多元插值的海上风电固定式基础结构设计方法,包括以下步骤:1)基于灵敏度分析确定各环境参数与风电基础结构耗钢量的相关系数,2)采用风电基础结构分析软件对相关程度较大的环境参数进行严格正交参数分析,3)选取实际工程的环境参数开展设计参数的多元插值计算,该方法可用于快速确定海上风电基础结构的设计参数。
Description
技术领域
本发明属于海上风电基础结构设计领域,涉及一种基于多元插值的海上风电固定式基础结构设计方法。
背景技术
随着我国工业的快速发展,对能源的需求日益增大。风能资源因利用成本低,分布区域广成为国家缓解能源供需矛盾、减轻环境污染的重要途径。2019年,我国风电 并网容量达到2.09亿千瓦,同比增长14%,风电发电量在全国发电量的比重为5.5%。 我国的风电场建设初期主要集中在陆上传统风资源优势区域(三北地区),但面临电力 消纳能力不足、弃风限电等问题。近年来,由于陆上风电市场趋于饱和,海上风力发 电成为风能资源开发新的增长点。我国具有较长的海岸线,海上风资源不仅平稳丰富, 而且区域平坦广阔,在开发过程中不会面临征地和噪音等问题。海上风电场在运行期 间基础不仅承受水动力载荷,还会受到上部风机气动载荷的长期作用,结构内部频繁 的交变应力计算会严重影响结构的设计进度。本发明提出了基于参数插值的海上风电 基础结构设计参数计算方法,通过提前开展正交参数分析建立海上风电基础结构的数 据库,在得到实际工程的环境勘测报告后可通过插值的方法快速得到结构的设计参数, 有利于海上风电场建设的“降本增效”。
发明内容
技术问题:本发明提供一种简单准确、快速合理的海上风电固定式基础结构设计方法。
技术方案:一种基于多元插值的海上风电固定式基础结构设计方法,包括以下步骤:
1)根据已有工程建立海上风电基础结构的基准有限元模型,利用蒙特卡罗模拟技术对海水深度、风机载荷、土体参数、波浪高度、波浪周期、海水流速、淤泥层厚 度等所有环境参数进行灵敏度分析,在选定环境参数取值范围和概率分布函数后,进 行N次以上随机抽样计算,得到各环境参数{EV}T和风电基础结构耗钢量SC之间的 线性相关系数矩阵,{EV}T为环境参数全集;
2)选取相关系数大于Sum的q个环境参数开展严格正交参数计算,每次计算时 仅改变一个环境参数,其他环境参数保持不变,计算目标为泥面处桩变形和侧摩阻力 满足设计限值条件下基础结构耗钢量最小。每个环境参数EVx共变换三次取值,分别 为上限值EVxu、基准值EVxb和下限值EVxl,EVx代表第x个环境参数,x为环境参数 序号,通过对每个工况进行设计计算,得到海上风电基础结构的设计参数值 {DV}T=(DV1,DV2…DVy…DVp),DVy代表第y个设计参数值,y为设计参数序号,p 为设计参数总数,数据库中共包含3q个海上风电结构设计工况;
3)对实际工程环境参数{EV}T=(EV1,EV2…EVx…EVq)进行多元插值计算,起 始设计组包含3q个工况,可以取出3q-1个除了EV1不同其余环境参数均相同的插值组, 分别采用每个插值组EV1u、EV1b与EV1l形成的曲线插值确定EV1对应的设计参数值, 经过第一次插值后设计组包含3q-1个工况,在这些工况中取出3q-2个除了EV2不同其 余环境参数均相同的插值组,分别采用每个插值组EV2u、EV2b与EV2l形成的曲线插 值确定EV2对应的设计参数值,经过第二次插值后设计组包含3q-2个工况,以此类推, 当进行第q次插值时,仅剩余一个插值组,采用EVqu、EVqb与EVql形成的曲线插值 确定EVq对应的设计参数值后计算结束。
有益效果:本发明与现有技术相比,具有以下优点:
(1)本发明提出的设计方法在预先建立海上风电基础结构设计标准库的基础上,通过插值的计算方法来快速实现实际工程的参数设计,便于在招投标过程中快速得出 结构的合理设计方案和工程造价。而传统设计方法是在明确结构的环境参数后再进行 设计,载荷分析和内力校核等会涉及多次疲劳计算和迭代过程,使其难以满足实际工 程设计的时间要求;
(2)本发明提出的数据库大小是由选取的环境参数数量决定的,多元插值的方 法仅需将每个环境参数的取值进行三次变换,总样本的计算量相对较小。而通过传统 回归分析的思路则需要对每个环境参数进行数百次变换,以反映各个环境参数的变化 规律,拟合的回归公式才能较为准确,总样本的计算量往往在5000次以上,大大增 加了数据库的计算成本。
附图说明
图1为本发明方法的流程图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步的说明。
1)根据已有工程建立海上风电基础结构的基准有限元模型,利用蒙特卡罗模拟技术对海水深度、风机载荷、土体参数、波浪高度、波浪周期、海水流速、淤泥层厚 度等所有环境参数进行灵敏度分析,每个环境参数的上下限范围和概率分布函数可基 于大量实际工程的地质勘探、水文资料和经验取值等进行确定,常用的概率分布函数 有均匀分布、指数分布、正态分布、对数正态分布和截断高斯分布等,如风机载荷的 取值应覆盖海上风机最小兆瓦数和最大兆瓦数对应的弯矩值,一般服从截断高斯分布, 共进行N次以上随机抽样计算,建议采用拉丁超立方抽样技术,可以避免直接抽样法 数据点集中导致计算效率低的问题。在得到各环境参数{EV}T和风电基础结构耗钢量 SC之间的线性相关系数矩阵后,便可以判断环境参数与结构耗钢量相关程度排名, {EV}T为环境参数全集;
2)选取相关系数大于Sum的q个环境参数开展严格正交参数计算,以排除影响 不大的环境参数,每次计算时仅改变一个环境参数,其他环境参数保持不变,计算目 标为泥面处桩变形和侧摩阻力满足设计限值条件下基础结构耗钢量最小。每个环境参 数EVx共变换三次取值,分别为上限值EVxu、基准值EVxb和下限值EVxl,这样可以 确保后续插值采用的是内插方法,而不是计算精度较低的外插方法,EVx代表第x个 环境参数,x为环境参数序号,通过Bladed和Sesam等数值分析软件对每个工况仅改 变一个环境参数进行设计计算,得到各个工况下海上风电基础结构设计参数值 {DV}T=(DV1,DV2…DVy…DVp),DVy代表第y个设计参数值,y为设计参数序号,p 为设计参数总数,设计参数一般为基础结构各个组件的长度、直径和径厚等,数据库 中共包含3q个海上风电结构设计工况;
3)对实际工程环境参数{EV}T=(EV1,EV2…EVx…EVq)进行多元参数插值计算, 起始设计组包含3q个工况,由于数据库是通过严格正交参数分析计算得到的,可以取 出3q-1个除了EV1不同其余环境参数均相同的插值组,分别采用每个插值组EV1u、EV1b与EV1l形成的曲线插值确定EV1对应的设计参数值,经过第一次插值后设计组包含 3q-1个工况,在这些工况中取出3q-2个除了EV2不同其余环境参数均相同的插值组, 分别采用每个插值组EV2u、EV2b与EV2l形成的曲线插值确定EV2对应的设计参数值, 经过第二次插值后设计组包含3q-2个工况,以此类推,当进行第q次插值时,仅剩余 一个插值组,采用EVqu、EVqb与EVql形成的曲线插值确定EVq对应的设计参数值后 计算结束,经过上述过程插值得到的最后一组设计参数即为实际环境参数下海上风电 基础结构的计算结果。
Claims (1)
1.一种基于多元插值的海上风电固定式基础结构设计方法,其特征在于,该方法包括以下步骤:
1)根据已有工程建立海上风电基础结构的基准有限元模型,利用蒙特卡罗模拟技术对海水深度、风机载荷、土体参数、波浪高度、波浪周期、海水流速、淤泥层厚度等所有环境参数进行灵敏度分析,在选定环境参数取值范围和概率分布函数后,进行N次以上随机抽样计算,得到各环境参数{EV}T和风电基础结构耗钢量SC之间的线性相关系数矩阵,{EV}T为环境参数全集;
2)选取相关系数大于Sum的q个环境参数开展严格正交参数计算,每次计算时仅改变一个环境参数,其他环境参数保持不变,计算目标为泥面处桩变形和侧摩阻力满足设计限值条件下基础结构耗钢量最小。每个环境参数EVx共变换三次取值,分别为上限值EVxu、基准值EVxb和下限值EVxl,EVx代表第x个环境参数,x为环境参数序号,通过对每个工况进行设计计算,得到海上风电基础结构的设计参数值{DV}T=(DV1,DV2…DVy…DVp),DVy代表第y个设计参数值,y为设计参数序号,p为设计参数总数,数据库中共包含3q个海上风电结构设计工况;
3)对实际工程环境参数{EV}T=(EV1,EV2…EVx…EVq)进行多元插值计算,起始设计组包含3q个工况,可以取出3q-1个除了EV1不同其余环境参数均相同的插值组,分别采用每个插值组EV1u、EV1b与EV1l形成的曲线插值确定EV1对应的设计参数值,经过第一次插值后设计组包含3q-1个工况,在这些工况中取出3q-2个除了EV2不同其余环境参数均相同的插值组,分别采用每个插值组EV2u、EV2b与EV2l形成的曲线插值确定EV2对应的设计参数值,经过第二次插值后设计组包含3q-2个工况,以此类推,当进行第q次插值时,仅剩余一个插值组,采用EVqu、EVqb与EVql形成的曲线插值确定EVq对应的设计参数值后计算结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011386827.2A CN112749429B (zh) | 2020-12-01 | 2020-12-01 | 一种基于多元插值的海上风电固定式基础结构设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011386827.2A CN112749429B (zh) | 2020-12-01 | 2020-12-01 | 一种基于多元插值的海上风电固定式基础结构设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112749429A true CN112749429A (zh) | 2021-05-04 |
CN112749429B CN112749429B (zh) | 2023-09-01 |
Family
ID=75648939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011386827.2A Active CN112749429B (zh) | 2020-12-01 | 2020-12-01 | 一种基于多元插值的海上风电固定式基础结构设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112749429B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114004436A (zh) * | 2021-06-10 | 2022-02-01 | 重庆大学 | 一种基于性能反应谱的自复位海上风电升压站抗震性能评估方法 |
CN116306109A (zh) * | 2023-02-03 | 2023-06-23 | 哈尔滨工业大学(深圳) | 基于时域模型修正的海上风机土结相互作用状态识别方法 |
CN114004436B (zh) * | 2021-06-10 | 2024-05-28 | 重庆大学 | 一种基于性能反应谱的自复位海上风电升压站抗震性能评估方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7630852B1 (en) * | 2005-11-04 | 2009-12-08 | Arizona Board Of Regents | Method of evaluating integrated circuit system performance using orthogonal polynomials |
CN103279651A (zh) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | 一种储层参数预测方法及装置 |
CN103645065A (zh) * | 2013-12-25 | 2014-03-19 | 国电联合动力技术有限公司 | 一种海上风机基础全时程耦合疲劳分析方法及系统 |
CN104239662A (zh) * | 2013-06-14 | 2014-12-24 | 中国科学院城市环境研究所 | 一种定量评价区域风能资源开发潜力的方法 |
CN105260953A (zh) * | 2015-12-02 | 2016-01-20 | 甘肃省电力公司风电技术中心 | 一种综合风电场理论出力评估仿真系统和评估方法 |
CN105814574A (zh) * | 2013-10-04 | 2016-07-27 | 塞昆纳姆股份有限公司 | 遗传变异的非侵入性评估的方法和过程 |
CN106202804A (zh) * | 2016-07-22 | 2016-12-07 | 北京临近空间飞行器系统工程研究所 | 基于数据库的复杂外形飞行器分布式热环境参数预测方法 |
CN109635494A (zh) * | 2018-12-28 | 2019-04-16 | 中国航天空气动力技术研究院 | 一种飞行试验与地面仿真气动力数据综合建模方法 |
CN109787250A (zh) * | 2018-12-30 | 2019-05-21 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于多元经验模态分解的电力系统低频振荡模式辨识方法 |
CN109918786A (zh) * | 2019-03-07 | 2019-06-21 | 龙源(北京)风电工程设计咨询有限公司 | 一种海上风机基础结构自动分析系统及方法 |
CN110593849A (zh) * | 2019-08-22 | 2019-12-20 | 魏峰 | 一种考虑相渗曲线动态变化的油藏数值模拟方法 |
CN111400849A (zh) * | 2018-12-28 | 2020-07-10 | 北京金风科创风电设备有限公司 | 确定风电机组的时序载荷的方法和装置 |
CN111488974A (zh) * | 2020-04-14 | 2020-08-04 | 四川北控清洁能源工程有限公司 | 基于深度学习神经网络的海洋风能降尺度方法 |
-
2020
- 2020-12-01 CN CN202011386827.2A patent/CN112749429B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7630852B1 (en) * | 2005-11-04 | 2009-12-08 | Arizona Board Of Regents | Method of evaluating integrated circuit system performance using orthogonal polynomials |
CN103279651A (zh) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | 一种储层参数预测方法及装置 |
CN104239662A (zh) * | 2013-06-14 | 2014-12-24 | 中国科学院城市环境研究所 | 一种定量评价区域风能资源开发潜力的方法 |
CN105814574A (zh) * | 2013-10-04 | 2016-07-27 | 塞昆纳姆股份有限公司 | 遗传变异的非侵入性评估的方法和过程 |
CN103645065A (zh) * | 2013-12-25 | 2014-03-19 | 国电联合动力技术有限公司 | 一种海上风机基础全时程耦合疲劳分析方法及系统 |
CN105260953A (zh) * | 2015-12-02 | 2016-01-20 | 甘肃省电力公司风电技术中心 | 一种综合风电场理论出力评估仿真系统和评估方法 |
CN106202804A (zh) * | 2016-07-22 | 2016-12-07 | 北京临近空间飞行器系统工程研究所 | 基于数据库的复杂外形飞行器分布式热环境参数预测方法 |
CN109635494A (zh) * | 2018-12-28 | 2019-04-16 | 中国航天空气动力技术研究院 | 一种飞行试验与地面仿真气动力数据综合建模方法 |
CN111400849A (zh) * | 2018-12-28 | 2020-07-10 | 北京金风科创风电设备有限公司 | 确定风电机组的时序载荷的方法和装置 |
CN109787250A (zh) * | 2018-12-30 | 2019-05-21 | 国网辽宁省电力有限公司电力科学研究院 | 一种基于多元经验模态分解的电力系统低频振荡模式辨识方法 |
CN109918786A (zh) * | 2019-03-07 | 2019-06-21 | 龙源(北京)风电工程设计咨询有限公司 | 一种海上风机基础结构自动分析系统及方法 |
CN110593849A (zh) * | 2019-08-22 | 2019-12-20 | 魏峰 | 一种考虑相渗曲线动态变化的油藏数值模拟方法 |
CN111488974A (zh) * | 2020-04-14 | 2020-08-04 | 四川北控清洁能源工程有限公司 | 基于深度学习神经网络的海洋风能降尺度方法 |
Non-Patent Citations (2)
Title |
---|
李维勃等: "基于径向基函数的冷却塔风场重构", vol. 36, no. 36, pages 226 - 234 * |
邹亮等: "海上风电机组基础结构选型程序开发", vol. 42, no. 42, pages 282 - 285 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114004436A (zh) * | 2021-06-10 | 2022-02-01 | 重庆大学 | 一种基于性能反应谱的自复位海上风电升压站抗震性能评估方法 |
CN114004436B (zh) * | 2021-06-10 | 2024-05-28 | 重庆大学 | 一种基于性能反应谱的自复位海上风电升压站抗震性能评估方法 |
CN116306109A (zh) * | 2023-02-03 | 2023-06-23 | 哈尔滨工业大学(深圳) | 基于时域模型修正的海上风机土结相互作用状态识别方法 |
CN116306109B (zh) * | 2023-02-03 | 2023-10-20 | 哈尔滨工业大学(深圳) | 基于时域模型修正的海上风机土结相互作用状态识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112749429B (zh) | 2023-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | A review of the impact of hydropower reservoirs on global climate change | |
Wang et al. | Assessment of the impact of climate change on hydropower potential in the Nanliujiang river basin of China | |
Zaharim et al. | Analyzing Malaysian wind speed data using statistical distribution | |
CN109740286A (zh) | 一种混合优化bp神经网络的湖泊水质预测模型构建方法 | |
CN117390894B (zh) | 一种预测极端风暴潮位的方法 | |
CN112766593B (zh) | 一种水利水电工程生态调度方案的优化方法 | |
CN111723505A (zh) | 一种流域水质水量监测系统 | |
CN114580152A (zh) | 一种基于多体耦合分析的浮式风电结构基础局部应力时域分析方法 | |
CN112749429A (zh) | 一种基于多元插值的海上风电固定式基础结构设计方法 | |
CN112818437A (zh) | 海上风电单桩基础设计计算优化出图的集成化分析方法 | |
Abdraman et al. | Wind Resource Assessment in the City of N'djamena in Chad | |
CN104933469A (zh) | 一种基于灰色生成扰动模型的短期风速预测方法 | |
CN112287301B (zh) | 一种基于全球水文模型的全球尺度水力发电能力的评估方法 | |
CN116341952A (zh) | 一种海洋风浪资源联合开发评估方法 | |
CN104331748B (zh) | 一种预测规划目标年风电场群持续功率曲线的方法 | |
Paturi et al. | Artificial neural networks modelling for power coefficient of Archimedes screw turbine for hydropower applications | |
CN112670997B (zh) | 考虑光伏不确定性的电热能源系统时序概率潮流计算方法 | |
Khayatnezhad et al. | Multivariate analysis of wind characteristics for optimal irrigation planning in Miandoab Plain, Urmia Lake | |
CN103065226B (zh) | 水电站水库长期优化调度决策水位的确定方法 | |
Huang et al. | Solar Radiation Forecasting based on Neural Network in Guangzhou | |
A Francis et al. | Statistical analysis of wind speed and evaluation of wind power density for Colaba, Mumbai | |
CN106326540A (zh) | 水电站水能动态模拟算法 | |
Lin et al. | On spatial distribution and determinants of urban photovoltaic utilization in China | |
Yuan et al. | Prediction of China's energy consumption structure | |
Wang et al. | Technological economic study for ocean energy development in China |
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 |