CN113536609B - 面向海上风电工程的雷达遮蔽区域模型构建方法 - Google Patents

面向海上风电工程的雷达遮蔽区域模型构建方法 Download PDF

Info

Publication number
CN113536609B
CN113536609B CN202111084672.1A CN202111084672A CN113536609B CN 113536609 B CN113536609 B CN 113536609B CN 202111084672 A CN202111084672 A CN 202111084672A CN 113536609 B CN113536609 B CN 113536609B
Authority
CN
China
Prior art keywords
radar
driven generator
wind driven
wind
area
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
Application number
CN202111084672.1A
Other languages
English (en)
Other versions
CN113536609A (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.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
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 Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN202111084672.1A priority Critical patent/CN113536609B/zh
Publication of CN113536609A publication Critical patent/CN113536609A/zh
Application granted granted Critical
Publication of CN113536609B publication Critical patent/CN113536609B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

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)
  • Wind Motors (AREA)

Abstract

本发明公开了一种面向海上风电工程的雷达遮蔽区域模型构建方法，建立了海上风力发电机叶片旋转遮挡模型，确定了风力发电机叶片旋转时雷达探测遮蔽区域面积的计算方法；然后根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，对雷达的主要探测特性能及可能产生的影响进行计算和分析，有助于掌握海上风电场对雷达探测信号的影响程度，可为海上风电场选址、布置及优化提供依据，同时也有利于海事监管部门精准识别雷达盲区及重点监测通航水域，可有效减小雷达探测盲区的通航安全隐患，为海上风电场水域的多形式协同安全监管提供参考。

Description

