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

    Insights into the radial water jet drilling technology-Application in a quarry

    2018-04-24 00:54:58ThomasReinshBoPaapSimonHahnVolkerWittigSineyvanenBerg

    Thomas Reinsh,Bo Paap,Simon Hahn,Volker Wittig,Siney van en Berg

    aGFZ German Research Centre for Geosciences,Telegrafenberg,14473 Potsdam,Germany

    bTNO,Sustainable Geo Energy,Utrecht,The Netherlands

    cInternational Geothermal Centre-GZB,Lennershofstr.140,44801 Bochum,Germany

    dWell Services Group,Phileas Foggstraat 65,NL-7825 AL Emmen,The Netherlands

    1.Introduction

    Applying high pressure water jets to penetrate rocks is extensively studied since the middle of the 20th century(e.g.Farmer and Attewell,1965).Applications range from cutting/carving applications(e.g.Harris and Mellor,1974;Summers and McGroarty,1982;Hagan,1992)to downhole drilling(e.g.Maurer et al.,1973;Deily,1977;Pols,1977a,b).In recent years,the application of the radial water jet drilling(RJD)technology is increasingly investigated to perforate and stimulate low performing wells(e.g.Buset et al.,2001;Bruni et al.,2007;Cirigliano and Talavera Blacutt,2007;Seywald and Marschall,2009;Abdel-Ghany et al.,2011;Elliott,2011;Cinelli and Kamel,2013).For the RJD technology,a bottom hole assembly,generally referred to as ‘def l ector shoe’,is connected to a tubing and lowered to the target depth.For cased hole intervals,a coiled tubing conveyed milling assembly is lowered through the tubing.At the bottom,the def l ector shoe def l ects the milling bit towards the casing.After milling a hole into the casing,a jetting assembly is lowered through the tubing.The coiled tubing conveyed jetting assembly consists of a self-propelled jetting nozzle attached to a flexible hose.This assembly is capable of jetting up to 100 m into the formation.

    Recently,the RJD technology has gained considerable interest for stimulating low performing geothermal wells,motivating numerical investigations to estimate the benefit from applying this technology(e.g.Peters et al.,2015).For enhanced or engineered geothermal systems,it is considered as a viable alternative to conventional hydraulic stimulation technologies and has triggered several research projects(e.g.Reinsch and Bruhn,2016).As one of the first applications in a geothermal environment in Europe,RJD stimulation was performed in a low performing injection well in Klaipeda,Lithuania.RJD was applied to jetting 12 laterals with a length up to 40 m each,leading to an increase of injectivity of about 14%although a maximum increase of 57%was suggested by numerical modeling(Nair et al.,2017).Analysing the sensitivity to different model parameters like lateral direction and length,which influences the predicted increase by up to 10%for individual parameters,underlines the importance of proper monitoring equipment to measure the laterals’geometry(Nair et al.,2017).

    While rock penetration experiments can be performed in the laboratory with good access to the process parameters as well as the jetted holes afterwards,RJD operations are performed downhole with limited access to process parameters and no option to investigate the jetted holes(fora review on different laboratoryand field experiments,please see Bl?cher et al.,2016).In order to bridge this gap,we performed a jetting experiment in a quarry with the state-of-the-art RJD technology.Monitoring and in fl uencing process parameters like fluid pressure and fl ow rate could easily be implemented.In addition,the jetted holes were accessible afterwards for geometry measurements and could be surveyed visually.Although jetting individual holes at ambient surface conditions is not comparable to downhole applications with regard to the reservoir pressure and temperature conditions as well as saturation and stress conditions within the rock,various aspects of the technology can be investigated:

    (1)To control jetting parameters(pressure and f l ow rate)and monitor their effects on process parameters such as rate of penetration(ROP),fluid return,and cutting return.

    (2)To monitor the direction of a jetted hole by acoustic/seismic measurements at the surface during the jetting process.

    (3)To inspect the inside of the jetted hole,thereby determining the shape of the jetted hole as well as its trajectory.

    (4)To observe the influence of local geological or structural interfaces on the jetting process as well as the hole geometry.

    (5)To determine the change in injectivity over the length of the jetted hole,especially when intersecting permeable features,e.g.a fracture.

    (6)To collect eroded rock particles/cuttings for further investigation.

    (7)To analyse if the jetting process parameters can be used to infer details about the jetted formation.

    (8)To test different jetting nozzles.

    1.1.Selection of quarry

    To ensure sufficient penetration depth and rate during the drilling experiment in a quarry,pre-quarry jetting tests were performed on different rock samples in yard tests at ambient surface conditions.Three types of sandstones were selected based on rock properties,mainly permeability and porosity,and the accessibility of the respective quarry(see Table 1).In addition,the pre-quarry tests were used to perform acoustic measurements during jetting to determine the frequency spectrum of the acoustic signal resulting from the jet-rock interaction.

    Either a static or a rotating RJD nozzle was attached to a flexible high pressure hose and connected to a high pressure triplex pump.To control the stand-off distance between the nozzle andthe rock surface,the nozzle was clamped to a hydraulic piston capable of steering it towards the rock surface.The rock samples were neither saturated nor immersed in water during the tests.Rock samples of approximately 50 cm×50 cm×50 cm were positioned in front of the nozzle and equipped with three types of acoustic sensors:

    Table 1Rock samples and respective quarries tested in the yard tests.The last column gives the International Geo Sample Number(IGSN).

    (1)Accelerometers for frequencies up to 800 Hz(Sercel DSU-3,200-800 Hz).

    (2)Geophones for frequencies up to 5000 Hz(GS-14).

    (3)Piezo-elements for frequencies up to 10,000 Hz(BK 4514-B).

    Due to the full overlap in frequency response of the geophones with the other two instruments,the geophones were considered redundant and not used for monitoring in the quarry experiment.Table 1 lists the maximum ROP achieved during the yard tests.It is shown that the highest ROP was obtained for the Cretaceous Sandstone from Gildehaus,Germany.For a static nozzle,significantly lower ROPs were observed compared to a nozzle with a selfpropelled rotating head,referred to as rotation nozzle.Therefore,it was decided to perform the jetting experiment mainly with a rotating nozzle in the Gildehaus quarry,Germany.

    Based on previous studies by Hardy Jr.(2003),Kovacevic et al.(1998),and Rabani et al.(2012)and on the pre-quarry rock jet tests in Bochum,the signal of the jetting activity in the quarry experiment was expected to lie in the frequency range of 500-10,000 Hz.Based on this,we decided to deploy three-component accelerometers combined with piezo-elements in the quarry,since they can together properly span the expected frequency band of the acoustic activity.

    1.2.Description of quarry and geology

    The quarry Gildehaus is located in North-West Germany,3.5 km west of Bad Bentheim and close to the border with the Netherlands(52°18′8.42′N,7°6′18.76′E).It is located about 50 m above mean sea level between the Ems-and Münsterland,within the Westf?lische Bucht.The stratigraphyof the sandstone canbe classified as Lower-Valanginian,Lower Cretatious(Haack,2007).

    At the onset of the Cretaceous time,a continental landmass developed in the northern part of Germany.Therefore,siliciclastic sediments were deposited.During the subsequent subsidence of the area,the Nieders?chsische Transtensionsbecken as well as the Münsterl?nder Kreidebucht was formed and a shelf sea,the chalk sea,was able to transgress.Within the Lower Cretaceous,diagenetic processes formed the up to 70 m thick Bentheimer sandstone.Onto the Bentheimer sandstone,marl and limestones were deposited within a depth up to 250 m in the sea.Climate changes at the end of the Cretatious led to a regression of the Cretaceous sea.Due to the periodical trans-and regression of the North Sea,successive sedimentary layers were deposited during the Tertiary age.Since the Quaternary,Tertiary and Cretaceous sediments were almost entirely eroded,in some locations such as Bad Bentheim,Cretaceous sediments are outcropping within sediments of mainly glacial and f l uviatile origin(Meschede,2015).Due to the uplift of a salt-diapir beneath Bad Bentheim,an anticline formed.The layering of the Bad Bentheim sandstone therefore dips to the south.

    The quarry Gildehaus covers an area of about 18,000 m2with a length of about 300 m.It can be subdivided into five excavation sites.At its lowest point,the quarry is about 30 m deep.Our jetting experiments were performed in the central part of the quarryat the southern side.

    Fig. 1.View to the east onto the quarry wall used for the jetting experiments.Subvertical red lines indicate major fracture zones.The sub-horizontal yellow line indicates the direction of the layering,whereas the thick sandstone layer can be seen below and the interlayered section can be seen above.

    Fig. 2.Fracture face in Gildehaus quarry.Reddish-brown colours originate from ironmanganese minerals.

    The outcropping sandstone consists of a homogeneous,weakly weathered,and non-graded sandstone.The lower part of the wall is made of a massive sandstone layer with a thickness of up to 10 m which is currently quarried,mainly in the western part of the quarry.It is overlain by a unit of about 5 m in thickness,consisting of 5-20 cm thick sandstone layers separated by thin clay layers(Fig. 1).On top,a massive sandstone with little layering can be found.

    Overall,the colour of the Bad Bentheim sandstone is grey-beige to brownish if weathered and light-grey to beige on fresh surfaces.Individual clay layers can clearly be identified by their darker appearance.On fracture faces,iron-manganese minerals are abundant and can be identified by their reddish-brown colour(see Fig. 2).

    On average,individual sandstone layers have a strike and dip of about 183°/28°.Fractures instead show a strike and dip of 25°/76°.Due to the anticline structure,individual fracture faces can dip to the south/south-west,rotated by 180°.Individual larger fractures show an offset of about 1-3 cm.The aperture of the individual fractures varies from<0.1 mm to 1 cm at the quarry wall.A single major fracture zone with a width of 17 cm was observed,which is divided into multiple smaller fractures.Individual fractures have a terraced,undulating,and smooth surface.The mainly well sorted fine-to-medium size grained sandstone is free of lime and consists of 95%quartz grains which are moderately rounded and cemented with a siliciclasic cement.The sandstone has a feldspar content of about 2%and can be classi fi ed as Arenite(DIN,2003;Stow,2008).A porosity up to 22.8%is reported(Klein et al.,2001).

    2.Experimental set-up

    2.1.Hydraulic set-up

    Before the jetting experiment,three ‘wellheads’were prepared at the quarry wall.A coring machine was used to drill 67 mm diameter holes up to 30 cm into the quarry wall.Into these holes,2′′(≈5.08 cm)WECO 1502 pup-joints with a length of about 0.6 m were cemented with composite mortar.After drying,a hand-held drilling machine was used to penetrate the mortar at the bottom of the WECO tubes.During jetting,a hydraulic system was attached to the pup-joints.Borrowing terminology fromoilfield applications,the necessaryset-up is herereferred toas X-mas tree.Fig. 3 shows a schematic overview of the set-up together with different sensors installed in the hydraulic system.

    Fig. 3.Wellhead set-up during jetting operation.For details on the geometry,please see Table 2.

    Table 2Length information(m)for the different jetted holes.See Fig. 3 for the nomenclature of different components.The depth(m)of the holes(e)is evaluated using information from drilling and subsequent logging campaigns.

    The X-mas tree started with a T-piece.The downward facing part redirected the back- fl ow including the cuttings from jetting.This fluid went through a valve and a crossover to a 1′hose.It was fi ltered to retrieve the cuttings and directed into a 1 m3Intermediate Bulk Container(IBC).Here,the return fluid volumes could be measured.On the collinear side of the T-piece,another T-piece was connected.The downward facing part was connected to a blind-cap with a 1/2′NPT thread,allowing for injection testing through this inlet.The collinear connection was connected to a valve.Through this valve,the jetting hose was fed into hole.During injection testing,both valves at the end of the X-mas tree as well as at the fluid return were closed.

    For jetting and injection testing,two pumps were used.For the injection testing,we used a ‘low pressure’pump.This pump is capable of delivering 60-200 L/min at a maximum pressure of 600 bar(1 bar=100 kPa).The maximum pressurewas set to 50 bar.For the jetting,a ‘high pressure’pump capable of delivering 11-30 L/min at a maximum pressure of 1000 bar was used.Due to the pressure rating of the equipment,the maximum operating pump pressure was set to 700 bar.

    The ‘low pressure’pump was connected to the wellhead via a 100 barrated 3/4′ID hose.The ‘high pressure’pump was connected to the jetting hose via a 700 bar rated flexible hose with a pressure gauge before the jetting hose.The jetting hose itself was a 700 bar rated flexible hose.For the first jetting experiments,we used a length of 50 m.Later,we shortened it to 15 m to increase the available fluid pressure at the nozzle.For RJD applications,pressure friction loss within the CT as well as the flexible hose determines the available pressure at the nozzle.In downhole applications,the pressure at the nozzle can be increased by decreasing the length of the CTor the flexible hose or by increasing diameters.In our jetting experiment,wedid notaim fora specific pressure at the nozzle.The experiment,therefore,can only give indications on performance parameters,e.g.for the ROP,when compared to downhole applications.Parameters like hole geometry and hole survey,however,should be comparable.

    During jetting and testing at one wellhead,the remaining two wellheads were plugged and pressure gauges were used to monitor any pressure changes.As none of the other wellheads was fluid fi lled and the formation was not saturated,no pressure change was measured.Fig. 4 shows the wellhead installation on the second jetted hole.High pressure hoses as well as pressure gauges and seismic receivers attached to the wall are shown.

    As we were not able to jet the planned 25 m in holes 1-3,it was decided to jet two more holes,i.e.holes 4 and 6.For holes 4 and 6,no wellhead with X-mas tree was installed.For holes 4 and 6,pilot holes were drilled with a 20 mm diameter rotary bit down to 15 cm and a 16 mm rotary bit down to 30 cm,respectively.Due to the missing wellhead,no injection tests could be performed.Cuttings were collected using a bucket placed along the vertical rock surface beneath the drilled hole.The position and direction of the pilot holes 4 and 6 were chosen not to intersect with any interfaces(bedding or fractures)within the surrounding rock mass.Hole 5 was marked but not drilled.

    2.2.Acoustic monitoring

    The acquisition layout is shown in Fig. 5.The monitoring instrumentation consisted of three-component accelerometers and unidirectional piezo-elements.Sensors were deployed with 2 m spacing along two lines termed L1 and L2 that were approximately parallel to ground level.Along line L1,both accelerometers and piezo-elements were placed to compare the performance of the two types and test the localisation method.Line L1 was used as a proxy for a downhole application where sensors can only be placed vertically in a well.In the quarry,we used the possibility to sample the acoustic wavefield in another plane,by placing accelerometers along line L2.

    Fig. 4.Wellhead of hole 2 installed in quarry wall.Wellhead 1 was equipped with a pressure sensor.Three-component geophones were cemented into the quarry wall(red)together with piezo-elements for high frequency recording.Hole numbers are indicated.

    Fig. 5.Geometry of sensors relative to individual jetted holes.The sensors are subdivided into a L1-and a L2-line.For the x′z and y′z projections,the coordinate system was rotated by-45°.

    During the jetting process,we expected a more or less continuous acoustic signal,and therefore needed to adopt techniques that are capable of processing continuous signals into finite signals to provide information on source location.The localisation technique thatweusedhererelied ondrill-bit interferometryas applied in the oil and gas industry,which was thoroughly explained and demonstrated byRectorand Marion(1991),and Polettoet al.(2004,2010,2014).

    We implemented the following workf l ow to estimate the location of the radial path:

    (1)Recorded data are conditioned by suppressing harmonic or cultural noise and isolating signal that comes directly from the interaction between nozzle jet and rock.Therefore,we compared data acquired during jetting periods to data from periods when there was no activity.We used a Butterworth bandpass f i lter to reject undesired signal.

    (2)Synthetic travel times of the P-wave arrival are calculated for all stations using a P-wave velocity model derived from active hammer shot data acquired at the quarry wall.Synthetic travel time data are calculated for all relevant source-receiver combinations.

    (3)The continuous data are divided into distinct files of specific duration.The duration is chosen such that it represents a time window within which we can assume that the nozzle source is at a stationary location.

    (4)Data of all station pairs are cross-correlated.Subsequently,the envelope of the resulting cross-correlated data is calculated,along which energy is stacked along synthetic delay times per synthetic source location.

    (5)This yields a maximum value for the synthetic source location for each time window that best explains the observed travel time data.

    2.3.Jetting program

    It was planned to jet three holes with different angles towards the local fracture network as well as the bedding with a length up to 25 m.After every 4 m,it was planned to perform a short-terminjection test.We assumed a decreasing injection pressure for a specific f l ow-rate over the length of the jetted hole.In case a fracture was intersected,this pressure was assumed to drop drastically.Data were supposed to show if any permeable features within the quarry wall were intersected.Between individual jetting intervals,the holes were inspected with a camera.During the tests,static and rotating nozzles were tested.Table 3 gives an overview of the different nozzle geometries tested during the experiment.Table 4 displays an overview about the different activities during the quarry experiment.

    Table 3Nozzles to be tested.

    3.Results

    We were able to jet individual holes with a length up to 12.3 m during the experiment.As we were not able to reach the anticipated 25 m,we decided to jet the additional holes 4 and 6.Injection tests into the unsaturated large open fractures,connected to the surface of the quarry wall,did not produce an evaluable dataset.Injection tests were therefore only performed in holes 1 and 2.In holes 4 and 6,total circulation loss was observed at a depth of 10.25 m and 2.9 m,respectively,indicating the connection to a high permeable fracture.For holes 1-3,the return f l ow data were inconclusive and a total circulation loss was not detected.Table 2 lists the geometric information for different holes that were jetted during the quarry experiment.

    3.1.Jetting performance

    During jetting,we measured the ROP in relation to jetting pressure and nozzle geometry.Using static nozzles(Table 3),we observed no ROP in the quarry.Individual ROP measurements were averaged over 25-50 cm intervals.For the rotating nozzles with 3 and 4 forward directed orif i ces,no conclusive results were obtained that favour one geometry over the other(Fig. 6).In addition,the length of the jetting hose was reduced from 50 m to 15 m during the experiment.Reducing the length of the hose decreased the pressure friction loss by about 100 bar at the maximum f l ow rate.For the additional pressure of 100 bar at the nozzle,a significant and sustainably increased ROP was not observed.On average,the ROP was a couple of cm/min.Anincreased ROP was observed in the first section of hole 1(0-3.5 m depth),the second interval in hole 2(1-2 m depth)and the single interval in hole 3.Holes 4 and 6,intentionally inclined not to intersect any structural feature(bedding or fracture),show a more stable ROP.

    Table 4List of experiments performed in different holes(Y=yes,N=no).

    Fig. 6.Rate of penetration for the different holes calculated for individual sections along each hole.The number associated to each bar indicates the number of forward directed orif i ces for the nozzle used in the test.‘l’and ‘s’indicate the length of the high pressure hose(long:50 m,and short:15 m).

    3.2.Inspection of jetted holes

    3.2.1.Visual inspection

    Each hole was inspected with an endoscope made from a USB camera fastened to a f i bre glass rod here referred to as endoscope(Fig. 7).The endoscope was pushed into each hole until the total depth or until a narrow spot prohibited the endoscope to progress deeper into the hole.Overa largeinterval of each hole,the diameter was a little larger than 20 mm(Fig. 8).However,the drift is significantly smaller.We were not able to push an 18 mm rod for more than 17 cm in any of the holes.On the other hand,a geometry sensor of 10 cm length and 12 mm diameter could be pushed along the holes.In individual intervals,although not very prominent,a star-shape profile was observed,indicating erosion caused by the backward directed nozzles.

    Fig. 9 shows the endoscopic view of three fractures that were intersected during jetting.Each fracture was intersected under a different angle varying from an estimated 25°to almost 90°.On the fracture faces,iron-manganese minerals can be observed.Before the 90°fracture,a larger excavated volume can be observed indicating that the fracture face was harder to penetrate.Without ROP in front of the fracture face,the weaker sandstone was eroded.

    Fig. 10 showsthe bottomof hole 2,wherea larger void spacewas hit.In front of the endoscope,the sandstone surface was slightly eroded.However,the nozzle did not further penetrate into the formation at this depth.From the visual inspection,it is hard to determine if a fracture or a bedding feature was intersected.Although,for individual fracture faces,the brownish colour is a good indicator for a fracture interface,thin layers of less resistant mudstone were washed out during jetting.If an offset corresponds to a fracture or a washed out lithological layer,it is difficult to distinguish.

    Fig. 7.Home-made endoscope used for inspecting holes.A USB camera was attached to a f i bre glass rod and pushed into all holes.

    3.2.2.Trajectory

    The geometry was measured with a three-dimensional(3D)accelerometer combined with a 3D magnetic field sensor.The sensor was attached to a f i bre glass rod and pushed into the hole.While pulling out of hole,a measurement of azimuth and inclination was performed every 25 cm.

    None of the measured trajectories is completely straight.All trajectories either have a rather constant change in direction over the entire length of each hole,or short radius changes at specific depths.An overview about the individual geometries is given in Fig. 11,where different two-dimensional(2D)projections are shown together with the dogleg severity,to indicate the change in direction per depth interval.For holes 1-3,the metal pipe was magnetized and in fl uenced the trajectory measurement.For the fi rst depth interval,where elevated magnetic field data indicated a strongly magnetized pipe,the initial direction was estimated to be normal to the quarry wall.When leaving the pipe,all holes show a high dogleg severity indicating that the jet hole was initiated slightly deviated from the axis of the wellhead.

    Fig. 8.Endoscopic view into the jetted hole.The red lines indicated the hole geometry.It has a diameter of about 20 mm.

    Fig. 9.Selected endoscopic view into two holes.(a)A fracture was intersected at an angle of about 45°.The fracture face is characterized by brownish colours.(b)A fracture was intersected at an angle close to 90°.A large volume was excavated before the fracture.The fracture face is characterized by brownish colours.And(c)A fracture was intersected at an estimated angle of 25°.The fracture is indicated by the dashed red line.

    Fig. 10.Endoscopic view into the second jetted hole.The jetting nozzle hit a larger void space.Here,ROP decreased to zero and the hole was terminated.

    Overall,the trajectories seemed to be diverted upwards over the length of each hole.Even when initially directed downwards,the curvature changed towards the top of the quarry.In addition to large radius direction changes towards the top of the quarry as well as within the horizontal plane,short radius changes can be observed in the dogleg severity plot.Apart from the end of the wellheads in holes 1-3,strong change in direction was observed for hole 1,where two distinct depths at 1-2 m as well as 4-5 m were located.From the visual inspection,apart from the onset of the jetted hole at the end of the pipe,a very heterogeneous hole geometry with larger cuttings remaining inside the hole was found in the first interval.This coincides with the depth of both thin bedding as well as intersection with a fracture.From the visual inspection,none of both interfaces can be excluded to have caused the change in direction.For the change in direction observed at greater depth,a larger oval-shaped caving was observed with a diameter of more than about 5 cm.At this depth,an interface was intersected.In addition,the hole split into two arms.To remove cuttings from the hole during jetting,the jetting hose was pulled out of hole from time to time.In this depth,a‘sidetrack’was initiated due to hole cleaning.Which one of both holes was entered with the geometry sensor cannot be determined from the visual inspection.The reason for the larger excavation causing the deviation might be the intersected interface.The deviations in hole 6,although initially directed such that it does not intersect any interface,show an almost constant change in direction at a depth larger than 2 m.Visual inspection showed that several interfaces including larger excavations could be observed from that depth.Neither structural features nor any bedding was observed before that depth.In hole 4,the geometry tool was stuck in about 1 m below surface.

    Fig. 11.Survey of jetted holes projected to a single reference point.

    3.2.3.Cutting analysis

    Fig. 12 shows the grain size distribution of cuttings retrieved from holes 1-3 where a reliable sampling was guaranteed by the hydraulic set-up.For all three holes,the grain size distribution is similarly corresponding to a medium-grained sand.From the grain size distribution,permeability values were calculated according to empirical formulas derived by Hazen(1893),Seelheim(1880),Beyer(1964),Bialas and Kleczkowski(1970)for unconsolidated sediments as reported in H?lting(1996).For each of the three holes,the calculated permeability values are in the range of 3-17 Darcy(Table 5).This value is very close to permeability values derived from laboratory experiments on intact rock samples of 0.8 Darcy at 50 bar of confining pressure(IGSN:GFTRE0066).

    3.3.Acoustic monitoring

    During the jetting tests at the Gildehaus quarry,acoustic signal resulting from jetting activity was successfully recorded using accelerometer and piezo-element instruments.Here,we demonstrate the localisation technique for the piezo-element and accelerometer measurements acquired while jetting in hole 1,where a penetration of 5.4 m was reached.Fig. 13 shows the data recorded bypiezo-elementsand accelerometerswhilethe maximum penetration at hole 1 was reached.Both datasets were frequency f i ltered in the same bandpass of 350-950 Hz for comparison.

    Fig. 12.Grain size distribution of cuttings sampled during the RJD operation from holes 1-3.

    Table 5List of formation permeability in Darcy calculated from grain size distribution of the cuttings as an average and standard deviation from holes 1-3.

    Fig. 14 shows the results of the localisation of the acoustic signal received by the piezo-elements and accelerometers while the radial drill at hole 1 was performed.The localisation result in Fig. 14 is shown at the start and termination of the jet,illustrating the source location with highest probability indicated in yellow.The figure shows increased probability(yellow)at the location where the radial is placed,and the estimated location remains close to the quarry wall during drilling.These results are based on a specific data selection with source-receiver offsets of 1-5 m and data inversion is based on a uniform velocity model of 1500 m/s.This velocity was determined from velocity analysis of a hammercalibration shotline conducted along the quarry wall.The location estimations found with the piezo-element and the accelerometer data agree with each other.Their estimated locations have an offset of 5 m compared to the final location as measured with the in situ accelerometer tool.

    4.Discussion

    4.1.The ROP

    Although the ROP was measured during jetting in different lithologicalsettings,with differentnozzle geometries and changing fluid pressures,a conclusive result could not be obtained.Neither jetting in the homogeneous rock formation(e.g.in the first 2 m of hole 6),nor varying nozzle pressure(e.g.in hole 2,where 50 m and 15 m high fluid pressures were used),nor changing the nozzle geometry from a nozzle with 4 to a nozzle with 3 forward directed orif i ces changed the ROP significantly.Individual high ROP values cannot be attributed to any specific downhole situation.One reason might be the spatial and temporal resolution of the ROP measurements.As individual measurements were taken every 25-50 cm,small-scale features like the effect of bedding or fractures on the ROP are obscured.From the geometry,it is known that all holes tend to be directed upwards into the part of the quarry with inter-bedded clay layers.An effect of the lithology on the ROP was not observed with the current set-up.Compared to a real field scale application,ROP values in the quarry were significantly lower than expected but corresponded to the ROPs observed during laboratory experiments at surface conditions.

    For the jetting experiments,rotating nozzles were applied.The nozzles created a mainly roundish hole geometry.Before individual fracture faces,larger excavations were observed,indicating a very slow ROP through the interface.Individual clay layers could not be observed visually.The soft clay was likely washed out,making it very difficult to distinguish between a fracture or several closely spaced fractures with a certain aperture and individual clay layers.From the visual inspection,however,it was observed that distinct interfaces can be intersected under different angles without changing the direction of the jetting nozzle.From the visual inspection combined with the geometry information,it can be concluded that all holes tend to orient themselves towards the bedded quarry interval,wherepotentiallysofterformations prevail.A clear indication from the ROP,however,was not observed.In larger cavities,either washed out by the jetting operation or naturally existent within the rock,the orientation seems to be influenced,as the backwards directed nozzles cannot orient the nozzle sufficiently.

    Fig. 13.Acoustic data acquired during jetting hole 1 along L1.(a)Acoustic waveform data recorded at 16 piezo-element stations.The vertical component data are plotted having the same orientation as the piezo-elements.(b)Acoustic waveform data recorded at 16 accelerometer stations.Both datasets were passed within the frequency band of 350-950 Hz.

    Fig. 14.Results of localization of acoustic data from hole 1.(a,b)Estimated location based respectively on the piezo-element and the accelerometer data during the start of the jetting.(c,d)Estimated location based respectively on the piezo-element and the accelerometer data when maximum penetration was reached.

    From the cutting analysis,assuming an unconsolidated sand,permeability values could be calculated.For the very homogenous,high porosity sandstone in Gildehaus,the calculated permeability values correspond very well to permeability data measured on intact rock samples under laboratory conditions.Although a sandstone as homogeneous as the Gildehaus sandstone will likely not be found in typical downhole applications,valuable data might be gathered from the cuttings.

    4.2.Acoustic monitoring

    Careful inspection of the acoustic data proved that only a limited part of the data was suited to perform inversion of continuous signal data to location estimation.Careful selection of data turned out to be especially important.Data should contain a minimum amount of undesired noise,such as man-made activity along the wall.Additionally,the strength and signature of acoustic data showed significant variations during individual jetting tests,which might be related to the various speeds at which the nozzle moved as well as the inhomogeneities within the rock.These signal variations appeared to have a significant effect on the robustness of the location estimates.This is partially related to the simpleone-dimensional(1D)velocitymodelusedherefor inverting travel times to source location estimations.The actual structure of the quarry wall is more complex,containing a number of fractures that will affect transmission and reflection behaviours of propagating waves.This will result in biases of measured travel times obtained from recorded station data,with respect to simulated travel times of the best fitting source location.Furthermore,reflections occurring along the free surface of the quarry wall may complicate an accurate estimation of the nozzle location and reconstruction of the jet path.Additionally,by using higher frequency information for small source-receiver offsets,the source estimation might be improved.A challenge here is that the accuracy of the geometry(source,receivers and geology)becomesincreasinglyimportantwhen considering higher frequencies to correctly predict time differences in travel time arrivals between stations.The accuracy of estimation of the jetted path is significantly smaller compared to the in situ location measurements.However,the acoustic localization can provide a mean of identifying zones where the nozzle has passed,and provides a clue in which direction the path has been set.It is encouraging that the inversion results of part of the piezoelement and accelerometer data agree with each other.The offset in location between acoustic inversion and accelerometer tool probably results from the simplif i ed homogeneous velocity model that we assumed.

    5.Conclusions

    (1)Monitoring the acoustic activity generated by the jetting process can in principle provide location estimations of the jetting nozzle within the rock mass.The success of this location method at the quarry site mainly depends on careful data selection,suppression of undesired noise,and the degree of realism accommodated in the forward velocity model used to calculate synthetic travel times that are in turn used to estimate the nozzle locations.A downhole application of geophones can provide insight in the approximate spatial distance towards the nozzle,provided that a sufficient aperture and sensor distribution can be reached within the well.At this stage,however,the accuracy of the jet geometry estimated from downhole acoustic measurements with the current data analysis procedure will be limited compared to in situ measurements.

    (2)The jetted holes have a round geometry but a limited drift close to the diameter of the nozzle itself.Intersecting various geological boundaries,the trajectory of the hole is far from straight.Although data suggest that interfaces can be intersected under different angles,it cannot be excluded that structural features caused the deviation of the hole.For this experiment,the spatial resolution of the geometry measurement was to low(25 cm),and the geometry information reasonably gained for individual interfaces(plural)from geological mapping at the rock surface was not accurate enough when projected to depth.

    (3)The pressure information from injection experiments alternating with the jetting could not be used to evaluate downhole permeabilities due to large void spaces and multiple high conductive features intersecting the holes.It cannot be excluded that pressure information during the jetting process downhole can be used to have an online information about permeability variations within the formation.Performance parameters like ROP,jetting pressure,fluid velocity,and return fluid volume,however,could not be used to assess the hydraulic variability downhole within the quarry.

    (4)As already observed in the pre-quarry tests,rotating nozzles proved to generate the highest ROP values in the quarry.

    (5)Collected rock cuttings could be used to accurately assess the permeability of the penetrated rock mass.It remains to be proven that this can be applied to the downhole situation as well.

    (6)To reduce the curvature of the jetted hole,the jetting nozzle should be stabilized.

    Conflict of interest

    The authors wish to confirm that there are no known Conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.

    Acknowledgements

    Abdel-Ghany M,Siso S,Hassan A,Pierpaolo P,Roberto C.New technology application,radial drilling petrobel, fi rst well in Egypt.In:Offshore Mediterranean conference and exhibition,23-25 March,Ravenna,Italy.Offshore Mediterranean conference;2011.OMC-2011-163.

    Beyer W.Zur Bestimmung der Wasserdurchl?ssigkeit von Kiesen und Sanden aus der Kornverteilungskurve.Wasserwirtschaft-Wassertechnik (WWT)1964;14(6):165-9(in German).

    Bialas Z,Kleczkowski A.On practical use of certain empirical formulas to determine coefficient of permeability k.In:Archivum Hydrotechniki,Warschau;1970.

    Bl?cher G,Peters E,Reinsch T,Petrauskas S,Valickas R,van den Berg S.Report on radial jet-drilling(RJD)stimulation technology.GFZ German Research Centre for Geosciences;2016.https://doi.org/10.2312/GFZ.6.2.2018.001.

    Bruni MA,Biasotti JH,Salomone GD.Radial drilling in Argentina.In:Latin American&Caribbean petroleum engineering conference,15-18 April,Buenos Aires,Argentina.Society of Petroleum Engineers;2007.https://doi.org/10.2118/107382-MS.SPE-107382-MS.

    Buset P,Riiber M,Eek A.Jet drilling tool:cost-effective lateral drilling technology for enhanced oil recovery.In:SPE/ICoTA coiled tubing Roundtable,7-8 March,Houston,Texas.Society of Petroleum Engineers;2001.https://doi.org/10.2118/68504-MS.SPE-68504-MS.

    Cinelli SD,Kamel AH.Novel technique to drill horizontal laterals revitalizes aging field.In:SPE/IADC drilling conference society of petroleum engineers,5-7 March,Amsterdam.The Netherlands;2013.https://doi.org/10.2118/163405-MS.SPE-163405-MS.

    Cirigliano RA,Talavera Blacutt JF.First experience in the application of radial perforation technology in deep wells.In:Latin American&Caribbean petroleum engineering conference,15-18 April,Buenos Aires,Argentina.Society of Petroleum Engineers;2007.https://doi.org/10.2118/107182-MS.SPE-107182-MS.

    Deily F.Five wells test high-pressure drilling.Oil and Gas Journal 1977;75(27):74-81.

    DIN.Geotechnische Erkundung und Untersuchung-Benennung,Beschreibung und Klassi fi zierung von Fels-Teil 1:Benennung und Beschreibung(Technical report).DIN Deutsches Institut für Normung e.V;2003(ISO 14689-1:2003);Deutsche Fassung EN ISO 14689-1:2003(in German).

    Elliott S.Coiled-tubing method drills radial laterals to improve oil production from a depleted reservoir.World Oil Magazine 2011.

    Farmer I,Attewell P.Rock penetration by high velocity water jet:a review of the general problem and an experimental study.International Journal of Rock Mechanics and Mining Sciences&Geomechanics Abstracts 1965;2(2):135-53.

    Haack.Haack-Weltatlas.Klett Ernst/Schulbuch.2007(in German).

    Hagan P.The cuttability of rock using a high pressure water jet.In:Western Australian conference on mining Geomechanics;1992.p.12.

    Hardy Jr HR.Acoustic emission/microseismic activity.Principles,techniques and geotechnical applications,vol.1.CRC Press;2003.

    Harris H,Mellor M.Cutting rock with water jets.International Journal of Rock Mechanics and Mining Sciences&Geomechanics Abstracts 1974;11(9):343-58.

    Hazen A.Some physical properties of sands and gravels with special reference to their use in fi ltration.1893.p.541-56.

    H?lting B.Hydrogeologie Einführung in die Allgemeine und Angewandte Hydrogeologie.5th ed.Ferdinand Enke Verlag;1996(in German).

    Klein E,Baud P,Reuschlé T,Wong TF.Mechanical behaviour and failure mode of Bentheim sandstone under triaxial compression.Physics and Chemistry of the Earth,Part A:Solid Earth and Geodesy 2001;26(1):21-5.

    Kovacevic R,Kwak HS,Mohan RS.Acoustic emission sensing as a tool for understanding the mechanisms of abrasive water jet drilling of difficult-to-machine materials.In:Proceedings of the Institution of mechanical engineers,Part B:Journal of engineering Manufacture,vol.212(1);1998.p.45-58.https://doi.org/10.1243/0954405981515491.

    Maurer WC,Heilhecker JK,Love WW.High-pressure drilling.Journal of Petroleum Technology 1973.https://doi.org/10.2118/3988-PA.

    Meschede M.Geologie Deutschlands.1st ed.Springer;2015.https://doi.org/10.1007/978-3-662-45298-1(in German).

    Nair R,Peters E,Sliaupa S,Valickas R,Petrauskas S.A case study of ra-dial jetting technology for enhancing geothermal energy systems at klaip?da geothermal demonstration plant.In:Proceedings,42nd Workshop on geothermal reservoir engineering.Stanford,California,USA:Stanford University;2017.p.11.

    Peters E,Veldkamp J,Pluymaekers M,Wilschut F.Radial drilling for Dutch geothermal applications.Technical Report.TNO Netherlands;2015.

    Poletto F,Malusa M,Miranda F,Tinivella U.Seismic-while-drilling by using dual sensors in drill strings.Geophysics 2004;69(5):1261-71.

    Poletto F,Corubolo P,Comelli P.Drill-bit seismic interferometry with and without pilot signals.Geophysical Prospecting 2010;58(2):257-65.

    Poletto F,Miranda F,Corubolo P,Schleifer A,Comelli P.Drill-bit seismic monitoring while drilling by downhole wired-pipe telemetry.Geophysical Prospecting 2014;62(4):702-18.

    Pols A.Hard-rock jetting-1:tests show jet-drilling hard-rock potential.Oil and Gas Journal 1977a;75(5):134-44.

    Pols A.Hard-rock jetting-conclusion:rock-type decided jetting economics.Oil and Gas Journal 1977b;75(6):71-5.

    Rabani A,Marinescu I,Axinte D.Acoustic emission energy transfer rate:a method for monitoring abrasive water jet milling.International Journal of Machine Tools and Manufacture 2012;61(Supp.C):80-9.

    Rector JW,Marion BP.The use of drill-bit energy as a downhole seismic source.Geophysics 1991;56(5):628-34.

    Reinsch T,Bruhn D.The SURE consortium.Novel productivity enhancement concept for a sustainable utilization of a geothermal resource-the SURE project.In:Proceedings,European geothermal Congress 2016,Strasbourg;2016.p.4.

    SeelheimF.MethodenzurBestimmungderDurchl?ssigkeitdesBodens.Zeitschriftfür Analytische Chemie 1880;19(1):387-418.https://doi.org/10.1007/BF01341054(in German).

    Seywald M,Marschall C.Usage of radial drilling by RAG.Oil Gas European Magazine 2009;3:114-6.

    Stow D.Sedimentgesteine im Gel?nde:Ein illustrierter Leitfaden(German Edition).Spektrum Akademischer Verlag;2008(in German).

    Summers DA,McGroarty S.The mechanism,and some parameters controlling,the water jet cutting of rock.In:The 23rd U.S.Symposium on rock mechanics(USRMS),25-27 August,Berkeley,California.American Rock Mechanics Association;1982.p.531-41.

    一级毛片aaaaaa免费看小| 久久久久久久精品精品| 亚州av有码| 在线观看三级黄色| 少妇人妻久久综合中文| 精品一区二区免费观看| 丰满人妻一区二区三区视频av| 久久久久久久久大av| 黄片wwwwww| 少妇丰满av| 国产国拍精品亚洲av在线观看| 婷婷色综合www| 国产人妻一区二区三区在| 2021少妇久久久久久久久久久| 男女免费视频国产| 黄色一级大片看看| 国产欧美日韩精品一区二区| 日韩欧美精品免费久久| 中文精品一卡2卡3卡4更新| 精品国产露脸久久av麻豆| 极品教师在线视频| 国产v大片淫在线免费观看| 欧美极品一区二区三区四区| 蜜桃久久精品国产亚洲av| 99热这里只有精品一区| 欧美日韩亚洲高清精品| 狠狠精品人妻久久久久久综合| 97在线人人人人妻| 欧美97在线视频| 成人综合一区亚洲| 伦精品一区二区三区| 欧美激情国产日韩精品一区| 久久精品人妻少妇| 成人影院久久| 少妇 在线观看| 男女无遮挡免费网站观看| 五月开心婷婷网| 中文欧美无线码| 一边亲一边摸免费视频| 人妻系列 视频| 少妇的逼水好多| 男女边摸边吃奶| 五月玫瑰六月丁香| 水蜜桃什么品种好| 中文资源天堂在线| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 99热这里只有是精品50| 纵有疾风起免费观看全集完整版| 国产男人的电影天堂91| 女性被躁到高潮视频| 又黄又爽又刺激的免费视频.| 亚洲成人av在线免费| 26uuu在线亚洲综合色| 国产综合精华液| 青春草国产在线视频| 久久久欧美国产精品| 午夜福利高清视频| 国产精品欧美亚洲77777| av女优亚洲男人天堂| 少妇精品久久久久久久| 免费观看av网站的网址| 少妇的逼好多水| 成年人午夜在线观看视频| 国产成人精品一,二区| 男女免费视频国产| 激情 狠狠 欧美| 国产中年淑女户外野战色| 80岁老熟妇乱子伦牲交| 中文天堂在线官网| 久久人人爽av亚洲精品天堂 | 国产精品不卡视频一区二区| 国产av码专区亚洲av| 久久久久久伊人网av| 少妇猛男粗大的猛烈进出视频| 熟妇人妻不卡中文字幕| 欧美区成人在线视频| 日韩制服骚丝袜av| 国产免费视频播放在线视频| 高清毛片免费看| 网址你懂的国产日韩在线| 免费观看a级毛片全部| 人人妻人人添人人爽欧美一区卜 | 黑丝袜美女国产一区| 一区在线观看完整版| 自拍欧美九色日韩亚洲蝌蚪91 | 国产色爽女视频免费观看| 国产av码专区亚洲av| 国产欧美另类精品又又久久亚洲欧美| 成年人午夜在线观看视频| 午夜免费观看性视频| 国产免费一级a男人的天堂| 亚洲av福利一区| 婷婷色综合www| 精品一区二区三卡| 91aial.com中文字幕在线观看| 国产精品女同一区二区软件| 国产精品蜜桃在线观看| 久久人妻熟女aⅴ| xxx大片免费视频| 国产精品一区www在线观看| 久久久a久久爽久久v久久| 国产免费福利视频在线观看| 成人影院久久| 丝瓜视频免费看黄片| 国产高清不卡午夜福利| 日韩av免费高清视频| 最近的中文字幕免费完整| 欧美97在线视频| 亚洲精品亚洲一区二区| 卡戴珊不雅视频在线播放| 高清在线视频一区二区三区| 精品人妻偷拍中文字幕| 国产亚洲最大av| 最近的中文字幕免费完整| 亚洲无线观看免费| 免费观看在线日韩| 最近中文字幕高清免费大全6| 一个人看的www免费观看视频| 国产日韩欧美在线精品| 全区人妻精品视频| 黄片wwwwww| 少妇猛男粗大的猛烈进出视频| 美女高潮的动态| 免费少妇av软件| 熟女av电影| 1000部很黄的大片| 亚洲av.av天堂| 亚洲高清免费不卡视频| 建设人人有责人人尽责人人享有的 | 91久久精品国产一区二区三区| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲网站| 精品一区二区三卡| 大香蕉久久网| 波野结衣二区三区在线| 午夜日本视频在线| 国产欧美另类精品又又久久亚洲欧美| 国产高清有码在线观看视频| 国产人妻一区二区三区在| 99久久人妻综合| 99re6热这里在线精品视频| 在线看a的网站| 少妇丰满av| 综合色丁香网| 乱系列少妇在线播放| 又爽又黄a免费视频| 成人免费观看视频高清| 天天躁夜夜躁狠狠久久av| 纵有疾风起免费观看全集完整版| 国产精品福利在线免费观看| 日韩一本色道免费dvd| 国产精品国产三级专区第一集| 又大又黄又爽视频免费| 国产亚洲最大av| 国产精品久久久久成人av| 99热这里只有精品一区| 欧美激情极品国产一区二区三区 | 国产av码专区亚洲av| 日日啪夜夜撸| 国产欧美另类精品又又久久亚洲欧美| 在线观看国产h片| 国产精品久久久久久精品电影小说 | 80岁老熟妇乱子伦牲交| 人体艺术视频欧美日本| 久久国产乱子免费精品| 国国产精品蜜臀av免费| 国语对白做爰xxxⅹ性视频网站| 国产视频内射| 中文字幕久久专区| 久久久亚洲精品成人影院| 国产精品精品国产色婷婷| 欧美 日韩 精品 国产| av福利片在线观看| 国产高清国产精品国产三级 | 成人国产av品久久久| 99九九线精品视频在线观看视频| 中文字幕制服av| 亚洲国产高清在线一区二区三| 麻豆国产97在线/欧美| 香蕉精品网在线| 插逼视频在线观看| 亚洲精品久久午夜乱码| 国产免费又黄又爽又色| 精品国产乱码久久久久久小说| 18禁裸乳无遮挡免费网站照片| 国产伦精品一区二区三区四那| 亚洲熟女精品中文字幕| 亚洲国产精品成人久久小说| 国产欧美亚洲国产| 噜噜噜噜噜久久久久久91| 狠狠精品人妻久久久久久综合| 高清黄色对白视频在线免费看 | 国产视频内射| av在线蜜桃| 国产精品女同一区二区软件| 一级毛片aaaaaa免费看小| 欧美日韩精品成人综合77777| 国产精品熟女久久久久浪| 国产伦理片在线播放av一区| 亚洲自偷自拍三级| 亚洲图色成人| 中文欧美无线码| 午夜福利影视在线免费观看| 另类亚洲欧美激情| 亚洲av.av天堂| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 久久久久人妻精品一区果冻| av在线老鸭窝| 在线观看一区二区三区| 久久精品国产亚洲网站| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| 高清毛片免费看| 亚洲激情五月婷婷啪啪| 久久 成人 亚洲| 大又大粗又爽又黄少妇毛片口| 一区二区三区精品91| 欧美成人一区二区免费高清观看| 国产精品99久久久久久久久| 高清av免费在线| 国产v大片淫在线免费观看| 人妻夜夜爽99麻豆av| 啦啦啦视频在线资源免费观看| 中文字幕久久专区| 精品国产乱码久久久久久小说| 亚洲精品久久午夜乱码| 性色av一级| 亚洲人成网站在线观看播放| 亚洲第一区二区三区不卡| 久久久久精品久久久久真实原创| 美女内射精品一级片tv| 国产爽快片一区二区三区| 岛国毛片在线播放| 亚洲美女视频黄频| 激情 狠狠 欧美| 美女高潮的动态| 高清av免费在线| 18禁在线无遮挡免费观看视频| 中文字幕免费在线视频6| 成人毛片60女人毛片免费| 日本爱情动作片www.在线观看| av黄色大香蕉| 精品久久久久久久久亚洲| 国产淫片久久久久久久久| 亚洲经典国产精华液单| 人人妻人人添人人爽欧美一区卜 | 又大又黄又爽视频免费| 婷婷色av中文字幕| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 女性生殖器流出的白浆| 天美传媒精品一区二区| 国产成人午夜福利电影在线观看| 精品久久久久久久久av| 高清在线视频一区二区三区| 边亲边吃奶的免费视频| 精品一区二区免费观看| 亚洲精品aⅴ在线观看| a级一级毛片免费在线观看| 日韩成人伦理影院| 欧美日韩在线观看h| 久久国产亚洲av麻豆专区| 国产欧美日韩精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 内地一区二区视频在线| 哪个播放器可以免费观看大片| 国产成人aa在线观看| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| 亚洲美女搞黄在线观看| 久久99精品国语久久久| 成年人午夜在线观看视频| 国产男女超爽视频在线观看| 丰满迷人的少妇在线观看| 亚洲三级黄色毛片| 中国三级夫妇交换| 欧美成人午夜免费资源| 亚洲欧美日韩无卡精品| 色视频www国产| av网站免费在线观看视频| 日韩欧美精品免费久久| 久久国产精品大桥未久av | 在现免费观看毛片| 美女脱内裤让男人舔精品视频| 超碰av人人做人人爽久久| 青春草视频在线免费观看| 男女下面进入的视频免费午夜| 乱系列少妇在线播放| 亚洲电影在线观看av| 久久 成人 亚洲| 亚洲无线观看免费| 成人漫画全彩无遮挡| 夜夜爽夜夜爽视频| 亚洲国产欧美在线一区| 观看免费一级毛片| 日本vs欧美在线观看视频 | 波野结衣二区三区在线| 毛片一级片免费看久久久久| 成人免费观看视频高清| 卡戴珊不雅视频在线播放| 99re6热这里在线精品视频| a级毛片免费高清观看在线播放| 久久久a久久爽久久v久久| 日韩不卡一区二区三区视频在线| 精品99又大又爽又粗少妇毛片| 国产精品成人在线| 欧美3d第一页| 国产精品国产三级专区第一集| av女优亚洲男人天堂| av.在线天堂| 18禁在线无遮挡免费观看视频| 亚洲欧美日韩东京热| 成人漫画全彩无遮挡| 人妻 亚洲 视频| 成人黄色视频免费在线看| 在线观看美女被高潮喷水网站| 99热网站在线观看| 欧美性感艳星| 国产乱来视频区| 国产欧美亚洲国产| 99热网站在线观看| 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 亚洲第一区二区三区不卡| 少妇的逼好多水| 美女高潮的动态| av天堂中文字幕网| 亚洲av男天堂| 水蜜桃什么品种好| 欧美+日韩+精品| 国产成人免费观看mmmm| 男人和女人高潮做爰伦理| 久久亚洲国产成人精品v| 视频区图区小说| 久久久色成人| 日日摸夜夜添夜夜添av毛片| 国产精品av视频在线免费观看| 色视频www国产| 汤姆久久久久久久影院中文字幕| 少妇 在线观看| 日韩一本色道免费dvd| 久久久久精品性色| 22中文网久久字幕| 少妇熟女欧美另类| 只有这里有精品99| 成人二区视频| 99热国产这里只有精品6| 国产精品蜜桃在线观看| 99热国产这里只有精品6| 欧美成人精品欧美一级黄| 国产 一区精品| a级毛色黄片| 免费人妻精品一区二区三区视频| 免费看av在线观看网站| 久久国产乱子免费精品| 中文字幕免费在线视频6| 欧美bdsm另类| 成年av动漫网址| 视频区图区小说| 欧美高清性xxxxhd video| 日韩精品有码人妻一区| 久久精品国产亚洲av天美| 久久精品国产自在天天线| 亚洲国产成人一精品久久久| 久久久久精品久久久久真实原创| 国产在线一区二区三区精| 日韩制服骚丝袜av| tube8黄色片| 国产精品一及| 黄色欧美视频在线观看| 免费av中文字幕在线| 午夜视频国产福利| 国产视频首页在线观看| 欧美老熟妇乱子伦牲交| 日本爱情动作片www.在线观看| 不卡视频在线观看欧美| 女人十人毛片免费观看3o分钟| 国产成人精品福利久久| 男男h啪啪无遮挡| 亚洲国产日韩一区二区| 亚洲成人手机| 成人亚洲欧美一区二区av| 国语对白做爰xxxⅹ性视频网站| 国产有黄有色有爽视频| 久久国内精品自在自线图片| 男女免费视频国产| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 最近2019中文字幕mv第一页| 午夜福利在线在线| 色视频www国产| 精品人妻视频免费看| 特大巨黑吊av在线直播| www.色视频.com| 欧美成人午夜免费资源| 97超视频在线观看视频| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 国产成人freesex在线| 舔av片在线| 国产精品99久久久久久久久| 在线观看免费高清a一片| 国产成人精品福利久久| 午夜福利影视在线免费观看| 亚洲欧美精品自产自拍| a 毛片基地| 国产永久视频网站| 麻豆乱淫一区二区| 3wmmmm亚洲av在线观看| 涩涩av久久男人的天堂| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 成年美女黄网站色视频大全免费 | 国产精品三级大全| 国产精品不卡视频一区二区| 国精品久久久久久国模美| 亚洲成色77777| 久久午夜福利片| 精品国产露脸久久av麻豆| 99久国产av精品国产电影| 边亲边吃奶的免费视频| 成年免费大片在线观看| 视频区图区小说| 乱码一卡2卡4卡精品| 国产亚洲91精品色在线| 久久毛片免费看一区二区三区| 日韩强制内射视频| 九草在线视频观看| 美女cb高潮喷水在线观看| 交换朋友夫妻互换小说| 热re99久久精品国产66热6| 3wmmmm亚洲av在线观看| 观看免费一级毛片| 激情五月婷婷亚洲| 日韩,欧美,国产一区二区三区| 日日摸夜夜添夜夜爱| 亚洲久久久国产精品| 男女国产视频网站| 亚洲国产成人一精品久久久| 特大巨黑吊av在线直播| 国产免费福利视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲性久久影院| 欧美精品国产亚洲| 大码成人一级视频| 国产日韩欧美在线精品| 亚洲精品国产成人久久av| 十八禁网站网址无遮挡 | 男人爽女人下面视频在线观看| 久久久久久九九精品二区国产| 亚洲欧美日韩另类电影网站 | 黄色视频在线播放观看不卡| 国产午夜精品久久久久久一区二区三区| 嘟嘟电影网在线观看| 菩萨蛮人人尽说江南好唐韦庄| 观看免费一级毛片| 男女边吃奶边做爰视频| 伦精品一区二区三区| 大片免费播放器 马上看| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区国产| 精品人妻视频免费看| 直男gayav资源| 黑人高潮一二区| 人妻 亚洲 视频| 黄色怎么调成土黄色| 久久久久久久久久人人人人人人| 精品人妻熟女av久视频| 美女主播在线视频| 观看免费一级毛片| 激情 狠狠 欧美| 精品酒店卫生间| 国内精品宾馆在线| 日产精品乱码卡一卡2卡三| 国产精品久久久久久久久免| 国产精品国产三级国产av玫瑰| 午夜视频国产福利| 91久久精品国产一区二区成人| 亚洲国产高清在线一区二区三| 久久精品夜色国产| 久久精品国产自在天天线| tube8黄色片| 搡老乐熟女国产| 久久久久久久久久久免费av| 国产精品嫩草影院av在线观看| 亚洲av中文av极速乱| 精品99又大又爽又粗少妇毛片| 大码成人一级视频| 欧美人与善性xxx| 精品久久久久久久末码| 日韩欧美 国产精品| 免费看日本二区| 18禁裸乳无遮挡动漫免费视频| 国产成人午夜福利电影在线观看| 国产免费视频播放在线视频| 交换朋友夫妻互换小说| 18禁动态无遮挡网站| 国产一区亚洲一区在线观看| 麻豆成人av视频| 一级毛片我不卡| 亚洲av中文av极速乱| 欧美成人午夜免费资源| 国产精品三级大全| 伦理电影大哥的女人| 亚洲国产色片| 在线观看免费日韩欧美大片 | 欧美极品一区二区三区四区| 午夜激情久久久久久久| 国产精品99久久久久久久久| 国产亚洲精品久久久com| 自拍偷自拍亚洲精品老妇| 亚洲图色成人| av一本久久久久| www.av在线官网国产| 久久久久网色| 看非洲黑人一级黄片| 又黄又爽又刺激的免费视频.| 亚洲国产精品一区三区| 成人亚洲欧美一区二区av| 免费黄色在线免费观看| www.色视频.com| 欧美高清成人免费视频www| 国产免费视频播放在线视频| 国产永久视频网站| 亚洲不卡免费看| 麻豆国产97在线/欧美| 黄色一级大片看看| 乱系列少妇在线播放| 26uuu在线亚洲综合色| 久久精品熟女亚洲av麻豆精品| 偷拍熟女少妇极品色| 国产精品久久久久久久电影| 热re99久久精品国产66热6| 一本—道久久a久久精品蜜桃钙片| 另类亚洲欧美激情| 能在线免费看毛片的网站| 亚洲真实伦在线观看| videos熟女内射| 日韩强制内射视频| 高清午夜精品一区二区三区| 欧美3d第一页| 啦啦啦在线观看免费高清www| 久久久久久久久久成人| 五月天丁香电影| 看非洲黑人一级黄片| 一级毛片aaaaaa免费看小| 啦啦啦视频在线资源免费观看| 国产爱豆传媒在线观看| 国产爽快片一区二区三区| 欧美国产精品一级二级三级 | 丰满人妻一区二区三区视频av| 国产毛片在线视频| 婷婷色综合www| 国产精品成人在线| 久久国产精品男人的天堂亚洲 | 制服丝袜香蕉在线| 18禁动态无遮挡网站| 色哟哟·www| 亚洲成色77777| 麻豆乱淫一区二区| 免费不卡的大黄色大毛片视频在线观看| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 又大又黄又爽视频免费| 国产又色又爽无遮挡免| 亚洲国产最新在线播放| 深夜a级毛片| 成人一区二区视频在线观看| 日日摸夜夜添夜夜爱| 久久精品人妻少妇| 精品亚洲成国产av| 亚洲精品久久午夜乱码| 国产永久视频网站| 久久精品国产a三级三级三级| 亚州av有码| 国产在线男女| 啦啦啦中文免费视频观看日本| 亚洲欧美一区二区三区黑人 | 免费观看的影片在线观看| 免费观看的影片在线观看| 亚洲精品一二三| 亚洲经典国产精华液单| 永久免费av网站大全| 国精品久久久久久国模美| 亚洲精华国产精华液的使用体验| 97超视频在线观看视频| 直男gayav资源| 观看av在线不卡| 日韩一区二区三区影片| 亚洲电影在线观看av| 亚洲精品国产av成人精品| 国产高清不卡午夜福利| 青青草视频在线视频观看| 欧美97在线视频| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产a三级三级三级| 国产精品三级大全| 老司机影院毛片| 水蜜桃什么品种好| 视频区图区小说| 国产精品人妻久久久久久| 亚洲国产最新在线播放| 91午夜精品亚洲一区二区三区| 中文字幕免费在线视频6| 一级a做视频免费观看| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 日日撸夜夜添| 伊人久久精品亚洲午夜| 丝袜脚勾引网站| 精品久久久久久电影网| 亚洲丝袜综合中文字幕| 好男人视频免费观看在线|