Shale parasequence analysis is an important part of sequence stratigraphy sudies. This paper proposed a systematic research method for analyzing shale parasequences including their delineation, division, characteristics and origins. The division method is established on the basis of lithofacies. Multi-method analysis and mutual verification were implemented by using auxiliary indicators(such as mineral compositions, geochemical indicators and wavelet values). A typical shale parasequence comprises a lower interval of deepening water-depth and an upper interval of shallowing water-depth(e.g., a shale parasequence including a high-total organic carbon(TOC) shale-low-TOC limy shale). Abrupt increases in pyrite content, TOC value, relative hydrocarbon generation potential((S_1+S_2)/TOC), and wavelet values are indicative of parasequence boundaries. The proposed research method was applied to study the upper fourth member of the Shahejie Formation in the Dongying depression, Bohai Bay Basin. Results show that there were seven types of parasequences developed. A singular and a dual structured parasequences were identified. Three factors controlling the development of the shale parasequences were identified including relative lake level change, terrestrial input and transgression. The development of high-TOC(>2%) shale parasequences was mainly controlled by biological and chemical sedimentation. The low-TOC(<2%) shale parasequences were mainly deposited by chemical sedimentation. The diversities of shale parasequences were caused by four major controlling factors including climate, relative lake level change, terrestrial input and emergency(e.g., transgression). Shale parasequence analysis is an important part of sequence stratigraphy sudies. This paper proposes a systematic research method for analyzing shale parasequences including their delineation, division, characteristics and origins. The division method is established on the basis of lithofacies. Multi-method analysis and mutual A typical shale parasequence includes a lower interval of deepening water-depth and an upper interval of shallowing water-depth (eg, a shale parasequence including a high-total organic carbon (TOC) shale-low-TOC limy shale). Abrupt increases in pyrite content, TOC value, relative hydrocarbon generation potential ((S_1 + S_2) / TOC), and wavelet values ​​are indicative of parasequence boundaries. proposed research method was applied to study the upper fourth member of the Shahejie Formation in the Dongying depression, Bohai Bay Basin Three factors controlling the development of the shale parasequences were identified including relative lake level change, terrestrial input and transgression. The development of high-TOC (> 2%) shale parasequences was mainly controlled by biological and chemical sedimentation. The low-TOC (<2%) shale parasequences were mainly deposited by chemical sedimentation. The diversities of shale parasequences were caused by the four major controlling factors including climate, relative lake level change, terrestrial input and emergency (eg, transgression).