面向海上风电工程的雷达遮蔽区域模型构建方法
技术领域
本发明涉及到海上风电场技术及雷达技术领域，具体涉及到面向海上风电工程的雷达遮蔽区域模型构建方法。
背景技术
我国大陆海岸线长达1.8万km，海上风能资源储备丰富，海上风电产业发展优势显著。据统计，我国近海水深在5～25m范围内50m高度风电开发潜力约200GW，水深在25～50m范围内的风电资源可开发量则达320GW。海上风电年平均利用小时数约2500h，比陆上风电高出约500h。虽然我国海上风电产业起步较晚，但凭借海上资源丰富、稳定等优势，近年来已进入发展快车道。自2010年我国首座、也是亚洲首座大型海上风电场——东海大桥海上风电场并网发电以来，经10多年的探索与快速发展，如今我国已成为名副其实的海上风电大国。2020年，我国海上风电新增装机量占全球的50.45%，位居世界第一。截至2020年底，我国海上风电累计装机容量占全球的28.12%，超过德国，仅次于英国，成为全球第二大海上风电市场。
根据《海上风电开发建设管理办法》的相关规定，海上风电场应当按照生态文明建设要求，统筹考虑开发强度和资源环境承载能力，原则上应在离岸距离不少于10km、滩涂宽度超过10km时海域水深不得少于10m的海域布局。通常，海上风电工程选址水深条件较好，其需做好自身保障与防护的同时，也需配备有效的监管、警示设备，防止船舶误入海上风电场。由于海上风电场通常占地范围广、风力发电机尺寸大且转速高，海上风电场建设有可能影响附近雷达探测海上目标的能力，对于海事安全监管造成一定的影响。近年来，由于在海上风电工程规划、设计阶段对雷达探测影响分析不足，而在被迫在后期调整海上风电场布置方案的工程案例频频出现，这将对海上风电工程建设实施进度造成巨大影响，也会造成前期投入资源的浪费。因此，在科学、合理布置海上风电场的过程中，精准评估海上风电场工程对雷达探测性能的影响程度显得至关重要。
中国发明专利申请（公开号：CN111796244A）在2020年公开了一种海上风电场对地波超视距雷达探测效能影响的评估方法，根据地波超视距雷达参数以及风电场的预选区域展开场景分析，筛选出需重点分析的对象；根据提供的风电机组参数建立风机电磁模型，实现单台风力发电机电磁模型的建立；按照拟建风电场风机经纬度坐标与雷达相对位置信息，在FEKO电磁计算软件中建阵列、建场景模型，对遮挡影响进行分析；计算单台风力发电机设备在不同状态下的RCS变化；多普勒频率的影响进行分析；推导多径效应理论公式；场景实测等效，电磁环境噪声数据，评估出风电场影响范围。该发明通过电磁模型的建立能够评价风电场的对电磁波的影响，但是对于风力发电机叶片旋转、风力发电机本身对雷达遮蔽的影响分析较少，而且其评价结果难以作为海上风电工程布置及航道布置的全面参考。
《雷达与对抗》2016年第36卷第1期公开了论文“海上风电场对岸基警戒雷达遮挡影响计算分析”，文中对风电场对雷达遮挡影响的理论计算分析主要利用雷达波绕射模型计算雷达的遮挡衰减值，然后从对岸基雷达距离、方位和高度等3个方面计算分析遮挡影响的结果，从而评估风电场工程对雷达系统的遮挡影响效果。这种方式主要用于对已知某风电场的评估，没有考虑风力发电机叶片旋转的影响及最大影响距离计算，不能够在海上风电场规划阶段优化海上风力发电机的布置。
发明内容
本发明的目的是针对现有技术存在的问题，提供面向海上风电工程的雷达遮蔽区域模型构建方法，利用几何法分析海上风电场对雷达遮挡的影响。
为实现上述目的，本发明采用的技术方案是：
面向海上风电工程的雷达遮蔽区域模型构建方法，包括设置在海上风电场中的风力发电机，以及设置在所述海上风电场附近的航道、雷达站；所述模型构建方法包括如下步骤：
步骤一：获取所述雷达站中雷达的类型及基本参数，获取所述海上风电场中风力发电机的设计参数；
步骤二：评估所述风力发电机叶片旋转对所述雷达电磁波遮挡的影响，首先确定雷达扫掠单台风力发电机时叶片产生的旋转角度
Figure 393396DEST_PATH_IMAGE001
，其次结合雷达旋转角度
Figure 257447DEST_PATH_IMAGE002
确定风力发电机叶片旋转遮挡雷达电磁波的面积S，并得到雷达扫描单台风力发电机时叶片所在平面内一定不会产生遮挡区域的百分比
Figure 275081DEST_PATH_IMAGE003
；
步骤三：评估步骤二中风力发电机叶片旋转对所述雷达电磁波遮挡的影响后，根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，包括确定雷达由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离d和阴影面积
Figure 97544DEST_PATH_IMAGE004
，得到风力发电机对航道影响的盲区距离e和风力发电机对雷达探测影响的角度α；
Figure 830489DEST_PATH_IMAGE005
，
其中，
Figure 498230DEST_PATH_IMAGE006
为风力发电机塔筒的半径，d1为雷达与风力发电机间的水平距离，d2为风力发电机与航道间的水平距离。
步骤四：根据雷达的基本参数，以及可接受的雷达探测盲区大小和影响角度，反推风力发电机与雷达站之间的距离和风力发电机与航道之间的距离，确定风力发电机和航道设置的位置是否合理，以及是否需要增设其他形式协同监管设备。
对于已规划或拟建的海上风电场，根据步骤二和步骤三中的计算模型分析风力发电机叶片和风力发电机对雷达遮蔽以及航道盲区的影响，判断是否需要优化风力发电机位置或是否需要补充协同监管设施。
本构建方法通过对雷达和风力发电机参数的分析，建立了海上风力发电机叶片旋转遮挡模型，确定了风力发电机叶片旋转时雷达探测遮蔽区域面积的计算方法；然后根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，对雷达的主要探测特性能及可能产生的影响进行计算和分析，有助于掌握海上风电场对雷达探测信号的影响程度，可为海上风电场选址、布置及优化提供依据，同时也有利于海事监管部门精准识别雷达盲区及重点监测通航水域，可有效减小雷达探测盲区的通航安全隐患，为海上风电场水域的多形式协同安全监管提供参考。
由于风力发电机叶片的尺寸较大而且通常会持续旋转，因此先行判断风力发电机叶片旋转对雷达遮蔽的影响，再进行风力发电机塔筒对雷达遮蔽的影响分析，在确定出风力发电机对航道影响的盲区距离e、风力发电机对雷达探测影响的角度α以及雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 635951DEST_PATH_IMAGE004
后，进而根据雷达实际扫掠情况，优化风力发电机布置或提出雷达探测盲区水域的监管设施补充方案。
进一步的，所述雷达的参数包括工作波长、脉冲重复频率、脉冲宽度、波束宽度和接收机通频带，并计算得出距离分辨力
Figure 363735DEST_PATH_IMAGE007
和方位分辨力
Figure 852485DEST_PATH_IMAGE008
。
进一步的，步骤二中叶片旋转角度
Figure 323918DEST_PATH_IMAGE001
的确定方法如下：
假设单台风力发电机单个叶片遮挡角度为
Figure 316145DEST_PATH_IMAGE009
，当雷达扫过单台风力发电机时叶片会产生的旋转角度
Figure 214831DEST_PATH_IMAGE001
，则会受到遮挡影响区域的角度为
Figure 925298DEST_PATH_IMAGE010
，一定不会受到遮挡影响区域的角度为
Figure 934842DEST_PATH_IMAGE011
，假设
Figure 312734DEST_PATH_IMAGE002
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 116742DEST_PATH_IMAGE012
，
由于
Figure 580084DEST_PATH_IMAGE013
，故
Figure 127740DEST_PATH_IMAGE014
，则
Figure 360138DEST_PATH_IMAGE015
，
其中
Figure 335048DEST_PATH_IMAGE016
为风力发电机叶片旋转角速度，
Figure 285686DEST_PATH_IMAGE017
为雷达扫掠角速度，r为风力发电机叶片长度，d为雷达与风力发电机间距离。
进一步的，在步骤二中风力发电机叶片旋转遮挡雷达电磁波的面积的确定方法如下：
首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 637033DEST_PATH_IMAGE001
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E；
求得单台风力发电机单个叶片遮挡雷达电磁波的面积SABCD为：
Figure 723938DEST_PATH_IMAGE018
，
则风力发电机叶片旋转时遮挡雷达电磁波的面积S表示为：
Figure 135327DEST_PATH_IMAGE019
；
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线
Figure 307683DEST_PATH_IMAGE020
与
Figure 217649DEST_PATH_IMAGE021
构建函数表达式：
Figure 159060DEST_PATH_IMAGE022
，
联立可解得E的位置为：
Figure 741351DEST_PATH_IMAGE023
，
SBCDE可由下式获得：
Figure 135424DEST_PATH_IMAGE024
，
其中，SDEG及SBCG可根据交点E及辅助点G的位置计算可得：
Figure 94152DEST_PATH_IMAGE025
。
进一步的，根据雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，计算雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 890070DEST_PATH_IMAGE003
；
Figure 908842DEST_PATH_IMAGE026
,
S为风力发电机叶片旋转时遮挡雷达电磁波的面积，将获取的w、r、h、n，以及计算得到的
Figure 790210DEST_PATH_IMAGE027
等数据代入SDEG和SBCG的计算公式中，得到SBCDE，进一步计算得到S，从而可以得到雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 21471DEST_PATH_IMAGE003
的数值；通过
Figure 937475DEST_PATH_IMAGE003
的数值进行分析，可以得到风力发电机叶片旋转对雷达电磁波的遮蔽影响，初步判断风力发电机对雷达探测的影响程度。
进一步的，分别构建d1（雷达与风力发电机间的距离）与
Figure 861568DEST_PATH_IMAGE001
、d1与
Figure 495812DEST_PATH_IMAGE003
的变化关系曲线图，以便直观明了的获取随距离变化风力发电机叶片旋转对雷达遮蔽影响程度。
与现有技术相比，本发明的有益效果是：1、建立风力发电机叶片旋转遮挡模型和雷达扫描风力发电机塔筒时形成盲区的计算模型，有助于掌握海上风电场对雷达探测信号的影响程度，可为海上风电场选址、布置及优化提供依据；2、为雷达站、航道等的布设提供参考，有效减小海上风电场遮蔽区域对雷达探测的影响；3、有助于海事监管部门精准识别雷达盲区及重点监测通航水域，可有效减小雷达探测盲区的通航安全隐患，为海上风电场水域的多形式协同安全监管提供参考。
附图说明
图1为本发明风力发电机旋转示意图；
图2为本发明风力发电机抽象模型旋转示意图；
图3为本发明雷达扫掠过程中风力发电机叶片旋转角度与雷达-风力发电机间距离变化关系；
图4为本发明风力发电机叶片平面一定不受遮挡影响区域百分比随雷达-风力发电机间距离变化关系；
图5为本发明海上风电场对雷达的直线遮挡侧视示意图；
图6为本发明海上风电场对雷达的直线遮挡俯视示意图。
具体实施方式
下面将结合本发明中的附图，对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动条件下所获得的所有其它实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，术语“中间”、“上”、“下”、“左”、“右”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
一、雷达主要参数的确定
选取某VTS雷达作为分析，主要参数如下：
（1）工作波长λ
雷达的工作波长即发射机产生的超高频脉冲波的波长。VTS雷达站雷达的工作频段为9375MHz，根据波长的计算公式可计算VTS雷达站的工作波长λ为3.2cm；船载雷达使用的工作波长允许范围为：S波段7.5～15cm，X波段2.5～3.75cm。
（2）脉冲重复频率PRF
雷达信号每秒内脉冲重复出现的次数，称为脉冲重复频率，即为PRF，与其对应的是相邻两个脉冲间的时间间隔，即为脉冲重复周期PRP，其对应关系为PRF=1/PRP。
（3）脉冲宽度
Figure 530764DEST_PATH_IMAGE028
脉冲宽度就是射频脉冲振荡持续的时间，一般用
Figure 35695DEST_PATH_IMAGE028
表示。一般选择0.05～2μs之内。雷达通常可按照近、中、远距离档以改换发射电磁波的脉冲宽度。
（4）波束宽度
在雷达波束水平面内的波束宽度成为水平波束宽度，用
Figure 396269DEST_PATH_IMAGE029
表示，它与雷达天线的结构、尺寸及使用波长有关。根据《某VTS用户指南》可知基站雷达的水平波束宽度≤0.45°，其后计算参照的最大值为0.45°。
（5）接收机通频带
接收机通频带表示接收机能有效放大的信号频率范围，一般指接收机对输入信号的电压放大倍数从中心频率
Figure 252230DEST_PATH_IMAGE030
的最大相对值1下降到0.707时的两个频率的差
Figure 90873DEST_PATH_IMAGE031
。现代船用雷达接收机通频带为3~25MHz。
（6）距离分辨力
雷达的距离分辨力是指雷达分辨同方向两个相邻目标的能力。同方向的两个目标逐渐靠拢时，两个目标的回波也会逐渐接近而至相切，此时两个目标间的实际距离即距离分辨力。当同方位的两个物标逐渐靠拢时，雷达屏上两个物标的回波亮点也将逐渐接近，当两个回波亮点相切时，两物标间的实际距离即为雷达的距离分辨力
Figure 715889DEST_PATH_IMAGE032
。
Figure 247364DEST_PATH_IMAGE033
，
其中，C为电磁波的传播速度，
Figure 325042DEST_PATH_IMAGE028
为发射脉冲宽度，
Figure 230025DEST_PATH_IMAGE034
为接收机通频带，
Figure 709548DEST_PATH_IMAGE035
为光点直径，D为屏幕直径，RD为所用量程距离。
雷达的距离分辨力主要取决于发射脉冲宽度、接收机通频带、量程以及屏幕光点尺寸大小等技术参数，通过查询获取某VTS雷达在不同量程下的距离分辨力如表1所示。
表1 雷达选取不同量程时的距离分辨力
Figure DEST_PATH_IMAGE036
（7）方位分辨力
雷达的方位分辨力是指雷达区分距离相同而方位相邻的两个目标的能力，以能区分开的两个目标间的最小夹角来表示。此夹角越小，表示雷达的方位分辨力越高，雷达图像的清晰度越好，以能分辨的两个物标间的最小方位夹角
Figure 880767DEST_PATH_IMAGE037
来表示，
Figure 711319DEST_PATH_IMAGE038
越小，表示雷达方位分辨力越高。
影响方位分辨力的主要因素是天线水平波束宽度
Figure 626186DEST_PATH_IMAGE039
、光点角尺寸
Figure 225794DEST_PATH_IMAGE040
（光点直径对屏中心的张角）和回波在屏幕扫描线上的位置，具体雷达的方位分辨力可表示为：
Figure 99072DEST_PATH_IMAGE041
，
基站雷达的水平波束宽度
Figure 151342DEST_PATH_IMAGE039
≤0.45°，光点尺寸
Figure 135478DEST_PATH_IMAGE040
约0.06°。
二、风力发电机参数的确定
某风力发电机底座半径为2.5m，轮毂高度90m，叶片长度为63.4m。
三、风力发电机叶片旋转对雷达电磁波的遮蔽模型构建
实际情况中，风力发电机叶片会对雷达电磁波产生遮挡影响，并且由于风力发电机叶片是处于旋转状态，在雷达扫掠单台风力发电机过程中，风力发电机叶片平面内只有风力发电机叶片对雷达探测存在遮挡影响，其他区域无遮挡。
如图1所示，单台风力发电机单个叶片遮挡角度为
Figure 589594DEST_PATH_IMAGE009
，当雷达扫过单台风力发电机时叶片会产生一定的旋转角度
Figure 633773DEST_PATH_IMAGE042
，则会受到遮挡影响区域的角度为
Figure 173339DEST_PATH_IMAGE010
，一定不会受到遮挡影响区域的角度为
Figure 961166DEST_PATH_IMAGE011
。由于单台风力发电机的3个叶片形状大小完全相同，因此首先单独考虑1个风力发电机叶片对雷达的影响情况：
假设
Figure 269788DEST_PATH_IMAGE043
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 484868DEST_PATH_IMAGE044
（1）
由于
Figure 246151DEST_PATH_IMAGE045
（2）
故
Figure 572090DEST_PATH_IMAGE046
（3）
Figure 797DEST_PATH_IMAGE047
（4）
其中
Figure 386779DEST_PATH_IMAGE048
为风力发电机叶片旋转角速度90（°/s），
Figure 900937DEST_PATH_IMAGE049
为雷达扫掠角速度
Figure 762058DEST_PATH_IMAGE049
=120（°/s），r为风力发电机叶片长度，d为雷达与风力发电机间距离。
进一步分析风力发电机叶片旋转情况下遮挡雷达电磁波的面积，如图2所示，首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 45272DEST_PATH_IMAGE042
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E。
因此，可求得单个风力发电机叶片遮挡的影响区域S ABCD 为：
Figure 602156DEST_PATH_IMAGE050
（5）
由于雷达扫描时风力发电机叶片的旋转角度为
Figure 603610DEST_PATH_IMAGE042
，则风力发电机叶片旋转的遮挡影响区域可表示为：
Figure 5772DEST_PATH_IMAGE051
（6）
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线LA′D与LAB构建函数表达式：
Figure 409072DEST_PATH_IMAGE052
（7）
Figure 136856DEST_PATH_IMAGE053
（8）
将式（6）与（7）联立可解得E的位置为：
Figure 360027DEST_PATH_IMAGE054
根据图2所示，SBCDE可由下式获得：
Figure 831460DEST_PATH_IMAGE055
（9）
其中，SDEG及SBCG可根据交点E及辅助点G的位置计算可得：
Figure 89266DEST_PATH_IMAGE056
（10）
Figure 722372DEST_PATH_IMAGE057
（11）
将公式（10）和（11）代入公式（9），并将公式（9）代入公式（5），同时将雷达扫掠单台风力发电机的旋转角度
Figure 698419DEST_PATH_IMAGE042
也代入公式（5）可计算雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，同时可得雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域面积百分比
Figure 707963DEST_PATH_IMAGE058
为：
Figure 820275DEST_PATH_IMAGE059
（12）
将获取的w、r、h、n，以及计算得到的
Figure 889863DEST_PATH_IMAGE042
等数据代入SDEG和SBCG的计算公式中，得到SBCDE的数值，进一步计算得到S的值，从而可以得到雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 822047DEST_PATH_IMAGE003
。从而根据计算可以得到雷达与海上风电场某边界点距离d1，并得到
Figure 900861DEST_PATH_IMAGE042
与
Figure 867680DEST_PATH_IMAGE003
。
Figure 108168DEST_PATH_IMAGE042
随d的变化关系见图3，
Figure 793228DEST_PATH_IMAGE003
随d的变化关系见图4。
由图3可知，雷达扫掠某台风力发电机时间内，叶片旋转角度与雷达-风力发电机间的距离成反比，由于雷达与风力发电机间距离从18～30km不等，雷达扫掠过程中风力发电机叶片旋转角度从0.18°～0.30°不等，在本实施例中雷达扫掠风力发电机时，风力发电机旋转角度非常小，对雷达探测影响不大。
由图4可知，随着雷达与风力发电机间距离变大，风力发电机叶片产生遮挡区域百分比减少，由于雷达与风力发电机间距离变化范围为18～30km，则风力发电机叶片所在平面一定不会产生遮挡区域百分比变化范围为95.23%～95.33%。可知在本实施例范围内单台风力发电机叶片相对其整个平面而言遮挡效应有限。但是，实际情况中，应结合海上风电场中各风力发电机布置，对某通航区域的遮挡影响进行叠加分析，进而结合叶片遮挡区域可接受范围，综合评估海上风电场中对雷达的影响。
四、雷达扫描盲区分析
以上述的雷达站为例，基于雷达波在空间匀速直线传播的特性，利用几何法分析海上风电场对雷达遮挡的影响。
通常，雷达与风力发电机间存在一定的距离，且雷达高度高于风力发电机高度。根据雷达波直线传递特性，雷达波对航道的覆盖范围会受到风力发电机的塔筒遮挡，形成雷达电磁波的阴影径向距离，部分航道区域会处于阴影范围，雷达与风力发电机间的位置关系示意图如图5和6所示。
根据雷达、海上风电场设计参数及布置方案，假设雷达高度为h1，风力发电机轮毂高度为h2，雷达与风力发电机间的水平距离为d1，则可推出雷达由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离d（最大遮挡影响距离）为：
Figure 410154DEST_PATH_IMAGE060
（13）
风力发电机会对雷达的扫描范围产生遮挡阴影区域，在已分析风力发电机叶片产生的遮挡影响后，再分析风力发电机塔筒对雷达的遮挡影响，以风力发电机轮毂直径来计算风力发电机对雷达电磁波的遮挡范围。设风力发电机塔筒半径为
Figure 228550DEST_PATH_IMAGE061
，为风力发电机与航道间的水平距离为d2，则雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 639939DEST_PATH_IMAGE062
（投射在海面上的阴影面积）为：
Figure 546716DEST_PATH_IMAGE063
（14）
Figure 967333DEST_PATH_IMAGE064
（15）
Figure 643165DEST_PATH_IMAGE065
（16）
其中，x1、x2表示计算中间变量，风力发电机对航道影响的盲区（宽度）距离e为：
Figure 491035DEST_PATH_IMAGE066
（17）
风力发电机对雷达探测影响的角度为：
Figure 885107DEST_PATH_IMAGE067
（18）
根据雷达的方位分辨力性能，不超过方位分辨力的方位遮挡对雷达的方位探测并无影响。对航道而言，海上风力发电机对雷达信号的多次反射、阻挡等会不同程度影响雷达性能和探测范围，影响雷达观测物标的参数主要有雷达的方位遮挡程度、形成遮蔽区域的大小等，可根据雷达电磁波阴影径向距离和盲区距离计算对海域通航的影响。
实施例一：
以某海域设立的海上风电场为例，进行雷达探测盲区计算结果及分析，风力发电机间距离雷达站20~30km。
利用直线传播模型分析海上风电工程对某航路的遮挡，风力发电机底座半径为2.5m，风力发电机轮毂高度90m，风力发电机塔筒半径2.5m左右，风力发电机叶片长度为63.4m，某海域附近的多个岸基雷达站的雷达天线高程分别为150m、190m、220m、180m。针对海上风电场与不同海事监管（VTS）雷达的相对位置关系以及对航道的遮挡情况不同，对于不同VTS雷达，我们选取影响相对较大的风力发电机进行遮挡影响分析。分析了1#、2#、3#、4#号风力发电机对雷达天线高程150m的雷达的遮挡影响，影响结果见表2。
表2 20~30km处风力发电机对某雷达站的影响
Figure 578257DEST_PATH_IMAGE068
根据表2可知，由于某雷达站与海上风电场距离较远，风力发电机对雷达的最大影响距离都超过39km，再将风力发电机与雷达的距离考虑在内，风力发电机对雷达的遮挡影响基本都达到雷达最大探测距离25海里（约46km），但由于风力发电机底座直径较小，对雷达的影响角度在0.02°左右，换算为对航道内产生的遮挡盲区宽度，最大结果为5.8m左右，因此理论计算的阴影盲区较小。但是，实际情况中，应结合海上风电场中各风力发电机布置，对某通航区域的遮挡影响进行叠加分析，进而结合塔筒遮挡区域可接受范围，综合评估海上风电场中对雷达的影响。
实施例二：
为了更好理解不同距离下风力发电机产生盲区的影响，本实施例与实施例一的区别在于将风力发电机设置在距离雷达站8~10km处，其编号分别为5#、6#、7#和8#。
表3 8~10km处风力发电机对某雷达站的影响
Figure 374174DEST_PATH_IMAGE069
根据表3知，由于该雷达站与海上风电场距离较近，风力发电机对雷达的影响较为明显，航道盲区和影响角度均大于实施例一。同理，对于更近距离的其他雷达站，或者距离其他航道更近的风力发电机可类比按照此方法计算分析。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (7)

