CN114063006A - 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法 - Google Patents

一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法 Download PDF

Info

Publication number
CN114063006A
CN114063006A CN202111347812.XA CN202111347812A CN114063006A CN 114063006 A CN114063006 A CN 114063006A CN 202111347812 A CN202111347812 A CN 202111347812A CN 114063006 A CN114063006 A CN 114063006A
Authority
CN
China
Prior art keywords
steering
array
direction vector
angle
equal
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.)
Pending
Application number
CN202111347812.XA
Other languages
English (en)
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.)
PLA Dalian Naval Academy
Original Assignee
PLA Dalian Naval Academy
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 PLA Dalian Naval Academy filed Critical PLA Dalian Naval Academy
Priority to CN202111347812.XA priority Critical patent/CN114063006A/zh
Publication of CN114063006A publication Critical patent/CN114063006A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/80Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
    • G01S3/802Systems for determining direction or deviation from predetermined direction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明提供一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法,属于水声探测与声纳信号处理技术领域。本发明通过建立阵列绕圆弧轨迹转向模型,针对大角度转向和小角度转向、左转向和右转向条件,利用阵列位圆弧所对应的圆心角将转向过程分解为5个转向阶段,给出5个阶段相邻4个时刻的时间点数值,并对每个转向阶段,以第1阵元为基准,给出当前时刻的波束形成方向向量,对于提高转向过程中目标方位估计性能具有较大的应用价值。

Description

