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

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

Info

Publication number
CN113536609A
CN113536609A CN202111084672.1A CN202111084672A CN113536609A CN 113536609 A CN113536609 A CN 113536609A CN 202111084672 A CN202111084672 A CN 202111084672A CN 113536609 A CN113536609 A CN 113536609A
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.)
Granted
Application number
CN202111084672.1A
Other languages
English (en)
Other versions
CN113536609B (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 354366DEST_PATH_IMAGE001
，其次结合雷达旋转角度
Figure 599403DEST_PATH_IMAGE002
确定风力发电机叶片旋转遮挡雷达电磁波的面积S，并得到雷达扫描单台风力发电机时叶片所在平面内一定不会产生遮挡区域的百分比
Figure 175878DEST_PATH_IMAGE003
；
步骤三：评估步骤二中风力发电机叶片旋转对所述雷达电磁波遮挡的影响后，根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，包括确定雷达由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离d和阴影面积
Figure 633404DEST_PATH_IMAGE004
，得到风力发电机对航道影响的盲区距离e和风力发电机对雷达探测影响的角度α；
Figure 572410DEST_PATH_IMAGE005
，
其中，
Figure 988348DEST_PATH_IMAGE006
为风力发电机塔筒的半径，d1为雷达与风力发电机间的水平距离，d2为风力发电机与航道间的水平距离。
步骤四：根据雷达的基本参数，以及可接受的雷达探测盲区大小和影响角度，反推风力发电机与雷达站之间的距离和风力发电机与航道之间的距离，确定风力发电机和航道设置的位置是否合理，以及是否需要增设其他形式协同监管设备。
对于已规划或拟建的海上风电场，根据步骤二和步骤三中的计算模型分析风力发电机叶片和风力发电机对雷达遮蔽以及航道盲区的影响，判断是否需要优化风力发电机位置或是否需要补充协同监管设施。
本构建方法通过对雷达和风力发电机参数的分析，建立了海上风力发电机叶片旋转遮挡模型，确定了风力发电机叶片旋转时雷达探测遮蔽区域面积的计算方法；然后根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，对雷达的主要探测特性能及可能产生的影响进行计算和分析，有助于掌握海上风电场对雷达探测信号的影响程度，可为海上风电场选址、布置及优化提供依据，同时也有利于海事监管部门精准识别雷达盲区及重点监测通航水域，可有效减小雷达探测盲区的通航安全隐患，为海上风电场水域的多形式协同安全监管提供参考。
由于风力发电机叶片的尺寸较大而且通常会持续旋转，因此先行判断风力发电机叶片旋转对雷达遮蔽的影响，再进行风力发电机塔筒对雷达遮蔽的影响分析，在确定出风力发电机对航道影响的盲区距离e、风力发电机对雷达探测影响的角度α以及雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 458643DEST_PATH_IMAGE004
后，进而根据雷达实际扫掠情况，优化风力发电机布置或提出雷达探测盲区水域的监管设施补充方案。
进一步的，所述雷达的参数包括工作波长、脉冲重复频率、脉冲宽度、波束宽度和接收机通频带，并计算得出距离分辨力
Figure 985440DEST_PATH_IMAGE007
和方位分辨力
Figure 716635DEST_PATH_IMAGE008
。
进一步的，步骤二中叶片旋转角度
Figure 303475DEST_PATH_IMAGE001
的确定方法如下：
假设单台风力发电机单个叶片遮挡角度为
Figure 261066DEST_PATH_IMAGE009
，当雷达扫过单台风力发电机时叶片会产生的旋转角度
Figure 325974DEST_PATH_IMAGE001
，则会受到遮挡影响区域的角度为
Figure 442835DEST_PATH_IMAGE010
，一定不会受到遮挡影响区域的角度为
Figure 934996DEST_PATH_IMAGE011
，假设
Figure 379884DEST_PATH_IMAGE002
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 254342DEST_PATH_IMAGE012
，
由于
Figure 366655DEST_PATH_IMAGE013
，故
Figure 560876DEST_PATH_IMAGE014
，则
Figure 758639DEST_PATH_IMAGE015
，
其中
Figure 430928DEST_PATH_IMAGE016
为风力发电机叶片旋转角速度，
Figure 663327DEST_PATH_IMAGE017
为雷达扫掠角速度，r为风力发电机叶片长度，d为雷达与风力发电机间距离。
进一步的，在步骤二中风力发电机叶片旋转遮挡雷达电磁波的面积的确定方法如下：
首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 903815DEST_PATH_IMAGE001
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E；
求得单台风力发电机单个叶片遮挡雷达电磁波的面积SABCD为：
Figure 447929DEST_PATH_IMAGE018
，
则风力发电机叶片旋转时遮挡雷达电磁波的面积S表示为：
Figure 64855DEST_PATH_IMAGE019
；
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线
Figure 10814DEST_PATH_IMAGE020
与
Figure 422204DEST_PATH_IMAGE021
构建函数表达式：
Figure 594559DEST_PATH_IMAGE022
，
联立可解得E的位置为：
Figure 608652DEST_PATH_IMAGE023
，
SBCDE可由下式获得：
Figure 815642DEST_PATH_IMAGE024
，
其中，SDEG及SBCG可根据交点E及辅助点G的位置计算可得：
Figure 522567DEST_PATH_IMAGE025
。
进一步的，根据雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，计算雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 916639DEST_PATH_IMAGE003
；
Figure 875368DEST_PATH_IMAGE026
,
S为风力发电机叶片旋转时遮挡雷达电磁波的面积，将获取的w、r、h、n，以及计算得到的
Figure 530340DEST_PATH_IMAGE027
等数据代入SDEG和SBCG的计算公式中，得到SBCDE，进一步计算得到S，从而可以得到雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 549112DEST_PATH_IMAGE003
的数值；通过
Figure 555114DEST_PATH_IMAGE003
的数值进行分析，可以得到风力发电机叶片旋转对雷达电磁波的遮蔽影响，初步判断风力发电机对雷达探测的影响程度。
进一步的，分别构建d1（雷达与风力发电机间的距离）与
Figure 51955DEST_PATH_IMAGE001
、d1与
Figure 967958DEST_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 751106DEST_PATH_IMAGE028
脉冲宽度就是射频脉冲振荡持续的时间，一般用
Figure 119771DEST_PATH_IMAGE028
表示。一般选择0.05～2μs之内。雷达通常可按照近、中、远距离档以改换发射电磁波的脉冲宽度。
（4）波束宽度
在雷达波束水平面内的波束宽度成为水平波束宽度，用
Figure 279356DEST_PATH_IMAGE029
表示，它与雷达天线的结构、尺寸及使用波长有关。根据《某VTS用户指南》可知基站雷达的水平波束宽度≤0.45°，其后计算参照的最大值为0.45°。
（5）接收机通频带
接收机通频带表示接收机能有效放大的信号频率范围，一般指接收机对输入信号的电压放大倍数从中心频率
Figure 315446DEST_PATH_IMAGE030
的最大相对值1下降到0.707时的两个频率的差
Figure 269495DEST_PATH_IMAGE031
。现代船用雷达接收机通频带为3~25MHz。
（6）距离分辨力
雷达的距离分辨力是指雷达分辨同方向两个相邻目标的能力。同方向的两个目标逐渐靠拢时，两个目标的回波也会逐渐接近而至相切，此时两个目标间的实际距离即距离分辨力。当同方位的两个物标逐渐靠拢时，雷达屏上两个物标的回波亮点也将逐渐接近，当两个回波亮点相切时，两物标间的实际距离即为雷达的距离分辨力
Figure 125456DEST_PATH_IMAGE032
。
Figure 229678DEST_PATH_IMAGE033
，
其中，C为电磁波的传播速度，
Figure 713749DEST_PATH_IMAGE028
为发射脉冲宽度，
Figure 245224DEST_PATH_IMAGE034
为接收机通频带，
Figure 447536DEST_PATH_IMAGE035
为光点直径，D为屏幕直径，RD为所用量程距离。
雷达的距离分辨力主要取决于发射脉冲宽度、接收机通频带、量程以及屏幕光点尺寸大小等技术参数，通过查询获取某VTS雷达在不同量程下的距离分辨力如表1所示。
表1 雷达选取不同量程时的距离分辨力
Figure 89869DEST_PATH_IMAGE037
（7）方位分辨力
雷达的方位分辨力是指雷达区分距离相同而方位相邻的两个目标的能力，以能区分开的两个目标间的最小夹角来表示。此夹角越小，表示雷达的方位分辨力越高，雷达图像的清晰度越好，以能分辨的两个物标间的最小方位夹角
Figure 834972DEST_PATH_IMAGE038
来表示，
Figure 396403DEST_PATH_IMAGE039
越小，表示雷达方位分辨力越高。
影响方位分辨力的主要因素是天线水平波束宽度
Figure 226956DEST_PATH_IMAGE040
、光点角尺寸
Figure 532035DEST_PATH_IMAGE041
（光点直径对屏中心的张角）和回波在屏幕扫描线上的位置，具体雷达的方位分辨力可表示为：
Figure 131644DEST_PATH_IMAGE042
，
基站雷达的水平波束宽度
Figure 739343DEST_PATH_IMAGE040
≤0.45°，光点尺寸
Figure 910387DEST_PATH_IMAGE041
约0.06°。
二、风力发电机参数的确定
某风力发电机底座半径为2.5m，轮毂高度90m，叶片长度为63.4m。
三、风力发电机叶片旋转对雷达电磁波的遮蔽模型构建
实际情况中，风力发电机叶片会对雷达电磁波产生遮挡影响，并且由于风力发电机叶片是处于旋转状态，在雷达扫掠单台风力发电机过程中，风力发电机叶片平面内只有风力发电机叶片对雷达探测存在遮挡影响，其他区域无遮挡。
如图1所示，单台风力发电机单个叶片遮挡角度为
Figure 894523DEST_PATH_IMAGE009
，当雷达扫过单台风力发电机时叶片会产生一定的旋转角度
Figure 473272DEST_PATH_IMAGE043
，则会受到遮挡影响区域的角度为
Figure 986293DEST_PATH_IMAGE010
，一定不会受到遮挡影响区域的角度为
Figure 384913DEST_PATH_IMAGE011
。由于单台风力发电机的3个叶片形状大小完全相同，因此首先单独考虑1个风力发电机叶片对雷达的影响情况：
假设
Figure 31795DEST_PATH_IMAGE044
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 340417DEST_PATH_IMAGE045
（1）
由于
Figure 414552DEST_PATH_IMAGE046
（2）
故
Figure 175835DEST_PATH_IMAGE047
（3）
Figure 626408DEST_PATH_IMAGE048
（4）
其中
Figure 789536DEST_PATH_IMAGE049
为风力发电机叶片旋转角速度90（°/s），
Figure 34572DEST_PATH_IMAGE050
为雷达扫掠角速度
Figure 814309DEST_PATH_IMAGE050
=120（°/s），r为风力发电机叶片长度，d为雷达与风力发电机间距离。
进一步分析风力发电机叶片旋转情况下遮挡雷达电磁波的面积，如图2所示，首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 678360DEST_PATH_IMAGE043
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E。
因此，可求得单个风力发电机叶片遮挡的影响区域S ABCD 为：
Figure 86208DEST_PATH_IMAGE051
（5）
由于雷达扫描时风力发电机叶片的旋转角度为
Figure 643091DEST_PATH_IMAGE043
，则风力发电机叶片旋转的遮挡影响区域可表示为：
Figure 503600DEST_PATH_IMAGE052
（6）
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线LA′D与LAB构建函数表达式：
Figure 436921DEST_PATH_IMAGE053
（7）
Figure 309062DEST_PATH_IMAGE054
（8）
将式（6）与（7）联立可解得E的位置为：
Figure 161480DEST_PATH_IMAGE055
根据上图，SBCDE可由下式获得：
Figure 650230DEST_PATH_IMAGE056
（9）
其中，SDEG及SBCG可根据交点E及辅助点G的位置计算可得：
Figure 980717DEST_PATH_IMAGE057
（10）
Figure 972944DEST_PATH_IMAGE058
（11）
将公式（10）和（11）代入公式（9），并将公式（9）代入公式（5），同时将雷达扫掠单台风力发电机的旋转角度
Figure 137209DEST_PATH_IMAGE043
也代入公式（5）可计算雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，同时可得雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域面积百分比
Figure 706731DEST_PATH_IMAGE059
为：
Figure 716275DEST_PATH_IMAGE060
（12）
将获取的w、r、h、n，以及计算得到的
Figure 953221DEST_PATH_IMAGE043
等数据代入SDEG和SBCG的计算公式中，得到SBCDE的数值，进一步计算得到S的值，从而可以得到雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 22809DEST_PATH_IMAGE003
。从而根据计算可以得到雷达与海上风电场某边界点距离d1，并得到
Figure 345206DEST_PATH_IMAGE043
与
Figure 892862DEST_PATH_IMAGE003
。
Figure 125260DEST_PATH_IMAGE043
随d的变化关系见图3，
Figure 224803DEST_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 909862DEST_PATH_IMAGE061
（13）
风力发电机会对雷达的扫描范围产生遮挡阴影区域，在已分析风力发电机叶片产生的遮挡影响后，再分析风力发电机塔筒对雷达的遮挡影响，以风力发电机轮毂直径来计算风力发电机对雷达电磁波的遮挡范围。设风力发电机塔筒半径为
Figure 651422DEST_PATH_IMAGE062
，为风力发电机与航道间的水平距离为d2，则雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 472748DEST_PATH_IMAGE063
（投射在海面上的阴影面积）为：
Figure 149716DEST_PATH_IMAGE064
（14）
Figure 186986DEST_PATH_IMAGE065
（15）
Figure 76444DEST_PATH_IMAGE066
（16）
其中，x1、x2表示计算中间变量，风力发电机对航道影响的盲区（宽度）距离e为：
Figure 142489DEST_PATH_IMAGE067
（17）
风力发电机对雷达探测影响的角度为：
Figure 724780DEST_PATH_IMAGE068
（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 243486DEST_PATH_IMAGE069
根据表2可知，由于某雷达站与海上风电场距离较远，风力发电机对雷达的最大影响距离都超过39km，再将风力发电机与雷达的距离考虑在内，风力发电机对雷达的遮挡影响基本都达到雷达最大探测距离25海里（约46km），但由于风力发电机底座直径较小，对雷达的影响角度在0.02°左右，换算为对航道内产生的遮挡盲区宽度，最大结果为5.8m左右，因此理论计算的阴影盲区较小。但是，实际情况中，应结合海上风电场中各风力发电机布置，对某通航区域的遮挡影响进行叠加分析，进而结合塔筒遮挡区域可接受范围，综合评估海上风电场中对雷达的影响。
实施例二：
为了更好理解不同距离下风力发电机产生盲区的影响，本实施例与实施例一的区别在于将风力发电机设置在距离雷达站8~10km处，其编号分别为5#、6#、7#和8#。
表3 8~10km处风力发电机对某雷达站的影响
Figure 936636DEST_PATH_IMAGE070
根据表3知，由于该雷达站与海上风电场距离较近，风力发电机对雷达的影响较为明显，航道盲区和影响角度均大于实施例一。同理，对于更近距离的其他雷达站，或者距离其他航道更近的风力发电机可类比按照此方法计算分析。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (7)

1.面向海上风电工程的雷达遮蔽区域模型构建方法，包括设置在海上风电场中的风力发电机，以及设置在所述海上风电场附近的航道、雷达站；其特征在于，所述模型构建方法包括如下步骤：
步骤一：获取所述雷达站中雷达的类型及基本参数，获取所述海上风电场中风力发电机的设计参数；
步骤二：评估所述风力发电机叶片旋转对所述雷达电磁波遮挡的影响，首先确定雷达扫掠单台风力发电机时叶片产生的旋转角度
Figure 798245DEST_PATH_IMAGE001
，其次结合雷达旋转角度
Figure 918647DEST_PATH_IMAGE002
确定单台风力发电机单个叶片旋转遮挡雷达电磁波的面积S，并得到雷达扫描单台风力发电机时叶片所在平面内一定不会产生遮挡区域的百分比
Figure 538198DEST_PATH_IMAGE003
；
步骤三：评估步骤二中风力发电机叶片旋转对所述雷达电磁波遮挡的影响后，根据风力发电机与雷达及航道之间的距离，建立雷达扫描风力发电机塔筒时形成盲区的计算模型，包括确定由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离d和雷达扫掠风力发电机时产生的遮挡阴影面积
Figure 402249DEST_PATH_IMAGE004
得到风力发电机对航道影响的盲区距离e和风力发电机对雷达探测影响的角度α；
Figure 482200DEST_PATH_IMAGE006
，
其中，
Figure 898138DEST_PATH_IMAGE007
为风力发电机塔筒的半径，d1为雷达与风力发电机间的水平距离，d2为风力发电机与航道间的水平距离；
步骤四：根据雷达的基本参数，以及可接受的雷达探测盲区大小和影响角度，反推风力发电机与雷达站之间的距离和风力发电机与航道之间的距离，确定风力发电机和航道设置的位置是否合理，以及是否需要增设其他形式协同监管设备。
2.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，所述雷达的参数包括工作波长、脉冲重复频率、脉冲宽度、波束宽度和接收机通频带，并计算得出距离分辨力
Figure DEST_PATH_IMAGE009
和方位分辨力
Figure 712641DEST_PATH_IMAGE011
。
3.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，步骤二中叶片旋转角度
Figure 177121DEST_PATH_IMAGE001
的确定方法如下：
假设单台风力发电机单个叶片遮挡角度为
Figure 314841DEST_PATH_IMAGE012
，当雷达扫过单台风力发电机时叶片会产生的旋转角度
Figure 636101DEST_PATH_IMAGE001
，则会受到遮挡影响区域的角度为
Figure 859272DEST_PATH_IMAGE013
，一定不会受到遮挡影响区域的角度为
Figure 861863DEST_PATH_IMAGE014
，假设
Figure 962412DEST_PATH_IMAGE002
表示雷达扫掠单台风力发电机的旋转角度，雷达扫掠单台风力发电机的时间为t，则
Figure 861098DEST_PATH_IMAGE016
；
由于
Figure 368302DEST_PATH_IMAGE018
，故
Figure 236901DEST_PATH_IMAGE020
，则
Figure 83635DEST_PATH_IMAGE022
；
其中
Figure 763009DEST_PATH_IMAGE023
为风力发电机叶片旋转角速度，
Figure 757510DEST_PATH_IMAGE024
为雷达扫掠角速度，r为风力发电机叶片长度，d为雷达与风力发电机间距离。
4.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，在步骤二中风力发电机叶片旋转遮挡雷达电磁波的面积S的确定方法如下：
首先将风力发电机叶片抽象为可计算模型，将风力发电机叶片的形状近似为一个三角形ABF和一个直角梯形BCDF的结合，则风力发电机叶片的初始位置为ABCD，旋转
Figure 39586DEST_PATH_IMAGE001
角度后位置为A′B′CD，风力发电机叶片的长度r即为AD之间的距离，B点到F点之间的距离设为w，直角梯形BCDF的短边CD的长度设为n，B′点到CD边之间的垂直距离设为h，线段A′D与线段AB的交点为E；
求得单台风力发电机单个叶片遮挡雷达电磁波的面积SABCD为：
Figure 131039DEST_PATH_IMAGE026
，
则风力发电机叶片旋转时遮挡雷达电磁波的面积S表示为：
Figure 168265DEST_PATH_IMAGE028
；
为计算风力发电机叶片旋转时的重叠面积SBCDE，对直线
Figure 587745DEST_PATH_IMAGE030
与LAB构建函数表达式：
Figure 50344DEST_PATH_IMAGE032
，
联立可解得E点的位置为：
Figure 871669DEST_PATH_IMAGE034
，
SBCDE可由下式获得：
Figure 79797DEST_PATH_IMAGE036
，
其中，SDEG及SBCG可根据交点E及辅助点据G的位置计算可得：
Figure 845628DEST_PATH_IMAGE038
。
5.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，根据雷达扫描单台风力发电机单个叶片时产生的遮挡区域面积S，计算雷达扫描单台风力发电机时叶片所在平面一定不会产生遮挡区域百分比
Figure 735086DEST_PATH_IMAGE003
；
Figure 473235DEST_PATH_IMAGE040
。
6.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，根据步骤二中获取的雷达旋转角度
Figure 930892DEST_PATH_IMAGE001
和一定不会产生遮挡区域百分比
Figure 59385DEST_PATH_IMAGE003
,分别构建
Figure 877169DEST_PATH_IMAGE001
、
Figure 204245DEST_PATH_IMAGE003
随雷达与风力发电机间距离d1的变化关系曲线。
7.根据权利要求1所述的面向海上风电工程的雷达遮蔽区域模型构建方法，其特征在于，假设雷达高度为h1，风力发电机轮毂高度为h2，则可推出由于风力发电机塔筒遮挡对雷达探测产生的阴影径向距离
Figure 957437DEST_PATH_IMAGE042
，并可推出雷达扫掠风力发电机塔筒时遮挡的阴影面积
Figure 947128DEST_PATH_IMAGE044
。
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 true CN113536609A (zh) 2021-10-22
CN113536609B 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)