1.面向海上风电工程的雷达遮蔽区域模型构建方法，包括设置在海上风电场中的风力发电机，以及设置在所述海上风电场附近的航道、雷达站；其特征在于，所述模型构建方法包括如下步骤：
步骤一：获取所述雷达站中雷达的类型及基本参数，获取所述海上风电场中风力发电机的设计参数；
步骤二：评估所述风力发电机叶片旋转对所述雷达电磁波遮挡的影响，首先确定雷达扫掠单台风力发电机时叶片产生的旋转角度
Figure 44163DEST_PATH_IMAGE001
，其次结合雷达旋转角度
Figure 637955DEST_PATH_IMAGE002
确定单台风力发电机单个叶片旋转遮挡雷达电磁波的面积S，并得到雷达扫描单台风力发电机时叶片所在平面内一定不会产生遮挡区域的百分比
Figure 650386DEST_PATH_IMAGE003
；
步骤三：评估步骤二中风力发电机叶片旋转对所述雷达电磁波遮挡的影响后，根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，包括确定由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离d和雷达扫掠风力发电机时产生的遮挡阴影面积
Figure 355037DEST_PATH_IMAGE004
，得到风力发电机对航道影响的盲区距离e和风力发电机对雷达探测影响的角度α；
Figure DEST_PATH_IMAGE005
，
其中，
Figure 772243DEST_PATH_IMAGE006
为风力发电机塔筒的半径，d1为雷达与风力发电机间的水平距离，d2为风力发电机与航道间的水平距离；
步骤四：根据雷达的基本参数，以及可接受的雷达探测盲区大小和影响角度，反推风力发电机与雷达站之间的距离和风力发电机与航道之间的距离，确定风力发电机和航道设置的位置是否合理，以及是否需要增设其他形式协同监管设备。
2.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，所述雷达的参数包括工作波长、脉冲重复频率、脉冲宽度、波束宽度和接收机通频带，并计算得出距离分辨力
Figure 271357DEST_PATH_IMAGE007
和方位分辨力
Figure 101910DEST_PATH_IMAGE008
。
3.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，步骤二中叶片旋转角度
Figure 220039DEST_PATH_IMAGE001
的确定方法如下：
假设单台风力发电机单个叶片遮挡角度为
Figure 881964DEST_PATH_IMAGE009
，当雷达扫过单台风力发电机时叶片会产生的旋转角度
Figure 427346DEST_PATH_IMAGE010
，则会受到遮挡影响区域的角度为
Figure 541933DEST_PATH_IMAGE011
，一定不会受到遮挡影响区域的角度为
Figure 526069DEST_PATH_IMAGE012
，假设
Figure 917867DEST_PATH_IMAGE002
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 24364DEST_PATH_IMAGE013
；
由于
Figure 501612DEST_PATH_IMAGE014
，故
Figure 289440DEST_PATH_IMAGE015
，则
Figure 660378DEST_PATH_IMAGE017
；
其中
Figure 547563DEST_PATH_IMAGE018
为风力发电机叶片旋转角速度，
Figure 636742DEST_PATH_IMAGE019
为雷达扫掠角速度，r为风力发电机叶片长度，d为雷达与风力发电机间距离。
4.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，在步骤二中风力发电机叶片旋转遮挡雷达电磁波的面积S的确定方法如下：
首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 900364DEST_PATH_IMAGE001
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E；
求得单台风力发电机单个叶片遮挡雷达电磁波的面积SABCD为：
Figure 329071DEST_PATH_IMAGE020
，
则风力发电机叶片旋转时遮挡雷达电磁波的面积S表示为：
Figure 511791DEST_PATH_IMAGE021
；
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线LA′D与LAB构建函数表达式：
Figure 226281DEST_PATH_IMAGE022
，
联立可解得E点的位置为：
Figure 152649DEST_PATH_IMAGE023
，
SBCDE可由下式获得：
Figure 107967DEST_PATH_IMAGE024
，
其中，SDEG及SBCG可根据交点E及辅助点据G的位置计算可得：
Figure 930429DEST_PATH_IMAGE025
。
5.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，根据雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，计算雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 728621DEST_PATH_IMAGE003
；
Figure 599625DEST_PATH_IMAGE026
；其中r为风力发电机叶片长度。
6.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，根据步骤二中获取的雷达旋转角度
Figure 534083DEST_PATH_IMAGE001
和一定不会产生遮挡区域百分比
Figure 527447DEST_PATH_IMAGE003
,分别构建
Figure 688301DEST_PATH_IMAGE001
、
Figure 222050DEST_PATH_IMAGE003
随雷达与风力发电机间距离d1的变化关系曲线。
7.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，假设雷达高度为h1，风力发电机轮毂高度为h2，风力发电机塔筒的半径
Figure 151960DEST_PATH_IMAGE006
，则可推出由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离
Figure 316225DEST_PATH_IMAGE027
，并可推出雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 823430DEST_PATH_IMAGE028
。
CN202111084672.1A 2021-09-16 2021-09-16 面向海上风电工程的雷达遮蔽区域模型构建方法 Active CN113536609B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111084672.1A CN113536609B (zh) 2021-09-16 2021-09-16 面向海上风电工程的雷达遮蔽区域模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111084672.1A CN113536609B (zh) 2021-09-16 2021-09-16 面向海上风电工程的雷达遮蔽区域模型构建方法