一种拖线阵声纳平台转向过程中波束形成方向向量的测算 方法
技术领域
本发明属于水声探测与声纳信号处理技术领域,涉及一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法。
背景技术
拖线阵声纳是探测低频、隐蔽性目标的重要工具,拖线阵声纳在探测过程中,需要使用波束形成方位谱估计技术,其核心是利用了阵列的方向向量。拖线阵声纳的平台在探测过程中,经常需要改变航向,此时柔性拖线阵列也随之转向。受平台转向的影响,阵列中各阵元的位置与理想线列阵阵元位置存在位置偏移,接收的阵列数据与理想线列阵的接收阵列数据存在较大的差异。这就导致使用理想条件下的方向向量,对有阵元误差的实际阵列数据进行目标方位估计,目标定位结果不准确。在平台转向过程中,准确评估阵列的实时方向向量,对于提高目标方位估计性能具有较大的价值。
发明内容
为了解决上述问题,本发明提供了一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法,在平台转向情况下,实时给出拖线阵声纳波束形成方向向量。
本发明的技术方案是:
一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法,步骤如下:
步骤1:建立阵列绕圆弧轨迹转向模型:假设阵列的第一阵元为1,最后一个阵元为N,d为阵元间距,θ为远场平面波入射方向,λ为声波波长,a(θ)为方向向量。
假设平台转向角度为α,转弯半径为r,阵列孔径为L=(N-1)d。定义αr<L为小角度转向,αr≥L为大角度转向。阵列转向过程中,阵列在圆环上的长度所对应的圆心角度数β随时间变化。根据圆心角度数β的取值,将阵列转向分成5个时间阶段,各时间阶段交界点依次定义为时刻t0、t1、t2、t3
步骤2:分别给出小角度和大角度、向左和向右转向情况下,5个时间阶段的方向向量。具体如下:
(一)小角度转向
小角度转向时,位转向圆环上阵列所对应的圆心角β的取值由图1所示,根据β的取值可以将转向过程分解为5个阶段。
情况1:小角度右转向
当平台和阵列小角度右转向时,各阶段临界时刻的阵列阵形示意图如图2所示。由图2可知,t0是基准时刻,对应于第1阵元到达转向圆环的时刻。
Figure BDA0003354848860000021
Figure BDA0003354848860000022
其中v是平台(阵列)速度。
(1)当0≤t≤t0时,阵列还未到达转向的圆环,呈标准的线列阵排列;此时方向向量为:
Figure BDA0003354848860000023
(2)当t0<t≤t1时,阵列的一部分位于转向圆环上,另一部分呈直线与转向圆环相切。令
Figure BDA0003354848860000024
Figure BDA0003354848860000025
表示下取整,k表示阵列在转向圆环上的最后一个阵元数;
Figure BDA0003354848860000026
m表示阵列位于第1和第k阵元之间的阵元数,βm表示第1至第m阵元的阵列所在圆弧对应的圆心角;阵列在圆环上对应的圆心角
Figure BDA0003354848860000027
此时方向向量:
Figure BDA0003354848860000031
(3)当t1<t≤t2时,阵列的第一个阵元出圆环,最后一个阵元未入圆环。令
Figure BDA0003354848860000032
k表示阵列前段出圆环的最后一个阵元数;
Figure BDA0003354848860000033
q表示阵列在转向圆环上的最后一个阵元数;
Figure BDA0003354848860000034
m表示阵列位于转向圆环上的第1个阵元数k+1和第q阵元之间的阵元数,βm表示阵列在转向圆环上的第一个切点至第m阵元的阵列所在圆弧对应的圆心角。此时方向向量:
Figure BDA0003354848860000035
(4)当t2<t≤t3时,第一个阵元位于圆环的切线上,最后一个阵元在转向圆环上。令
Figure BDA0003354848860000036
k表示阵列前段出圆环的最后一个阵元数;
Figure BDA0003354848860000041
m表示阵列位于转向圆环上的第1个阵元数k+1和第N阵元之间的阵元数,βm表示阵列在转向圆环上的切点至第m阵元的阵列所在圆弧对应的圆心角。此时方向向量:
Figure BDA0003354848860000042
(5)当t>t3时,最后一个阵元已出圆环,所有阵元都位于圆环切线上,呈一条均匀线列阵排列;此时方向向量为:
Figure BDA0003354848860000043
情况2:小角度左转向
此时,转向的5个阶段与小角度右转向相同,如图1所示,且各阶段的字母含义与情况1一一对应。当平台和阵列小角度左转向时,各阶段临界时刻的阵列阵形示意图如图3所示。t0是基准时刻,
Figure BDA0003354848860000044
Figure BDA0003354848860000045
其中v是平台(阵列)速度。
在小角度右转向各时间阶段的方向向量结果上,将变量r、β、βm都分别变为原值的相反数(即乘-1),即可得到小角度左转向各时间阶段的方向向量。具体为:
(1)当0≤t≤t0时,阵列还未到达转向的圆环,呈标准的线列阵排列;此时方向向量为:
Figure BDA0003354848860000051
(2)当t0<t≤t1时,阵列的一部分位于转向圆环上,另一部分呈直线与转向圆环相切。令
Figure BDA0003354848860000052
此时方向向量为:
Figure BDA0003354848860000053
(3)当t1<t≤t2时,阵列的第一个阵元出圆环,最后一个阵元未入圆环。
Figure BDA0003354848860000054
此时方向向量为:
Figure BDA0003354848860000055
(4)当t2<t≤t3时,第一个阵元位于圆环的切线上,最后一个阵元在转向圆环上。令
Figure BDA0003354848860000061
此时方向向量为:
Figure BDA0003354848860000062
(5)当t>t3时,最后一个阵元已驶出圆环,所有阵元都位于圆环切线上,呈一条均匀线列阵排列。此时方向向量为:
Figure BDA0003354848860000063
(二)大角度转向
大角度转向时,位转向圆环上阵列所对应的圆心角β的取值由图4所示,根据β的取值可以将转向过程分解为5个阶段。
情况3:大角度右转向
当平台和阵列大角度右转向时,各阶段临界时刻的阵列阵形示意图如图5所示。由图5可知,t0是基准时刻,
Figure BDA0003354848860000064
(1)当0≤t≤t0时,阵列还未到达转向的圆环,呈标准的线列阵排列。此时方向向量为:
Figure BDA0003354848860000065
(2)当t0<t≤t1时,阵列的一部分位于转向圆环上,另一部分呈直线与转向圆环相切。令
Figure BDA0003354848860000071
k表示阵列在转向圆环上的最后一个阵元数;
Figure BDA0003354848860000072
m表示阵列位于第1和第k阵元之间的阵元数,βm表示第1阵元至第m阵元的阵列所在圆弧对应的圆心角;阵列在圆环上对应的圆心角
Figure BDA0003354848860000073
此时方向向量为:
Figure BDA0003354848860000074
(3)当t1<t≤t2时,阵列的第一个阵元和最后一个阵元都在圆环内。令
Figure BDA0003354848860000075
m表示阵列位于第1阵元和第N阵元之间的阵元数,βm表示第1阵元至第m阵元的阵列所在圆弧对应的圆心角。此时方向向量为:
Figure BDA0003354848860000076
(4)当t2<t≤t3时,第一个阵元位于圆环的切线上,最后一个阵元在转向圆环上。令
Figure BDA0003354848860000077
m表示阵列位于转向圆环上的第1个阵元数k+1和第N阵元之间的阵元数,βm表示阵列在转向圆环上的切点至第m阵元的阵列所在圆弧对应的圆心角。此时方向向量为:
Figure BDA0003354848860000081
(5)当t>t3时,最后一个阵元已驶出圆环,所有阵元都位于圆环切线上,呈一条均匀线列阵排列。此时方向向量为:
Figure BDA0003354848860000082
情况4:大角度左转向
此时,转向的5个阶段与大角度右转向相同,如图4所示,且各阶段的字母含义与情况3一一对应。当平台和阵列大角度左转向时,各阶段临界时刻的阵列阵形示意图如图6所示,各临界时刻分别为:t0是基准时刻,
Figure BDA0003354848860000083
Figure BDA0003354848860000084
在大角度右转向各时间阶段的方向向量结果上,将变量r、β、βm都分别变为原值的相反数(即乘-1),即可得到大角度左转向各时间阶段的方向向量。具体为:
(1)当0≤t≤t0时,阵列还未到达转向的圆环,呈标准的线列阵排列。此时方向向量为:
Figure BDA0003354848860000085
(2)当t0<t≤t1时,阵列的一部分位于转向圆环上,另一部分呈直线与转向圆环相切。令
Figure BDA0003354848860000091
此时方向向量为:
Figure BDA0003354848860000092
(3)当t1<t≤t2时,阵列的第一个阵元和最后一个阵元都在圆环内。令
Figure BDA0003354848860000093
此时方向向量为:
Figure BDA0003354848860000094
(4)当t2<t≤t3时,第一个阵元位于圆环的切线上,最后一个阵元在转向圆环上。令
Figure BDA0003354848860000095
此时方向向量为:
Figure BDA0003354848860000096
(5)当t>t3时,最后一个阵元已驶出圆环,所有阵元都位于圆环切线上,呈一条均匀线列阵排列。此时方向向量为:
Figure BDA0003354848860000101
本发明的有益效果:本发明通过建立阵列绕圆弧轨迹转向模型,针对大角度转向和小角度转向、左转向和右转向条件,利用阵列位圆弧所对应的圆心角将转向过程分解为5个转向阶段,给出5个阶段相邻4个时刻的时间点数值,并对每个转向阶段,以第1阵元为基准,给出当前时刻的波束形成方向向量,对于提高转向过程中目标方位估计性能具有较大的应用价值。
附图说明
图1表示阵列小角度转向过程中,阵列位圆环上的长度所对应的圆心角度数β随时间的变化;β的取值将阵列转向分成5个时间阶段,各阶段交界点分别定义为时刻t0、t1、t2、t3,α为平台或阵列转向角度。
图2表示阵列小角度右转向时,在4个临界时刻点对应的阵列阵形和转向圆环的位置关系。其中,O0(0,0)是进入转向圆环的位置,并以此点建立坐标系x0y0,平台需要转向的角度为α。O1(0,0)、O2(0,0)、O3(0,0)分别对应于第1阵元当前的位置所建立的当前时刻坐标系x1y1、x2y2、x3y3。所建立的坐标系x0,x1,x2,x3方向,指向阵列第1阵元前进方向的左侧垂直方向。y0,y1,y2,y3方向指向第1阵元向第2阵元的切向方向。图3、图5、图6中的相应坐标原点和参数也以相同方式设置。
图3表示阵列小角度左转向时,在4个临界时刻点对应的阵列阵形和转向圆环的位置关系。
图4表示阵列大角度转向过程中,阵列在圆环上所对应的圆心角度数β随时间的变化;β的取值将阵列转向分成5个时间阶段,各阶段交界点分别定义为时刻t0、t1、t2、t3;其中βmax是阵列在转向圆弧上所对应的最大圆心角,图5、图6中相同。
图5表示阵列大角度右转向时,在4个临界时刻点对应的阵列阵形和转向圆环的位置关系。
图6表示阵列大角度左转向时,在4个临界时刻点对应的阵列阵形和转向圆环的位置关系。
具体实施方式
以下结合发明内容和附图详细叙述本发明的具体实施例子。
一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法,步骤如下:
步骤1:确定平台转弯方向、转弯角度α、转弯半径r、阵元数N和阵元间隔d。计算L=(N-1)d。若αr<L,定义为小角度转向,αr≥L定义为大角度转向。
步骤2:当平台和阵列属于小角度转向时,令:t0是第1阵元进入转向圆环的基准时刻,
Figure BDA0003354848860000111
当平台和阵列属于大角度转向时,令:t0是基准时刻,
Figure BDA0003354848860000112
步骤3:判断转向方向和转向角度的四种可能情况,分别给出每种情况转向过程中5个阶段的方向向量。
情况1.小角度右转向
(1)当0≤t≤t0时,
Figure BDA0003354848860000113
(2)当t0<t≤t1时,
Figure BDA0003354848860000121
其中,
Figure BDA0003354848860000122
Figure BDA0003354848860000127
表示下取整,
Figure BDA0003354848860000123
(3)当t1<t≤t2时,
Figure BDA0003354848860000124
其中,
Figure BDA0003354848860000125
(4)当t2<t≤t3时,
Figure BDA0003354848860000126
其中,
Figure BDA0003354848860000131
(5)当t>t3时,
Figure BDA0003354848860000132
情况2.小角度左转向
(1)当0≤t≤t0时,
Figure BDA0003354848860000133
(2)当t0<t≤t1时,
Figure BDA0003354848860000134
其中,
Figure BDA0003354848860000135
(3)当t1<t≤t2时,
Figure BDA0003354848860000141
其中,
Figure BDA0003354848860000142
(4)当t2<t≤t3时,
Figure BDA0003354848860000143
其中,
Figure BDA0003354848860000144
(5)当t>t3时,
Figure BDA0003354848860000145
情况3.大角度右转向
(1)当0≤t≤t0时,
Figure BDA0003354848860000146
(2)当t0<t≤t1时,
Figure BDA0003354848860000151
其中,
Figure BDA0003354848860000152
(3)当t1<t≤t2时,
Figure BDA0003354848860000153
其中
Figure BDA0003354848860000154
(4)当t2<t≤t3时,
Figure BDA0003354848860000155
其中,
Figure BDA0003354848860000156
(5)当t>t3时,
Figure BDA0003354848860000157
情况4.大角度左转向
(1)当0≤t≤t0时,
Figure BDA0003354848860000161
(2)当t0<t≤t1时,
Figure BDA0003354848860000162
其中,
Figure BDA0003354848860000163
(3)当t1<t≤t2时,
Figure BDA0003354848860000164
其中,
Figure BDA0003354848860000165
(4)当t2<t≤t3时,
Figure BDA0003354848860000166
其中,
Figure BDA0003354848860000167
(5)当t>t3时,
Figure BDA0003354848860000171

