• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Improvement of circuit oscillation generated by underwater high voltage pulse discharges based on pulse power thyristor

    2023-03-15 00:54:32YingboYU于營(yíng)波andZhongjianKANG康忠健
    Plasma Science and Technology 2023年3期

    Yingbo YU (于營(yíng)波) and Zhongjian KANG (康忠健)

    Institute of New Energy, China University of Petroleum Huadong-Qingdao Campus, Qingdao 266580,People’s Republic of China

    Abstract High voltage fracturing technology was widely used in the field of reservoir reconstruction due to its advantages of being clean, pollution-free, and high-efficiency.However, high-frequency circuit oscillation occurs during the underwater high voltage pulse discharge process,which brings security risks to the stability of the pulse fracturing system.In order to solve this problem,an underwater pulse power discharge system was established, the circuit oscillation generation conditions were analyzed and the circuit oscillation suppression method was proposed.Firstly, the system structure was introduced and the charging model of the energy storage capacitor was established by the state space average method.Next, the electrode high-voltage breakdown model was established through COMSOL software,the electrode breakdown process was analyzed according to the electron density distribution image, and the plasma channel impedance was estimated based on the conductivity simulation results.Then the underwater pulse power discharge process and the circuit oscillation generation condition were analyzed, and the circuit oscillation suppression strategy of using the thyristor to replace the gas spark switch was proposed.Finally, laboratory experiments were carried out to verify the precision of the theoretical model and the suppression effect of circuit oscillation.The experimental results show that the voltage variation of the energy storage capacitor, the impedance change of the pulse power discharge process,and the equivalent circuit in each discharge stage were consistent with the theoretical model.The proposed oscillation suppression strategy cannot only prevent the damage caused by circuit oscillation but also reduce the damping oscillation time by 77.1%,which can greatly improve the stability of the system.This research has potential application value in the field of underwater pulse power discharge for reservoir reconstruction.

    Keywords: underwater high voltage pulse discharge, circuit oscillation suppression, state space average method, pulse power thyristor

    1.Introduction

    Reservoir reconstruction plays an important role in enhancing oil and gas recovery, and the reservoir improvement techniques consist mainly of multistage hydrofracturing [1-3], high energy gas fracturing [4-6], the ultrasonic method [7, 8], and electrical explosion method[9-11],etc.High voltage pulse fracturing plays a more and more important role in the process of reservoir reconstruction due to its advantages of high efficiency, being environmentally friendly, and easy operation.Therefore, it was widely used in the applications of oil plug removal [12], coal fracturing[13],permeability enhancement[14,15],and shale gas stimulation [16], and it would have a better effect if combined with frequency resonance technology [17].

    However, there generates high-frequency oscillation during the high voltage discharge process, which brings security risks to the stability of the system.The high-power oscillation current reduces the service life of the energy storage capacitor, and the peak oscillation voltage also puts forward higher requirements for the withstand voltage of the energy storage circuit,but we often neglect the importance of oscillation suppression in improving the characteristics of pulse discharge.

    For example,the pulse power switch is very important in pulse fracturing technology, and the characteristics and performance of the pulse power switch directly affect the efficiency of pulse power discharge and the reservoir fracturing effect.The gas spark switch was the most widely used [18]due to its simple structure and high-power breakdown capability, though it has the disadvantages of low repetition rate and serious switch corrosion.Scholars solve its weakness by improving the recovery rate of hold-off voltage through adding airflow between electrodes [19], decreasing electrode space [20], injecting different gases between electrodes[21, 22] and changing switch geometry structure [23].These improvements increase switching power and repetition rate,which is of positive significance for repeated pulse fracturing,but the problem of electrode corrosion and high-frequency circuit oscillation has not been solved.Electrode corrosion leads to inconsistent breakdown voltage[24],which makes it difficult to accurately control the discharge time interval,resulting in a poor reservoir resonance effect.The peak high voltage and the pulse current generated by high-frequency circuit oscillation can cause the system over-voltage and overtemperature, and the lifetime of the high voltage film capacitor decreases sharply under such circumstances [25], leading to irreversible damage to the energy storage capacitor and peripheral circuits.

    In order to suppress the potential threat caused by highfrequency circuit oscillation,a pulsed power discharge system based on a gas spark switch was established in this work.The equivalent charge-discharge circuit model was built, the impedance change tendency of the discharge electrode during the breakdown process was simulated by using COMSOL software, and the current-voltage characteristics of underwater high voltage pulse discharges were analyzed.The generation conditions of electrical pulse discharge oscillation were analyzed through the RLC equivalent circuit, and the strategy of using pulse power thyristors to suppress reverse current was proposed.

    2.Simulation model analysis and modeling of energy storage system

    In order to analyze the high voltage pulse discharge oscillation process, the equivalent circuit of a pulse power system based on a gas spark switch was studied,and the time-varying energy storage circuit was simplified based on the state space averaging method.

    2.1.Equivalent circuits of each part of the system

    The high-voltage pulse fracturing system includes a control module, rectifier inverter module, high-voltage armored cable,and energy transducer.The control module and rectifier inverter module were arranged on the ground, and the transducer was placed in the underground mine.The high-voltage cable was connected with the transducer,and the cable length was up to 3 km, enabling the transducer to perform highvoltage pulse fracturing work at any perforation position in the mine within 3 km depth.Each module of the pulse fracturing system and its equivalent circuit are shown in figure 1.In figure 1, the three-phase AC voltage source is the power input of the system, and the rectifier inverter module receives control signals from the control system.The length of a high-voltage armored cable is more than 3 km,which can be equivalent to a series circuit of resistance and inductance,where the resistance and inductance are 93 Ω and 250 μH,respectively.The transducer is mainly composed of an energy storage structure,pulse power switch,and discharge structure,and it works in an underground mine full of water.The energy storage structure was a fourfold voltage rectifier circuit consisting of four 10 μF high-voltage polymer film capacitors and four high-voltage silicon stacks.The pulse power switch was two-electrode self-breakdown gas spark switch, and the discharge structure was a threaded cylindrical electrode,whose material was copper tungsten alloy resistant to high temperature and oxidation corrosion.The shell was made of high tensile strength metal with good sealing characteristics,which can isolate the energy storage circuit and power switch in the shell to avoid external water erosion and discharge pulse impact.Besides, the shell was connected to the ground and slotted at the tail to form an energy transmission channel.

    2.2.Equivalent model of nonlinear time-varying circuit

    The voltage on the energy storage capacitor in the voltage doubling rectifier circuit presents a nonlinear time-varying relationship, therefore, it is necessary to model the timevarying circuit for better analysis.The input voltage is a threephase AC voltage source, which can be expressed as

    where ua, ub, and ucare the three items of three-phase voltages,and the controller adjusts the rectifier voltage amplitude by controlling the conduction angle α of the rectifier bridge,and the relationship between the control angle α and the output rectifier voltage Ureccan be expressed as

    Figure 1.Equivalent circuit of high-voltage pulse fracturing system.

    Figure 2.System equivalent circuit simplification.(a) The complex AC time-varying circuit and (b) the simplified circuit.

    The inverter module converts the rectifier voltage to AC square wave voltage with a duty cycle of 50%,then the boost transformer enlarges the AC voltage by nine times.Therefore,at the input of the high-voltage cable, the relationship between the input voltage Uacand rectifier control angle α was obtained

    The control module and the rectifier inverter module were equivalent to the controlled AC source, and the system equivalent circuit is shown in figure 2(a).The circuit in figure 2(a)was still a complex AC time-varying circuit,and it was necessary to be further equivalent to the circuit displayed in figure 2(b)for better analysis of the pulse power discharge process.

    In figure 2, figure 2(a) is a voltage doubling rectifier circuit, which belongs to a complex AC time-varying circuit,and it is necessary to convert figure 2(a) to the DC circuit shown in figure 2(b) for better analyzing the circuit characteristics.

    The state-space averaging method has been widely applied for nonlinear time-varying circuits analysis [26, 27]due to its ability to convert time-varying circuits into linearity circuits [28], so this method was adopted for figure 2(a) circuit simplification.The AC power supply in figure 2(a) was divided into positive and negative stages according to voltage direction, and the simplified circuit for each stage was illustrated in figure 3.

    In figure 3,upand unare the amplitudes of AC source in positive and negative voltage directions, Ipand Inare the positive and negative flowing currents, respectively.In addition, the RL series circuit was used to represent the armored cable impedance.In the circuit corresponding to the positive voltage stage, the loop voltage relations of meshes 1 and 2 can be listed as

    Figure 3.Positive voltage stage (a) and negative voltage stage (b)simplified circuit.

    where equations (4) and (5) are the voltage relationship of meshes 1 and 2 in figure 3, respectively.Besides, uc1-uc4denote the voltages of capacitors C1-C4, respectively.As an energy storage element,the capacitor C1 keeps absorbing the energy of the power supply upfor the first half cycle.In the meanwhile, the potential on C3 was moved to C2 to realize potential balance.The voltage on inductor L would increases linearly in a very short time according to the theory of the small signal model [29], so, if the rising voltage on L was noted as vg(t), the following expression can be obtained

    where vg(t)can be equivalent to the voltage increase value on the inductance at every moment, which was far less than the value of Ripin equation (4), so the value of L(dip/dt) can be neglected, and equation (4) can be overwritten as

    Therefore, the voltage variation trend on capacitor C1 can be expressed as

    The same analysis can also be applied to the circuit corresponding to the negative voltage stage, and the loop voltage relations of meshes 3 and 4 can be listed as

    where equations (9) and (10) are the voltage relationship of meshes 3 and 4 in figure 3, respectively.Different from the previous stage, capacitor C1 plays the role of the voltage source in this stage, and the capacitor C2 was charged by unand C1 for another half cycle.During this period, the potential on C4 was moved to C3 to realize potential balance.Noting the voltage on C1 as u0, the voltage on C1 and C2 as u1and the relationship between u0and u1satisfies

    where u0can be replaced with the following expression(referring to equation (8))

    Since there exists the relationship up=un(referring to equation (3)), the voltage on capacitor C1 was noted asuc1′,and its voltage expression in a single cycle can be expressed as

    where T and τ are the single period time and time constant,respectively, and the full cycle voltage change tendency on C1 can be estimated by using the state space average method according to circuit time constant τ.The value of T in this design is 0.5 ms, and substituting the value into equations (11)-(13), then we obtain the time constant value τ=0.0233.As is known to all,the value of the time constant τ is equivalent to the value of Ceqmultiplied by Reqin figure 2(b), and the Ceqcan be calculated according to the principle of energy conservation, the storage energy and discharge energy satisfy

    where the storage energy value of the voltage doubling circuit is on the left of the equation and the discharge energy is on the right,besides,the value of C is 10 μf,and there exists the relationship uc4=uc3=uc2=2uc1.Finally, we obtain the circuit parameters in figure 2(b)

    Furthermore, the theoretical expressions of discharge voltage and energy storage voltage are obtained

    where udand ucrepresent the discharge voltage and energy storage voltage, respectively, and ufis the reverse voltage after the discharge process.

    3.Analysis of electrode breakdown process

    The impedance changes dramatically when the electrode gap was broken down by high voltage, and it is difficult to observe the plasma channel formation process as well as impedance changing tendency since the electrode is wrapped by an insulation medium.The structure of the gas spark switch and discharge electrode was introduced in this section,and the high voltage discharge simulation was carried out by COMSOL software, then the plasma channel formation process and impedance changing tendency were analyzed based on simulation modeling results.

    3.1.Gas spark switch and discharge electrode structure

    The main components of the gas spark switch and discharge electrode are made of copper tungsten alloy resistant to high temperature and oxidation corrosion, and their structures are depicted in figure 4.

    Figure 4.The structure of gas spark switch (a) and discharge electrode (b).

    In figure 4, there is a ceramic wall wrapped around the gas switch, which plays a key role in heat isolation and high voltage insulation.Besides, the electrode plate in the gas switch is made of copper, which has good conductivity and heat dissipation ability, and there are small threaded holes distributed around the middle big threaded hole in the plate,the middle threaded hole is used to fix the electrode and the around smaller holes are designed to fix ceramic tubes for electrode plate supporting function.The electrodes in gas switch as well as in discharge structure are made of copper tungsten alloy with high strength and oxidation corrosion resistance.In addition, the electrode structure is designed as cylindrical with thread,which is not only convenient to adjust the electrode gap but also conducive to replacement and maintenance.The working environment of the discharge electrode is a water-filled underground mine, and the outer wall of the discharge electrode is a metal shell with high tensile strength and good sealing characteristics,the shell was connected to ground and slotted at the tail to form an energy transmission channel.

    3.2.High voltage breakdown process of electrode

    The main components of gas spark switch and discharge structure are copper tungsten alloy electrodes, and there are some similarities between gas spark discharge and pulsed arc electrohydraulic discharge.When the electrode discharges in gas,the high voltage directly breaks through the gas to form a plasma channel.In the process of pulsed arc electrohydraulic discharge, there are both the collision ionization of charged particles and the expansion and development of hot bubbles,and the electron avalanche develops along the bubble holes[30-32], and finally completes the whole process of pulsed electrohydraulic discharge.

    In this work, COMSOL 6.0 multi-physical field simulation software was used to simulate the electrode breakdown process with argon as the gas between the electrode gaps.In the simulation software, the electrode material was set as copper tungsten alloy, and the electrode diameter and electrode gap were 32 mm and 10 mm, respectively.To simplify the simulation model, the closed space inside the bubble is chosen as the simulation domain,and the hydrate interference factors around the bubbles are removed.The electrode breakdown process can be observed through the development of electron density on the electrode surface.Two-dimensional axisymmetric modeling was adopted, and its simulation circuit and modeling graphics are exhibited in figure 5.

    In figure 5,Vinis the external DC voltage,which was set to 15 000 V, and in the external RC series circuit, R and C were set to 120 Ω and 8.1 μF, respectively.In order to make the initial environment universal,the electron density value in the simulation domain is set to 1013m-3with reference to many studies [17, 33, 34], and the electron density variation during the electrode breakdown process was shown in figure 6.

    In figure 6, figure 6(a) was the original electron density,the space electrons gradually move to the anode under the influence of the applied electric field and finally converge to the anode tip, and the electron migration process can be seen from figures 6(b)and(c),corresponding to the time of 0-1.58 μs.As electrons gather at the anode tip, the electrons were emitted from anode to cathode to form a plasma channel during the period of 1.67-2.18 μs, and the electrode breakdown process can be seen from the development of electron density referring to figures 6(d)-(f).The electron density increases sharply at this stage,and the whole breakdown time was within 0.6 μs.

    Figure 5.The structure of simulation circuit.

    Figure 6.Electron density distribution during electrode breakdown process.

    Figure 7.Conductivity variation under different currents.

    Figure 8.The RLC equivalent discharge circuit.

    Figure 9.Capacitor voltage variation characteristics in different RLC parameters.

    3.3.Analysis of impedance characteristics of plasma channel

    In order to study the impedance variation characteristics in the process of electrode breakdown, the plasma channel current was set based on the variation trend of the actual oscillation current from 10 to 610 A,and the conductivity changes under different currents as displayed in figure 7.

    As shown in figure 7, the conductivity shows an increasing tendency with the increase of current, and the maximum conductivity value increases from 5.65 × 103to 9.45 × 103S m-1.In addition, the plasma channel radius would expand with the increase of conduction current,and the radius of the plasma channel was distributed from 1.5 to 3 mm when the current varies from 10 to -610 A.The equivalent resistance Rplsof plasma channel can be calculated based on the conductivity S using the following equation

    where L represents the gap length of the electrode and r represents the cross-section radius of the plasma channel.The equivalent resistance was distributed in the range of 0.037-0.35 Ω after substituting data into the calculation.

    It can be seen from the simulation results that the high voltage breakdown time of the electrode is very short(within 1 μs), and the equivalent resistance gradually decreases with the increase of the plasma channel current,and the equivalent resistance changes in a certain range to form a relatively stable state.The circuit oscillation conditions can be analyzed in detail through the impedance change trend.

    4.Underwater pulse power discharge process analysis

    4.1.Discharge equivalent circuit analysis

    The discharge circuit of the energy storage capacitor in the process of pulse power discharge can be equivalent to the RLC series circuit, and the schematic illustration is in figure 8.

    In figure 8, R, L, and C represent plasma channel equivalent resistor,circuit inductance,and a storage capacitor,respectively, and the circuit relationship can be described as follows.

    Figure 10.Equivalent circuit and corresponding voltage and current on electrode.

    where UL, UC, and i represent inductance voltage, discharge voltage, and circuit current, respectively.The relationship between voltage variation characteristics of the energy storage capacitor and R, L, C parameters was shown in figure 9.

    It can be seen from figure 9 that, when R < 2(L/C)1/2,the system exhibits underdamped oscillation, and the peak oscillation will cause great damage to the system.When the electrode breaks down in the discharge process, the system meets the damping oscillation condition according to the simulation results, which were also consistent with the relevant research [35, 36].When circuit oscillation occurs, the capacitor voltage and current satisfy

    where A and K are constants,φ and β are angles,the capacitor voltage and the current show oscillation attenuation trend.

    4.2.Circuit equivalence in different discharge stages

    The gas spark switch and the discharge structure meet the circuit oscillation condition during the high voltage breakdown process according to previous study.Therefore, the discharge process was divided into four stages according to the circuit oscillation characteristics, namely, the gas spark switch breakdown stage, the discharge electrode pre-breakdown stage, the electrode discharge stage, and the damping oscillation stage.The equivalent circuit of each stage and corresponding discharge voltage and current on the electrode are exhibited in figure 10.

    In figure 10, the gas spark switch breakdown stage(marked 1) and electrode discharge stage (marked 3) have something in common,the impedance of both stages changes dramatically in a very short time, resulting in a high-quality factor to generate an oscillation peak.In the discharge electrode pre-breakdown stage (marked 2), due to the fact that hydrate filled between two discharge electrodes and it took some time for liquid gasification and ionization to produce an electro-hydraulic effect, during this stage, the discharge electrode presents a resistance character and shows a zero input response in equivalent discharge circuit.In the damping oscillation stage (marked 4), the impedance of plasma channels formed in previous stages (marked 1 and 3) varies with the current in a small range showing a relatively stable characteristic, resulting in an oscillation attenuation trend.

    4.3.Improvement strategy

    The gas spark switch and the discharge structure can generate high-frequency circuit oscillation during the high voltage breakdown process,and the values of peak oscillation voltage and current are much larger than the average discharge voltage and current,which brings great potential safety hazards to the circuit insulation and components voltage endurance.This kind of oscillation harm is particularly obvious to the energy storage capacitor,in the case of high-frequency discharge,the multiple cycle oscillation peak pulse current can make the capacitor heat up, resulting in the service life decrease of the energy storage capacitor due to over-voltage and over-temperature working conditions.In order to improve the discharge process without loss of discharge strength, the most effective method is to suppress the reverse current.

    The thyristor has a high switching power characteristic and the ability to pass through ultra-high surge current,which has been widely used in many high-power applications [37],in addition, its reverse current automatic shutdown characteristic was very suitable for reverse oscillation current suppress in underwater pulse power discharges.If the gas spark switch was replaced by a pulse power thyristor, the oscillation of electrode discharge voltage and current waveform can be well improved, and the waveform theoretical improvement can be displayed in figure 11.

    In figure 11,the original high-frequency oscillation in the high-voltage breakdown process of the gas spark switch would completely disappear since the gas spark switch was replaced by thyristors.In the electrode discharge stage, the voltage energy on the energy storage capacitor was so high that the reverse oscillation voltage cannot make the thyristor current flow backward,and the current flow direction was still dominated by the energy storage voltage.Therefore, in this stage, the reverse oscillation current can be suppressed without turning off the thyristor.

    In the damping oscillation stage, the voltage of the energy storage capacitor decreases,and the reverse oscillation voltage would be higher than that of the energy storage capacitor, so, at this time, the current flow direction was dominated by the oscillation voltage.However, when the current reverse oscillation occurs,the thyristor would be shut down due to the reverse current automatic shutdown characteristic, and the damping oscillation process will be terminated in advance after only half a cycle,which greatly reduces the discharge time.

    Figure 11.Electrode discharge waveform theoretical improvement.(a)The original discharge waveform and(b)the discharge waveform after oscillation suppression.

    5.Underwater pulse power discharge process analysis

    5.1.Introduction to experimental device

    In order to study the underwater pulse power discharge process of a gas spark switch as the pulse power switch and verify the thyristor suppression effect on pulse power discharge oscillation,the experimental platform was designed as shown in figure 12.

    In figure 12, the control unit (marked 1) can adjust the voltage of the energy storage capacitor (marked 2) by controlling the rectifier angle, and the storage capacitor and discharge electrode (marked 5) were connected at both ends of the two-electrode gas spark switch(marked 3)to form a series connection.In addition, the electrode was immersed in the cuboid water cylinder and fixed on the metal frame connected to the ground, and the pulse power thyristor (marked 4) was the comparison test component.

    Figure 12.Experimental platform for underwater pulse power discharge.

    5.2.Energy storage experiment

    The charging voltage is adjusted by the control unit,when the voltage of the energy storage capacitor reaches the gas spark switch self-breakdown voltage and the pulse power discharge condition is reached.However, there will be reverse voltage ufafter the pulse power discharge, and the actual charging waveform was compared with the theoretical charging waveform as depicted in figure 13.

    It can be seen from figure 13 that the theoretical modeling curve almost overlaps with the simulation curve, and it also shows a good consistency compared with the actual charging curve.The error between theoretical calculation and actual charging waveform was calculated, and the results were that the average error, absolute error and mean square deviation were 2.69%, 4.02%, and 7.27%, respectively.The error is small, which indicates that the theoretical modeling can accurately reflect the circuit characteristics of the charging circuit.

    Figure 13.Comparison between actual charging curve and theoretical modeling.

    5.3.Underwater pulse power discharge experiment

    5.3.1.Typical waveform of pulse power discharge.When the voltage of the energy storage capacitor reaches the discharge voltage, the electrode produces discharge action after the gas spark switch breakdown by high voltage, the voltage-current characteristic curves on the discharge electrode and energy storage capacitor were recorded by oscilloscope DL350,three sets of typical curves were selected,as exhibited in figure 14.In figure 14, it can be seen from figures 14(a1) and (b1)that the gas spark switch breakdown stage, discharge electrode pre-breakdown stage, electrode discharge stage,and damping oscillation stage were distinguished clearly, the oscillation peak voltage was more than two times higher than the discharge voltage,and the damping oscillation attenuation time is about 175 ms for 5-6 cycles.It seems that the gas spark switch breakdown stage was missing in figures 14(a2)and(b2),this is because the digital oscilloscope sampling time is 1 μs (DL350 sample rate is 1 M), which is larger than the electrode breakdown time (0.6 μs in simulation), and based on this reason, the oscillation process of electrode discharge stage was often missing in the recording curve, which was also suitable for gas spark switch breakdown stage.Therefore,the record maximum oscillation peak value was often smaller than the actual value.When the time of discharge electrode pre-breakdown stage is very short, there will be two highfrequency oscillation waveforms suspected to be connected together, as shown in figures 14(a3) and (b3).

    5.3.2.Parameter identification of damping oscillation.In the discharge electrode pre-breakdown stage, the circuit is equivalent to the RC series circuit shown in figure 10, and the discharge electrode voltage satisfies zero input response,the relationship can be expressed as

    where Udrepresents the electrode voltage, U represents the voltage before attenuation, R and C represent the values of equivalent resistance of the electrode and energy storage capacitor, respectively.The equivalent resistance can be calculated by substituting a two-point value in the prebreakdown stage, and after averaging multiple calculations,the resistance value satisfies R=220 Ω.

    In the damping oscillation stage, the voltage and current relation satisfies equation (21).For the oscillation current,when the condition satisfiesωt+β=and there will besin (ωt+β)= 1,these points are the vertexes of the current damping oscillation curve, which decay exponentially, and can be expressed as

    The time difference between two adjacent vertices is exactly a period T, when the value of current damping oscillation between two adjacent vertices and the time between two points are known, the values of R and L can be obtained by the following formula

    The average period of damped oscillation is 35 μs measured in figure 14, and the values of R and L were calculated as R=0.1515 Ω,L=3.787 μH.The peak voltage of circuit oscillation is related to the quality factor Q, the higher the quality value Q, the higher the peak oscillation voltage, and the expression of quality factor Q satisfies

    In the gas spark switch breakdown stage and electrode discharge stage, because of the large breakdown current, the equivalent resistor satisfies R < 0.1515.Assuming that the inductance value remains unchanged,and the quality factor Q was calculated as expressed Q > 4.5, therefore, high peak circuit oscillation will be formed in the two stages.

    The electrode discharge voltage waveform in the same test was selected, and the calculation parameter was substituted into the equivalent circuit, the comparison between the actual waveform and equivalent circuit waveform is illustrated in figure 15.

    Figure 14.Typical voltage-current characteristic curves on discharge electrode(a1),(a2)and(a3)and energy storage capacitor(b1),(b2)and(b3).(a1) and (b1) is a set of the typical voltage-current characteristic curve on the discharge electrode and energy storage capacitor,respectively, (a2) and (b2), (a3) and (c3) are another two sets.

    In figure 15, the actual discharge waveform in the prebreakdown stage and damping oscillation stage is consistent with the equivalent circuit model of the theoretical analysis,besides, the calculated plasma channel resistance and inductance values are within the theoretical calculation range of the simulation model, indicating that the theoretical modeling build by COMSOL software can accurately reflect the impedance variation characteristics during the discharge process.

    5.4.Improvement strategy based on thyristor

    In order to verify the oscillation suppression effect of using thyristors as a pulse power switch, the pulse power thyristor 5STP6500 was selected to replace the gas spark switch for the comparison test.The pulse power discharge comparison waveform in a single cycle was shown in figure 16.

    In figure 16,there are oblivious peak oscillation voltages on the inserted figure in the original discharge waveform,while such oscillation is not found in the insert figure in the waveform after oscillation suppression,indicating that reverse turn-off characteristics of the thyristor can effectively suppress the peak voltage caused by high-frequency oscillation.In addition, the voltage and current of the energy storage capacitor will experience the damping oscillation attenuation process in the damping oscillation stage on the original discharge waveform, which lasts for more than 175 ms, and during this period, the energy storage capacitor will experience 5 to 6 times charging and discharging process.The life of a high-voltage film capacitor is often measured by the number of charging and discharging cycles, therefore, multiple charging and discharging in a short time will make the capacitor life decay faster [25].In the waveform after oscillation suppression, the current damping oscillation process only maintains half a cycle, and the oscillation time was reduced to 40 ms, which is 77.1% lower than that before suppression.In addition, the energy storage capacitor only charges and discharges one cycle in a single discharge process,which improves the life of the energy storage capacitor.

    Through comparative tests, it is concluded that using thyristors as a pulse power switch to replace the traditional gas spark switch cannot only prevent the damage of continuous high-frequency oscillation current to the energy storage capacitor, but also greatly reduce the voltage requirements of components, and reduce the discharge time,which can significantly improve the oscillation characteristics of underwater pulse discharge circuit.

    Figure 15.Comparison between actual waveform and equivalent circuit waveform.

    Figure 16.Comparison between original discharge waveform(a)and the waveform after oscillation suppression(b).The insert picture was repeated discharge waveform on storage capacitor.

    6.Conclusions

    An underwater pulsed power discharge system was established,the conditions of circuit oscillation were analyzed based on the theoretical model,and the circuit oscillation suppression method by using pulse power thyristors to replace the gas spark switch was proposed.The method was verified by experiment, and the results show that the high-frequency oscillation was significantly suppressed, and the damping oscillation time was reduced by 71%.Thyristor has the advantages of low cost, small volume,large conduction power, and high surge current, which was widely used in high voltage rectification and high voltage direct current transmission (HVDC), however, it is rarely used in the field of underwater pulse power discharge.As a semi-controlled device, the thyristor has the characteristic of automatic turn-off under reverse current, therefore, it can suppress the circuit oscillation process by blocking the reverse current and reduce the oscillation time by cutting off the circuit,which was suitable for the application of underwater pulsed power discharge.The series diode in the discharge circuit can also suppress the reverse current, which seems to be a simple method.In fact, the ultra-high voltage high-current diode is expensive and has a large volume and uncontrollable factors, it would not turn off when experiencing the reverse current,and the circuit would still experience a long period of oscillation.Compared with fully controlled devices such as IGBT and MOSFET,thyristor has the advantages of higher withstand voltage,higher switching power,and a stronger anti-interference ability for current signal control, and the automatic turn-off characteristics under reverse current make it easier to control the discharge process by focusing only on the driving signal.This study has a potential application value in the field of stability improvement of pulse power discharge equipment.

    Acknowledgments

    This work was financially supported by the National Science and Technology Major Project (No.2016ZX05034004).

    大陆偷拍与自拍| 免费观看人在逋| 亚洲国产欧美日韩在线播放| 日韩免费高清中文字幕av| 777米奇影视久久| 久99久视频精品免费| 1024视频免费在线观看| 黄色视频,在线免费观看| 人妻一区二区av| 免费黄频网站在线观看国产| 国产精品久久久久久精品古装| 水蜜桃什么品种好| 悠悠久久av| 亚洲精品在线美女| 午夜两性在线视频| 国产精品98久久久久久宅男小说| 亚洲自偷自拍图片 自拍| 中文字幕精品免费在线观看视频| 国产成人精品久久二区二区免费| 成人18禁高潮啪啪吃奶动态图| 久久ye,这里只有精品| 老司机午夜福利在线观看视频| 热99re8久久精品国产| 免费不卡黄色视频| 视频区图区小说| 在线观看免费午夜福利视频| 无人区码免费观看不卡| 99久久精品国产亚洲精品| e午夜精品久久久久久久| 国产精品秋霞免费鲁丝片| 免费不卡黄色视频| 一级作爱视频免费观看| 嫁个100分男人电影在线观看| 中文字幕人妻丝袜制服| 久久中文字幕人妻熟女| 国产三级黄色录像| 亚洲久久久国产精品| 日韩欧美一区二区三区在线观看 | 一本大道久久a久久精品| 91在线观看av| 久久精品亚洲精品国产色婷小说| 亚洲色图综合在线观看| 国产在线精品亚洲第一网站| 多毛熟女@视频| 亚洲成a人片在线一区二区| 欧美日韩福利视频一区二区| cao死你这个sao货| 岛国毛片在线播放| 国产亚洲欧美在线一区二区| 亚洲片人在线观看| 日韩欧美三级三区| 国产一区二区激情短视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人精品无人区| 成在线人永久免费视频| av有码第一页| 国内毛片毛片毛片毛片毛片| 大陆偷拍与自拍| 欧美日韩成人在线一区二区| 亚洲免费av在线视频| 精品第一国产精品| 激情在线观看视频在线高清 | 麻豆av在线久日| 日韩大码丰满熟妇| 成在线人永久免费视频| 亚洲少妇的诱惑av| 十八禁高潮呻吟视频| 久久人人97超碰香蕉20202| 国产高清视频在线播放一区| 久热爱精品视频在线9| 美女国产高潮福利片在线看| 亚洲精品国产区一区二| 亚洲成人免费电影在线观看| 窝窝影院91人妻| 欧美乱码精品一区二区三区| 成人精品一区二区免费| 美女福利国产在线| 窝窝影院91人妻| 国产成+人综合+亚洲专区| 中文字幕精品免费在线观看视频| 曰老女人黄片| 麻豆成人av在线观看| 精品国产超薄肉色丝袜足j| 亚洲专区中文字幕在线| 久久精品国产a三级三级三级| 无遮挡黄片免费观看| 亚洲第一av免费看| 久久精品aⅴ一区二区三区四区| 一本一本久久a久久精品综合妖精| 亚洲情色 制服丝袜| 国产激情欧美一区二区| 丝袜美足系列| 免费观看精品视频网站| 久久精品国产综合久久久| 9191精品国产免费久久| 一级毛片精品| 看免费av毛片| 成人精品一区二区免费| 成人永久免费在线观看视频| 欧美最黄视频在线播放免费 | 午夜亚洲福利在线播放| 免费观看人在逋| 香蕉丝袜av| 欧美日韩乱码在线| 久久人妻福利社区极品人妻图片| 久久人妻熟女aⅴ| 香蕉丝袜av| 免费高清在线观看日韩| 视频区欧美日本亚洲| 国产男女超爽视频在线观看| 免费女性裸体啪啪无遮挡网站| 男女之事视频高清在线观看| 天堂中文最新版在线下载| 国产亚洲精品久久久久5区| 男女免费视频国产| 天堂中文最新版在线下载| 日韩人妻精品一区2区三区| 国产成人影院久久av| 一区二区三区精品91| 亚洲人成电影观看| 久久精品亚洲av国产电影网| 大香蕉久久成人网| 一级作爱视频免费观看| 99香蕉大伊视频| 后天国语完整版免费观看| 亚洲精品自拍成人| 中亚洲国语对白在线视频| 最近最新中文字幕大全免费视频| 久久人妻福利社区极品人妻图片| 91精品三级在线观看| 深夜精品福利| 国产欧美日韩精品亚洲av| 人人妻人人爽人人添夜夜欢视频| 另类亚洲欧美激情| 黄色女人牲交| 国产91精品成人一区二区三区| 人妻一区二区av| 99国产精品免费福利视频| 久久天堂一区二区三区四区| 国产精品二区激情视频| 亚洲精品美女久久av网站| 亚洲av成人av| 久久国产精品大桥未久av| 18禁美女被吸乳视频| 精品福利永久在线观看| 久久精品91无色码中文字幕| 久久久国产成人免费| 成人三级做爰电影| 久久久久视频综合| 老司机福利观看| av天堂久久9| 99久久99久久久精品蜜桃| 欧美国产精品va在线观看不卡| 国产精品免费视频内射| 精品国产超薄肉色丝袜足j| 日韩制服丝袜自拍偷拍| 美女午夜性视频免费| 亚洲视频免费观看视频| 99精国产麻豆久久婷婷| 丝袜美腿诱惑在线| 性色av乱码一区二区三区2| 麻豆av在线久日| 亚洲人成电影免费在线| 国产区一区二久久| 校园春色视频在线观看| 国产aⅴ精品一区二区三区波| 久久热在线av| 国精品久久久久久国模美| 真人做人爱边吃奶动态| 久久精品亚洲精品国产色婷小说| 欧美 日韩 精品 国产| 18禁国产床啪视频网站| 日韩视频一区二区在线观看| a在线观看视频网站| e午夜精品久久久久久久| 国产高清videossex| 亚洲中文日韩欧美视频| 国产精品99久久99久久久不卡| 亚洲专区中文字幕在线| 黄色女人牲交| 欧美日本中文国产一区发布| 久久久久久久精品吃奶| 视频在线观看一区二区三区| 三级毛片av免费| 多毛熟女@视频| 91老司机精品| 国产免费男女视频| 国产无遮挡羞羞视频在线观看| 色在线成人网| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯 | 久久人人爽av亚洲精品天堂| 人成视频在线观看免费观看| 亚洲精华国产精华精| 两性夫妻黄色片| 校园春色视频在线观看| 精品第一国产精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品成人av观看孕妇| 飞空精品影院首页| 国产野战对白在线观看| 久久久精品免费免费高清| 一二三四社区在线视频社区8| 午夜福利在线观看吧| 老司机午夜十八禁免费视频| 婷婷丁香在线五月| 国精品久久久久久国模美| 欧美不卡视频在线免费观看 | 高清毛片免费观看视频网站 | 久久天堂一区二区三区四区| 天堂√8在线中文| 久久中文字幕人妻熟女| 欧美乱妇无乱码| 成人黄色视频免费在线看| 最新在线观看一区二区三区| 国产高清国产精品国产三级| 欧美国产精品va在线观看不卡| 午夜影院日韩av| 狠狠狠狠99中文字幕| 欧美日韩福利视频一区二区| 91麻豆av在线| 人人澡人人妻人| 欧美精品啪啪一区二区三区| 日韩免费av在线播放| 天天影视国产精品| 亚洲熟妇中文字幕五十中出 | 成在线人永久免费视频| 久99久视频精品免费| 亚洲欧美日韩高清在线视频| 国产免费男女视频| 国产亚洲精品第一综合不卡| 国产欧美日韩精品亚洲av| 成年人黄色毛片网站| 脱女人内裤的视频| 在线观看www视频免费| 午夜福利乱码中文字幕| 成在线人永久免费视频| 欧美成狂野欧美在线观看| 久久久国产成人精品二区 | 黄色a级毛片大全视频| 老司机在亚洲福利影院| 久久精品国产a三级三级三级| 十八禁网站免费在线| 夜夜爽天天搞| 很黄的视频免费| 免费久久久久久久精品成人欧美视频| 黄色怎么调成土黄色| 免费日韩欧美在线观看| 欧美 日韩 精品 国产| 在线免费观看的www视频| 国产精品久久视频播放| bbb黄色大片| 久久久久久久久免费视频了| 久久久久国产精品人妻aⅴ院 | 宅男免费午夜| 满18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 国产午夜精品久久久久久| 亚洲精品中文字幕一二三四区| 精品久久久久久电影网| 人妻丰满熟妇av一区二区三区 | 又紧又爽又黄一区二区| 一级a爱片免费观看的视频| 国产免费现黄频在线看| 美女午夜性视频免费| 露出奶头的视频| 欧美日韩福利视频一区二区| 两性夫妻黄色片| 一区二区三区国产精品乱码| 久久人人97超碰香蕉20202| 亚洲,欧美精品.| av网站免费在线观看视频| aaaaa片日本免费| 夜夜躁狠狠躁天天躁| 黄色成人免费大全| 国产日韩欧美亚洲二区| 十八禁网站免费在线| 亚洲国产欧美一区二区综合| av片东京热男人的天堂| 精品一区二区三区av网在线观看| 欧美国产精品va在线观看不卡| 99久久人妻综合| 天堂动漫精品| 国产成人精品久久二区二区91| 操出白浆在线播放| 国产高清videossex| 免费日韩欧美在线观看| 久久精品国产a三级三级三级| 欧美黄色片欧美黄色片| 亚洲av日韩在线播放| 精品久久蜜臀av无| 黄色毛片三级朝国网站| 亚洲五月色婷婷综合| 亚洲欧美色中文字幕在线| av有码第一页| 十八禁人妻一区二区| 婷婷丁香在线五月| 王馨瑶露胸无遮挡在线观看| 亚洲精品在线美女| 飞空精品影院首页| 亚洲欧美色中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美一区视频在线观看| 亚洲精品国产精品久久久不卡| 日韩人妻精品一区2区三区| 成人国语在线视频| 成人精品一区二区免费| 精品国产亚洲在线| 国产精品久久久人人做人人爽| 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 免费观看精品视频网站| 国产亚洲一区二区精品| 丰满的人妻完整版| 亚洲一码二码三码区别大吗| 欧美丝袜亚洲另类 | 俄罗斯特黄特色一大片| 操出白浆在线播放| 怎么达到女性高潮| 一区二区三区激情视频| 精品一品国产午夜福利视频| 免费观看人在逋| 99久久精品国产亚洲精品| 亚洲av片天天在线观看| 成人av一区二区三区在线看| 欧美黄色片欧美黄色片| 男女免费视频国产| 在线观看免费午夜福利视频| 欧美久久黑人一区二区| 国产精品98久久久久久宅男小说| 丁香欧美五月| 99国产综合亚洲精品| 日本vs欧美在线观看视频| 男男h啪啪无遮挡| 欧美日韩亚洲高清精品| 捣出白浆h1v1| tocl精华| 精品少妇久久久久久888优播| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 色精品久久人妻99蜜桃| 伦理电影免费视频| 成年版毛片免费区| 精品一品国产午夜福利视频| cao死你这个sao货| 校园春色视频在线观看| 色婷婷久久久亚洲欧美| tube8黄色片| 狠狠狠狠99中文字幕| 精品午夜福利视频在线观看一区| 国产在视频线精品| 嫩草影视91久久| 大香蕉久久网| 女警被强在线播放| 99国产精品99久久久久| av线在线观看网站| 久久久国产精品麻豆| 日韩一卡2卡3卡4卡2021年| 午夜日韩欧美国产| 99国产极品粉嫩在线观看| 精品国产乱码久久久久久男人| 两人在一起打扑克的视频| 亚洲av片天天在线观看| 国产亚洲精品第一综合不卡| 满18在线观看网站| 黄色怎么调成土黄色| 久久国产精品影院| 在线永久观看黄色视频| 国产黄色免费在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美网| 在线观看www视频免费| 一二三四社区在线视频社区8| 两性夫妻黄色片| 亚洲avbb在线观看| 免费观看精品视频网站| 国产亚洲一区二区精品| 乱人伦中国视频| 国产有黄有色有爽视频| 欧美成人免费av一区二区三区 | 国产日韩欧美亚洲二区| 狠狠婷婷综合久久久久久88av| 午夜成年电影在线免费观看| 建设人人有责人人尽责人人享有的| 捣出白浆h1v1| 久久久久久亚洲精品国产蜜桃av| 757午夜福利合集在线观看| 久久精品aⅴ一区二区三区四区| 成人手机av| 色综合婷婷激情| 亚洲国产看品久久| 99riav亚洲国产免费| 免费一级毛片在线播放高清视频 | 国产成人免费观看mmmm| 国产精品国产高清国产av | 欧美成人免费av一区二区三区 | 午夜福利乱码中文字幕| 久久性视频一级片| 精品高清国产在线一区| 色在线成人网| 亚洲精品自拍成人| 国产1区2区3区精品| 日韩成人在线观看一区二区三区| 国产成人精品久久二区二区91| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三区在线| 叶爱在线成人免费视频播放| 国产成人啪精品午夜网站| 身体一侧抽搐| 亚洲av电影在线进入| 成人三级做爰电影| 一区二区三区激情视频| 俄罗斯特黄特色一大片| av一本久久久久| 久久久久久人人人人人| 亚洲欧美一区二区三区黑人| 制服人妻中文乱码| 日韩免费av在线播放| 亚洲国产精品一区二区三区在线| 大香蕉久久成人网| 免费在线观看日本一区| 国产一区有黄有色的免费视频| 丝袜人妻中文字幕| 99国产精品99久久久久| 国产真人三级小视频在线观看| 超色免费av| 精品高清国产在线一区| 久久狼人影院| 欧美人与性动交α欧美精品济南到| 欧美性长视频在线观看| 亚洲成av片中文字幕在线观看| 很黄的视频免费| 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 18禁裸乳无遮挡动漫免费视频| 在线观看66精品国产| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 欧美乱色亚洲激情| 午夜福利,免费看| 国产亚洲欧美98| 亚洲精品在线美女| 在线看a的网站| 一级毛片精品| 精品人妻在线不人妻| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 一二三四社区在线视频社区8| 91麻豆av在线| 国产一区二区三区视频了| 性色av乱码一区二区三区2| 成人国语在线视频| 亚洲三区欧美一区| 午夜福利乱码中文字幕| 国产亚洲精品第一综合不卡| 日韩大码丰满熟妇| 91九色精品人成在线观看| 80岁老熟妇乱子伦牲交| 午夜老司机福利片| 亚洲精品成人av观看孕妇| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 国产真人三级小视频在线观看| 女警被强在线播放| 女同久久另类99精品国产91| 久久国产精品大桥未久av| 最近最新中文字幕大全电影3 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲久久久国产精品| 777米奇影视久久| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 一边摸一边抽搐一进一出视频| 亚洲av成人av| 校园春色视频在线观看| 久久人人97超碰香蕉20202| 亚洲精品自拍成人| 韩国av一区二区三区四区| 高清视频免费观看一区二区| 老司机午夜福利在线观看视频| 欧美 日韩 精品 国产| av超薄肉色丝袜交足视频| 天天躁日日躁夜夜躁夜夜| 人妻一区二区av| 少妇的丰满在线观看| 亚洲国产欧美日韩在线播放| 最近最新中文字幕大全免费视频| 热re99久久国产66热| www.999成人在线观看| 夫妻午夜视频| 免费黄频网站在线观看国产| 99香蕉大伊视频| 国产精品1区2区在线观看. | 国产色视频综合| 欧美精品啪啪一区二区三区| 19禁男女啪啪无遮挡网站| 黑人猛操日本美女一级片| 悠悠久久av| 1024视频免费在线观看| 在线观看舔阴道视频| 精品人妻熟女毛片av久久网站| 欧美另类亚洲清纯唯美| 精品电影一区二区在线| 免费看十八禁软件| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕 | 久久国产乱子伦精品免费另类| 后天国语完整版免费观看| 水蜜桃什么品种好| 叶爱在线成人免费视频播放| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 亚洲人成伊人成综合网2020| 999久久久精品免费观看国产| 免费人成视频x8x8入口观看| 王馨瑶露胸无遮挡在线观看| 中文字幕人妻丝袜一区二区| 日韩欧美在线二视频 | bbb黄色大片| 999久久久精品免费观看国产| 亚洲精品自拍成人| 热99国产精品久久久久久7| 亚洲av日韩精品久久久久久密| 久久国产精品人妻蜜桃| 国产不卡一卡二| 性少妇av在线| 少妇粗大呻吟视频| 久久久久国内视频| 国产精品国产高清国产av | 怎么达到女性高潮| 日韩欧美一区二区三区在线观看 | 在线观看免费高清a一片| 久久中文看片网| 欧美中文综合在线视频| 色在线成人网| 久久这里只有精品19| 亚洲精品国产区一区二| 亚洲精华国产精华精| 国产一卡二卡三卡精品| 精品一区二区三区四区五区乱码| 精品一区二区三区视频在线观看免费 | 中文字幕色久视频| 久久精品国产a三级三级三级| 久久久国产成人免费| 久久香蕉精品热| 亚洲av欧美aⅴ国产| 女性被躁到高潮视频| 欧美日本中文国产一区发布| 久久午夜综合久久蜜桃| 国产欧美日韩一区二区三| 自拍欧美九色日韩亚洲蝌蚪91| 丝瓜视频免费看黄片| 99香蕉大伊视频| 亚洲精品中文字幕在线视频| 欧美亚洲 丝袜 人妻 在线| 水蜜桃什么品种好| 久久精品国产亚洲av高清一级| 天天躁日日躁夜夜躁夜夜| 国产三级黄色录像| 午夜免费鲁丝| 国产精品免费一区二区三区在线 | 一级毛片高清免费大全| 国产不卡一卡二| 亚洲成a人片在线一区二区| 韩国精品一区二区三区| 在线观看免费午夜福利视频| 999精品在线视频| 大香蕉久久成人网| 建设人人有责人人尽责人人享有的| 久久精品91无色码中文字幕| 9热在线视频观看99| 男女下面插进去视频免费观看| 亚洲自偷自拍图片 自拍| 极品教师在线免费播放| 搡老熟女国产l中国老女人| 91成年电影在线观看| 99久久99久久久精品蜜桃| 国产免费现黄频在线看| 啦啦啦视频在线资源免费观看| 国产欧美亚洲国产| 窝窝影院91人妻| 国产精品九九99| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 日日夜夜操网爽| 成人免费观看视频高清| 日日夜夜操网爽| 大香蕉久久网| 美女福利国产在线| 久久国产精品大桥未久av| 香蕉久久夜色| 精品免费久久久久久久清纯 | 国产亚洲欧美在线一区二区| 搡老熟女国产l中国老女人| 久久久久久久精品吃奶| 久久久久国产精品人妻aⅴ院 | 欧美人与性动交α欧美软件| 91字幕亚洲| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 露出奶头的视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美在线一区亚洲| 在线观看66精品国产| 国产亚洲精品久久久久5区| 9热在线视频观看99| 日韩大码丰满熟妇| 飞空精品影院首页| 99热网站在线观看| 亚洲国产看品久久| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲欧美在线一区二区| 天天躁日日躁夜夜躁夜夜| 大型av网站在线播放| 一区二区三区激情视频|