Publications (2)

Publication Number Publication Date
CN113536609A CN113536609A (zh) 2021-10-22
CN113536609B true CN113536609B (zh) 2021-11-30

Family

ID=78123225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111084672.1A Active CN113536609B (zh) 2021-09-16 2021-09-16 面向海上风电工程的雷达遮蔽区域模型构建方法

Country Status (1)

Country Link
CN (1) CN113536609B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114510846B (zh) * 2022-04-18 2022-07-22 天津航大天元航空技术有限公司 一种风力发电场的安全评估方法、装置及电子设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106121913A (zh) * 2016-06-20 2016-11-16 广东科诺勘测工程有限公司 一种海上风电场的风机排布方法
CN111796244A (zh) * 2020-06-09 2020-10-20 哈尔滨工业大学 一种海上风电场对地波超视距雷达探测效能影响的评估方法
CN112502911A (zh) * 2020-11-30 2021-03-16 东方电气风电有限公司 一种实时预测叶片通过塔筒时扫塔风险的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DK2877741T3 (da) * 2012-07-27 2019-05-20 Univ Texas Tech System System og fremgangsmåde til at evaluere vindstrømningsfelter ved anvendelse af fjernregistreringsindretninger
CN111610494B (zh) * 2020-05-27 2022-08-19 武汉理工大学 一种vts雷达配置信号覆盖优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106121913A (zh) * 2016-06-20 2016-11-16 广东科诺勘测工程有限公司 一种海上风电场的风机排布方法
CN111796244A (zh) * 2020-06-09 2020-10-20 哈尔滨工业大学 一种海上风电场对地波超视距雷达探测效能影响的评估方法
CN112502911A (zh) * 2020-11-30 2021-03-16 东方电气风电有限公司 一种实时预测叶片通过塔筒时扫塔风险的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Raza M B 等.Impact of off-shore wind turbine on forward scattering of marine radar signals.《2017 IEEE MTT-S International Microwave Symposium (IMS)》.2017, *
刘克中 等.海上风电场对航海雷达探测性能影响研究.《武汉理工大学学报(交通科学与工程版)》.2010,第34卷(第3期), *