Claims (1)

1.一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法,其特征在于,步骤如下:
步骤1:建立阵列绕圆弧轨迹转向模型:设阵列的第一阵元为1,最后一个阵元为N,d为阵元间距,θ为远场平面波入射方向,λ为声波波长,a(θ)为方向向量;
设拖线阵声纳平台转向角度为α,转弯半径为r,阵列孔径为L=(N-1)d;定义αr<L为小角度转向,αr≥L为大角度转向;
阵列转向过程中,阵列在圆环上的长度所对应的圆心角度数β随时间变化;根据圆心角度数β的取值,将阵列转向分成5个时间阶段,各时间阶段交界点依次定义为时刻t0、t1、t2、t3
步骤2:当拖线阵声纳平台和阵列属于小角度转向时,令t0为第1阵元进入转向圆环的基准时刻,
Figure FDA0003354848850000011
当平台和阵列属于大角度转向时,令t0为第1阵元进入转向圆环的基准时刻,
Figure FDA0003354848850000012
步骤3:判断转向方向和转向角度的四种情况,分别给出每种情况下转向过程中5个阶段的方向向量,具体如下:
情况1.小角度右转向
(1)当0≤t≤t0时,方向向量为:
Figure FDA0003354848850000013
(2)当t0<t≤t1时,方向向量为:
Figure FDA0003354848850000021
其中,
Figure FDA0003354848850000022
Figure FDA0003354848850000023
表示下取整;
Figure FDA0003354848850000024
1≤m≤k;
Figure FDA0003354848850000025
(3)当t1<t≤t2时,方向向量为:
Figure FDA0003354848850000026
其中,
Figure FDA0003354848850000027
k+1≤m≤q;
(4)当t2<t≤t3时,方向向量为:
Figure FDA0003354848850000028
其中,
Figure FDA0003354848850000031
k+1≤m≤N;
(5)当t>t3时,方向向量为:
Figure FDA0003354848850000032
情况2.小角度左转向
(1)当0≤t≤t0时,方向向量为:
Figure FDA0003354848850000033
(2)当t0<t≤t1时,方向向量为:
Figure FDA0003354848850000034
其中,
Figure FDA0003354848850000035
1≤m≤k,
Figure FDA0003354848850000036
(3)当t1<t≤t2时,方向向量为:
Figure FDA0003354848850000041
其中,
Figure FDA0003354848850000042
k+1≤m≤q;
(4)当t2<t≤t3时,方向向量为:
Figure FDA0003354848850000043
其中,
Figure FDA0003354848850000044
k+1≤m≤N;
(5)当t>t3时,方向向量为:
Figure FDA0003354848850000045
情况3.大角度右转向
(1)当0≤t≤t0时,方向向量为:
Figure FDA0003354848850000046
(2)当t0<t≤t1时,方向向量为:
Figure FDA0003354848850000051
其中,
Figure FDA0003354848850000052
1≤m≤k,
Figure FDA0003354848850000053
(3)当t1<t≤t2时,方向向量为:
Figure FDA0003354848850000054
其中
Figure FDA0003354848850000055
1≤m≤N;
(4)当t2<t≤t3时,方向向量为:
Figure FDA0003354848850000056
其中,
Figure FDA0003354848850000057
k+1≤m≤N;
(5)当t>t3时,方向向量为:
Figure FDA0003354848850000058
情况4.大角度左转向
(1)当0≤t≤t0时,方向向量为:
Figure FDA0003354848850000061
(2)当t0<t≤t1时,方向向量为:
Figure FDA0003354848850000062
其中,
Figure FDA0003354848850000063
1≤m≤k,
Figure FDA0003354848850000064
(3)当t1<t≤t2时,方向向量为:
Figure FDA0003354848850000065
其中,
Figure FDA0003354848850000066
1≤m≤N;
(4)当t2<t≤t3时,方向向量为:
Figure FDA0003354848850000067
其中,
Figure FDA0003354848850000068
k+1≤m≤N;
(5)当t>t3时,方向向量为:
Figure FDA0003354848850000071
CN202111347812.XA 2021-11-15 2021-11-15 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法 Pending CN114063006A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111347812.XA CN114063006A (zh) 2021-11-15 2021-11-15 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111347812.XA CN114063006A (zh) 2021-11-15 2021-11-15 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法