Cited By (1)

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

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140028496A1 (en) * 2012-07-27 2014-01-30 Texas Tech University System Apparatus and method for using radar to evaluate wind flow fields for wind energy applications
CN106121913A (zh) * 2016-06-20 2016-11-16 广东科诺勘测工程有限公司 一种海上风电场的风机排布方法
CN111610494A (zh) * 2020-05-27 2020-09-01 武汉理工大学 一种vts雷达配置信号覆盖优化方法
CN111796244A (zh) * 2020-06-09 2020-10-20 哈尔滨工业大学 一种海上风电场对地波超视距雷达探测效能影响的评估方法
CN112502911A (zh) * 2020-11-30 2021-03-16 东方电气风电有限公司 一种实时预测叶片通过塔筒时扫塔风险的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140028496A1 (en) * 2012-07-27 2014-01-30 Texas Tech University System Apparatus and method for using radar to evaluate wind flow fields for wind energy applications
CN106121913A (zh) * 2016-06-20 2016-11-16 广东科诺勘测工程有限公司 一种海上风电场的风机排布方法
CN111610494A (zh) * 2020-05-27 2020-09-01 武汉理工大学 一种vts雷达配置信号覆盖优化方法
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)》 *
刘克中 等: "海上风电场对航海雷达探测性能影响研究", 《武汉理工大学学报(交通科学与工程版)》 *

Cited By (2)

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

Also Published As

Publication number Publication date
CN113536609B (zh) 2021-11-30

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
CN106019286B (zh) 一种机载气象雷达多扫描多重频全空域气象目标探测方法
Norin et al. Doppler weather radars and wind turbines
JP5238531B2 (ja) レーダ装置、海洋レーダ観測装置およびドップラ周波数データ算出方法
De la Vega et al. Software tool for the analysis of potential impact of wind farms on radiocommunication services
CN113536609B (zh) 面向海上风电工程的雷达遮蔽区域模型构建方法
CN109001846A (zh) 一种复杂地形条件下s波段与x波段雷达组网测雨方法
CN110726980A (zh) 一种机载相控阵天气雷达对地表杂波的分析方法
CN115032605A (zh) 一种基于目标能力画像的雷达目标体系分析方法
Angulo et al. Estimating reflectivity values from wind turbines for analyzing the potential impact on weather radar services
CN111736150A (zh) 一种面向远距离低功耗探鸟雷达的探测方法
CN104199002B (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
CN114738207A (zh) 一种全天候非接触式风电叶片运行状态监测方法
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