Also Published As

Publication number Publication date
CN113536609A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN111796244B (zh) 一种海上风电场对地波超视距雷达探测效能影响的评估方法
Poupart Wind farms impact on radar aviation interests-final report
EP2877741B1 (en) System and method for evaluating wind flow fields using remote sensing devices
CN104215946B (zh) 一种天地波混合雷达回波谱仿真方法
CN106845018B (zh) 风电场对气象雷达降雨量影响的分析与定量化评估方法
CN102112892A (zh) 雷达系统和方法
Norin A quantitative analysis of the impact of wind turbines on operational Doppler weather radar data
Norin et al. Doppler weather radars and wind turbines
CN106019286B (zh) 一种机载气象雷达多扫描多重频全空域气象目标探测方法
CN113536609B (zh) 面向海上风电工程的雷达遮蔽区域模型构建方法
CN110988884B (zh) 一种基于高频地波雷达的中纬度电离层探测方法
CN109001846A (zh) 一种复杂地形条件下s波段与x波段雷达组网测雨方法
Lorente et al. Evaluating the surface circulation in the Ebro delta (northeastern Spain) with quality-controlled high-frequency radar measurements
CN110726980A (zh) 一种机载相控阵天气雷达对地表杂波的分析方法
Angulo et al. Estimating reflectivity values from wind turbines for analyzing the potential impact on weather radar services
JP5238531B2 (ja) レーダ装置、海洋レーダ観測装置およびドップラ周波数データ算出方法
CN104199002B (zh) 抑制风力发电旋叶虚警目标的雷达信息处理方法
CN111736150A (zh) 一种面向远距离低功耗探鸟雷达的探测方法
CN115032605A (zh) 一种基于目标能力画像的雷达目标体系分析方法
JP2005351853A (ja) 光波レーダ装置
CN108594220B (zh) 近海风电场水域的海事雷达图像模拟系统及方法
Heathershaw et al. Wave direction estimates in coastal waters using radar
CN114114253A (zh) 一种基于人工智能算法的海浪监测阵列雷达
Sozen et al. Scatter and Doppler effect of wind power plants to land radars
CN110261836B (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