Publications (1)

Publication Number Publication Date
CN114063006A true CN114063006A (zh) 2022-02-18

Family

ID=80271995

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111347812.XA Pending CN114063006A (zh) 2021-11-15 2021-11-15 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法

Country Status (1)

Country Link
CN (1) CN114063006A (zh)

Similar Documents

Publication Publication Date Title
CN111537982B (zh) 一种畸变拖曳阵线谱特征增强方法及系统
CN114397643B (zh) 一种基于超短基线水声定位系统的声线修正方法
CN108919199A (zh) 多波束成像声纳稀疏阵列的旁瓣抑制方法、阵列稀疏方法
CN104931919B (zh) 一种二维平面圆环阵列的测向方法
CN104730513A (zh) 一种分级子阵聚焦mvdr波束形成方法
CN112987004A (zh) 一种浅海环境下基于水平阵列的水面水下目标分类方法
CN104977565A (zh) 一种基于分布式多阵列实现的水下目标三维定位方法
CN114063006A (zh) 一种拖线阵声纳平台转向过程中波束形成方向向量的测算方法
CN107202975B (zh) 一种二维矢量阵阵元姿态误差校正方法
CN109541526A (zh) 一种利用矩阵变换的圆环阵方位估计方法
CN109254265A (zh) 一种基于麦克风阵列的鸣笛车辆定位方法
CN112216298A (zh) 双麦克风阵列声源定向方法、装置及设备
CN113848556B (zh) 基于多波束测深声呐波束图像的水深范围快速提取方法
CN115656994A (zh) 双基地有源探测拖曳阵阵形实时校准方法
CN114200432A (zh) 一种拖线阵声纳平台转向过程中阵列接收信号的测算方法
CN113075645B (zh) 一种基于主成分分析-密度聚类的畸变阵形线谱增强方法
CN115469314A (zh) 一种均匀圆环阵稳健水下目标方位跟踪方法及系统
CN111551943B (zh) 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统
CN111722178A (zh) 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法
CN113138365A (zh) 一种基于深度学习的单矢量水听器方位估计方法
CN110824484A (zh) 一种基于恒模算法的阵元位置估计方法
CN113589299B (zh) 一种基于最优化模型的拖曳线列阵阵形估计方法
CN115015831B (zh) 一种脉冲噪声和非均匀噪声联合影响下的大规模阵列目标方位估计方法
CN113009419B (zh) 一种基于频域互相关匹配的目标深度估计方法
CN103792528B (zh) 一种基于对角减载的水声阵列Bartlett波束形成的方法

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